--- ray/src/rt/m_brdf.c	1991/06/17 08:53:58	1.10
+++ ray/src/rt/m_brdf.c	1991/08/12 08:20:53	1.16
@@ -191,7 +191,7 @@ register RAY  *r;
 	BRDFDAT  nd;
 	double  transtest, transdist;
 	COLOR  ctmp;
-	double  dtmp;
+	double  dtmp, tspect, rspecr;
 	register int  i;
 						/* check arguments */
 	switch (m->otype) {
@@ -246,32 +246,31 @@ register RAY  *r;
 				break;
 		if (i-3 != nd.dp->nd)
 			objerror(m, USER, "dimension error");
-		if (!fundefined(m->oargs.sarg[3]))
-			loadfunc(m->oargs.sarg[2]);
+		funcfile(m->oargs.sarg[2]);
 	} else if (m->otype == MAT_BRTDF) {
 		nd.dp = NULL;
-		if (!fundefined(m->oargs.sarg[7]))
-			loadfunc(m->oargs.sarg[9]);
+		funcfile(m->oargs.sarg[9]);
 	} else {
 		nd.dp = NULL;
-		if (!fundefined(m->oargs.sarg[0]))
-			loadfunc(m->oargs.sarg[1]);
+		funcfile(m->oargs.sarg[1]);
 	}
 						/* set special variables */
 	setbrdfunc(&nd);
 						/* compute transmitted ray */
+	tspect = 0.;
 	if (m->otype == MAT_BRTDF && nd.tspec > FTINY) {
 		RAY  sr;
 		errno = 0;
-		setcolor(ctmp, varvalue(m->oargs.sarg[0]),
-				varvalue(m->oargs.sarg[1]),
-				varvalue(m->oargs.sarg[2]));
-		scalecolor(ctmp, nd.tspec);
+		setcolor(ctmp, varvalue(m->oargs.sarg[3]),
+				varvalue(m->oargs.sarg[4]),
+				varvalue(m->oargs.sarg[5]));
+		scalecolor(ctmp, nd.trans);
 		if (errno)
 			objerror(m, WARNING, "compute error");
-		else if ((dtmp = bright(ctmp)) > FTINY &&
-				rayorigin(&sr, r, TRANS, dtmp) == 0) {
-			if (DOT(r->pert,r->pert) > FTINY*FTINY) {
+		else if ((tspect = bright(ctmp)) > FTINY &&
+				rayorigin(&sr, r, TRANS, tspect) == 0) {
+			if (!(r->crtype & SHADOW) &&
+					DOT(r->pert,r->pert) > FTINY*FTINY) {
 				for (i = 0; i < 3; i++)	/* perturb direction */
 					sr.rdir[i] = r->rdir[i] -
 							.75*r->pert[i];
@@ -290,17 +289,17 @@ register RAY  *r;
 	if (r->crtype & SHADOW)			/* the rest is shadow */
 		return;
 						/* compute reflected ray */
+	rspecr = 0.;
 	if (m->otype == MAT_BRTDF && nd.rspec > FTINY) {
 		RAY  sr;
 		errno = 0;
-		setcolor(ctmp, varvalue(m->oargs.sarg[3]),
-				varvalue(m->oargs.sarg[4]),
-				varvalue(m->oargs.sarg[5]));
-		scalecolor(ctmp, nd.rspec);
+		setcolor(ctmp, varvalue(m->oargs.sarg[0]),
+				varvalue(m->oargs.sarg[1]),
+				varvalue(m->oargs.sarg[2]));
 		if (errno)
 			objerror(m, WARNING, "compute error");
-		else if ((dtmp = bright(ctmp)) > FTINY &&
-				rayorigin(&sr, r, REFLECTED, dtmp) == 0) {
+		else if ((rspecr = bright(ctmp)) > FTINY &&
+				rayorigin(&sr, r, REFLECTED, rspecr) == 0) {
 			for (i = 0; i < 3; i++)
 				sr.rdir[i] = r->rdir[i] +
 						2.0*nd.pdot*nd.pnorm[i];
@@ -310,22 +309,16 @@ register RAY  *r;
 		}
 	}
 						/* compute ambient */
-	if (nd.rdiff > FTINY) {
+	if ((dtmp = 1.0-nd.trans-rspecr) > FTINY) {
 		ambient(ctmp, r);
-		if (m->otype == MAT_BRTDF)
-			scalecolor(ctmp, nd.rdiff);
-		else
-			scalecolor(ctmp, 1.0-nd.trans);
+		scalecolor(ctmp, dtmp);
 		multcolor(ctmp, nd.mcolor);	/* modified by material color */
 		addcolor(r->rcol, ctmp);	/* add to returned color */
 	}
-	if (nd.tdiff > FTINY) {			/* from other side */
+	if ((dtmp = nd.trans-tspect) > FTINY) {	/* from other side */
 		flipsurface(r);
 		ambient(ctmp, r);
-		if (m->otype == MAT_BRTDF)
-			scalecolor(ctmp, nd.tdiff);
-		else
-			scalecolor(ctmp, nd.trans);
+		scalecolor(ctmp, dtmp);
 		multcolor(ctmp, nd.mcolor);
 		addcolor(r->rcol, ctmp);
 		flipsurface(r);
@@ -350,7 +343,8 @@ register BRDFDAT  *np;
 	varset("NxP", '=', vec[0]/funcxf.sca);
 	varset("NyP", '=', vec[1]/funcxf.sca);
 	varset("NzP", '=', vec[2]/funcxf.sca);
-	varset("RdotP", '=', np->pdot);
+	varset("RdotP", '=', np->pdot <= -1.0 ? -1.0 :
+			np->pdot >= 1.0 ? 1.0 : np->pdot);
 	varset("CrP", '=', colval(np->mcolor,RED));
 	varset("CgP", '=', colval(np->mcolor,GRN));
 	varset("CbP", '=', colval(np->mcolor,BLU));