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.31 by greg, Wed Nov 15 18:02:52 2023 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;
74 >        COLOR  pcol, ctrans, talb;
75          int  hastexture;
76 +        int     flatsurface;
77          double  refl, trans;
78          FVECT  dnorm;
79          double  d1, d2;
80          RAY  p;
81 <        register int  i;
81 >        int  i;
82  
83 +        /* PMAP: skip refracted shadow or ambient ray if accounted for in
84 +           photon map */
85 +        if (shadowRayInPmap(r) || ambRayInPmap(r))
86 +                return(1);
87 +        
88          if (m->oargs.nfargs != (m->otype==MAT_DIELECTRIC ? 5 : 8))
89                  objerror(m, USER, "bad arguments");
90  
# Line 91 | Line 96 | m_dielectric(  /* color a ray which hit a dielectric in
96                  VCOPY(dnorm, r->ron);
97                  cos1 = r->rod;
98          }
99 +        flatsurface = r->ro != NULL && isflat(r->ro->otype) &&
100 +                        !hastexture | (r->crtype & AMBIENT);
101 +
102                                                  /* index of refraction */
103          if (m->otype == MAT_DIELECTRIC)
104                  nratio = m->oargs.farg[3] + m->oargs.farg[4]/MLAMBDA;
105          else
106                  nratio = m->oargs.farg[3] / m->oargs.farg[7];
107 <        
107 >
108 >        scolor_rgb(pcol, r->pcol);
109          if (cos1 < 0.0) {                       /* inside */
110                  hastexture = -hastexture;
111                  cos1 = -cos1;
112                  dnorm[0] = -dnorm[0];
113                  dnorm[1] = -dnorm[1];
114                  dnorm[2] = -dnorm[2];
115 <                setcolor(r->cext, -mylog(m->oargs.farg[0]*colval(r->pcol,RED)),
116 <                                 -mylog(m->oargs.farg[1]*colval(r->pcol,GRN)),
117 <                                 -mylog(m->oargs.farg[2]*colval(r->pcol,BLU)));
115 >                setcolor(r->cext, -mylog(m->oargs.farg[0]*colval(pcol,RED)),
116 >                                 -mylog(m->oargs.farg[1]*colval(pcol,GRN)),
117 >                                 -mylog(m->oargs.farg[2]*colval(pcol,BLU)));
118                  setcolor(r->albedo, 0., 0., 0.);
119                  r->gecc = 0.;
120                  if (m->otype == MAT_INTERFACE) {
121                          setcolor(ctrans,
122 <                                -mylog(m->oargs.farg[4]*colval(r->pcol,RED)),
123 <                                -mylog(m->oargs.farg[5]*colval(r->pcol,GRN)),
124 <                                -mylog(m->oargs.farg[6]*colval(r->pcol,BLU)));
122 >                                -mylog(m->oargs.farg[4]*colval(pcol,RED)),
123 >                                -mylog(m->oargs.farg[5]*colval(pcol,GRN)),
124 >                                -mylog(m->oargs.farg[6]*colval(pcol,BLU)));
125                          setcolor(talb, 0., 0., 0.);
126                  } else {
127                          copycolor(ctrans, cextinction);
# Line 121 | Line 130 | m_dielectric(  /* color a ray which hit a dielectric in
130          } else {                                /* outside */
131                  nratio = 1.0 / nratio;
132  
133 <                setcolor(ctrans, -mylog(m->oargs.farg[0]*colval(r->pcol,RED)),
134 <                                 -mylog(m->oargs.farg[1]*colval(r->pcol,GRN)),
135 <                                 -mylog(m->oargs.farg[2]*colval(r->pcol,BLU)));
133 >                setcolor(ctrans, -mylog(m->oargs.farg[0]*colval(pcol,RED)),
134 >                                 -mylog(m->oargs.farg[1]*colval(pcol,GRN)),
135 >                                 -mylog(m->oargs.farg[2]*colval(pcol,BLU)));
136                  setcolor(talb, 0., 0., 0.);
137                  if (m->otype == MAT_INTERFACE) {
138                          setcolor(r->cext,
139 <                                -mylog(m->oargs.farg[4]*colval(r->pcol,RED)),
140 <                                -mylog(m->oargs.farg[5]*colval(r->pcol,GRN)),
141 <                                -mylog(m->oargs.farg[6]*colval(r->pcol,BLU)));
139 >                                -mylog(m->oargs.farg[4]*colval(pcol,RED)),
140 >                                -mylog(m->oargs.farg[5]*colval(pcol,GRN)),
141 >                                -mylog(m->oargs.farg[6]*colval(pcol,BLU)));
142                          setcolor(r->albedo, 0., 0., 0.);
143                          r->gecc = 0.;
144                  }
# Line 159 | Line 168 | m_dielectric(  /* color a ray which hit a dielectric in
168  
169                  trans *= nratio*nratio;         /* solid angle ratio */
170  
171 <                if (rayorigin(&p, r, REFRACTED, trans) == 0) {
171 >                setscolor(p.rcoef, trans, trans, trans);
172  
173 +                if (rayorigin(&p, REFRACTED, r, p.rcoef) == 0) {
174 +
175                                                  /* compute refracted ray */
176                          d1 = nratio*cos1 - cos2;
177                          for (i = 0; i < 3; i++)
# Line 173 | Line 184 | m_dielectric(  /* color a ray which hit a dielectric in
184                                          p.rdir[i] = nratio*r->rdir[i] +
185                                                          d1*r->ron[i];
186                                  normalize(p.rdir);      /* not exact */
187 <                        }
187 >                        } else
188 >                                checknorm(p.rdir);
189   #ifdef  DISPERSE
190                          if (m->otype != MAT_DIELECTRIC
191                                          || r->rod > 0.0
# Line 187 | Line 199 | m_dielectric(  /* color a ray which hit a dielectric in
199                                  copycolor(p.cext, ctrans);
200                                  copycolor(p.albedo, talb);
201                                  rayvalue(&p);
202 <                                scalecolor(p.rcol, trans);
203 <                                addcolor(r->rcol, p.rcol);
204 <                                if (nratio >= 1.0-FTINY && nratio <= 1.0+FTINY)
205 <                                        r->rt = r->rot + p.rt;
202 >                                smultscolor(p.rcol, p.rcoef);
203 >                                saddscolor(r->rcol, p.rcol);
204 >                                                /* virtual distance */
205 >                                if (flatsurface ||
206 >                                        (1.-FTINY <= nratio) &
207 >                                                (nratio <= 1.+FTINY))
208 >                                        r->rxt = r->rot + raydistance(&p);
209                          }
210                  }
211          }
212 <                
212 >        setscolor(p.rcoef, refl, refl, refl);
213 >
214          if (!(r->crtype & SHADOW) &&
215 <                        rayorigin(&p, r, REFLECTED, refl) == 0) {
215 >                        rayorigin(&p, REFLECTED, r, p.rcoef) == 0) {
216  
217                                          /* compute reflected ray */
218 <                for (i = 0; i < 3; i++)
203 <                        p.rdir[i] = r->rdir[i] + 2.0*cos1*dnorm[i];
218 >                VSUM(p.rdir, r->rdir, dnorm, 2.*cos1);
219                                          /* accidental penetration? */
220                  if (hastexture && DOT(p.rdir,r->ron)*hastexture <= FTINY)
221 <                        for (i = 0; i < 3; i++)         /* ignore texture */
222 <                                p.rdir[i] = r->rdir[i] + 2.0*r->rod*r->ron[i];
208 <
221 >                        VSUM(p.rdir, r->rdir, r->ron, 2.*r->rod);
222 >                checknorm(p.rdir);
223                  rayvalue(&p);                   /* reflected ray value */
224  
225 <                scalecolor(p.rcol, refl);       /* color contribution */
226 <                addcolor(r->rcol, p.rcol);
225 >                smultscolor(p.rcol, p.rcoef);   /* color contribution */
226 >                copyscolor(r->mcol, p.rcol);
227 >                saddscolor(r->rcol, p.rcol);
228 >                                                /* virtual distance */
229 >                r->rmt = r->rot;
230 >                if (flatsurface)
231 >                        r->rmt += raydistance(&p);
232          }
233                                  /* rayvalue() computes absorption */
234          return(1);
# Line 228 | Line 247 | disperse(  /* check light sources for dispersion */
247          COLOR  abt
248   )
249   {
250 <        RAY  sray, *entray;
250 >        RAY  sray;
251 >        const RAY  *entray;
252          FVECT  v1, v2, n1, n2;
253          FVECT  dv, v2Xdv;
254          double  v2Xdvv2Xdv;
# Line 283 | Line 303 | disperse(  /* check light sources for dispersion */
303          VCOPY(n2, r->ron);
304  
305                                          /* first order dispersion approx. */
306 <        dtmp1 = DOT(n1, v1);
307 <        dtmp2 = DOT(n2, v2);
306 >        dtmp1 = 1./DOT(n1, v1);
307 >        dtmp2 = 1./DOT(n2, v2);
308          for (i = 0; i < 3; i++)
309 <                dv[i] = v1[i] + v2[i] - n1[i]/dtmp1 - n2[i]/dtmp2;
309 >                dv[i] = v1[i] + v2[i] - n1[i]*dtmp1 - n2[i]*dtmp2;
310                  
311          if (DOT(dv, dv) <= FTINY)       /* null effect */
312                  return(0);
# Line 329 | Line 349 | disperse(  /* check light sources for dispersion */
349                  dtmp1 = sqrt(si.dom  / v2Xdvv2Xdv / PI);
350  
351                                                          /* compute first ray */
352 <                for (i = 0; i < 3; i++)
333 <                        vtmp2[i] = sray.rdir[i] + dtmp1*vtmp1[i];
352 >                VSUM(vtmp2, sray.rdir, vtmp1, dtmp1);
353  
354                  l1 = lambda(m, v2, dv, vtmp2);          /* first lambda */
355                  if (l1 < 0)
356                          continue;
357                                                          /* compute second ray */
358 <                for (i = 0; i < 3; i++)
340 <                        vtmp2[i] = sray.rdir[i] - dtmp1*vtmp1[i];
358 >                VSUM(vtmp2, sray.rdir, vtmp1, -dtmp1);
359  
360                  l2 = lambda(m, v2, dv, vtmp2);          /* second lambda */
361                  if (l2 < 0)
362                          continue;
363                                          /* compute color from spectrum */
364 <                if (l1 < l2)
364 >                if (l1 < l2)            /* XXX should use direct spectral xfer */
365                          spec_rgb(ctmp, l1, l2);
366                  else
367                          spec_rgb(ctmp, l2, l1);
368 <                multcolor(ctmp, sray.rcol);
368 >                multscolor(ctmp, sray.rcol);
369                  scalecolor(ctmp, tr);
370 <                addcolor(r->rcol, ctmp);
370 >                saddcolor(r->rcol, ctmp);
371                  success++;
372          }
373          return(success);
# Line 358 | Line 376 | disperse(  /* check light sources for dispersion */
376  
377   static int
378   lambda(                 /* compute lambda for material */
379 <        register OBJREC  *m,
379 >        OBJREC  *m,
380          FVECT  v2,
381          FVECT  dv,
382          FVECT  lr
# Line 370 | Line 388 | lambda(                        /* compute lambda for material */
388  
389          fcross(lrXdv, lr, dv);
390          for (i = 0; i < 3; i++)
391 <                if (lrXdv[i] > FTINY || lrXdv[i] < -FTINY)
391 >                if ((lrXdv[i] > FTINY) | (lrXdv[i] < -FTINY))
392                          break;
393          if (i >= 3)
394                  return(-1);

Diff Legend

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