--- ray/src/rt/rv2.c	1994/12/20 20:18:37	2.26
+++ ray/src/rt/rv2.c	2003/02/25 02:47:23	2.39
@@ -1,27 +1,25 @@
-/* Copyright (c) 1994 Regents of the University of California */
-
 #ifndef lint
-static char SCCSid[] = "$SunId$ LBL";
+static const char	RCSid[] = "$Id: rv2.c,v 2.39 2003/02/25 02:47:23 greg Exp $";
 #endif
-
 /*
  *  rv2.c - command routines used in tracing a view.
  *
- *     3/24/87
+ *  External symbols declared in rpaint.h
  */
 
+#include "copyright.h"
+
 #include  "ray.h"
 
-#include  "octree.h"
-
 #include  "otypes.h"
 
 #include  "rpaint.h"
 
-#include  "resolu.h"
-
 #include  <ctype.h>
 
+extern int  psample;			/* pixel sample size */
+extern double  maxdiff;			/* max. sample difference */
+
 #define  CTRL(c)	((c)-'@')
 
 #ifdef  SMLFLT
@@ -30,17 +28,13 @@ static char SCCSid[] = "$SunId$ LBL";
 #define  sscanvec(s,v)	(sscanf(s,"%lf %lf %lf",v,v+1,v+2)==3)
 #endif
 
-extern char  *atos();
-
-extern FILE  *popen();
-
 extern char  rifname[128];		/* rad input file name */
 
-extern char  VersionID[];
 extern char  *progname;
 extern char  *octname;
 
 
+void
 getframe(s)				/* get a new frame */
 char  *s;
 {
@@ -50,6 +44,7 @@ char  *s;
 }
 
 
+void
 getrepaint(s)				/* get area and repaint */
 char  *s;
 {
@@ -61,6 +56,7 @@ char  *s;
 }
 
 
+void
 getview(s)				/* get/show view parameters */
 char  *s;
 {
@@ -165,6 +161,7 @@ char  *s;
 }
 
 
+void
 lastview(s)				/* return to a previous view */
 char  *s;
 {
@@ -198,6 +195,7 @@ char  *s;
 }
 
 
+void
 saveview(s)				/* save view to rad file */
 char  *s;
 {
@@ -234,6 +232,7 @@ char  *s;
 }
 
 
+void
 loadview(s)				/* load view from rad file */
 char  *s;
 {
@@ -275,6 +274,7 @@ char  *s;
 }
 
 
+void
 getaim(s)				/* aim camera */
 char  *s;
 {
@@ -286,13 +286,15 @@ char  *s;
 	nv.type = ourview.type;
 	VCOPY(nv.vp, ourview.vp);
 	VCOPY(nv.vup, ourview.vup);
-	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
 	nv.horiz = ourview.horiz; nv.vert = ourview.vert;
+	nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
+	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
 	zoomview(&nv, zfact);
 	newview(&nv);
 }
 
 
+void
 getmove(s)				/* move camera */
 char  *s;
 {
@@ -305,6 +307,7 @@ char  *s;
 }
 
 
+void
 getrotate(s)				/* rotate camera */
 char  *s;
 {
@@ -321,6 +324,7 @@ char  *s;
 	VCOPY(nv.vp, ourview.vp);
 	VCOPY(nv.vup, ourview.vup);
 	nv.hoff = ourview.hoff; nv.voff = ourview.voff;
+	nv.vfore = ourview.vfore; nv.vaft = ourview.vaft;
 	spinvector(nv.vdir, ourview.vdir, ourview.vup, angle*(PI/180.));
 	if (elev != 0.0) {
 		fcross(v1, nv.vdir, ourview.vup);
@@ -333,6 +337,7 @@ char  *s;
 }
 
 
+void
 getpivot(s)				/* pivot viewpoint */
 register char  *s;
 {
@@ -344,12 +349,13 @@ register char  *s;
 		error(COMMAND, "missing angle");
 		return;
 	}
-	if (getinterest(sskip(sskip(s)), 0, vc, &mag) < 0)
+	if (getinterest(sskip2(s,2), 0, vc, &mag) < 0)
 		return;
 	moveview(angle, elev, mag, vc);
 }
 
 
+void
 getexposure(s)				/* get new exposure */
 char  *s;
 {
@@ -415,13 +421,15 @@ char  *s;
 	redraw();
 }
 
+typedef union {int i; double d; COLOR C;}	*MyUptr;
 
-getparam(str, dsc, typ, ptr)		/* get variable from user */
+int
+getparam(str, dsc, typ, p)		/* get variable from user */
 char  *str, *dsc;
 int  typ;
-register union {int i; double d; COLOR C;}  *ptr;
+void  *p;
 {
-	extern char  *index();
+	register MyUptr  ptr = (MyUptr)p;
 	int  i0;
 	double  d0, d1, d2;
 	char  buf[48];
@@ -479,32 +487,15 @@ register union {int i; double d; COLOR C;}  *ptr;
 }
 
 
+void
 setparam(s)				/* get/set program parameter */
 register char  *s;
 {
-	extern int  psample;
-	extern double  maxdiff;
-	extern double  minweight;
-	extern int  maxdepth;
-	extern double  dstrsrc;
-	extern double  shadthresh;
-	extern double  shadcert;
-	extern COLOR  ambval;
-	extern double  ambacc;
-	extern int  ambres;
-	extern int  ambdiv;
-	extern int  ambssamp;
-	extern int  ambounce;
-	extern int  directvis;
-	extern double  srcsizerat;
-	extern int  do_irrad;
-	extern double  specjitter;
-	extern double  specthresh;
 	char  buf[128];
 	
 	if (s[0] == '\0') {
 		(*dev->comout)(
-		"aa ab ad ar as av b dc di dj ds dt i lr lw ps pt sj st: ");
+		"aa ab ad ar as av aw b dc dv dj ds dt i lr lw me ma mg ms ps pt sj st bv: ");
 		(*dev->comin)(buf, NULL);
 		s = buf;
 	}
@@ -512,10 +503,12 @@ register char  *s;
 	case 'l':			/* limit */
 		switch (s[1]) {
 		case 'w':			/* weight */
-			getparam(s+2, "limit weight", 'r', &minweight);
+			getparam(s+2, "limit weight", 'r',
+					(void *)&minweight);
 			break;
 		case 'r':			/* reflection */
-			getparam(s+2, "limit reflection", 'i', &maxdepth);
+			getparam(s+2, "limit reflection", 'i',
+					(void *)&maxdepth);
 			break;
 		default:
 			goto badparam;
@@ -524,65 +517,118 @@ register char  *s;
 	case 'd':			/* direct */
 		switch (s[1]) {
 		case 'j':			/* jitter */
-			getparam(s+2, "direct jitter", 'r', &dstrsrc);
+			getparam(s+2, "direct jitter", 'r',
+					(void *)&dstrsrc);
 			break;
 		case 'c':			/* certainty */
-			getparam(s+2, "direct certainty", 'r', &shadcert);
+			getparam(s+2, "direct certainty", 'r',
+					(void *)&shadcert);
 			break;
 		case 't':			/* threshold */
-			getparam(s+2, "direct threshold", 'r', &shadthresh);
+			getparam(s+2, "direct threshold", 'r',
+					(void *)&shadthresh);
 			break;
 		case 'v':			/* visibility */
-			getparam(s+2, "direct visibility",
-					'b', &directvis);
+			getparam(s+2, "direct visibility", 'b',
+					(void *)&directvis);
 			break;
 		case 's':			/* sampling */
-			getparam(s+2, "direct sampling", 'r', &srcsizerat);
+			getparam(s+2, "direct sampling", 'r',
+					(void *)&srcsizerat);
 			break;
 		default:
 			goto badparam;
 		}
 		break;
-	case 'b':			/* black and white */
-		getparam(s+1, "black and white", 'b', &greyscale);
+	case 'b':			/* back faces or black and white */
+		switch (s[1]) {
+		case 'v':			/* back face visibility */
+			getparam(s+2, "back face visibility", 'b',
+					(void *)&backvis);
+			break;
+		case '\0':			/* black and white */
+		case ' ':
+		case 'y': case 'Y': case 't': case 'T': case '1': case '+':
+		case 'n': case 'N': case 'f': case 'F': case '0': case '-':
+			getparam(s+1, "black and white", 'b',
+					(void *)&greyscale);
+			break;
+		default:
+			goto badparam;
+		}
 		break;
 	case 'i':			/* irradiance */
-		getparam(s+1, "irradiance", 'b', &do_irrad);
+		getparam(s+1, "irradiance", 'b',
+				(void *)&do_irrad);
 		break;
 	case 'a':			/* ambient */
 		switch (s[1]) {
 		case 'v':			/* value */
-			getparam(s+2, "ambient value", 'C', (COLOR *)ambval);
+			getparam(s+2, "ambient value", 'C',
+					(void *)ambval);
 			break;
+		case 'w':			/* weight */
+			getparam(s+2, "ambient value weight", 'i',
+					(void *)&ambvwt);
+			break;
 		case 'a':			/* accuracy */
-			if (getparam(s+2, "ambient accuracy", 'r', &ambacc))
+			if (getparam(s+2, "ambient accuracy", 'r',
+					(void *)&ambacc))
 				setambacc(ambacc);
 			break;
 		case 'd':			/* divisions */
-			getparam(s+2, "ambient divisions", 'i', &ambdiv);
+			getparam(s+2, "ambient divisions", 'i',
+					(void *)&ambdiv);
 			break;
 		case 's':			/* samples */
-			getparam(s+2, "ambient super-samples", 'i', &ambssamp);
+			getparam(s+2, "ambient super-samples", 'i',
+					(void *)&ambssamp);
 			break;
 		case 'b':			/* bounces */
-			getparam(s+2, "ambient bounces", 'i', &ambounce);
+			getparam(s+2, "ambient bounces", 'i',
+					(void *)&ambounce);
 			break;
 		case 'r':
-			if (getparam(s+2, "ambient resolution", 'i', &ambres))
+			if (getparam(s+2, "ambient resolution", 'i',
+					(void *)&ambres))
 				setambres(ambres);
 			break;
 		default:
 			goto badparam;
 		}
 		break;
+	case 'm':			/* medium */
+		switch (s[1]) {
+		case 'e':			/* extinction coefficient */
+			getparam(s+2, "extinction coefficient", 'C',
+					(void *)cextinction);
+			break;
+		case 'a':			/* scattering albedo */
+			getparam(s+2, "scattering albedo", 'C',
+					(void *)salbedo);
+			break;
+		case 'g':			/* scattering eccentricity */
+			getparam(s+2, "scattering eccentricity", 'r',
+					(void *)&seccg);
+			break;
+		case 's':			/* sampling distance */
+			getparam(s+2, "mist sampling distance", 'r',
+					(void *)&ssampdist);
+			break;
+		default:
+			goto badparam;
+		}
+		break;
 	case 'p':			/* pixel */
 		switch (s[1]) {
 		case 's':			/* sample */
-			if (getparam(s+2, "pixel sample", 'i', &psample))
+			if (getparam(s+2, "pixel sample", 'i',
+					(void *)&psample))
 				pdepth = 0;
 			break;
 		case 't':			/* threshold */
-			if (getparam(s+2, "pixel threshold", 'r', &maxdiff))
+			if (getparam(s+2, "pixel threshold", 'r',
+					(void *)&maxdiff))
 				pdepth = 0;
 			break;
 		default:
@@ -592,10 +638,12 @@ register char  *s;
 	case 's':			/* specular */
 		switch (s[1]) {
 		case 'j':			/* jitter */
-			getparam(s+2, "specular jitter", 'r', &specjitter);
+			getparam(s+2, "specular jitter", 'r',
+					(void *)&specjitter);
 			break;
 		case 't':			/* threshold */
-			getparam(s+2, "specular threshold", 'r', &specthresh);
+			getparam(s+2, "specular threshold", 'r',
+					(void *)&specthresh);
 			break;
 		default:
 			goto badparam;
@@ -613,17 +661,19 @@ badparam:
 }
 
 
+void
 traceray(s)				/* trace a single ray */
 char  *s;
 {
 	char  buf[128];
 	int  x, y;
+	OBJREC	*ino;
 	RAY  thisray;
 
 	thisray.rmax = 0.0;
 
 	if (!sscanvec(s, thisray.rorg) ||
-			!sscanvec(sskip(sskip(sskip(s))), thisray.rdir)) {
+			!sscanvec(sskip2(s,3), thisray.rdir)) {
 
 		if (dev->getcur == NULL)
 			return;
@@ -655,17 +705,21 @@ char  *s;
 					objptr(thisray.ro->omod)->oname,
 				ofun[thisray.ro->otype].funame,
 				thisray.ro->oname);
+		if ((ino = objptr(thisray.robj)) != thisray.ro)
+			sprintf(buf+strlen(buf), " in %s \"%s\"",
+					ofun[ino->otype].funame, ino->oname);
 		(*dev->comout)(buf);
 		(*dev->comin)(buf, NULL);
 		if (thisray.rot >= FHUGE)
 			(*dev->comout)("at infinity");
 		else {
-			sprintf(buf, "at (%.6g %.6g %.6g)", thisray.rop[0],
-					thisray.rop[1], thisray.rop[2]);
+			sprintf(buf, "at (%.6g %.6g %.6g) (%.6g)",
+					thisray.rop[0], thisray.rop[1],
+					thisray.rop[2], thisray.rt);
 			(*dev->comout)(buf);
 		}
 		(*dev->comin)(buf, NULL);
-		sprintf(buf, "value (%.5g %.5g %.5g) (%.1fL)",
+		sprintf(buf, "value (%.5g %.5g %.5g) (%.3gL)",
 				colval(thisray.rcol,RED),
 				colval(thisray.rcol,GRN),
 				colval(thisray.rcol,BLU),
@@ -676,6 +730,7 @@ char  *s;
 }
 
 
+void
 writepict(s)				/* write the picture to a file */
 char  *s;
 {
@@ -714,6 +769,7 @@ char  *s;
 	else
 		putc('\n', fp);
 	fprintf(fp, "SOFTWARE= %s\n", VersionID);
+	fputnow(fp);
 	if (exposure != 1.0)
 		fputexpos(exposure, fp);
 	if (dev->pixaspect != 1.0)
@@ -734,7 +790,7 @@ char  *s;
 		if (fwritecolrs(scanline, hresolu, fp) < 0)
 			break;
 	}
-	free((char *)scanline);
+	free((void *)scanline);
 	if (fclose(fp) < 0)
 		error(COMMAND, "write error");
 }