--- ray/src/cv/bsdf2klems.c	2013/09/26 17:05:00	2.10
+++ ray/src/cv/bsdf2klems.c	2014/03/12 22:24:59	2.14
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: bsdf2klems.c,v 2.10 2013/09/26 17:05:00 greg Exp $";
+static const char RCSid[] = "$Id: bsdf2klems.c,v 2.14 2014/03/12 22:24:59 greg Exp $";
 #endif
 /*
  * Load measured BSDF interpolant and write out as XML file with Klems matrix.
@@ -27,9 +27,48 @@ static const char	*kbasis = "LBNL/Klems Full";
 static int		npsamps = 256;
 				/* limit on number of RBF lobes */
 static int		lobe_lim = 15000;
+				/* progress bar length */
+static int		do_prog = 79;
 
+
+/* Start new progress bar */
+#define prog_start(s)	if (do_prog) fprintf(stderr, "%s: %s...\n", progname, s); else
+
+/* Draw progress bar of the appropriate length */
+static void
+prog_show(double frac)
+{
+	char	pbar[256];
+	int	nchars;
+
+	if (do_prog <= 1) return;
+	if (do_prog > sizeof(pbar)-2)
+		do_prog = sizeof(pbar)-2;
+	if (frac < 0) frac = 0;
+	else if (frac > 1) frac = 1;
+	nchars = do_prog*frac + .5;
+	pbar[0] = '\r';
+	memset(pbar+1, '*', nchars);
+	memset(pbar+1+nchars, '-', do_prog-nchars);
+	pbar[do_prog+1] = '\0';
+	fputs(pbar, stderr);
+}
+
+/* Finish progress bar */
+static void
+prog_done(void)
+{
+	int	n = do_prog;
+
+	if (n <= 1) return;
+	fputc('\r', stderr);
+	while (n--)
+		fputc(' ', stderr);
+	fputc('\r', stderr);
+}
+
 /* Return angle basis corresponding to the given name */
-ANGLE_BASIS *
+static ANGLE_BASIS *
 get_basis(const char *bn)
 {
 	int	n = nabases;
@@ -58,8 +97,10 @@ xml_header(int ac, char *av[])
 static void
 xml_prologue(const SDData *sd)
 {
-	const char	*matn = (sd && sd->matn[0]) ? sd->matn : "Name";
-	const char	*makr = (sd && sd->makr[0]) ? sd->makr : "Manufacturer";
+	const char	*matn = (sd && sd->matn[0]) ? sd->matn :
+				bsdf_name[0] ? bsdf_name : "Unknown";
+	const char	*makr = (sd && sd->makr[0]) ? sd->makr :
+				bsdf_manuf[0] ? bsdf_manuf : "Unknown";
 	ANGLE_BASIS	*abp = get_basis(kbasis);
 	int		i;
 
@@ -295,8 +336,10 @@ eval_function(char *funame)
 		printf("\t%.3e\n", sum/npsamps);
 	    }
 	    putchar('\n');
+	    prog_show((j+1.)/abp->nangles);
 	}
 	data_epilogue();			/* finish output */
+	prog_done();
 }
 
 /* Interpolate and output a radial basis function BSDF representation */
@@ -331,12 +374,14 @@ eval_rbf(void)
 		    else
 			bo_getvec(vout, j+(n+frandom())/npsamps, abp);
 
-		    sum += eval_rbfrep(rbf, vout) / vout[2];
+		    sum += eval_rbfrep(rbf, vout);
 		}
-		bsdfarr[j*abp->nangles + i] = sum*output_orient/npsamps;
+		fo_getvec(vout, j+.5, abp);	/* use centered secant */
+		bsdfarr[j*abp->nangles + i] = sum / (npsamps*vout[2]);
 	    }
 	    if (rbf != NULL)
 		free(rbf);
+	    prog_show((i+1.)/abp->nangles);
 	}
 	n = 0;					/* write out our matrix */
 	for (j = 0; j < abp->nangles; j++) {
@@ -345,6 +390,7 @@ eval_rbf(void)
 	    putchar('\n');
 	}
 	data_epilogue();			/* finish output */
+	prog_done();
 }
 
 /* Read in BSDF and interpolate as Klems matrix representation */
@@ -390,6 +436,9 @@ main(int argc, char *argv[])
 		case 'l':
 			lobe_lim = atoi(argv[++i]);
 			break;
+		case 'p':
+			do_prog = atoi(argv[i]+2);
+			break;
 		default:
 			goto userr;
 		}
@@ -398,8 +447,7 @@ main(int argc, char *argv[])
 			fprintf(stderr,
 	"%s: need single function with 6 arguments: bsdf(ix,iy,iz,ox,oy,oz)\n",
 					progname);
-			fprintf(stderr, "\tor 3 arguments using Dx,Dy,Dz: bsdf(ix,iy,iz)\n",
-					progname);
+			fprintf(stderr, "\tor 3 arguments using Dx,Dy,Dz: bsdf(ix,iy,iz)\n");
 			goto userr;
 		}
 		++eclock;
@@ -408,16 +456,20 @@ main(int argc, char *argv[])
 		if (dofwd) {
 			input_orient = -1;
 			output_orient = -1;
-			eval_function(argv[i]);		/* outside reflectance */
+			prog_start("Evaluating outside reflectance");
+			eval_function(argv[i]);
 			output_orient = 1;
-			eval_function(argv[i]);		/* outside -> inside */
+			prog_start("Evaluating outside->inside transmission");
+			eval_function(argv[i]);
 		}
 		if (dobwd) {
 			input_orient = 1;
 			output_orient = 1;
-			eval_function(argv[i]);		/* inside reflectance */
+			prog_start("Evaluating inside reflectance");
+			eval_function(argv[i]);
 			output_orient = -1;
-			eval_function(argv[i]);		/* inside -> outside */
+			prog_start("Evaluating inside->outside transmission");
+			eval_function(argv[i]);
 		}
 		xml_epilogue();			/* finish XML output & exit */
 		return(0);
@@ -433,6 +485,7 @@ main(int argc, char *argv[])
 	if (i < argc) {				/* open input files if given */
 		int	nbsdf = 0;
 		for ( ; i < argc; i++) {	/* interpolate each component */
+			char	pbuf[256];
 			FILE	*fpin = fopen(argv[i], "rb");
 			if (fpin == NULL) {
 				fprintf(stderr, "%s: cannot open BSDF interpolant '%s'\n",
@@ -446,6 +499,8 @@ main(int argc, char *argv[])
 				xml_header(argc, argv);
 				xml_prologue(NULL);
 			}
+			sprintf(pbuf, "Interpolating component '%s'", argv[i]);
+			prog_start(pbuf);
 			eval_rbf();
 		}
 		xml_epilogue();			/* finish XML output & exit */
@@ -456,6 +511,7 @@ main(int argc, char *argv[])
 		return(1);
 	xml_header(argc, argv);			/* start XML output */
 	xml_prologue(NULL);
+	prog_start("Interpolating from standard input");
 	eval_rbf();				/* resample dist. */
 	xml_epilogue();				/* finish XML output & exit */
 	return(0);