--- ray/src/rt/rpict.c	1990/11/01 12:06:01	1.23
+++ ray/src/rt/rpict.c	1991/05/21 17:41:35	1.25
@@ -1,4 +1,4 @@
-/* Copyright (c) 1986 Regents of the University of California */
+/* Copyright (c) 1991 Regents of the University of California */
 
 #ifndef lint
 static char SCCSid[] = "$SunId$ LBL";
@@ -24,6 +24,10 @@ static char SCCSid[] = "$SunId$ LBL";
 
 #include  "random.h"
 
+int  dimlist[MAXDIM];			/* sampling dimensions */
+int  ndims = 0;				/* number of sampling dimensions */
+int  samplendx;				/* sample index number */
+
 VIEW  ourview = STDVIEW;		/* view parameters */
 int  hresolu = 512;			/* horizontal resolution */
 int  vresolu = 512;			/* vertical resolution */
@@ -32,6 +36,7 @@ double  pixaspect = 1.0;		/* pixel aspect ratio */
 int  psample = 4;			/* pixel sample size */
 double  maxdiff = .05;			/* max. difference for interpolation */
 double  dstrpix = 0.67;			/* square pixel distribution */
+int  psuper = 2;			/* pixel super-sampling rate */
 
 double  dstrsrc = 0.0;			/* square source distribution */
 double  shadthresh = .05;		/* shadow threshold */
@@ -53,6 +58,11 @@ int  ralrm = 0;				/* seconds between reports */
 
 double  pctdone = 0.0;			/* percentage done */
 
+long  tlastrept = 0L;			/* time at last report */
+
+extern long  time();
+extern long  tstart;			/* starting time */
+
 extern long  nrays;			/* number of rays traced */
 
 #define  MAXDIV		32		/* maximum sample size */
@@ -72,9 +82,9 @@ int  code;
 }
 
 
+#ifdef BSD
 report()		/* report progress */
 {
-#ifdef BSD
 	struct rusage  rubuf;
 	double  t;
 
@@ -87,15 +97,19 @@ report()		/* report progress */
 
 	sprintf(errmsg, "%ld rays, %4.2f%% done after %5.4f CPU hours\n",
 			nrays, pctdone, t/3600.0);
+	eputs(errmsg);
+	tlastrept = time((long *)0);
+}
 #else
+report()		/* report progress */
+{
 	signal(SIGALRM, report);
-	sprintf(errmsg, "%ld rays, %4.2f%% done\n", nrays, pctdone);
-#endif
+	tlastrept = time((long *)0);
+	sprintf(errmsg, "%ld rays, %4.2f%% done after %5.4f hours\n",
+			nrays, pctdone, (tlastrept-tstart)/3600.0);
 	eputs(errmsg);
-
-	if (ralrm > 0)
-		alarm(ralrm);
 }
+#endif
 
 
 render(zfile, oldfile)				/* render the scene */
@@ -144,13 +158,14 @@ char  *zfile, *oldfile;
 	if (zfd != -1 && i > 0 &&
 			lseek(zfd, (long)i*hresolu*sizeof(float), 0) == -1)
 		error(SYSTEM, "z file seek error in render");
+	pctdone = 100.0*i/vresolu;
+	if (ralrm > 0)			/* report init stats */
+		report();
 	ypos = vresolu-1 - i;
 	fillscanline(scanbar[0], zbar[0], hresolu, ypos, psample);
 	ystep = psample;
 						/* compute scanlines */
 	for (ypos -= ystep; ypos > -ystep; ypos -= ystep) {
-							/* record progress */
-		pctdone = 100.0*(vresolu-1-ypos-ystep)/vresolu;
 							/* bottom adjust? */
 		if (ypos < 0) {
 			ystep += ypos;
@@ -177,6 +192,10 @@ char  *zfile, *oldfile;
 		}
 		if (fflush(stdout) == EOF)
 			goto writerr;
+							/* record progress */
+		pctdone = 100.0*(vresolu-1-ypos)/vresolu;
+		if (ralrm > 0 && time((long *)0) >= tlastrept+ralrm)
+			report();
 	}
 						/* clean up */
 	if (zfd != -1) {
@@ -319,7 +338,7 @@ pixvalue(col, x, y)		/* compute pixel value */
 COLOR  col;			/* returned color */
 int  x, y;			/* pixel position */
 {
-	static RAY  thisray;	/* our ray for this pixel */
+	static RAY  thisray;
 
 	if (viewray(thisray.rorg, thisray.rdir, &ourview,
 			(x+pixjitter())/hresolu, (y+pixjitter())/vresolu) < 0) {
@@ -328,7 +347,9 @@ int  x, y;			/* pixel position */
 	}
 
 	rayorigin(&thisray, NULL, PRIMARY, 1.0);
-	
+
+	samplendx = 3*y + x;			/* set pixel index */
+
 	rayvalue(&thisray);			/* trace ray */
 
 	copycolor(col, thisray.rcol);		/* return color */