--- ray/src/rt/rpict.c 1990/01/08 13:38:01 1.12 +++ ray/src/rt/rpict.c 1990/11/01 12:06:01 1.23 @@ -15,7 +15,10 @@ static char SCCSid[] = "$SunId$ LBL"; #ifdef BSD #include #include +#else +#include #endif +#include #include "view.h" @@ -24,6 +27,7 @@ static char SCCSid[] = "$SunId$ LBL"; VIEW ourview = STDVIEW; /* view parameters */ int hresolu = 512; /* horizontal resolution */ int vresolu = 512; /* vertical resolution */ +double pixaspect = 1.0; /* pixel aspect ratio */ int psample = 4; /* pixel sample size */ double maxdiff = .05; /* max. difference for interpolation */ @@ -84,6 +88,7 @@ report() /* report progress */ sprintf(errmsg, "%ld rays, %4.2f%% done after %5.4f CPU hours\n", nrays, pctdone, t/3600.0); #else + signal(SIGALRM, report); sprintf(errmsg, "%ld rays, %4.2f%% done\n", nrays, pctdone); #endif eputs(errmsg); @@ -96,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; @@ -116,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); } @@ -126,7 +133,7 @@ char *zfile, *oldfile; goto memerr; } } else { - zfp = NULL; + zfd = -1; for (i = 0; i <= psample; i++) zbar[i] = NULL; } @@ -134,61 +141,54 @@ 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 && i > 0 && + 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; + if (close(zfd) == -1) + goto writerr; 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++) @@ -214,20 +214,17 @@ 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); - } } @@ -239,7 +236,6 @@ int xres, y, ysize; COLOR vline[MAXDIV+1]; float zline[MAXDIV+1]; int b = ysize; - double z; register int i, j; for (i = 0; i < xres; i++) { @@ -325,8 +321,11 @@ int x, y; /* pixel position */ { static RAY thisray; /* our ray for this pixel */ - viewray(thisray.rorg, thisray.rdir, &ourview, - (x+pixjitter())/hresolu, (y+pixjitter())/vresolu); + if (viewray(thisray.rorg, thisray.rdir, &ourview, + (x+pixjitter())/hresolu, (y+pixjitter())/vresolu) < 0) { + setcolor(col, 0.0, 0.0, 0.0); + return(0.0); + } rayorigin(&thisray, NULL, PRIMARY, 1.0); @@ -334,7 +333,7 @@ int x, y; /* pixel position */ copycolor(col, thisray.rcol); /* return color */ - return(thisray.rot); /* return distance */ + return(thisray.rt); /* return distance */ }