--- ray/src/rt/rtrace.c	1995/12/08 18:49:12	2.20
+++ ray/src/rt/rtrace.c	2003/06/05 19:29:34	2.30
@@ -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.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,24 +21,31 @@ static char SCCSid[] = "$SunId$ LBL";
  *  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 */
 
 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 */
+
+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 */
@@ -59,7 +63,7 @@ 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 */
 
@@ -71,7 +75,9 @@ 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 */
 int  ambres = 128;			/* ambient resolution */
 int  ambdiv = 512;			/* ambient divisions */
@@ -80,22 +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 */
 
 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;
 {
@@ -121,12 +129,14 @@ int  f;
 }
 
 
+void
 rtrace(fname)				/* 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 +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);
@@ -161,20 +171,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");
@@ -192,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)
@@ -218,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;
@@ -245,21 +263,39 @@ register char  *vs;
 }
 
 
-rad(org, dir)			/* compute and print ray value(s) */
+void
+bogusray()			/* 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);
+}
+
+
+void
+rad(org, dir, dmax)		/* compute and print ray value(s) */
 FVECT  org, 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);
 }
 
 
+void
 irrad(org, dir)			/* compute immediate irradiance value */
 FVECT  org, dir;
 {
@@ -282,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;
@@ -296,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];
@@ -333,6 +370,7 @@ FILE  *fp;
 }
 
 
+void
 tranotify(obj)			/* record new modifier */
 OBJECT	obj;
 {
@@ -340,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++)
@@ -355,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;
@@ -375,7 +418,7 @@ RAY  *r;
 }
 
 
-static
+static void
 tabin(r)				/* tab in appropriate amount */
 RAY  *r;
 {
@@ -386,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]);
@@ -396,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]);
@@ -406,9 +449,9 @@ register RAY  *r;
 }
 
 
-static
+static void
 oputv(r)				/* print value */
-register RAY  *r;
+RAY  *r;
 {
 	COLR  cout;
 	
@@ -425,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]);
@@ -457,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]);
@@ -473,7 +525,7 @@ register RAY  *r;
 }
 
 
-static
+static void
 oputn(r)				/* print perturbed normal */
 RAY  *r;
 {
@@ -492,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);
@@ -504,27 +556,30 @@ 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)
-		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
+static void
 puta(v)				/* print ascii value */
 double  v;
 {
@@ -532,7 +587,7 @@ double  v;
 }
 
 
-static
+static void
 putd(v)				/* print binary double */
 double  v;
 {
@@ -540,7 +595,7 @@ double  v;
 }
 
 
-static
+static void
 putf(v)				/* print binary float */
 double  v;
 {