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.23 by greg, Wed Aug 7 05:10:09 2013 UTC vs.
Revision 2.29 by greg, Tue Nov 13 19:58:33 2018 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=0, transtest=0;
78 –        double  mirdist=0, mirtest=0;
77          int     flatsurface;
78          double  refl, trans;
79          FVECT  dnorm;
80          double  d1, d2;
81          RAY  p;
82 <        register int  i;
82 >        int  i;
83  
84 +        /* PMAP: skip refracted shadow or ambient ray if accounted for in
85 +           photon map */
86 +        if (shadowRayInPmap(r) || ambRayInPmap(r))
87 +                return(1);
88 +        
89          if (m->oargs.nfargs != (m->otype==MAT_DIELECTRIC ? 5 : 8))
90                  objerror(m, USER, "bad arguments");
91  
# Line 94 | Line 97 | m_dielectric(  /* color a ray which hit a dielectric in
97                  VCOPY(dnorm, r->ron);
98                  cos1 = r->rod;
99          }
100 <        flatsurface = !hastexture && r->ro != NULL && isflat(r->ro->otype);
100 >        flatsurface = r->ro != NULL && isflat(r->ro->otype) &&
101 >                        !hastexture | (r->crtype & AMBIENT);
102  
103                                                  /* index of refraction */
104          if (m->otype == MAT_DIELECTRIC)
# Line 199 | Line 203 | m_dielectric(  /* color a ray which hit a dielectric in
203                                  addcolor(r->rcol, p.rcol);
204                                                  /* virtual distance */
205                                  if (flatsurface ||
206 <                                        (1.-FTINY <= nratio &&
207 <                                                nratio <= 1.+FTINY)) {
208 <                                        transtest = 2*bright(p.rcol);
205 <                                        transdist = r->rot + p.rt;
206 <                                }
206 >                                        (1.-FTINY <= nratio) &
207 >                                                (nratio <= 1.+FTINY))
208 >                                        r->rxt = r->rot + raydistance(&p);
209                          }
210                  }
211          }
# Line 221 | Line 223 | m_dielectric(  /* color a ray which hit a dielectric in
223                  rayvalue(&p);                   /* reflected ray value */
224  
225                  multcolor(p.rcol, p.rcoef);     /* color contribution */
226 +                copycolor(r->mcol, p.rcol);
227                  addcolor(r->rcol, p.rcol);
228                                                  /* virtual distance */
229 <                if (flatsurface) {
230 <                        mirtest = 2*bright(p.rcol);
228 <                        mirdist = r->rot + p.rt;
229 <                }
229 >                if (flatsurface)
230 >                        r->rmt = r->rot + raydistance(&p);
231          }
231 –                                /* check distance to return */
232 –        d1 = bright(r->rcol);
233 –        if (transtest > d1)
234 –                r->rt = transdist;
235 –        else if (mirtest > d1)
236 –                r->rt = mirdist;
232                                  /* rayvalue() computes absorption */
233          return(1);
234   }
# Line 380 | Line 375 | disperse(  /* check light sources for dispersion */
375  
376   static int
377   lambda(                 /* compute lambda for material */
378 <        register OBJREC  *m,
378 >        OBJREC  *m,
379          FVECT  v2,
380          FVECT  dv,
381          FVECT  lr

Diff Legend

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