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

Comparing ray/src/rt/glass.c (file contents):
Revision 2.18 by greg, Sun Sep 26 15:51:15 2010 UTC vs.
Revision 2.23 by greg, Tue Feb 24 19:39:26 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   /*
16   *  This definition of glass provides for a quick calculation
# Line 41 | Line 42 | static const char RCSid[] = "$Id$";
42   #define  RINDEX         1.52            /* refractive index of glass */
43  
44  
45 < extern int
45 > int
46   m_glass(                /* color a ray which hit a thin glass surface */
47          OBJREC  *m,
48 <        register RAY  *r
48 >        RAY  *r
49   )
50   {
51          COLOR  mcolor;
52          double  pdot;
53          FVECT  pnorm;
54 <        double  rindex, cos2;
54 >        double  rindex=0, cos2;
55          COLOR  trans, refl;
56          int  hastexture, hastrans;
57          double  d, r1e, r1m;
58          double  transtest, transdist;
59          double  mirtest, mirdist;
60          RAY  p;
61 <        register int  i;
61 >        int  i;
62 >
63 >        /* PMAP: skip refracted shadow ray if accounted for in photon map */
64 >   if (shadowRayInPmap(r))
65 >                return(1);
66                                                  /* check arguments */
67          if (m->oargs.nfargs == 3)
68                  rindex = RINDEX;                /* default value of n */
# Line 65 | Line 70 | m_glass(               /* color a ray which hit a thin glass surfac
70                  rindex = m->oargs.farg[3];      /* use their value */
71          else
72                  objerror(m, USER, "bad arguments");
73 +                                                /* check back face visibility */
74 +        if (!backvis && r->rod <= 0.0) {
75 +                raytrans(r);
76 +                return(1);
77 +        }
78                                                  /* check transmission */
79          setcolor(mcolor, m->oargs.farg[0], m->oargs.farg[1], m->oargs.farg[2]);
80          if ((hastrans = (intens(mcolor) > 1e-15))) {
# Line 111 | Line 121 | m_glass(               /* color a ray which hit a thin glass surfac
121                  multcolor(trans, r->pcol);      /* modify by pattern */
122                                                  /* transmitted ray */
123                  if (rayorigin(&p, TRANS, r, trans) == 0) {
124 <                        if (!(r->crtype & SHADOW) && hastexture) {
125 <                                for (i = 0; i < 3; i++) /* perturb direction */
116 <                                        p.rdir[i] = r->rdir[i] +
117 <                                                2.*(1.-rindex)*r->pert[i];
124 >                        if (!(r->crtype & (SHADOW|AMBIENT)) && hastexture) {
125 >                                VSUM(p.rdir, r->rdir, r->pert, 2.*(1.-rindex));
126                                  if (normalize(p.rdir) == 0.0) {
127                                          objerror(m, WARNING, "bad perturbation");
128                                          VCOPY(p.rdir, r->rdir);
# Line 143 | Line 151 | m_glass(               /* color a ray which hit a thin glass surfac
151          }
152                                                  /* reflected ray */
153          if (rayorigin(&p, REFLECTED, r, refl) == 0) {
154 <                for (i = 0; i < 3; i++)
147 <                        p.rdir[i] = r->rdir[i] + 2.0*pdot*pnorm[i];
154 >                VSUM(p.rdir, r->rdir, pnorm, 2.*pdot);
155                  checknorm(p.rdir);
156                  rayvalue(&p);
157                  multcolor(p.rcol, p.rcoef);
158                  addcolor(r->rcol, p.rcol);
159 <                if (!hastexture && r->ro != NULL && isflat(r->ro->otype)) {
159 >                if (r->ro != NULL && isflat(r->ro->otype) &&
160 >                                !hastexture | (r->crtype & AMBIENT)) {
161                          mirtest = 2.0*bright(p.rcol);
162                          mirdist = r->rot + p.rt;
163                  }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines