--- ray/src/rt/rpict.c 1991/06/25 14:06:12 1.29 +++ ray/src/rt/rpict.c 1991/12/10 17:42:01 2.3 @@ -15,13 +15,15 @@ static char SCCSid[] = "$SunId$ LBL"; #ifdef BSD #include #include -#else -#include #endif + +#include #include #include "view.h" +#include "resolu.h" + #include "random.h" int dimlist[MAXDIM]; /* sampling dimensions */ @@ -40,8 +42,10 @@ 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 directrelay = 0; /* number of source relays */ +int directrelay = 1; /* number of source relays */ int vspretest = 512; /* virtual source pretest density */ +int directinvis = 0; /* sources invisible? */ +double srcsizerat = .25; /* maximum ratio source size/dist. */ int maxdepth = 6; /* maximum recursion depth */ double minweight = 5e-3; /* minimum ray weight */ @@ -66,7 +70,7 @@ extern long tstart; /* starting time */ extern long nrays; /* number of rays traced */ -#define MAXDIV 15 /* maximum sample size */ +#define MAXDIV 16 /* maximum sample size */ #define pixjitter() (.5+dstrpix*(.5-frandom())) @@ -119,8 +123,10 @@ char *zfile, *oldfile; extern long lseek(); COLOR *scanbar[MAXDIV+1]; /* scanline arrays of pixel values */ float *zbar[MAXDIV+1]; /* z values */ + char *sampdens; /* previous sample density */ int ypos; /* current scanline */ int ystep; /* current y step size */ + int hstep; /* h step size */ int zfd; COLOR *colptr; float *zptr; @@ -140,6 +146,16 @@ char *zfile, *oldfile; if (scanbar[i] == NULL) goto memerr; } + hstep = (psample*140+49)/99; /* quincunx sampling */ + ystep = (psample*99+70)/140; + if (hstep > 2) { + i = hresolu/hstep + 2; + if ((sampdens = malloc(i)) == NULL) + goto memerr; + while (i--) + sampdens[i] = hstep; + } else + sampdens = NULL; /* open z file */ if (zfile != NULL) { if ((zfd = open(zfile, O_WRONLY|O_CREAT, 0666)) == -1) { @@ -157,7 +173,7 @@ char *zfile, *oldfile; zbar[i] = NULL; } /* write out boundaries */ - fputresolu(YMAJOR|YDECR, hresolu, vresolu, stdout); + fprtresolu(hresolu, vresolu, stdout); /* recover file and compute first */ i = salvage(oldfile); if (zfd != -1 && i > 0 && @@ -166,9 +182,12 @@ char *zfile, *oldfile; pctdone = 100.0*i/vresolu; if (ralrm > 0) /* report init stats */ report(); +#ifndef BSD + else +#endif + signal(SIGALRM, report); ypos = vresolu-1 - i; - fillscanline(scanbar[0], zbar[0], hresolu, ypos, psample); - ystep = psample; + fillscanline(scanbar[0], zbar[0], sampdens, hresolu, ypos, hstep); /* compute scanlines */ for (ypos -= ystep; ypos > -ystep; ypos -= ystep) { /* bottom adjust? */ @@ -183,10 +202,14 @@ char *zfile, *oldfile; zbar[ystep] = zbar[0]; zbar[0] = zptr; /* fill base line */ - fillscanline(scanbar[0], zbar[0], hresolu, ypos, psample); + fillscanline(scanbar[0], zbar[0], sampdens, + hresolu, ypos, hstep); /* fill bar */ fillscanbar(scanbar, zbar, hresolu, ypos, ystep); /* write it out */ +#ifndef BSD + signal(SIGALRM, SIG_IGN); /* don't interrupt writes */ +#endif for (i = ystep; i > 0; i--) { if (zfd != -1 && write(zfd, (char *)zbar[i], hresolu*sizeof(float)) @@ -201,6 +224,10 @@ char *zfile, *oldfile; pctdone = 100.0*(vresolu-1-ypos)/vresolu; if (ralrm > 0 && time((long *)0) >= tlastrept+ralrm) report(); +#ifndef BSD + else + signal(SIGALRM, report); +#endif } /* clean up */ if (zfd != -1) { @@ -217,6 +244,8 @@ char *zfile, *oldfile; goto writerr; for (i = 0; i <= psample; i++) free((char *)scanbar[i]); + if (sampdens != NULL) + free(sampdens); pctdone = 100.0; return; writerr: @@ -226,29 +255,38 @@ memerr: } -fillscanline(scanline, zline, xres, y, xstep) /* fill scan line at y */ +fillscanline(scanline, zline, sd, xres, y, xstep) /* fill scan at y */ register COLOR *scanline; register float *zline; +register char *sd; int xres, y, xstep; { - int b = xstep; + static int nc = 0; /* number of calls */ + int bl = xstep, b = xstep; double z; register int i; z = pixvalue(scanline[0], 0, y); if (zline) zline[0] = z; - - for (i = xstep; i < xres-1+xstep; i += xstep) { + /* zig-zag start for quincunx pattern */ + for (i = ++nc & 1 ? xstep : xstep/2; 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 (sd) b = sd[0] > sd[1] ? sd[0] : sd[1]; + if (i <= xstep) + b = fillsample(scanline, zline, 0, y, i, 0, b/2); + else + b = fillsample(scanline+i-xstep, + zline ? zline+i-xstep : NULL, + i-xstep, y, xstep, 0, b/2); + if (sd) *sd++ = nc & 1 ? bl : b; + bl = b; } + if (sd && nc & 1) *sd = bl; } @@ -353,7 +391,7 @@ int x, y; /* pixel position */ rayorigin(&thisray, NULL, PRIMARY, 1.0); - samplendx = 3*y + x; /* set pixel index */ + samplendx = pixnumber(x,y,hresolu,vresolu); /* set pixel index */ rayvalue(&thisray); /* trace ray */ @@ -382,7 +420,7 @@ char *oldfile; /* discard header */ getheader(fp, NULL); /* get picture size */ - if (fgetresolu(&x, &y, fp) != (YMAJOR|YDECR)) { + if (!fscnresolu(&x, &y, fp)) { sprintf(errmsg, "bad recover file \"%s\"", oldfile); error(WARNING, errmsg); fclose(fp); @@ -412,4 +450,16 @@ char *oldfile; return(y); writerr: error(SYSTEM, "write error in salvage"); +} + + +int +pixnumber(x, y, xres, yres) /* compute pixel index (brushed) */ +register int x, y; +int xres, yres; +{ + x -= y; + while (x < 0) + x += xres; + return((((x>>2)*yres + y) << 2) + (x & 3)); }