--- ray/src/px/pcond.c	1996/10/11 10:46:57	3.5
+++ ray/src/px/pcond.c	1997/01/08 21:40:38	3.7
@@ -1,4 +1,4 @@
-/* Copyright (c) 1996 Regents of the University of California */
+/* Copyright (c) 1997 Regents of the University of California */
 
 #ifndef lint
 static char SCCSid[] = "$SunId$ LBL";
@@ -26,6 +26,7 @@ char	*progname;			/* global argv[0] */
 
 char	*infn;				/* input file name */
 FILE	*infp;				/* input stream */
+FILE	*mapfp = NULL;			/* tone-mapping function stream */
 VIEW	ourview = STDVIEW;		/* picture view */
 int	gotview = 0;			/* picture has view */
 double	pixaspect = 1.0;		/* pixel aspect ratio */
@@ -113,6 +114,15 @@ char	*argv[];
 			if (i+1 >= argc) goto userr;
 			ldmin = atof(argv[++i]);
 			break;
+		case 'm':
+			if (i+1 >= argc) goto userr;
+			if ((mapfp = fopen(argv[++i], "w")) == NULL) {
+				fprintf(stderr,
+					"%s: cannot open for writing\n",
+						argv[i]);
+				exit(1);
+			}
+			break;
 		default:
 			goto userr;
 		}
@@ -147,9 +157,11 @@ char	*argv[];
 	if (outprims != inprims)
 		fputprims(outprims, stdout);
 	mapimage();			/* map the picture */
+	if (mapfp != NULL)		/* write out basic mapping */
+		putmapping(mapfp);
 	exit(0);
 userr:
-	fprintf(stderr, "Usage: %s [-{h|a|v|s|c|l|w}[+-]][-e ev][-p xr yr xg yg xb yb xw yw|-f mbf.cal][-t Ldmax][-b Ldmin] inpic [outpic]\n",
+	fprintf(stderr, "Usage: %s [-{h|a|v|s|c|l|w}[+-]][-e ev][-p xr yr xg yg xb yb xw yw|-f mbf.cal][-t Ldmax][-b Ldmin][-m mapfile] inpic [outpic]\n",
 			progname);
 	exit(1);
 #undef bool
@@ -345,8 +357,8 @@ fovhist()			/* create foveal sampled image and histogr
 	if (lwmin < LMIN) lwmin = LMIN;
 	if (lwmax > LMAX) lwmax = LMAX;
 					/* compute histogram */
-	bwmin = Bl(lwmin)*(1. - .01/HISTRES);
-	bwmax = Bl(lwmax)*(1. + .01/HISTRES);
+	bwmin = Bl(lwmin);
+	bwmax = Bl(lwmax);
 	bwavg = 0.;
 	for (y = 0; y < fvyr; y++)
 		for (x = 0; x < fvxr; x++) {
@@ -354,9 +366,9 @@ fovhist()			/* create foveal sampled image and histogr
 					frandom() > centprob(x,y))
 				continue;
 			l = plum(fovscan(y)[x]);
+			if (l < lwmin) continue;
+			if (l > lwmax) continue;
 			b = Bl(l);
-			if (b < bwmin) continue;
-			if (b > bwmax) continue;
 			bwavg += b;
 			bwhc(b)++;
 			histot++;