ViewVC Help

View File | Revision Log | Show Annotations | Download File | Root Listing

root/Development/ray/src/rt/aniso.c

Comparing ray/src/rt/aniso.c (file contents):
Revision 2.22 by greg, Fri Oct 16 10:20:27 1992 UTC vs.
Revision 2.23 by greg, Fri Feb 12 10:41:05 1993 UTC

#	Line 109 | Line 109 | double  omega;                 /* light source size */
109		h[0] = ldir[0] - np->rp->rdir[0];
110		h[1] = ldir[1] - np->rp->rdir[1];
111		h[2] = ldir[2] - np->rp->rdir[2];
112	–	normalize(h);
112		/* ellipse */
113		dtmp1 = DOT(np->u, h);
114		dtmp1 *= dtmp1 / au2;
115		dtmp2 = DOT(np->v, h);
116		dtmp2 *= dtmp2 / av2;
117		/* gaussian */
118	<	dtmp = (dtmp1 + dtmp2) / (1.0 + DOT(np->pnorm, h));
119	<	dtmp = exp(-2.0*dtmp) * (1.0/4.0/PI)
118	>	dtmp = DOT(np->pnorm, h);
119	>	dtmp = (dtmp1 + dtmp2) / (dtmp*dtmp);
120	>	dtmp = exp(-dtmp) * (0.25/PI)
121		* sqrt(ldot/(np->pdot*au2*av2));
122		/* worth using? */
123		if (dtmp > FTINY) {
#	Line 155 | Line 155 | double  omega;                 /* light source size */
155		dtmp = 1.0 - dtmp1*dtmp1/dtmp;
156		if (dtmp > FTINY*FTINY) {
157		dtmp1 = DOT(h,np->u);
158	<	dtmp1 = dtmp1*dtmp1 / au2;
158	>	dtmp1 *= dtmp1 / au2;
159		dtmp2 = DOT(h,np->v);
160	<	dtmp2 = dtmp2*dtmp2 / av2;
160	>	dtmp2 *= dtmp2 / av2;
161		dtmp = (dtmp1 + dtmp2) / dtmp;
162		}
163		} else

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines (old)
>	Changed lines (new)