--- ray/src/rt/rtrace.c 1997/10/17 10:14:37 2.25 +++ ray/src/rt/rtrace.c 2003/06/05 19:29:34 2.30 @@ -1,15 +1,12 @@ -/* Copyright (c) 1997 Silicon Graphics, Inc. */ - #ifndef lint -static char SCCSid[] = "$SunId$ SGI"; +static const char RCSid[] = "$Id: rtrace.c,v 2.30 2003/06/05 19:29:34 schorsch Exp $"; #endif - /* * rtrace.c - program and variables for individual ray tracing. - * - * 6/11/86 */ +#include "copyright.h" + /* * Input is in the form: * @@ -24,14 +21,14 @@ static char SCCSid[] = "$SunId$ SGI"; * irradiance values are desired. */ +#include "platform.h" #include "ray.h" - -#include "octree.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 */ @@ -43,6 +40,12 @@ 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 */ + +extern void ambnotify(), tranotify(); +void (*addobjnotify[])() = {ambnotify, tranotify, NULL}; 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 */ @@ -72,6 +75,7 @@ int backvis = 1; /* back face visibility */ int maxdepth = 6; /* maximum recursion depth */ double minweight = 4e-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.2; /* ambient accuracy */ @@ -82,24 +86,24 @@ 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 */ -extern OBJREC Aftplane; /* aft clipping object */ - static RAY thisray; /* for our convenience */ -static int oputo(), oputd(), oputv(), oputl(), oputL(), +static void oputo(), oputd(), oputv(), oputl(), oputL(), oputc(), oputp(), oputn(), oputN(), oputs(), oputw(), oputm(); -static int ourtrace(), tabin(); -static int (*ray_out[16])(), (*every_out[16])(); +static void ourtrace(), tabin(); +static void (*ray_out[16])(), (*every_out[16])(); static int castonly = 0; -static int puta(), putf(), putd(); +static void puta(), putf(), putd(); -static int (*putreal)(); +static void (*putreal)(); +void bogusray(), rad(), irrad(), printvals(); + +void quit(code) /* quit program */ int code; { @@ -125,6 +129,7 @@ int f; } +void rtrace(fname) /* trace rays from file */ char *fname; { @@ -140,9 +145,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); @@ -169,7 +174,7 @@ char *fname; d = normalize(direc); if (d == 0.0) { /* zero ==> flush */ bogusray(); - if (--nextflush <= 0) { + if (--nextflush <= 0 || vcount <= 0) { fflush(stdout); nextflush = hresolu; } @@ -202,8 +207,8 @@ char *fname; setoutput(vs) /* set up output tables */ register char *vs; { - extern int (*trace)(); - register int (**table)() = ray_out; + extern void (*trace)(); + register void (**table)() = ray_out; castonly = 1; while (*vs) @@ -228,6 +233,9 @@ register char *vs; *table++ = oputl; castonly = 0; break; + case 'c': /* local coordinates */ + *table++ = oputc; + break; case 'L': /* single ray length */ *table++ = oputL; break; @@ -255,6 +263,7 @@ register char *vs; } +void bogusray() /* print out empty record */ { thisray.rorg[0] = thisray.rorg[1] = thisray.rorg[2] = @@ -264,6 +273,7 @@ bogusray() /* print out empty record */ } +void rad(org, dir, dmax) /* compute and print ray value(s) */ FVECT org, dir; double dmax; @@ -285,6 +295,7 @@ double dmax; } +void irrad(org, dir) /* compute immediate irradiance value */ FVECT org, dir; { @@ -307,10 +318,11 @@ FVECT org, dir; } +void printvals(r) /* print requested ray values */ RAY *r; { - register int (**tp)(); + register void (**tp)(); if (ray_out[0] == NULL) return; @@ -321,12 +333,12 @@ RAY *r; } +int getvec(vec, fmt, fp) /* 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]; @@ -358,6 +370,7 @@ FILE *fp; } +void tranotify(obj) /* record new modifier */ OBJECT obj; { @@ -365,6 +378,11 @@ OBJECT obj; 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++) @@ -380,11 +398,11 @@ OBJECT obj; } -static +static void ourtrace(r) /* print ray values */ RAY *r; { - register int (**tp)(); + register void (**tp)(); if (every_out[0] == NULL) return; @@ -400,7 +418,7 @@ RAY *r; } -static +static void tabin(r) /* tab in appropriate amount */ RAY *r; { @@ -411,9 +429,9 @@ RAY *r; } -static +static void oputo(r) /* print origin */ -register RAY *r; +RAY *r; { (*putreal)(r->rorg[0]); (*putreal)(r->rorg[1]); @@ -421,9 +439,9 @@ register RAY *r; } -static +static void oputd(r) /* print direction */ -register RAY *r; +RAY *r; { (*putreal)(r->rdir[0]); (*putreal)(r->rdir[1]); @@ -431,9 +449,9 @@ register RAY *r; } -static +static void oputv(r) /* print value */ -register RAY *r; +RAY *r; { COLR cout; @@ -450,25 +468,34 @@ register RAY *r; } -static +static void oputl(r) /* print effective distance */ -register RAY *r; +RAY *r; { (*putreal)(r->rt); } -static +static void oputL(r) /* print single ray length */ -register RAY *r; +RAY *r; { (*putreal)(r->rot); } -static +static void +oputc(r) /* print local coordinates */ +RAY *r; +{ + (*putreal)(r->uv[0]); + (*putreal)(r->uv[1]); +} + + +static void oputp(r) /* print point */ -register RAY *r; +RAY *r; { if (r->rot < FHUGE) { (*putreal)(r->rop[0]); @@ -482,9 +509,9 @@ register RAY *r; } -static +static void oputN(r) /* print unperturbed normal */ -register RAY *r; +RAY *r; { if (r->rot < FHUGE) { (*putreal)(r->ron[0]); @@ -498,7 +525,7 @@ register RAY *r; } -static +static void oputn(r) /* print perturbed normal */ RAY *r; { @@ -517,9 +544,9 @@ RAY *r; } -static +static void oputs(r) /* print name */ -register RAY *r; +RAY *r; { if (r->ro != NULL) fputs(r->ro->oname, stdout); @@ -529,17 +556,17 @@ register RAY *r; } -static +static void oputw(r) /* print weight */ -register RAY *r; +RAY *r; { (*putreal)(r->rweight); } -static +static void oputm(r) /* print modifier */ -register RAY *r; +RAY *r; { if (r->ro != NULL) if (r->ro->omod != OVOID) @@ -552,7 +579,7 @@ register RAY *r; } -static +static void puta(v) /* print ascii value */ double v; { @@ -560,7 +587,7 @@ double v; } -static +static void putd(v) /* print binary double */ double v; { @@ -568,7 +595,7 @@ double v; } -static +static void putf(v) /* print binary float */ double v; {