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.21 by greg, Sun Sep 26 15:51:15 2010 UTC vs.
Revision 2.22 by greg, Mon Oct 25 22:57:45 2010 UTC

# Line 213 | Line 213 | m_dielectric(  /* color a ray which hit a dielectric in
213                          rayorigin(&p, REFLECTED, r, p.rcoef) == 0) {
214  
215                                          /* compute reflected ray */
216 <                for (i = 0; i < 3; i++)
217 <                        p.rdir[i] = r->rdir[i] + 2.0*cos1*dnorm[i];
216 >                VSUM(p.rdir, r->rdir, dnorm, 2.*cos1);
217                                          /* accidental penetration? */
218                  if (hastexture && DOT(p.rdir,r->ron)*hastexture <= FTINY)
219 <                        for (i = 0; i < 3; i++)         /* ignore texture */
221 <                                p.rdir[i] = r->rdir[i] + 2.0*r->rod*r->ron[i];
219 >                        VSUM(p.rdir, r->rdir, r->ron, 2.*r->rod);
220                  checknorm(p.rdir);
221                  rayvalue(&p);                   /* reflected ray value */
222  
# Line 309 | Line 307 | disperse(  /* check light sources for dispersion */
307          VCOPY(n2, r->ron);
308  
309                                          /* first order dispersion approx. */
310 <        dtmp1 = DOT(n1, v1);
311 <        dtmp2 = DOT(n2, v2);
310 >        dtmp1 = 1./DOT(n1, v1);
311 >        dtmp2 = 1./DOT(n2, v2);
312          for (i = 0; i < 3; i++)
313 <                dv[i] = v1[i] + v2[i] - n1[i]/dtmp1 - n2[i]/dtmp2;
313 >                dv[i] = v1[i] + v2[i] - n1[i]*dtmp1 - n2[i]*dtmp2;
314                  
315          if (DOT(dv, dv) <= FTINY)       /* null effect */
316                  return(0);
# Line 355 | Line 353 | disperse(  /* check light sources for dispersion */
353                  dtmp1 = sqrt(si.dom  / v2Xdvv2Xdv / PI);
354  
355                                                          /* compute first ray */
356 <                for (i = 0; i < 3; i++)
359 <                        vtmp2[i] = sray.rdir[i] + dtmp1*vtmp1[i];
356 >                VSUM(vtmp2, sray.rdir, vtmp1, dtmp1);
357  
358                  l1 = lambda(m, v2, dv, vtmp2);          /* first lambda */
359                  if (l1 < 0)
360                          continue;
361                                                          /* compute second ray */
362 <                for (i = 0; i < 3; i++)
366 <                        vtmp2[i] = sray.rdir[i] - dtmp1*vtmp1[i];
362 >                VSUM(vtmp2, sray.rdir, vtmp1, -dtmp1);
363  
364                  l2 = lambda(m, v2, dv, vtmp2);          /* second lambda */
365                  if (l2 < 0)
# Line 396 | Line 392 | lambda(                        /* compute lambda for material */
392  
393          fcross(lrXdv, lr, dv);
394          for (i = 0; i < 3; i++)
395 <                if (lrXdv[i] > FTINY || lrXdv[i] < -FTINY)
395 >                if ((lrXdv[i] > FTINY) | (lrXdv[i] < -FTINY))
396                          break;
397          if (i >= 3)
398                  return(-1);

Diff Legend

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