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 2.16 by greg, Tue Apr 19 01:15:06 2005 UTC vs.
Revision 2.17 by greg, Tue Jun 21 15:06:50 2005 UTC

# Line 52 | Line 52 | m_glass(               /* color a ray which hit a thin glass surfac
52          FVECT  pnorm;
53          double  rindex, cos2;
54          COLOR  trans, refl;
55 <        int  hastexture;
55 >        int  hastexture, hastrans;
56          double  d, r1e, r1m;
57          double  transtest, transdist;
58          double  mirtest, mirdist;
# Line 65 | Line 65 | m_glass(               /* color a ray which hit a thin glass surfac
65                  rindex = m->oargs.farg[3];      /* use their value */
66          else
67                  objerror(m, USER, "bad arguments");
68 <
68 >                                                /* check transmission */
69          setcolor(mcolor, m->oargs.farg[0], m->oargs.farg[1], m->oargs.farg[2]);
70 <
70 >        if ((hastrans = (intens(mcolor) > 1e-15))) {
71 >                for (i = 0; i < 3; i++)
72 >                        if (colval(mcolor,i) < 1e-15)
73 >                                colval(mcolor,i) = 1e-15;
74 >        } else if (r->crtype & SHADOW)
75 >                return(1);
76                                                  /* get modifiers */
77          raytexture(r, m->omod);
78          if (r->rod < 0.0)                       /* reorient if necessary */
# Line 84 | Line 89 | m_glass(               /* color a ray which hit a thin glass surfac
89                                                  /* angular transmission */
90          cos2 = sqrt( (1.0-1.0/(rindex*rindex)) +
91                       pdot*pdot/(rindex*rindex) );
92 <        setcolor(mcolor, pow(colval(mcolor,RED), 1.0/cos2),
93 <                         pow(colval(mcolor,GRN), 1.0/cos2),
94 <                         pow(colval(mcolor,BLU), 1.0/cos2));
92 >        if (hastrans)
93 >                setcolor(mcolor, pow(colval(mcolor,RED), 1.0/cos2),
94 >                                 pow(colval(mcolor,GRN), 1.0/cos2),
95 >                                 pow(colval(mcolor,BLU), 1.0/cos2));
96  
97                                                  /* compute reflection */
98          r1e = (pdot - rindex*cos2) / (pdot + rindex*cos2);
99          r1e *= r1e;
100          r1m = (1.0/pdot - rindex/cos2) / (1.0/pdot + rindex/cos2);
101          r1m *= r1m;
102 <                                                /* compute transmittance */
103 <        for (i = 0; i < 3; i++) {
104 <                d = colval(mcolor, i);
105 <                colval(trans,i) = .5*(1.0-r1e)*(1.0-r1e)*d/(1.0-r1e*r1e*d*d);
106 <                colval(trans,i) += .5*(1.0-r1m)*(1.0-r1m)*d/(1.0-r1m*r1m*d*d);
107 <        }
108 <        multcolor(trans, r->pcol);              /* modify by pattern */
102 >                                                /* compute transmission */
103 >        if (hastrans) {
104 >                for (i = 0; i < 3; i++) {
105 >                        d = colval(mcolor, i);
106 >                        colval(trans,i) = .5*(1.0-r1e)*(1.0-r1e)*d /
107 >                                                        (1.0-r1e*r1e*d*d);
108 >                        colval(trans,i) += .5*(1.0-r1m)*(1.0-r1m)*d /
109 >                                                        (1.0-r1m*r1m*d*d);
110 >                }
111 >                multcolor(trans, r->pcol);      /* modify by pattern */
112                                                  /* transmitted ray */
113 <        if (rayorigin(&p, TRANS, r, trans) == 0) {
114 <                if (!(r->crtype & SHADOW) && hastexture) {
115 <                        for (i = 0; i < 3; i++)         /* perturb direction */
116 <                                p.rdir[i] = r->rdir[i] +
113 >                if (rayorigin(&p, TRANS, r, trans) == 0) {
114 >                        if (!(r->crtype & SHADOW) && hastexture) {
115 >                                for (i = 0; i < 3; i++) /* perturb direction */
116 >                                        p.rdir[i] = r->rdir[i] +
117                                                  2.*(1.-rindex)*r->pert[i];
118 <                        if (normalize(p.rdir) == 0.0) {
119 <                                objerror(m, WARNING, "bad perturbation");
118 >                                if (normalize(p.rdir) == 0.0) {
119 >                                        objerror(m, WARNING, "bad perturbation");
120 >                                        VCOPY(p.rdir, r->rdir);
121 >                                }
122 >                        } else {
123                                  VCOPY(p.rdir, r->rdir);
124 +                                transtest = 2;
125                          }
126 <                } else {
127 <                        VCOPY(p.rdir, r->rdir);
128 <                        transtest = 2;
126 >                        rayvalue(&p);
127 >                        multcolor(p.rcol, p.rcoef);
128 >                        addcolor(r->rcol, p.rcol);
129 >                        transtest *= bright(p.rcol);
130 >                        transdist = r->rot + p.rt;
131                  }
117 –                rayvalue(&p);
118 –                multcolor(p.rcol, p.rcoef);
119 –                addcolor(r->rcol, p.rcol);
120 –                transtest *= bright(p.rcol);
121 –                transdist = r->rot + p.rt;
132          }
123 –
133          if (r->crtype & SHADOW) {               /* skip reflected ray */
134                  r->rt = transdist;
135                  return(1);

Diff Legend

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