--- ray/src/rt/m_direct.c	1991/07/16 15:44:30	1.1
+++ ray/src/rt/m_direct.c	1995/09/07 12:35:21	2.6
@@ -1,4 +1,4 @@
-/* Copyright (c) 1991 Regents of the University of California */
+/* Copyright (c) 1994 Regents of the University of California */
 
 #ifndef lint
 static char SCCSid[] = "$SunId$ LBL";
@@ -15,35 +15,39 @@ static char SCCSid[] = "$SunId$ LBL";
 
 #include  "source.h"
 
+#include  "func.h"
+
 /*
  * The arguments for MAT_DIRECT1 are:
  *
- *	4+ coef1 dx1 dy1 dz1 transform..
+ *	5+ coef1 dx1 dy1 dz1 funcfile transform..
  *	0
  *	n A1 A2 .. An
  *
  * The arguments for MAT_DIRECT2 are:
  *
- *	8+ coef1 dx1 dy1 dz1 coef2 dx2 dy2 dz2 transform..
+ *	9+ coef1 dx1 dy1 dz1 coef2 dx2 dy2 dz2 funcfile transform..
  *	0
  *	n A1 A2 .. An
  */
 
 
-extern double  varvalue();
-
 int  dir_proj();
 VSMATERIAL  direct1_vs = {dir_proj, 1};
 VSMATERIAL  direct2_vs = {dir_proj, 2};
 
+#define getdfunc(m)	( (m)->otype == MAT_DIRECT1 ? \
+				getfunc(m, 4, 0xf, 1) : \
+				getfunc(m, 8, 0xff, 1) )
 
+
 m_direct(m, r)			/* shade redirected ray */
 register OBJREC  *m;
 register RAY  *r;
 {
 					/* check if source ray */
 	if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
-		return;				/* got the wrong guy */
+		return(1);			/* got the wrong guy */
 					/* compute first projection */
 	if (m->otype == MAT_DIRECT1 ||
 			(r->rsrc < 0 || source[r->rsrc].sa.sv.pn == 0))
@@ -52,6 +56,7 @@ register RAY  *r;
 	if (m->otype == MAT_DIRECT2 &&
 			(r->rsrc < 0 || source[r->rsrc].sa.sv.pn == 1))
 		redirect(m, r, 1);
+	return(1);
 }
 
 
@@ -60,27 +65,54 @@ OBJREC  *m;
 RAY  *r;
 int  n;
 {
-	register char  **sa;
+	MFUNC  *mf;
+	register EPNODE  **va;
+	FVECT  nsdir;
 	RAY  nr;
 	double  coef;
 	register int  j;
 					/* set up function */
+	mf = getdfunc(m);
 	setfunc(m, r);
-	if (m->oargs.nsargs < 4+4*n)
-		objerror(m, USER, "too few arguments");
-	sa = m->oargs.sarg + 4*n;
+					/* assign direction variable */
+	if (r->rsrc >= 0) {
+		register SRCREC  *sp = source + source[r->rsrc].sa.sv.sn;
+
+		if (sp->sflags & SDISTANT)
+			VCOPY(nsdir, sp->sloc);
+		else {
+			for (j = 0; j < 3; j++)
+				nsdir[j] = sp->sloc[j] - r->rop[j];
+			normalize(nsdir);
+		}
+		multv3(nsdir, nsdir, funcxf.xfm);
+		varset("DxA", '=', nsdir[0]/funcxf.sca);
+		varset("DyA", '=', nsdir[1]/funcxf.sca);
+		varset("DzA", '=', nsdir[2]/funcxf.sca);
+	} else {
+		varset("DxA", '=', 0.0);
+		varset("DyA", '=', 0.0);
+		varset("DzA", '=', 0.0);
+	}
 					/* compute coefficient */
 	errno = 0;
-	coef = varvalue(sa[0]);
+	va = mf->ep + 4*n;
+	coef = evalue(va[0]);
 	if (errno)
 		goto computerr;
 	if (coef <= FTINY || rayorigin(&nr, r, TRANS, coef) < 0)
 		return(0);
-					/* compute direction */
-	errno = 0;
-	for (j = 0; j < 3; j++)
-		nr.rdir[j] = varvalue(sa[j+1]);
-	if (errno || normalize(nr.rdir) == 0.0)
+	va++;				/* compute direction */
+	for (j = 0; j < 3; j++) {
+		nr.rdir[j] = evalue(va[j]);
+		if (errno)
+			goto computerr;
+	}
+	if (mf->f != &unitxf)
+		multv3(nr.rdir, nr.rdir, mf->f->xfm);
+	if (r->rox != NULL)
+		multv3(nr.rdir, nr.rdir, r->rox->f.xfm);
+	if (normalize(nr.rdir) == 0.0)
 		goto computerr;
 					/* compute value */
 	if (r->rsrc >= 0)
@@ -102,16 +134,12 @@ SRCREC  *s;
 int  n;
 {
 	RAY  tr;
-	register OBJREC  *m;
-	char  **sa;
+	OBJREC  *m;
+	MFUNC  *mf;
+	EPNODE  **va;
 	FVECT  cent, newdir, nv, h;
-	double  olddot, newdot, od;
+	double  coef, olddot, newdot, od;
 	register int  i, j;
-				/* get material arguments */
-	m = objptr(o->omod);
-	if (m->oargs.nsargs < 4+4*n)
-		objerror(m, USER, "too few arguments");
-	sa = m->oargs.sarg + 4*n;
 				/* initialize test ray */
 	getmaxdisk(cent, o);
 	if (s->sflags & SDISTANT)
@@ -122,7 +150,7 @@ int  n;
 	else {
 		for (i = 0; i < 3; i++) {
 			tr.rdir[i] = cent[i] - s->sloc[i];
-			tr.rorg[i] = cent[i] - .05*tr.rdir[i];
+			tr.rorg[i] = s->sloc[i];
 		}
 		if (normalize(tr.rdir) == 0.0)
 			return(0);		/* at source! */
@@ -131,26 +159,39 @@ int  n;
 	olddot = DOT(tr.rdir, nv);
 	if (olddot <= FTINY && olddot >= -FTINY)
 		return(0);		/* old dir parallels plane */
+	tr.rmax = 0.0;
 	rayorigin(&tr, NULL, PRIMARY, 1.0);
 	if (!(*ofun[o->otype].funp)(o, &tr))
 		return(0);		/* no intersection! */
 				/* compute redirection */
+	m = vsmaterial(o);
+	mf = getdfunc(m);
 	setfunc(m, &tr);
+	varset("DxA", '=', 0.0);
+	varset("DyA", '=', 0.0);
+	varset("DzA", '=', 0.0);
 	errno = 0;
-	if (varvalue(sa[0]) <= FTINY)
-		return(0);		/* insignificant */
+	va = mf->ep + 4*n;
+	coef = evalue(va[0]);
 	if (errno)
 		goto computerr;
-	for (i = 0; i < 3; i++)
-		newdir[i] = varvalue(sa[i+1]);
-	if (errno)
-		goto computerr;
+	if (coef <= FTINY)
+		return(0);		/* insignificant */
+	va++;
+	for (i = 0; i < 3; i++) {
+		newdir[i] = evalue(va[i]);
+		if (errno)
+			goto computerr;
+	}
+	if (mf->f != &unitxf)
+		multv3(newdir, newdir, mf->f->xfm);
+					/* normalization unnecessary */
 	newdot = DOT(newdir, nv);
 	if (newdot <= FTINY && newdot >= -FTINY)
 		return(0);		/* new dir parallels plane */
 				/* everything OK -- compute shear */
 	for (i = 0; i < 3; i++)
-		h[i] = tr.rdir[i]/olddot + newdir[i]/newdot;
+		h[i] = newdir[i]/newdot - tr.rdir[i]/olddot;
 	setident4(pm);
 	for (j = 0; j < 3; j++) {
 		for (i = 0; i < 3; i++)