--- ray/src/common/objutil.c	2021/04/07 16:40:28	2.16
+++ ray/src/common/objutil.c	2022/01/15 02:00:21	2.21
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: objutil.c,v 2.16 2021/04/07 16:40:28 greg Exp $";
+static const char RCSid[] = "$Id: objutil.c,v 2.21 2022/01/15 02:00:21 greg Exp $";
 #endif
 /*
  *  Basic .OBJ scene handling routines.
@@ -573,6 +573,7 @@ clearComments(Scene *sc)
 	while (sc->ndescr > 0)
 		freeqstr(sc->descr[--sc->ndescr]);
 	efree((char *)sc->descr);
+	sc->descr = NULL;
 	sc->ndescr = 0;
 }
 
@@ -700,6 +701,89 @@ addFace(Scene *sc, VNDX vid[], int nv)
 	return(f);
 }
 
+/* Get neighbor vertices: malloc array with valence in index[0] */
+int *
+getVertexNeighbors(Scene *sc, int vid)
+{
+	int	*varr;
+	Face	*f;
+	int	j=0;
+
+	if (sc == NULL || (vid < 0) | (vid >= sc->nverts) ||
+			(f = sc->vert[vid].vflist) == NULL)
+		return(NULL);
+
+	varr = (int *)emalloc(sizeof(int)*CHUNKSIZ);
+	varr[0] = 0;				/* add unique neighbors */
+	for ( ; f != NULL; f = f->v[j].fnext) {
+		int	i, cvi;
+		for (j = f->nv; j--; )		/* find ourself in poly */
+			if (f->v[j].vid == vid)
+				break;
+		if (j < 0)			/* this is an error */
+			break;
+		if (f->nv < 3)			/* also never happens */
+			continue;
+						/* check previous neighbor */
+		cvi = f->v[ (j > 0) ? j-1 : f->nv-1 ].vid;
+		for (i = varr[0]+1; --i; )	/* making sure not in list */
+			if (varr[i] == cvi)
+				break;
+		if (!i) {			/* new neighbor? */
+			varr = chunk_alloc(int, varr, varr[0]+1);
+			varr[++varr[0]] = cvi;
+		}
+						/* check next neighbor */
+		cvi = f->v[ (j < f->nv-1) ? j+1 : 0 ].vid;
+		for (i = varr[0]+1; --i; )
+			if (varr[i] == cvi)
+				break;
+		if (!i) {
+			varr = chunk_alloc(int, varr, varr[0]+1);
+			varr[++varr[0]] = cvi;
+		}
+	}
+	if (!varr[0]) {
+		efree((char *)varr);		/* something went awry */
+		return(NULL);
+	}
+						/* tighten array & return */
+	return((int *)erealloc(varr, sizeof(int)*(varr[0]+1)));
+}
+
+/* Callback for growBoundingBox() */
+static int
+addBBox(Scene *sc, Face *f, void *p)
+{
+	double	(*bbox)[3] = (double (*)[3])p;
+	int	i, j;
+
+	for (i = f->nv; i-- > 0; ) {
+		double	*p3 = sc->vert[f->v[i].vid].p;
+		for (j = 3; j--; ) {
+			if (p3[j] < bbox[0][j])
+				bbox[0][j] = p3[j];
+			if (p3[j] > bbox[1][j])
+				bbox[1][j] = p3[j];
+		}
+	}
+	return(1);
+}
+
+/* Expand bounding box min & max (initialize bbox to all zeroes) */
+int
+growBoundingBox(Scene *sc, double bbox[2][3], int flreq, int flexc)
+{
+	if (sc == NULL || sc->nfaces <= 0 || bbox == NULL)
+		return(0);
+
+	if (VABSEQ(bbox[0], bbox[1])) {		/* first run */
+		bbox[0][0] = bbox[0][1] = bbox[0][2] = FHUGE;
+		bbox[1][0] = bbox[1][1] = bbox[1][2] = -FHUGE;
+	}
+	return(foreachFace(sc, addBBox, flreq, flexc, bbox));
+}
+
 /* Allocate an empty scene */
 Scene *
 newScene(void)
@@ -730,14 +814,14 @@ dupScene(const Scene *osc, int flreq, int flexc)
 	for (i = 0; i < osc->ndescr; i++)
 		addComment(sc, osc->descr[i]);
 	if (osc->ngrps > 1) {
-		sc->grpname = (char **)erealloc((char *)sc->grpname,
+		sc->grpname = (char **)erealloc(sc->grpname,
 				sizeof(char *) * (osc->ngrps+(CHUNKSIZ-1)));
 		for (i = 1; i < osc->ngrps; i++)
 			sc->grpname[i] = savqstr(osc->grpname[i]);
 		sc->ngrps = osc->ngrps;
 	}
 	if (osc->nmats > 1) {
-		sc->matname = (char **)erealloc((char *)sc->matname,
+		sc->matname = (char **)erealloc(sc->matname,
 				sizeof(char *) * (osc->nmats+(CHUNKSIZ-1)));
 		for (i = 1; i < osc->nmats; i++)
 			sc->matname[i] = savqstr(osc->matname[i]);
@@ -799,7 +883,7 @@ addScene(Scene *scdst, const Scene *scsrc)
 	if (scsrc->nfaces <= 0)
 		return(0);
 					/* map vertices */
-	vert_map = (int *)emalloc(sizeof(int *)*scsrc->nverts);
+	vert_map = (int *)emalloc(sizeof(int)*scsrc->nverts);
 	for (i = 0; i < scsrc->nverts; i++) {
 		const Vertex	*v = scsrc->vert + i;
 		if (v->vflist == NULL) {
@@ -810,14 +894,14 @@ addScene(Scene *scdst, const Scene *scsrc)
 	}
 	tex_off = scdst->ntex;		/* append texture coords */
 	if (scsrc->ntex > 0) {
-		scdst->tex = (TexCoord *)erealloc((char *)scdst->tex,
+		scdst->tex = (TexCoord *)erealloc(scdst->tex,
 			sizeof(TexCoord)*(tex_off+scsrc->ntex+(CHUNKSIZ-1)));
 		memcpy(scdst->tex+tex_off, scsrc->tex,
 				sizeof(TexCoord)*scsrc->ntex);
 	}
 	norm_off = scdst->nnorms;	/* append normals */
 	if (scsrc->nnorms > 0) {
-		scdst->norm = (Normal *)erealloc((char *)scdst->norm,
+		scdst->norm = (Normal *)erealloc(scdst->norm,
 			sizeof(Normal)*(norm_off+scsrc->nnorms+(CHUNKSIZ-1)));
 		memcpy(scdst->norm+norm_off, scsrc->norm,
 				sizeof(Normal)*scsrc->nnorms);
@@ -830,12 +914,10 @@ addScene(Scene *scdst, const Scene *scsrc)
 		if (f->mat != cur_mat)
 			setMaterial(scdst, scsrc->matname[cur_mat = f->mat]);
 		if (f->nv > vllen) {
-			if (vlist == my_vlist)
-				vlist = (VNDX *)emalloc(
-						sizeof(VNDX)*(vllen = f->nv));
-			else
-				vlist = (VNDX *)erealloc((char *)vlist,
-						sizeof(VNDX)*(vllen = f->nv));
+			vlist = (VNDX *)( vlist == my_vlist ?
+				emalloc(sizeof(VNDX)*f->nv) :
+				erealloc(vlist, sizeof(VNDX)*f->nv) );
+			vllen = f->nv;
 		}
 		memset(vlist, 0xff, sizeof(VNDX)*f->nv);
 		for (i = f->nv; i-- > 0; ) {
@@ -849,8 +931,8 @@ addScene(Scene *scdst, const Scene *scsrc)
 		fcnt += (addFace(scdst, vlist, f->nv) != NULL);
 	}
 					/* clean up */
-	if (vlist != my_vlist) efree((char *)vlist);
-	efree((char *)vert_map);
+	if (vlist != my_vlist) efree(vlist);
+	efree(vert_map);
 	return(fcnt);
 }
 
@@ -910,7 +992,7 @@ xfmScene(Scene *sc, const char *xfm)
 	if (!*xfm)
 		return(0);
 					/* parse string into words */
-	xav[0] = strcpy((char *)emalloc(strlen(xfm)+1), xfm);
+	xav[0] = savqstr((char *)xfm);
 	xac = 1; i = 0;
 	for ( ; ; ) {
 		while (!isspace(xfm[++i]))
@@ -921,14 +1003,14 @@ xfmScene(Scene *sc, const char *xfm)
 		if (!xfm[i])
 			break;
 		if (xac >= MAXAC-1) {
-			free(xav[0]);
+			freeqstr(xav[0]);
 			return(0);
 		}
 		xav[xac++] = xav[0] + i;
 	}
 	xav[xac] = NULL;
 	i = xfScene(sc, xac, xav);
-	efree((char *)xav[0]);
+	freeqstr(xav[0]);
 	return(i);
 }
 #undef MAXAC
@@ -961,7 +1043,7 @@ deleteUnreferenced(Scene *sc)
 	}
 	if (nused == sc->nverts)
 		goto skip_pos;
-	sc->vert = (Vertex *)erealloc((char *)sc->vert,
+	sc->vert = (Vertex *)erealloc(sc->vert,
 				sizeof(Vertex)*(nused+(CHUNKSIZ-1)));
 	sc->nverts = nused;
 	for (f = sc->flist; f != NULL; f = f->next)
@@ -987,7 +1069,7 @@ skip_pos:
 	}
 	if (nused == sc->ntex)
 		goto skip_tex;
-	sc->tex = (TexCoord *)erealloc((char *)sc->tex,
+	sc->tex = (TexCoord *)erealloc(sc->tex,
 				sizeof(TexCoord)*(nused+(CHUNKSIZ-1)));
 	sc->ntex = nused;
 	for (f = sc->flist; f != NULL; f = f->next)
@@ -1012,7 +1094,7 @@ skip_tex:
 	}
 	if (nused == sc->nnorms)
 		goto skip_norms;
-	sc->norm = (Normal *)erealloc((char *)sc->norm,
+	sc->norm = (Normal *)erealloc(sc->norm,
 				sizeof(Normal)*(nused+(CHUNKSIZ-1)));
 	sc->nnorms = nused;
 	for (f = sc->flist; f != NULL; f = f->next)
@@ -1021,7 +1103,7 @@ skip_tex:
 				f->v[i].nid = vmap[f->v[i].nid];
 skip_norms:
 						/* clean up */
-	efree((char *)vmap);
+	efree(vmap);
 }
 
 /* Free a scene */
@@ -1036,16 +1118,16 @@ freeScene(Scene *sc)
 	clearComments(sc);
 	for (i = sc->ngrps; i-- > 0; )
 		freeqstr(sc->grpname[i]);
-	efree((char *)sc->grpname);
+	efree(sc->grpname);
 	for (i = sc->nmats; i-- > 0; )
 		freeqstr(sc->matname[i]);
-	efree((char *)sc->matname);
-	efree((char *)sc->vert);
-	efree((char *)sc->tex);
-	efree((char *)sc->norm);
+	efree(sc->matname);
+	efree(sc->vert);
+	efree(sc->tex);
+	efree(sc->norm);
 	while ((f = sc->flist) != NULL) {
 		sc->flist = f->next;
-		efree((char *)f);
+		efree(f);
 	}
-	efree((char *)sc);
+	efree(sc);
 }