--- ray/src/rt/rtrace.c 2004/11/05 17:36:56 2.38 +++ ray/src/rt/rtrace.c 2009/12/14 07:31:37 2.62 @@ -1,5 +1,5 @@ #ifndef lint -static const char RCSid[] = "$Id: rtrace.c,v 2.38 2004/11/05 17:36:56 greg Exp $"; +static const char RCSid[] = "$Id: rtrace.c,v 2.62 2009/12/14 07:31:37 greg Exp $"; #endif /* * rtrace.c - program and variables for individual ray tracing. @@ -29,106 +29,70 @@ static const char RCSid[] = "$Id: rtrace.c,v 2.38 2004 #include "source.h" #include "otypes.h" #include "resolu.h" +#include "random.h" -CUBE thescene; /* our scene */ -OBJECT nsceneobjs; /* number of objects in our scene */ +extern int inform; /* input format */ +extern int outform; /* output format */ +extern char *outvals; /* output values */ -int dimlist[MAXDIM]; /* sampling dimensions */ -int ndims = 0; /* number of sampling dimensions */ -int samplendx = 0; /* index for this sample */ +extern int imm_irrad; /* compute immediate irradiance? */ +extern int lim_dist; /* limit distance? */ -int imm_irrad = 0; /* compute immediate irradiance? */ -int lim_dist = 0; /* limit distance? */ +extern char *tralist[]; /* list of modifers to trace (or no) */ +extern int traincl; /* include == 1, exclude == 0 */ -int inform = 'a'; /* input format */ -int outform = 'a'; /* output format */ -char *outvals = "v"; /* output specification */ +extern int hresolu; /* horizontal resolution */ +extern int vresolu; /* vertical resolution */ -int do_irrad = 0; /* compute irradiance? */ +static int castonly = 0; -void (*trace)() = NULL; /* trace call */ - -char *tralist[128]; /* list of modifers to trace (or no) */ -int traincl = -1; /* include == 1, exclude == 0 */ -#define MAXTSET 511 /* maximum number in trace set */ +#ifndef MAXTSET +#define MAXTSET 8191 /* maximum number in trace set */ +#endif OBJECT traset[MAXTSET+1]={0}; /* trace include/exclude set */ -int hresolu = 0; /* horizontal (scan) size */ -int vresolu = 0; /* vertical resolution */ - -double dstrsrc = 0.0; /* square source distribution */ -double shadthresh = .03; /* shadow threshold */ -double shadcert = .75; /* shadow certainty */ -int directrelay = 2; /* number of source relays */ -int vspretest = 512; /* virtual source pretest density */ -int directvis = 1; /* sources visible? */ -double srcsizerat = .2; /* maximum ratio source size/dist. */ - -COLOR cextinction = BLKCOLOR; /* global extinction coefficient */ -COLOR salbedo = BLKCOLOR; /* global scattering albedo */ -double seccg = 0.; /* global scattering eccentricity */ -double ssampdist = 0.; /* scatter sampling distance */ - -double specthresh = .15; /* specular sampling threshold */ -double specjitter = 1.; /* specular sampling jitter */ - -int backvis = 1; /* back face visibility */ - -int maxdepth = 8; /* maximum recursion depth */ -double minweight = 2e-3; /* minimum ray weight */ - -char *ambfile = NULL; /* ambient file name */ -COLOR ambval = BLKCOLOR; /* ambient value */ -int ambvwt = 0; /* initial weight for ambient value */ -double ambacc = 0.15; /* ambient accuracy */ -int ambres = 256; /* ambient resolution */ -int ambdiv = 1024; /* ambient divisions */ -int ambssamp = 512; /* ambient super-samples */ -int ambounce = 0; /* ambient bounces */ -char *amblist[AMBLLEN]; /* ambient include/exclude list */ -int ambincl = -1; /* include == 1, exclude == 0 */ - -static int castonly = 0; - static RAY thisray; /* for our convenience */ typedef void putf_t(double v); static putf_t puta, putd, putf; typedef void oputf_t(RAY *r); -static oputf_t oputo, oputd, oputv, oputl, oputL, oputc, - oputp, oputn, oputN, oputs, oputw, oputm; +static oputf_t oputo, oputd, oputv, oputV, oputl, oputL, oputc, oputp, + oputn, oputN, oputs, oputw, oputW, oputm, oputM, oputtilde; static void setoutput(char *vs); -static void tranotify(OBJECT obj); +extern void tranotify(OBJECT obj); static void bogusray(void); -static void rad(FVECT org, FVECT dir, double dmax); -static void irrad(FVECT org, FVECT dir); -static void printvals(RAY *r); -static int getvec(FVECT vec, int fmt, FILE *fp); -static void tabin(RAY *r); -static void ourtrace(RAY *r); +static void raycast(RAY *r); +static void rayirrad(RAY *r); +static void rtcompute(FVECT org, FVECT dir, double dmax); +static int printvals(RAY *r); +static int getvec(FVECT vec, int fmt, FILE *fp); +static void tabin(RAY *r); +static void ourtrace(RAY *r); static oputf_t *ray_out[16], *every_out[16]; static putf_t *putreal; -void (*addobjnotify[])() = {ambnotify, tranotify, NULL}; - void quit( /* quit program */ int code ) { -#ifndef NON_POSIX /* XXX we don't clean up elsewhere? */ - headclean(); /* delete header file */ - pfclean(); /* clean up persist files */ + if (ray_pnprocs > 0) /* close children if any */ + ray_pclose(0); +#ifndef NON_POSIX + else if (!ray_pnprocs) { + headclean(); /* delete header file */ + pfclean(); /* clean up persist files */ + } #endif exit(code); } -extern char * +char * formstr( /* return format identifier */ int f ) @@ -145,10 +109,12 @@ formstr( /* return format identifier */ extern void rtrace( /* trace rays from file */ - char *fname + char *fname, + int nproc ) { - long vcount = hresolu>1 ? hresolu*vresolu : vresolu; + unsigned long vcount = (hresolu > 1) ? (unsigned long)hresolu*vresolu + : vresolu; long nextflush = hresolu; FILE *fp; double d; @@ -160,12 +126,14 @@ rtrace( /* trace rays from file */ sprintf(errmsg, "cannot open input file \"%s\"", fname); error(SYSTEM, errmsg); } -#ifdef _WIN32 if (inform != 'a') SET_FILE_BINARY(fp); -#endif /* set up output */ setoutput(outvals); + if (imm_irrad) + castonly = 0; + else if (castonly) + nproc = 1; /* don't bother multiprocessing */ switch (outform) { case 'a': putreal = puta; break; case 'f': putreal = putf; break; @@ -177,6 +145,10 @@ rtrace( /* trace rays from file */ default: error(CONSISTENCY, "botched output format"); } + if (nproc > 1) { /* start multiprocessing */ + ray_popen(nproc); + ray_fifo_out = printvals; + } if (hresolu > 0) { if (vresolu > 0) fprtresolu(hresolu, vresolu, stdout); @@ -188,48 +160,68 @@ rtrace( /* trace rays from file */ d = normalize(direc); if (d == 0.0) { /* zero ==> flush */ + if (nproc > 1 && ray_fifo_flush() < 0) + error(USER, "lost children"); bogusray(); - if (--nextflush <= 0 || vcount <= 0) { + if (--nextflush <= 0 || !vcount) { fflush(stdout); nextflush = hresolu; } - } else { - samplendx++; - /* compute and print */ - if (imm_irrad) - irrad(orig, direc); - else - rad(orig, direc, lim_dist ? d : 0.0); + } else { /* compute and print */ + rtcompute(orig, direc, lim_dist ? d : 0.0); /* flush if time */ - if (--nextflush == 0) { + if (!--nextflush) { + if (nproc > 1 && ray_fifo_flush() < 0) + error(USER, "lost children"); fflush(stdout); nextflush = hresolu; } } if (ferror(stdout)) error(SYSTEM, "write error"); - if (--vcount == 0) /* check for end */ + if (vcount && !--vcount) /* check for end */ break; } - fflush(stdout); - if (vcount > 0) - error(USER, "read error"); + if (nproc > 1) { /* clean up children */ + if (ray_fifo_flush() < 0) + error(USER, "unable to complete processing"); + ray_pclose(0); + } + if (fflush(stdout) < 0) + error(SYSTEM, "write error"); + if (vcount) + error(USER, "unexpected EOF on input"); if (fname != NULL) fclose(fp); } static void +trace_sources(void) /* trace rays to light sources, also */ +{ + int sn; + + for (sn = 0; sn < nsources; sn++) + source[sn].sflags |= SFOLLOW; +} + + +static void setoutput( /* set up output tables */ - register char *vs + char *vs ) { - register oputf_t **table = ray_out; + oputf_t **table = ray_out; castonly = 1; while (*vs) switch (*vs++) { + case 'T': /* trace sources */ + if (!*vs) break; + trace_sources(); + /* fall through */ case 't': /* trace */ + if (!*vs) break; *table = NULL; table = every_out; trace = ourtrace; @@ -245,6 +237,12 @@ setoutput( /* set up output tables */ *table++ = oputv; castonly = 0; break; + case 'V': /* contribution */ + *table++ = oputV; + if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0)) + error(WARNING, + "-otV accuracy depends on -aa 0 -as 0"); + break; case 'l': /* effective distance */ *table++ = oputl; castonly = 0; @@ -271,9 +269,21 @@ setoutput( /* set up output tables */ case 'w': /* weight */ *table++ = oputw; break; + case 'W': /* coefficient */ + *table++ = oputW; + if (ambounce > 0 && (ambacc > FTINY || ambssamp > 0)) + error(WARNING, + "-otW accuracy depends on -aa 0 -as 0"); + break; case 'm': /* modifier */ *table++ = oputm; break; + case 'M': /* material */ + *table++ = oputM; + break; + case '~': /* tilde */ + *table++ = oputtilde; + break; } *table = NULL; } @@ -284,80 +294,101 @@ bogusray(void) /* print out empty record */ { thisray.rorg[0] = thisray.rorg[1] = thisray.rorg[2] = thisray.rdir[0] = thisray.rdir[1] = thisray.rdir[2] = 0.0; - rayorigin(&thisray, NULL, PRIMARY, 1.0); + thisray.rmax = 0.0; + rayorigin(&thisray, PRIMARY, NULL, NULL); printvals(&thisray); } static void -rad( /* compute and print ray value(s) */ - FVECT org, - FVECT dir, - double dmax +raycast( /* compute first ray intersection only */ + RAY *r ) { - VCOPY(thisray.rorg, org); - VCOPY(thisray.rdir, dir); - thisray.rmax = dmax; - rayorigin(&thisray, NULL, PRIMARY, 1.0); - if (castonly) { - if (!localhit(&thisray, &thescene)) { - if (thisray.ro == &Aftplane) { /* clipped */ - thisray.ro = NULL; - thisray.rot = FHUGE; - } else - sourcehit(&thisray); - } - } else - rayvalue(&thisray); - printvals(&thisray); + if (!localhit(r, &thescene)) { + if (r->ro == &Aftplane) { /* clipped */ + r->ro = NULL; + r->rot = FHUGE; + } else + sourcehit(r); + } } static void -irrad( /* compute immediate irradiance value */ - FVECT org, - FVECT dir +rayirrad( /* compute irradiance rather than radiance */ + RAY *r ) { - register int i; + void (*old_revf)(RAY *) = r->revf; - for (i = 0; i < 3; i++) { - thisray.rorg[i] = org[i] + dir[i]; - thisray.rdir[i] = -dir[i]; + r->rot = 1e-5; /* pretend we hit surface */ + VSUM(r->rop, r->rorg, r->rdir, r->rot); + r->ron[0] = -r->rdir[0]; + r->ron[1] = -r->rdir[1]; + r->ron[2] = -r->rdir[2]; + r->rod = 1.0; + /* compute result */ + r->revf = raytrace; + (*ofun[Lamb.otype].funp)(&Lamb, r); + r->revf = old_revf; +} + + +static void +rtcompute( /* compute and print ray value(s) */ + FVECT org, + FVECT dir, + double dmax +) +{ + /* set up ray */ + rayorigin(&thisray, PRIMARY, NULL, NULL); + if (imm_irrad) { + VSUM(thisray.rorg, org, dir, 1.1e-4); + thisray.rdir[0] = -dir[0]; + thisray.rdir[1] = -dir[1]; + thisray.rdir[2] = -dir[2]; + thisray.rmax = 0.0; + thisray.revf = rayirrad; + } else { + VCOPY(thisray.rorg, org); + VCOPY(thisray.rdir, dir); + thisray.rmax = dmax; + if (castonly) + thisray.revf = raycast; } - rayorigin(&thisray, NULL, PRIMARY, 1.0); - /* pretend we hit surface */ - thisray.rot = 1.0-1e-4; - thisray.rod = 1.0; - VCOPY(thisray.ron, dir); - for (i = 0; i < 3; i++) /* fudge factor */ - thisray.rop[i] = org[i] + 1e-4*dir[i]; - /* compute and print */ - (*ofun[Lamb.otype].funp)(&Lamb, &thisray); + if (ray_pnprocs > 1) { /* multiprocessing FIFO? */ + if (ray_fifo_in(&thisray) < 0) + error(USER, "lost children"); + return; + } + samplendx++; /* else do it ourselves */ + rayvalue(&thisray); printvals(&thisray); } -static void +static int printvals( /* print requested ray values */ RAY *r ) { - register oputf_t **tp; + oputf_t **tp; if (ray_out[0] == NULL) - return; + return(0); for (tp = ray_out; *tp != NULL; tp++) (**tp)(r); if (outform == 'a') putchar('\n'); + return(1); } static int getvec( /* get a vector from fp */ - register FVECT vec, + FVECT vec, int fmt, FILE *fp ) @@ -365,7 +396,7 @@ getvec( /* get a vector from fp */ static float vf[3]; static double vd[3]; char buf[32]; - register int i; + int i; switch (fmt) { case 'a': /* ascii */ @@ -379,12 +410,12 @@ getvec( /* get a vector from fp */ case 'f': /* binary float */ if (fread((char *)vf, sizeof(float), 3, fp) != 3) return(-1); - vec[0] = vf[0]; vec[1] = vf[1]; vec[2] = vf[2]; + VCOPY(vec, vf); break; case 'd': /* binary double */ if (fread((char *)vd, sizeof(double), 3, fp) != 3) return(-1); - vec[0] = vd[0]; vec[1] = vd[1]; vec[2] = vd[2]; + VCOPY(vec, vd); break; default: error(CONSISTENCY, "botched input format"); @@ -393,14 +424,14 @@ getvec( /* get a vector from fp */ } -static void +void tranotify( /* record new modifier */ OBJECT obj ) { static int hitlimit = 0; - register OBJREC *o = objptr(obj); - register char **tralp; + OBJREC *o = objptr(obj); + char **tralp; if (obj == OVOID) { /* starting over */ traset[0] = 0; @@ -427,7 +458,7 @@ ourtrace( /* print ray values */ RAY *r ) { - register oputf_t **tp; + oputf_t **tp; if (every_out[0] == NULL) return; @@ -439,7 +470,8 @@ ourtrace( /* print ray values */ tabin(r); for (tp = every_out; *tp != NULL; tp++) (**tp)(r); - putchar('\n'); + if (outform == 'a') + putchar('\n'); } @@ -448,7 +480,7 @@ tabin( /* tab in appropriate amount */ RAY *r ) { - register RAY *rp; + const RAY *rp; for (rp = r->parent; rp != NULL; rp = rp->parent) putchar('\t'); @@ -482,9 +514,8 @@ oputv( /* print value */ RAY *r ) { - COLR cout; - if (outform == 'c') { + COLR cout; setcolr(cout, colval(r->rcol,RED), colval(r->rcol,GRN), colval(r->rcol,BLU)); @@ -498,6 +529,21 @@ oputv( /* print value */ static void +oputV( /* print value contribution */ + RAY *r +) +{ + double contr[3]; + + raycontrib(contr, r, PRIMARY); + multcolor(contr, r->rcol); + (*putreal)(contr[RED]); + (*putreal)(contr[GRN]); + (*putreal)(contr[BLU]); +} + + +static void oputl( /* print effective distance */ RAY *r ) @@ -602,6 +648,20 @@ oputw( /* print weight */ static void +oputW( /* print coefficient */ + RAY *r +) +{ + double contr[3]; + + raycontrib(contr, r, PRIMARY); + (*putreal)(contr[RED]); + (*putreal)(contr[GRN]); + (*putreal)(contr[BLU]); +} + + +static void oputm( /* print modifier */ RAY *r ) @@ -614,6 +674,33 @@ oputm( /* print modifier */ else putchar('*'); putchar('\t'); +} + + +static void +oputM( /* print material */ + RAY *r +) +{ + OBJREC *mat; + + if (r->ro != NULL) { + if ((mat = findmaterial(r->ro)) != NULL) + fputs(mat->oname, stdout); + else + fputs(VOIDID, stdout); + } else + putchar('*'); + putchar('\t'); +} + + +static void +oputtilde( /* output tilde (spacer) */ + RAY *r +) +{ + fputs("~\t", stdout); }