--- ray/src/cv/bsdf2ttree.c	2013/11/26 17:33:55	2.22
+++ ray/src/cv/bsdf2ttree.c	2014/03/11 19:37:45	2.25
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: bsdf2ttree.c,v 2.22 2013/11/26 17:33:55 greg Exp $";
+static const char RCSid[] = "$Id: bsdf2ttree.c,v 2.25 2014/03/11 19:37:45 greg Exp $";
 #endif
 /*
  * Load measured BSDF interpolant and write out as XML file with tensor tree.
@@ -11,6 +11,7 @@ static const char RCSid[] = "$Id: bsdf2ttree.c,v 2.22 
 #include <stdio.h>
 #include <stdlib.h>
 #include <math.h>
+#include "random.h"
 #include "platform.h"
 #include "rtprocess.h"
 #include "calcomp.h"
@@ -134,6 +135,9 @@ eval_isotropic(char *funame)
 			exit(1);
 		}
 		SET_FILE_BINARY(ofp);
+#ifdef getc_unlocked				/* avoid lock/unlock overhead */
+		flockfile(ofp);
+#endif
 	} else
 		fputs("{\n", stdout);
 						/* need to assign Dx, Dy, Dz? */
@@ -169,7 +173,7 @@ eval_isotropic(char *funame)
 			    if (abs_diff(bsdf, last_bsdf) > ssamp_thresh) {
 				sum = 0;	/* super-sample voxel */
 				for (ssi = nssamp; ssi--; ) {
-				    SDmultiSamp(ssa, 3, (ssi+drand48())/nssamp);
+				    SDmultiSamp(ssa, 3, (ssi+frandom())/nssamp);
 				    ssvec[0] = 2.*(ix+ssa[0])/sqres - 1.;
 				    ssvec[1] = .0;
 				    ssvec[2] = input_orient *
@@ -238,6 +242,10 @@ eval_anisotropic(char *funame)
 					progname);
 			exit(1);
 		}
+		SET_FILE_BINARY(ofp);
+#ifdef getc_unlocked				/* avoid lock/unlock overhead */
+		flockfile(ofp);
+#endif
 	} else
 		fputs("{\n", stdout);
 						/* need to assign Dx, Dy, Dz? */
@@ -274,7 +282,7 @@ eval_anisotropic(char *funame)
 			    if (abs_diff(bsdf, last_bsdf) > ssamp_thresh) {
 				sum = 0;	/* super-sample voxel */
 				for (ssi = nssamp; ssi--; ) {
-				    SDmultiSamp(ssa, 4, (ssi+drand48())/nssamp);
+				    SDmultiSamp(ssa, 4, (ssi+frandom())/nssamp);
 				    SDsquare2disk(ssvec, 1.-(ix+ssa[0])/sqres,
 						1.-(iy+ssa[1])/sqres);
 				    ssvec[2] = output_orient *
@@ -374,8 +382,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;