--- ray/src/rt/rpict.c	1990/01/11 08:30:42	1.14
+++ ray/src/rt/rpict.c	1990/03/01 08:19:01	1.19
@@ -18,6 +18,7 @@ static char SCCSid[] = "$SunId$ LBL";
 #else
 #include  <signal.h>
 #endif
+#include  <fcntl.h>
 
 #include  "view.h"
 
@@ -100,10 +101,12 @@ report()		/* report progress */
 render(zfile, oldfile)				/* render the scene */
 char  *zfile, *oldfile;
 {
+	extern long  lseek();
 	COLOR  *scanbar[MAXDIV+1];	/* scanline arrays of pixel values */
 	float  *zbar[MAXDIV+1];		/* z values */
 	int  ypos;			/* current scanline */
-	FILE  *zfp;
+	int  ystep;			/* current y step size */
+	int  zfd;
 	COLOR  *colptr;
 	float  *zptr;
 	register int  i;
@@ -120,7 +123,7 @@ char  *zfile, *oldfile;
 	}
 					/* open z file */
 	if (zfile != NULL) {
-		if ((zfp = fopen(zfile, "a+")) == NULL) {
+		if ((zfd = open(zfile, O_WRONLY|O_CREAT, 0666)) == -1) {
 			sprintf(errmsg, "cannot open z file \"%s\"", zfile);
 			error(SYSTEM, errmsg);
 		}
@@ -130,7 +133,7 @@ char  *zfile, *oldfile;
 				goto memerr;
 		}
 	} else {
-		zfp = NULL;
+		zfd = -1;
 		for (i = 0; i <= psample; i++)
 			zbar[i] = NULL;
 	}
@@ -138,61 +141,52 @@ char  *zfile, *oldfile;
 	fputresolu(YMAJOR|YDECR, hresolu, vresolu, stdout);
 					/* recover file and compute first */
 	i = salvage(oldfile);
-	if (zfp != NULL && fseek(zfp, (long)i*hresolu*sizeof(float), 0) == EOF)
+	if (zfd != -1 && lseek(zfd, (long)i*hresolu*sizeof(float), 0) == -1)
 		error(SYSTEM, "z file seek error in render");
 	ypos = vresolu-1 - i;
 	fillscanline(scanbar[0], zbar[0], hresolu, ypos, psample);
+	ystep = psample;
 						/* compute scanlines */
-	for (ypos -= psample; ypos >= 0; ypos -= psample) {
-	
-		pctdone = 100.0*(vresolu-ypos-psample)/vresolu;
-
-		colptr = scanbar[psample];		/* move base to top */
-		scanbar[psample] = scanbar[0];
+	for (ypos -= ystep; ypos > -ystep; ypos -= ystep) {
+							/* record progress */
+		pctdone = 100.0*(vresolu-1-ypos-ystep)/vresolu;
+							/* bottom adjust? */
+		if (ypos < 0) {
+			ystep += ypos;
+			ypos = 0;
+		}
+		colptr = scanbar[ystep];		/* move base to top */
+		scanbar[ystep] = scanbar[0];
 		scanbar[0] = colptr;
-		zptr = zbar[psample];
-		zbar[psample] = zbar[0];
+		zptr = zbar[ystep];
+		zbar[ystep] = zbar[0];
 		zbar[0] = zptr;
 							/* fill base line */
 		fillscanline(scanbar[0], zbar[0], hresolu, ypos, psample);
 							/* fill bar */
-		fillscanbar(scanbar, zbar, hresolu, ypos, psample);
+		fillscanbar(scanbar, zbar, hresolu, ypos, ystep);
 							/* write it out */
-		for (i = psample; i > 0; i--) {
-			if (zfp != NULL && fwrite(zbar[i],sizeof(float),hresolu,zfp) != hresolu)
+		for (i = ystep; i > 0; i--) {
+			if (zfd != -1 && write(zfd, (char *)zbar[i],
+					hresolu*sizeof(float))
+					< hresolu*sizeof(float))
 				goto writerr;
-			if (fwritescan(scanbar[i],hresolu,stdout) < 0)
+			if (fwritescan(scanbar[i], hresolu, stdout) < 0)
 				goto writerr;
 		}
-		if (zfp != NULL && fflush(zfp) == EOF)
-			goto writerr;
 		if (fflush(stdout) == EOF)
 			goto writerr;
 	}
-						/* compute residual */
-	colptr = scanbar[psample];
-	scanbar[psample] = scanbar[0];
-	scanbar[0] = colptr;
-	zptr = zbar[psample];
-	zbar[psample] = zbar[0];
-	zbar[0] = zptr;
-	if (ypos > -psample) {
-		fillscanline(scanbar[-ypos], zbar[-ypos], hresolu, 0, psample);
-		fillscanbar(scanbar-ypos, zbar-ypos, hresolu, 0, psample+ypos);
-	}
-	for (i = psample; i+ypos >= 0; i--) {
-		if (zfp != NULL && fwrite(zbar[i],sizeof(float),hresolu,zfp) != hresolu)
-			goto writerr;
-		if (fwritescan(scanbar[i], hresolu, stdout) < 0)
-			goto writerr;
-	}
 						/* clean up */
-	if (zfp != NULL) {
-		if (fclose(zfp) == EOF)
+	if (zfd != -1) {
+		if (write(zfd, (char *)zbar[0], hresolu*sizeof(float))
+				< hresolu*sizeof(float))
 			goto writerr;
+		close(zfd);
 		for (i = 0; i <= psample; i++)
 			free((char *)zbar[i]);
 	}
+	fwritescan(scanbar[0], hresolu, stdout);
 	if (fflush(stdout) == EOF)
 		goto writerr;
 	for (i = 0; i <= psample; i++)
@@ -218,19 +212,16 @@ int  xres, y, xstep;
 	z = pixvalue(scanline[0], 0, y);
 	if (zline) zline[0] = z;
 
-	for (i = xstep; i < xres; i += xstep) {
-	
+	for (i = xstep; i < xres-1+xstep; i += xstep) {
+		if (i >= xres) {
+			xstep += xres-1-i;
+			i = xres-1;
+		}
 		z = pixvalue(scanline[i], i, y);
 		if (zline) zline[i] = z;
 		
 		b = fillsample(scanline+i-xstep, zline ? zline+i-xstep : NULL,
 				i-xstep, y, xstep, 0, b/2);
-	}
-	if (i-xstep < xres-1) {
-		z = pixvalue(scanline[xres-1], xres-1, y);
-		if (zline) zline[xres-1] = z;
-		fillsample(scanline+i-xstep, zline ? zline+i-xstep : NULL,
-				i-xstep, y, xres-1-(i-xstep), 0, b/2);
 	}
 }