--- ray/src/rt/normal.c	2003/02/25 02:47:22	2.39
+++ ray/src/rt/normal.c	2005/01/05 19:34:11	2.49
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char	RCSid[] = "$Id: normal.c,v 2.39 2003/02/25 02:47:22 greg Exp $";
+static const char RCSid[] = "$Id: normal.c,v 2.49 2005/01/05 19:34:11 greg Exp $";
 #endif
 /*
  *  normal.c - shading function for normal materials.
@@ -14,18 +14,18 @@ static const char	RCSid[] = "$Id: normal.c,v 2.39 2003
 #include "copyright.h"
 
 #include  "ray.h"
-
+#include  "ambient.h"
+#include  "source.h"
 #include  "otypes.h"
-
+#include  "rtotypes.h"
 #include  "random.h"
 
 #ifndef  MAXITER
 #define  MAXITER	10		/* maximum # specular ray attempts */
 #endif
 					/* estimate of Fresnel function */
-#define  FRESNE(ci)	(exp(-6.0*(ci)) - 0.00247875217)
+#define  FRESNE(ci)	(exp(-5.85*(ci)) - 0.00287989916)
 
-static void  gaussamp();
 
 /*
  *	This routine implements the isotropic Gaussian
@@ -64,16 +64,21 @@ typedef struct {
 	double  pdot;		/* perturbed dot product */
 }  NORMDAT;		/* normal material data */
 
+static srcdirf_t dirnorm;
+static void gaussamp(RAY  *r, NORMDAT  *np);
 
+
 static void
-dirnorm(cval, np, ldir, omega)		/* compute source contribution */
-COLOR  cval;			/* returned coefficient */
-register NORMDAT  *np;		/* material data */
-FVECT  ldir;			/* light source direction */
-double  omega;			/* light source size */
+dirnorm(		/* compute source contribution */
+	COLOR  cval,			/* returned coefficient */
+	void  *nnp,		/* material data */
+	FVECT  ldir,			/* light source direction */
+	double  omega			/* light source size */
+)
 {
+	register NORMDAT *np = nnp;
 	double  ldot;
-	double  ldiff;
+	double  lrdiff, ltdiff;
 	double  dtmp, d2;
 	FVECT  vtmp;
 	COLOR  ctmp;
@@ -86,18 +91,23 @@ double  omega;			/* light source size */
 		return;		/* wrong side */
 
 				/* Fresnel estimate */
-	ldiff = np->rdiff;
-	if (np->specfl & SP_PURE && (np->rspec > FTINY & ldiff > FTINY))
-		ldiff *= 1. - FRESNE(fabs(ldot));
+	lrdiff = np->rdiff;
+	ltdiff = np->tdiff;
+	if (np->specfl & SP_PURE && np->rspec > FTINY &&
+			(lrdiff > FTINY) | (ltdiff > FTINY)) {
+		dtmp = 1. - FRESNE(fabs(ldot));
+		lrdiff *= dtmp;
+		ltdiff *= dtmp;
+	}
 
-	if (ldot > FTINY && ldiff > FTINY) {
+	if (ldot > FTINY && lrdiff > FTINY) {
 		/*
 		 *  Compute and add diffuse reflected component to returned
 		 *  color.  The diffuse reflected component will always be
 		 *  modified by the color of the material.
 		 */
 		copycolor(ctmp, np->mcolor);
-		dtmp = ldot * omega * ldiff / PI;
+		dtmp = ldot * omega * lrdiff * (1.0/PI);
 		scalecolor(ctmp, dtmp);
 		addcolor(cval, ctmp);
 	}
@@ -110,7 +120,7 @@ double  omega;			/* light source size */
 		dtmp = np->alpha2;
 						/* + source if flat */
 		if (np->specfl & SP_FLAT)
-			dtmp += omega/(4.0*PI);
+			dtmp += omega * (0.25/PI);
 						/* half vector */
 		vtmp[0] = ldir[0] - np->rp->rdir[0];
 		vtmp[1] = ldir[1] - np->rp->rdir[1];
@@ -119,21 +129,21 @@ double  omega;			/* light source size */
 		d2 *= d2;
 		d2 = (DOT(vtmp,vtmp) - d2) / d2;
 						/* gaussian */
-		dtmp = exp(-d2/dtmp)/(4.*PI*dtmp);
+		dtmp = exp(-d2/dtmp)/(4.*PI * np->pdot * dtmp);
 						/* worth using? */
 		if (dtmp > FTINY) {
 			copycolor(ctmp, np->scolor);
-			dtmp *= omega * sqrt(ldot/np->pdot);
+			dtmp *= omega;
 			scalecolor(ctmp, dtmp);
 			addcolor(cval, ctmp);
 		}
 	}
-	if (ldot < -FTINY && np->tdiff > FTINY) {
+	if (ldot < -FTINY && ltdiff > FTINY) {
 		/*
 		 *  Compute diffuse transmission.
 		 */
 		copycolor(ctmp, np->mcolor);
-		dtmp = -ldot * omega * np->tdiff / PI;
+		dtmp = -ldot * omega * ltdiff * (1.0/PI);
 		scalecolor(ctmp, dtmp);
 		addcolor(cval, ctmp);
 	}
@@ -143,13 +153,14 @@ double  omega;			/* light source size */
 		 *  is always modified by material color.
 		 */
 						/* roughness + source */
-		dtmp = np->alpha2 + omega/PI;
+		dtmp = np->alpha2 + omega*(1.0/PI);
 						/* gaussian */
-		dtmp = exp((2.*DOT(np->prdir,ldir)-2.)/dtmp)/(PI*dtmp);
+		dtmp = exp((2.*DOT(np->prdir,ldir)-2.)/dtmp) /
+					(PI*np->pdot*dtmp);
 						/* worth using? */
 		if (dtmp > FTINY) {
 			copycolor(ctmp, np->mcolor);
-			dtmp *= np->tspec * omega * sqrt(-ldot/np->pdot);
+			dtmp *= np->tspec * omega;
 			scalecolor(ctmp, dtmp);
 			addcolor(cval, ctmp);
 		}
@@ -157,10 +168,11 @@ double  omega;			/* light source size */
 }
 
 
-int
-m_normal(m, r)			/* color a ray that hit something normal */
-register OBJREC  *m;
-register RAY  *r;
+extern int
+m_normal(			/* color a ray that hit something normal */
+	register OBJREC  *m,
+	register RAY  *r
+)
 {
 	NORMDAT  nd;
 	double  fest;
@@ -182,8 +194,10 @@ register RAY  *r;
 			raytrans(r);
 			return(1);
 		}
+		raytexture(r, m->omod);
 		flipsurface(r);			/* reorient if backvis */
-	}
+	} else
+		raytexture(r, m->omod);
 	nd.mp = m;
 	nd.rp = r;
 						/* get material color */
@@ -195,16 +209,15 @@ register RAY  *r;
 	nd.alpha2 = m->oargs.farg[4];
 	if ((nd.alpha2 *= nd.alpha2) <= FTINY)
 		nd.specfl |= SP_PURE;
-	if (r->ro != NULL && isflat(r->ro->otype))
-		nd.specfl |= SP_FLAT;
-						/* get modifiers */
-	raytexture(r, m->omod);
-	if (hastexture = DOT(r->pert,r->pert) > FTINY*FTINY)
+
+	if ( (hastexture = (DOT(r->pert,r->pert) > FTINY*FTINY)) ) {
 		nd.pdot = raynormal(nd.pnorm, r);	/* perturb normal */
-	else {
+	} else {
 		VCOPY(nd.pnorm, r->ron);
 		nd.pdot = r->rod;
 	}
+	if (r->ro != NULL && isflat(r->ro->otype))
+		nd.specfl |= SP_FLAT;
 	if (nd.pdot < .001)
 		nd.pdot = .001;			/* non-zero for dirnorm() */
 	multcolor(nd.mcolor, r->pcol);		/* modify material color */
@@ -350,9 +363,10 @@ register RAY  *r;
 
 
 static void
-gaussamp(r, np)			/* sample gaussian specular */
-RAY  *r;
-register NORMDAT  *np;
+gaussamp(			/* sample gaussian specular */
+	RAY  *r,
+	register NORMDAT  *np
+)
 {
 	RAY  sr;
 	FVECT  u, v, h;