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

Comparing ray/src/rt/glass.c (file contents):
Revision 1.11 by greg, Tue Oct 29 16:30:46 1991 UTC vs.
Revision 2.3 by greg, Fri Oct 2 16:17:31 1992 UTC

# Line 32 | Line 32 | static char SCCSid[] = "$SunId$ LBL";
32   *      tn = (sqrt(.8402528435+.0072522239*Tn*Tn)-.9166530661)/.0036261119/Tn
33   *
34   *  The transmission of standard 88% transmissivity glass is 0.96.
35 + *  A refractive index other than the default can be used by giving
36 + *  it as the fourth real argument.  The above formula no longer applies.
37 + *
38   *  If we appear to hit the back side of the surface, then we
39   *  turn the normal around.
40   */
# Line 43 | Line 46 | m_glass(m, r)          /* color a ray which hit a thin glass s
46   OBJREC  *m;
47   register RAY  *r;
48   {
46 –        double  sqrt(), pow();
49          COLOR  mcolor;
50          double  pdot;
51          FVECT  pnorm;
52 <        double  cos2;
52 >        double  rindex, cos2;
53          COLOR  trans, refl;
54          double  d, r1e, r1m;
55          double  transtest, transdist;
56          RAY  p;
57          register int  i;
58 <
59 <        if (m->oargs.nfargs != 3)
58 >                                                /* check arguments */
59 >        if (m->oargs.nfargs == 3)
60 >                rindex = RINDEX;                /* default value of n */
61 >        else if (m->oargs.nfargs == 4)
62 >                rindex = m->oargs.farg[3];      /* use their value */
63 >        else
64                  objerror(m, USER, "bad arguments");
65  
66          setcolor(mcolor, m->oargs.farg[0], m->oargs.farg[1], m->oargs.farg[2]);
# Line 66 | Line 72 | register RAY  *r;
72          raytexture(r, m->omod);
73          pdot = raynormal(pnorm, r);
74                                                  /* angular transmission */
75 <        cos2 = sqrt( (1.0-1.0/RINDEX/RINDEX) +
76 <                     pdot*pdot/(RINDEX*RINDEX) );
75 >        cos2 = sqrt( (1.0-1.0/(rindex*rindex)) +
76 >                     pdot*pdot/(rindex*rindex) );
77          setcolor(mcolor, pow(colval(mcolor,RED), 1.0/cos2),
78                           pow(colval(mcolor,GRN), 1.0/cos2),
79                           pow(colval(mcolor,BLU), 1.0/cos2));
80  
81                                                  /* compute reflection */
82 <        r1e = (pdot - RINDEX*cos2) / (pdot + RINDEX*cos2);
82 >        r1e = (pdot - rindex*cos2) / (pdot + rindex*cos2);
83          r1e *= r1e;
84 <        r1m = (1.0/pdot - RINDEX/cos2) / (1.0/pdot + RINDEX/cos2);
84 >        r1m = (1.0/pdot - rindex/cos2) / (1.0/pdot + rindex/cos2);
85          r1m *= r1m;
86                                                  /* compute transmittance */
87          for (i = 0; i < 3; i++) {
# Line 88 | Line 94 | register RAY  *r;
94                  if (!(r->crtype & SHADOW) &&
95                                  DOT(r->pert,r->pert) > FTINY*FTINY) {
96                          for (i = 0; i < 3; i++)         /* perturb direction */
97 <                                p.rdir[i] = r->rdir[i] - r->pert[i]/RINDEX;
97 >                                p.rdir[i] = r->rdir[i] +
98 >                                                2.*(1.-rindex)*r->pert[i];
99                          normalize(p.rdir);
100                  } else {
101                          VCOPY(p.rdir, r->rdir);

Diff Legend

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