--- ray/src/rt/rtrace.c 1995/12/08 18:49:12 2.20 +++ ray/src/rt/rtrace.c 2004/03/30 16:13:01 2.37 @@ -1,15 +1,12 @@ -/* Copyright (c) 1995 Regents of the University of California */ - #ifndef lint -static char SCCSid[] = "$SunId$ LBL"; +static const char RCSid[] = "$Id: rtrace.c,v 2.37 2004/03/30 16:13:01 schorsch Exp $"; #endif - /* * rtrace.c - program and variables for individual ray tracing. - * - * 6/11/86 */ +#include "copyright.h" + /* * Input is in the form: * @@ -24,24 +21,33 @@ static char SCCSid[] = "$SunId$ LBL"; * irradiance values are desired. */ -#include "ray.h" +#include -#include "octree.h" - +#include "platform.h" +#include "ray.h" +#include "ambient.h" +#include "source.h" #include "otypes.h" - #include "resolu.h" +CUBE thescene; /* our scene */ +OBJECT nsceneobjs; /* number of objects in our scene */ + int dimlist[MAXDIM]; /* sampling dimensions */ int ndims = 0; /* number of sampling dimensions */ int samplendx = 0; /* index for this sample */ int imm_irrad = 0; /* compute immediate irradiance? */ +int lim_dist = 0; /* limit distance? */ int inform = 'a'; /* input format */ int outform = 'a'; /* output format */ char *outvals = "v"; /* output specification */ +int do_irrad = 0; /* compute irradiance? */ + +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 */ @@ -51,15 +57,15 @@ int hresolu = 0; /* horizontal (scan) size */ int vresolu = 0; /* vertical resolution */ double dstrsrc = 0.0; /* square source distribution */ -double shadthresh = .05; /* shadow threshold */ -double shadcert = .5; /* shadow certainty */ +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 */ -double salbedo = 0.; /* global scattering albedo */ +COLOR salbedo = BLKCOLOR; /* global scattering albedo */ double seccg = 0.; /* global scattering eccentricity */ double ssampdist = 0.; /* scatter sampling distance */ @@ -68,38 +74,53 @@ double specjitter = 1.; /* specular sampling jitter int backvis = 1; /* back face visibility */ -int maxdepth = 6; /* maximum recursion depth */ -double minweight = 4e-3; /* minimum ray weight */ +int maxdepth = 8; /* maximum recursion depth */ +double minweight = 2e-3; /* minimum ray weight */ +char *ambfile = NULL; /* ambient file name */ COLOR ambval = BLKCOLOR; /* ambient value */ -double ambacc = 0.2; /* ambient accuracy */ -int ambres = 128; /* ambient resolution */ -int ambdiv = 512; /* ambient divisions */ -int ambssamp = 0; /* ambient super-samples */ +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[128]; /* ambient include/exclude list */ int ambincl = -1; /* include == 1, exclude == 0 */ -extern OBJREC Lamb; /* a Lambertian surface */ +static int castonly = 0; static RAY thisray; /* for our convenience */ -static int oputo(), oputd(), oputv(), oputl(), oputL(), - oputp(), oputn(), oputN(), oputs(), oputw(), oputm(); +typedef void putf_t(double v); +static putf_t puta, putd, putf; -static int ourtrace(), tabin(); -static int (*ray_out[16])(), (*every_out[16])(); -static int castonly = 0; +typedef void oputf_t(RAY *r); +static oputf_t oputo, oputd, oputv, oputl, oputL, oputc, + oputp, oputn, oputN, oputs, oputw, oputm; -static int puta(), putf(), putd(); +static void setoutput(char *vs); +static 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 int (*putreal)(); +static oputf_t *ray_out[16], *every_out[16]; +static putf_t *putreal; +void (*addobjnotify[])() = {ambnotify, tranotify, NULL}; -quit(code) /* quit program */ -int code; + +void +quit( /* quit program */ + int code +) { -#ifndef NIX +#ifndef NON_POSIX /* XXX we don't clean up elsewhere? */ headclean(); /* delete header file */ pfclean(); /* clean up persist files */ #endif @@ -107,9 +128,10 @@ int code; } -char * -formstr(f) /* return format identifier */ -int f; +extern char * +formstr( /* return format identifier */ + int f +) { switch (f) { case 'a': return("ascii"); @@ -121,12 +143,15 @@ int f; } -rtrace(fname) /* trace rays from file */ -char *fname; +extern void +rtrace( /* trace rays from file */ + char *fname +) { long vcount = hresolu>1 ? hresolu*vresolu : vresolu; long nextflush = hresolu; FILE *fp; + double d; FVECT orig, direc; /* set up input */ if (fname == NULL) @@ -135,9 +160,9 @@ char *fname; sprintf(errmsg, "cannot open input file \"%s\"", fname); error(SYSTEM, errmsg); } -#ifdef MSDOS +#ifdef _WIN32 if (inform != 'a') - setmode(fileno(fp), O_BINARY); + SET_FILE_BINARY(fp); #endif /* set up output */ setoutput(outvals); @@ -161,20 +186,25 @@ char *fname; while (getvec(orig, inform, fp) == 0 && getvec(direc, inform, fp) == 0) { - if (normalize(direc) == 0.0) { /* zero ==> flush */ - fflush(stdout); - continue; - } - samplendx++; + d = normalize(direc); + if (d == 0.0) { /* zero ==> flush */ + bogusray(); + if (--nextflush <= 0 || vcount <= 0) { + fflush(stdout); + nextflush = hresolu; + } + } else { + samplendx++; /* compute and print */ - if (imm_irrad) - irrad(orig, direc); - else - rad(orig, direc); + if (imm_irrad) + irrad(orig, direc); + else + rad(orig, direc, lim_dist ? d : 0.0); /* flush if time */ - if (--nextflush == 0) { - fflush(stdout); - nextflush = hresolu; + if (--nextflush == 0) { + fflush(stdout); + nextflush = hresolu; + } } if (ferror(stdout)) error(SYSTEM, "write error"); @@ -189,11 +219,12 @@ char *fname; } -setoutput(vs) /* set up output tables */ -register char *vs; +static void +setoutput( /* set up output tables */ + register char *vs +) { - extern int (*trace)(); - register int (**table)() = ray_out; + register oputf_t **table = ray_out; castonly = 1; while (*vs) @@ -218,6 +249,9 @@ register char *vs; *table++ = oputl; castonly = 0; break; + case 'c': /* local coordinates */ + *table++ = oputc; + break; case 'L': /* single ray length */ *table++ = oputL; break; @@ -245,23 +279,46 @@ register char *vs; } -rad(org, dir) /* compute and print ray value(s) */ -FVECT org, dir; +static void +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); + printvals(&thisray); +} + + +static void +rad( /* compute and print ray value(s) */ + FVECT org, + FVECT dir, + double dmax +) +{ VCOPY(thisray.rorg, org); VCOPY(thisray.rdir, dir); - thisray.rmax = 0.0; + thisray.rmax = dmax; rayorigin(&thisray, NULL, PRIMARY, 1.0); - if (castonly) - localhit(&thisray, &thescene) || sourcehit(&thisray); - else + 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); } -irrad(org, dir) /* compute immediate irradiance value */ -FVECT org, dir; +static void +irrad( /* compute immediate irradiance value */ + FVECT org, + FVECT dir +) { register int i; @@ -282,10 +339,12 @@ FVECT org, dir; } -printvals(r) /* print requested ray values */ -RAY *r; +static void +printvals( /* print requested ray values */ + RAY *r +) { - register int (**tp)(); + register oputf_t **tp; if (ray_out[0] == NULL) return; @@ -296,12 +355,13 @@ RAY *r; } -getvec(vec, fmt, fp) /* get a vector from fp */ -register FVECT vec; -int fmt; -FILE *fp; +static int +getvec( /* get a vector from fp */ + register FVECT vec, + int fmt, + FILE *fp +) { - extern char *fgetword(); static float vf[3]; static double vd[3]; char buf[32]; @@ -333,13 +393,20 @@ FILE *fp; } -tranotify(obj) /* record new modifier */ -OBJECT obj; +static void +tranotify( /* record new modifier */ + OBJECT obj +) { static int hitlimit = 0; register OBJREC *o = objptr(obj); register char **tralp; + if (obj == OVOID) { /* starting over */ + traset[0] = 0; + hitlimit = 0; + return; + } if (hitlimit || !ismodifier(o->otype)) return; for (tralp = tralist; *tralp != NULL; tralp++) @@ -355,11 +422,12 @@ OBJECT obj; } -static -ourtrace(r) /* print ray values */ -RAY *r; +static void +ourtrace( /* print ray values */ + RAY *r +) { - register int (**tp)(); + register oputf_t **tp; if (every_out[0] == NULL) return; @@ -375,9 +443,10 @@ RAY *r; } -static -tabin(r) /* tab in appropriate amount */ -RAY *r; +static void +tabin( /* tab in appropriate amount */ + RAY *r +) { register RAY *rp; @@ -386,9 +455,10 @@ RAY *r; } -static -oputo(r) /* print origin */ -register RAY *r; +static void +oputo( /* print origin */ + RAY *r +) { (*putreal)(r->rorg[0]); (*putreal)(r->rorg[1]); @@ -396,9 +466,10 @@ register RAY *r; } -static -oputd(r) /* print direction */ -register RAY *r; +static void +oputd( /* print direction */ + RAY *r +) { (*putreal)(r->rdir[0]); (*putreal)(r->rdir[1]); @@ -406,9 +477,10 @@ register RAY *r; } -static -oputv(r) /* print value */ -register RAY *r; +static void +oputv( /* print value */ + RAY *r +) { COLR cout; @@ -425,26 +497,39 @@ register RAY *r; } -static -oputl(r) /* print effective distance */ -register RAY *r; +static void +oputl( /* print effective distance */ + RAY *r +) { (*putreal)(r->rt); } -static -oputL(r) /* print single ray length */ -register RAY *r; +static void +oputL( /* print single ray length */ + RAY *r +) { (*putreal)(r->rot); } -static -oputp(r) /* print point */ -register RAY *r; +static void +oputc( /* print local coordinates */ + RAY *r +) { + (*putreal)(r->uv[0]); + (*putreal)(r->uv[1]); +} + + +static void +oputp( /* print point */ + RAY *r +) +{ if (r->rot < FHUGE) { (*putreal)(r->rop[0]); (*putreal)(r->rop[1]); @@ -457,9 +542,10 @@ register RAY *r; } -static -oputN(r) /* print unperturbed normal */ -register RAY *r; +static void +oputN( /* print unperturbed normal */ + RAY *r +) { if (r->rot < FHUGE) { (*putreal)(r->ron[0]); @@ -473,9 +559,10 @@ register RAY *r; } -static -oputn(r) /* print perturbed normal */ -RAY *r; +static void +oputn( /* print perturbed normal */ + RAY *r +) { FVECT pnorm; @@ -492,9 +579,10 @@ RAY *r; } -static -oputs(r) /* print name */ -register RAY *r; +static void +oputs( /* print name */ + RAY *r +) { if (r->ro != NULL) fputs(r->ro->oname, stdout); @@ -504,35 +592,41 @@ register RAY *r; } -static -oputw(r) /* print weight */ -register RAY *r; +static void +oputw( /* print weight */ + RAY *r +) { (*putreal)(r->rweight); } -static -oputm(r) /* print modifier */ -register RAY *r; +static void +oputm( /* print modifier */ + RAY *r +) { if (r->ro != NULL) - fputs(objptr(r->ro->omod)->oname, stdout); + if (r->ro->omod != OVOID) + fputs(objptr(r->ro->omod)->oname, stdout); + else + fputs(VOIDID, stdout); else putchar('*'); putchar('\t'); } -static -puta(v) /* print ascii value */ -double v; +static void +puta( /* print ascii value */ + double v +) { printf("%e\t", v); } -static +static void putd(v) /* print binary double */ double v; { @@ -540,7 +634,7 @@ double v; } -static +static void putf(v) /* print binary float */ double v; {