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.27 by greg, Thu May 21 13:54:59 2015 UTC

# Line 10 | Line 10 | static const char      RCSid[] = "$Id$";
10   #include  "ray.h"
11   #include  "otypes.h"
12   #include  "rtotypes.h"
13 + #include  "pmapmat.h"
14  
15   #ifdef  DISPERSE
16   #include  "source.h"
# Line 50 | Line 51 | static double mylog(double  x);
51  
52   #define  MINCOS         0.997           /* minimum dot product for dispersion */
53  
53 –
54   static double
55   mylog(          /* special log for extinction coefficients */
56          double  x
# Line 64 | Line 64 | mylog(         /* special log for extinction coefficients */
64   }
65  
66  
67 < extern int
67 > int
68   m_dielectric(   /* color a ray which hit a dielectric interface */
69          OBJREC  *m,
70 <        register RAY  *r
70 >        RAY  *r
71   )
72   {
73          double  cos1, cos2, nratio;
74          COLOR  ctrans;
75          COLOR  talb;
76          int  hastexture;
77 <        double  transdist, transtest=0;
78 <        double  mirdist, mirtest=0;
77 >        double  transdist=0, transtest=0;
78 >        double  mirdist=0, mirtest=0;
79          int     flatsurface;
80          double  refl, trans;
81          FVECT  dnorm;
82          double  d1, d2;
83          RAY  p;
84 <        register int  i;
84 >        int  i;
85  
86 +        /* PMAP: skip refracted shadow or ambient ray if accounted for in
87 +           photon map */
88 +        if (shadowRayInPmap(r))
89 +                return(1);
90 +        
91          if (m->oargs.nfargs != (m->otype==MAT_DIELECTRIC ? 5 : 8))
92                  objerror(m, USER, "bad arguments");
93  
# Line 94 | Line 99 | m_dielectric(  /* color a ray which hit a dielectric in
99                  VCOPY(dnorm, r->ron);
100                  cos1 = r->rod;
101          }
102 <        flatsurface = !hastexture && r->ro != NULL && isflat(r->ro->otype);
102 >        flatsurface = r->ro != NULL && isflat(r->ro->otype) &&
103 >                        !hastexture | (r->crtype & AMBIENT);
104  
105                                                  /* index of refraction */
106          if (m->otype == MAT_DIELECTRIC)
# Line 199 | Line 205 | m_dielectric(  /* color a ray which hit a dielectric in
205                                  addcolor(r->rcol, p.rcol);
206                                                  /* virtual distance */
207                                  if (flatsurface ||
208 <                                        (1.-FTINY <= nratio &&
209 <                                                nratio <= 1.+FTINY)) {
208 >                                        (1.-FTINY <= nratio) &
209 >                                                (nratio <= 1.+FTINY)) {
210                                          transtest = 2*bright(p.rcol);
211                                          transdist = r->rot + p.rt;
212                                  }
# Line 213 | Line 219 | m_dielectric(  /* color a ray which hit a dielectric in
219                          rayorigin(&p, REFLECTED, r, p.rcoef) == 0) {
220  
221                                          /* compute reflected ray */
222 <                for (i = 0; i < 3; i++)
217 <                        p.rdir[i] = r->rdir[i] + 2.0*cos1*dnorm[i];
222 >                VSUM(p.rdir, r->rdir, dnorm, 2.*cos1);
223                                          /* accidental penetration? */
224                  if (hastexture && DOT(p.rdir,r->ron)*hastexture <= FTINY)
225 <                        for (i = 0; i < 3; i++)         /* ignore texture */
221 <                                p.rdir[i] = r->rdir[i] + 2.0*r->rod*r->ron[i];
225 >                        VSUM(p.rdir, r->rdir, r->ron, 2.*r->rod);
226                  checknorm(p.rdir);
227                  rayvalue(&p);                   /* reflected ray value */
228  
# Line 309 | Line 313 | disperse(  /* check light sources for dispersion */
313          VCOPY(n2, r->ron);
314  
315                                          /* first order dispersion approx. */
316 <        dtmp1 = DOT(n1, v1);
317 <        dtmp2 = DOT(n2, v2);
316 >        dtmp1 = 1./DOT(n1, v1);
317 >        dtmp2 = 1./DOT(n2, v2);
318          for (i = 0; i < 3; i++)
319 <                dv[i] = v1[i] + v2[i] - n1[i]/dtmp1 - n2[i]/dtmp2;
319 >                dv[i] = v1[i] + v2[i] - n1[i]*dtmp1 - n2[i]*dtmp2;
320                  
321          if (DOT(dv, dv) <= FTINY)       /* null effect */
322                  return(0);
# Line 355 | Line 359 | disperse(  /* check light sources for dispersion */
359                  dtmp1 = sqrt(si.dom  / v2Xdvv2Xdv / PI);
360  
361                                                          /* compute first ray */
362 <                for (i = 0; i < 3; i++)
359 <                        vtmp2[i] = sray.rdir[i] + dtmp1*vtmp1[i];
362 >                VSUM(vtmp2, sray.rdir, vtmp1, dtmp1);
363  
364                  l1 = lambda(m, v2, dv, vtmp2);          /* first lambda */
365                  if (l1 < 0)
366                          continue;
367                                                          /* compute second ray */
368 <                for (i = 0; i < 3; i++)
366 <                        vtmp2[i] = sray.rdir[i] - dtmp1*vtmp1[i];
368 >                VSUM(vtmp2, sray.rdir, vtmp1, -dtmp1);
369  
370                  l2 = lambda(m, v2, dv, vtmp2);          /* second lambda */
371                  if (l2 < 0)
# Line 384 | Line 386 | disperse(  /* check light sources for dispersion */
386  
387   static int
388   lambda(                 /* compute lambda for material */
389 <        register OBJREC  *m,
389 >        OBJREC  *m,
390          FVECT  v2,
391          FVECT  dv,
392          FVECT  lr
# Line 396 | Line 398 | lambda(                        /* compute lambda for material */
398  
399          fcross(lrXdv, lr, dv);
400          for (i = 0; i < 3; i++)
401 <                if (lrXdv[i] > FTINY || lrXdv[i] < -FTINY)
401 >                if ((lrXdv[i] > FTINY) | (lrXdv[i] < -FTINY))
402                          break;
403          if (i >= 3)
404                  return(-1);

Diff Legend

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