--- ray/src/rt/rpict.c	1989/06/07 21:00:57	1.4
+++ ray/src/rt/rpict.c	1989/10/03 12:17:09	1.10
@@ -29,13 +29,14 @@ double  dstrpix = 0.67;			/* square pixel distribution
 
 double  dstrsrc = 0.0;			/* square source distribution */
 double  shadthresh = .05;		/* shadow threshold */
+double  shadcert = .5;			/* shadow certainty */
 
 int  maxdepth = 6;			/* maximum recursion depth */
 double  minweight = 5e-3;		/* minimum ray weight */
 
 COLOR  ambval = BLKCOLOR;		/* ambient value */
 double  ambacc = 0.2;			/* ambient accuracy */
-int  ambres = 128;			/* ambient resolution */
+int  ambres = 32;			/* ambient resolution */
 int  ambdiv = 128;			/* ambient divisions */
 int  ambssamp = 0;			/* ambient super-samples */
 int  ambounce = 0;			/* ambient bounces */
@@ -95,45 +96,58 @@ char  *oldfile;
 	int  ypos;			/* current scanline */
 	COLOR  *colptr;
 	register int  i;
-
+					/* check sampling */
 	if (psample < 1)
 		psample = 1;
 	else if (psample > MAXDIV)
 		psample = MAXDIV;
-
-	ourview.hresolu -= ourview.hresolu % psample;
-	ourview.vresolu -= ourview.vresolu % psample;
-		
+					/* allocate scanlines */
 	for (i = 0; i <= psample; i++) {
-		scanbar[i] = (COLOR *)malloc((ourview.hresolu+1)*sizeof(COLOR));
+		scanbar[i] = (COLOR *)malloc(ourview.hresolu*sizeof(COLOR));
 		if (scanbar[i] == NULL)
 			error(SYSTEM, "out of memory in render");
 	}
-	
 					/* write out boundaries */
-	printf("-Y %d +X %d\n", ourview.vresolu, ourview.hresolu);
-
-	ypos = ourview.vresolu - salvage(oldfile);	/* find top line */
-	fillscanline(scanbar[0], ypos, psample);	/* top scan */
-
-	for (ypos -= psample; ypos > -psample; ypos -= psample) {
+	fputresolu(YMAJOR|YDECR, ourview.hresolu, ourview.vresolu, stdout);
+					/* recover file and compute first */
+	ypos = ourview.vresolu-1 - salvage(oldfile);
+	fillscanline(scanbar[0], ourview.hresolu, ypos, psample);
+						/* compute scanlines */
+	for (ypos -= psample; ypos >= 0; ypos -= psample) {
 	
-		colptr = scanbar[psample];		/* get last scanline */
+		pctdone = 100.0*(ourview.vresolu-ypos-psample)/ourview.vresolu;
+
+		colptr = scanbar[psample];		/* move base to top */
 		scanbar[psample] = scanbar[0];
 		scanbar[0] = colptr;
-
-		fillscanline(scanbar[0], ypos, psample);	/* base scan */
-	
-		fillscanbar(scanbar, ypos, psample);
-		
-		for (i = psample-1; i >= 0; i--)
+							/* fill base line */
+		fillscanline(scanbar[0], ourview.hresolu, ypos, psample);
+							/* fill bar */
+		fillscanbar(scanbar, ourview.hresolu, ypos, psample);
+							/* write it out */
+		for (i = psample; i > 0; i--)
 			if (fwritescan(scanbar[i], ourview.hresolu, stdout) < 0)
 				goto writerr;
 		if (fflush(stdout) == EOF)
 			goto writerr;
-		pctdone = 100.0*(ourview.vresolu-ypos)/ourview.vresolu;
 	}
-		
+						/* compute residual */
+	colptr = scanbar[psample];
+	scanbar[psample] = scanbar[0];
+	scanbar[0] = colptr;
+	if (ypos > -psample) {
+		fillscanline(scanbar[-ypos], ourview.hresolu,
+				0, psample);
+		fillscanbar(scanbar-ypos, ourview.hresolu,
+				0, psample+ypos);
+	}
+	for (i = psample; i+ypos >= 0; i--)
+		if (fwritescan(scanbar[i], ourview.hresolu, stdout) < 0)
+			goto writerr;
+	if (fflush(stdout) == EOF)
+		goto writerr;
+	pctdone = 100.0;
+						/* free scanlines */
 	for (i = 0; i <= psample; i++)
 		free((char *)scanbar[i]);
 	return;
@@ -142,33 +156,38 @@ writerr:
 }
 
 
-fillscanline(scanline, y, xstep)		/* fill scan line at y */
+fillscanline(scanline, xres, y, xstep)		/* fill scan line at y */
 register COLOR  *scanline;
-int  y, xstep;
+int  xres, y, xstep;
 {
 	int  b = xstep;
 	register int  i;
 	
 	pixvalue(scanline[0], 0, y);
 
-	for (i = xstep; i <= ourview.hresolu; i += xstep) {
+	for (i = xstep; i < xres; i += xstep) {
 	
 		pixvalue(scanline[i], i, y);
 		
 		b = fillsample(scanline+i-xstep, i-xstep, y, xstep, 0, b/2);
 	}
+	if (i-xstep < xres-1) {
+		pixvalue(scanline[xres-1], xres-1, y);
+		fillsample(scanline+i-xstep, i-xstep, y,
+				xres-1-(i-xstep), 0, b/2);
+	}
 }
 
 
-fillscanbar(scanbar, y, ysize)		/* fill interior */
+fillscanbar(scanbar, xres, y, ysize)		/* fill interior */
 register COLOR  *scanbar[];
-int  y, ysize;
+int  xres, y, ysize;
 {
 	COLOR  vline[MAXDIV+1];
 	int  b = ysize;
 	register int  i, j;
 	
-	for (i = 0; i < ourview.hresolu; i++) {
+	for (i = 0; i < xres; i++) {
 		
 		copycolor(vline[0], scanbar[0][i]);
 		copycolor(vline[ysize], scanbar[ysize][i]);
@@ -254,7 +273,8 @@ char  *oldfile;
 
 	if (oldfile == NULL)
 		return(0);
-	else if ((fp = fopen(oldfile, "r")) == NULL) {
+	
+	if ((fp = fopen(oldfile, "r")) == NULL) {
 		sprintf(errmsg, "cannot open recover file \"%s\"", oldfile);
 		error(WARNING, errmsg);
 		return(0);
@@ -262,7 +282,7 @@ char  *oldfile;
 				/* discard header */
 	getheader(fp, NULL);
 				/* get picture size */
-	if (fscanf(fp, "-Y %d +X %d\n", &y, &x) != 2) {
+	if (fgetresolu(&x, &y, fp) != (YMAJOR|YDECR)) {
 		sprintf(errmsg, "bad recover file \"%s\"", oldfile);
 		error(WARNING, errmsg);
 		fclose(fp);
@@ -273,7 +293,6 @@ char  *oldfile;
 		sprintf(errmsg, "resolution mismatch in recover file \"%s\"",
 				oldfile);
 		error(USER, errmsg);
-		return(0);
 	}
 
 	scanline = (COLR *)malloc(ourview.hresolu*sizeof(COLR));