--- ray/src/gen/gensurf.c	2007/10/18 03:59:54	2.18
+++ ray/src/gen/gensurf.c	2012/03/11 22:27:17	2.21
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: gensurf.c,v 2.18 2007/10/18 03:59:54 greg Exp $";
+static const char RCSid[] = "$Id: gensurf.c,v 2.21 2012/03/11 22:27:17 greg Exp $";
 #endif
 /*
  *  gensurf.c - program to generate functional surfaces
@@ -35,7 +35,7 @@ char  VNAME[] = 	"valid";		/* valid vertex name */
 #define  pvect(p)	printf(vformat, (p)[0], (p)[1], (p)[2])
 
 char  vformat[] = "%18.12g %18.12g %18.12g\n";
-char  tsargs[] = "4 surf_dx surf_dy surf_dz surf.cal\n";
+char  tsargs[] = "4 surf_dx surf_dy surf_dz surf.cal";
 char  texname[] = "Phong";
 
 int  smooth = 0;		/* apply smoothing? */
@@ -58,6 +58,7 @@ double  l_hermite(), l_bezier(), l_bspline(), l_datava
 
 typedef struct {
 	int  valid;	/* point is valid (vertex number) */
+	int  nvalid;	/* normal is valid */
 	FVECT  p;	/* vertex position */
 	FVECT  n;	/* average normal */
 	RREAL  uv[2];	/* (u,v) position */
@@ -229,7 +230,7 @@ loaddata(		/* load point data from file */
 		size--;
 	}
 	if (size == (m+n+1)*pointsize) {	/* no border after all */
-		dp = (RREAL *)realloc((void *)datarec.data,
+		dp = (RREAL *)realloc(datarec.data,
 				m*n*pointsize*sizeof(RREAL));
 		if (dp != NULL)
 			datarec.data = dp;
@@ -298,15 +299,19 @@ putobjrow(			/* output vertex row to .OBJ */
 )
 {
 	static int	nverts = 0;
+	static int	nnorms = 0;
 
 	for ( ; n-- >= 0; rp++) {
 		if (!rp->valid)
 			continue;
 		fputs("v ", stdout);
 		pvect(rp->p);
-		if (smooth && !ZEROVECT(rp->n))
+		if (smooth && !ZEROVECT(rp->n)) {
 			printf("\tvn %.9g %.9g %.9g\n",
 					rp->n[0], rp->n[1], rp->n[2]);
+			rp->nvalid = ++nnorms;
+		} else
+			rp->nvalid = 0;
 		printf("\tvt %.9g %.9g\n", rp->uv[0], rp->uv[1]);
 		rp->valid = ++nverts;
 	}
@@ -344,35 +349,24 @@ putsquare(		/* put out a square */
 	if (!(ok1 | ok2))
 		return;
 	if (objout) {			/* output .OBJ faces */
-		int	p0n=0, p1n=0, p2n=0, p3n=0;
-		if (smooth) {
-			if (!ZEROVECT(p0->n))
-				p0n = p0->valid;
-			if (!ZEROVECT(p1->n))
-				p1n = p1->valid;
-			if (!ZEROVECT(p2->n))
-				p2n = p2->valid;
-			if (!ZEROVECT(p3->n))
-				p3n = p3->valid;
-		}
 		if (ok1 & ok2 && fdot(vc1,vc2) >= 1.0-FTINY*FTINY) {
 			printf("f %d/%d/%d %d/%d/%d %d/%d/%d %d/%d/%d\n",
-					p0->valid, p0->valid, p0n,
-					p1->valid, p1->valid, p1n,
-					p3->valid, p3->valid, p3n,
-					p2->valid, p2->valid, p2n);
+					p0->valid, p0->valid, p0->nvalid,
+					p1->valid, p1->valid, p1->nvalid,
+					p3->valid, p3->valid, p3->nvalid,
+					p2->valid, p2->valid, p2->nvalid);
 			return;
 		}
 		if (ok1)
 			printf("f %d/%d/%d %d/%d/%d %d/%d/%d\n",
-					p0->valid, p0->valid, p0n,
-					p1->valid, p1->valid, p1n,
-					p2->valid, p2->valid, p2n);
+					p0->valid, p0->valid, p0->nvalid,
+					p1->valid, p1->valid, p1->nvalid,
+					p2->valid, p2->valid, p2->nvalid);
 		if (ok2)
 			printf("f %d/%d/%d %d/%d/%d %d/%d/%d\n",
-					p2->valid, p2->valid, p2n,
-					p1->valid, p1->valid, p1n,
-					p3->valid, p3->valid, p3n);
+					p2->valid, p2->valid, p2->nvalid,
+					p1->valid, p1->valid, p1->nvalid,
+					p3->valid, p3->valid, p3->nvalid);
 		return;
 	}
 					/* compute normal interpolation */
@@ -382,8 +376,7 @@ putsquare(		/* put out a square */
 	if (ok1 & ok2 && fdot(vc1,vc2) >= 1.0-FTINY*FTINY) {
 		printf("\n%s ", modname);
 		if (axis != -1) {
-			printf("texfunc %s\n", texname);
-			printf(tsargs);
+			printf("texfunc %s\n%s\n", texname, tsargs);
 			printf("0\n13\t%d\n", axis);
 			pvect(norm[0]);
 			pvect(norm[1]);
@@ -405,8 +398,7 @@ putsquare(		/* put out a square */
 	if (ok1) {
 		printf("\n%s ", modname);
 		if (axis != -1) {
-			printf("texfunc %s\n", texname);
-			printf(tsargs);
+			printf("texfunc %s\n%s\n", texname, tsargs);
 			printf("0\n13\t%d\n", axis);
 			pvect(norm[0]);
 			pvect(norm[1]);
@@ -424,8 +416,7 @@ putsquare(		/* put out a square */
 	if (ok2) {
 		printf("\n%s ", modname);
 		if (axis != -1) {
-			printf("texfunc %s\n", texname);
-			printf(tsargs);
+			printf("texfunc %s\n%s\n", texname, tsargs);
 			printf("0\n13\t%d\n", axis);
 			pvect(norm[0]);
 			pvect(norm[1]);