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

Comparing ray/src/rt/dielectric.c (file contents):
Revision 2.18 by schorsch, Tue Mar 30 16:13:01 2004 UTC vs.
Revision 2.21 by greg, Sun Sep 26 15:51:15 2010 UTC

# Line 74 | Line 74 | m_dielectric(  /* color a ray which hit a dielectric in
74          COLOR  ctrans;
75          COLOR  talb;
76          int  hastexture;
77 +        double  transdist, transtest=0;
78 +        double  mirdist, mirtest=0;
79 +        int     flatsurface;
80          double  refl, trans;
81          FVECT  dnorm;
82          double  d1, d2;
# Line 91 | Line 94 | m_dielectric(  /* color a ray which hit a dielectric in
94                  VCOPY(dnorm, r->ron);
95                  cos1 = r->rod;
96          }
97 +        flatsurface = !hastexture && r->ro != NULL && isflat(r->ro->otype);
98 +
99                                                  /* index of refraction */
100          if (m->otype == MAT_DIELECTRIC)
101                  nratio = m->oargs.farg[3] + m->oargs.farg[4]/MLAMBDA;
# Line 159 | Line 164 | m_dielectric(  /* color a ray which hit a dielectric in
164  
165                  trans *= nratio*nratio;         /* solid angle ratio */
166  
167 <                if (rayorigin(&p, r, REFRACTED, trans) == 0) {
167 >                setcolor(p.rcoef, trans, trans, trans);
168  
169 +                if (rayorigin(&p, REFRACTED, r, p.rcoef) == 0) {
170 +
171                                                  /* compute refracted ray */
172                          d1 = nratio*cos1 - cos2;
173                          for (i = 0; i < 3; i++)
# Line 173 | Line 180 | m_dielectric(  /* color a ray which hit a dielectric in
180                                          p.rdir[i] = nratio*r->rdir[i] +
181                                                          d1*r->ron[i];
182                                  normalize(p.rdir);      /* not exact */
183 <                        }
183 >                        } else
184 >                                checknorm(p.rdir);
185   #ifdef  DISPERSE
186                          if (m->otype != MAT_DIELECTRIC
187                                          || r->rod > 0.0
# Line 187 | Line 195 | m_dielectric(  /* color a ray which hit a dielectric in
195                                  copycolor(p.cext, ctrans);
196                                  copycolor(p.albedo, talb);
197                                  rayvalue(&p);
198 <                                scalecolor(p.rcol, trans);
198 >                                multcolor(p.rcol, p.rcoef);
199                                  addcolor(r->rcol, p.rcol);
200 <                                if (nratio >= 1.0-FTINY && nratio <= 1.0+FTINY)
201 <                                        r->rt = r->rot + p.rt;
200 >                                                /* virtual distance */
201 >                                if (flatsurface ||
202 >                                        (1.-FTINY <= nratio &&
203 >                                                nratio <= 1.+FTINY)) {
204 >                                        transtest = 2*bright(p.rcol);
205 >                                        transdist = r->rot + p.rt;
206 >                                }
207                          }
208                  }
209          }
210 <                
210 >        setcolor(p.rcoef, refl, refl, refl);
211 >
212          if (!(r->crtype & SHADOW) &&
213 <                        rayorigin(&p, r, REFLECTED, refl) == 0) {
213 >                        rayorigin(&p, REFLECTED, r, p.rcoef) == 0) {
214  
215                                          /* compute reflected ray */
216                  for (i = 0; i < 3; i++)
# Line 205 | Line 219 | m_dielectric(  /* color a ray which hit a dielectric in
219                  if (hastexture && DOT(p.rdir,r->ron)*hastexture <= FTINY)
220                          for (i = 0; i < 3; i++)         /* ignore texture */
221                                  p.rdir[i] = r->rdir[i] + 2.0*r->rod*r->ron[i];
222 <
222 >                checknorm(p.rdir);
223                  rayvalue(&p);                   /* reflected ray value */
224  
225 <                scalecolor(p.rcol, refl);       /* color contribution */
225 >                multcolor(p.rcol, p.rcoef);     /* color contribution */
226                  addcolor(r->rcol, p.rcol);
227 +                                                /* virtual distance */
228 +                if (flatsurface) {
229 +                        mirtest = 2*bright(p.rcol);
230 +                        mirdist = r->rot + p.rt;
231 +                }
232          }
233 +                                /* check distance to return */
234 +        d1 = bright(r->rcol);
235 +        if (transtest > d1)
236 +                r->rt = transdist;
237 +        else if (mirtest > d1)
238 +                r->rt = mirdist;
239                                  /* rayvalue() computes absorption */
240          return(1);
241   }
# Line 228 | Line 253 | disperse(  /* check light sources for dispersion */
253          COLOR  abt
254   )
255   {
256 <        RAY  sray, *entray;
256 >        RAY  sray;
257 >        const RAY  *entray;
258          FVECT  v1, v2, n1, n2;
259          FVECT  dv, v2Xdv;
260          double  v2Xdvv2Xdv;

Diff Legend

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