ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/rt/aniso.c
(Generate patch)

Comparing ray/src/rt/aniso.c (file contents):
Revision 2.32 by greg, Wed Apr 24 15:47:25 1996 UTC vs.
Revision 2.33 by gwlarson, Wed Dec 16 18:14:57 1998 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1996 Regents of the University of California */
1 > /* Copyright (c) 1998 Silicon Graphics, Inc. */
2  
3   #ifndef lint
4 < static char SCCSid[] = "$SunId$ LBL";
4 > static char SCCSid[] = "$SunId$ SGI";
5   #endif
6  
7   /*
# Line 359 | Line 359 | register ANISODAT  *np;
359                                  d = urand(ilhash(dimlist,ndims)+samplendx);
360                          multisamp(rv, 2, d);
361                          d = 2.0*PI * rv[0];
362 <                        cosp = cos(d) * np->u_alpha;
363 <                        sinp = sin(d) * np->v_alpha;
362 >                        cosp = tcos(d) * np->u_alpha;
363 >                        sinp = tsin(d) * np->v_alpha;
364                          d = sqrt(cosp*cosp + sinp*sinp);
365                          cosp /= d;
366                          sinp /= d;
# Line 397 | Line 397 | register ANISODAT  *np;
397                                  d = urand(ilhash(dimlist,ndims)+1823+samplendx);
398                          multisamp(rv, 2, d);
399                          d = 2.0*PI * rv[0];
400 <                        cosp = cos(d) * np->u_alpha;
401 <                        sinp = sin(d) * np->v_alpha;
400 >                        cosp = tcos(d) * np->u_alpha;
401 >                        sinp = tsin(d) * np->v_alpha;
402                          d = sqrt(cosp*cosp + sinp*sinp);
403                          cosp /= d;
404                          sinp /= d;
# Line 408 | Line 408 | register ANISODAT  *np;
408                          else
409                                  d = sqrt(-log(rv[1]) /
410                                          (cosp*cosp/(np->u_alpha*np->u_alpha) +
411 <                                         sinp*sinp/(np->v_alpha*np->u_alpha)));
411 >                                         sinp*sinp/(np->v_alpha*np->v_alpha)));
412                          for (i = 0; i < 3; i++)
413                                  sr.rdir[i] = np->prdir[i] +
414                                                  d*(cosp*np->u[i] + sinp*np->v[i]);

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines