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.13 by gwlarson, Wed Jun 17 12:53:07 1998 UTC vs.
Revision 2.15 by greg, Sat Feb 22 02:07:28 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1998 Silicon Graphics, Inc. */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ SGI";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   *  dielectric.c - shading function for transparent materials.
6 + */
7 +
8 + /* ====================================================================
9 + * The Radiance Software License, Version 1.0
10   *
11 < *     9/6/85
11 > * Copyright (c) 1990 - 2002 The Regents of the University of California,
12 > * through Lawrence Berkeley National Laboratory.   All rights reserved.
13 > *
14 > * Redistribution and use in source and binary forms, with or without
15 > * modification, are permitted provided that the following conditions
16 > * are met:
17 > *
18 > * 1. Redistributions of source code must retain the above copyright
19 > *         notice, this list of conditions and the following disclaimer.
20 > *
21 > * 2. Redistributions in binary form must reproduce the above copyright
22 > *       notice, this list of conditions and the following disclaimer in
23 > *       the documentation and/or other materials provided with the
24 > *       distribution.
25 > *
26 > * 3. The end-user documentation included with the redistribution,
27 > *           if any, must include the following acknowledgment:
28 > *             "This product includes Radiance software
29 > *                 (http://radsite.lbl.gov/)
30 > *                 developed by the Lawrence Berkeley National Laboratory
31 > *               (http://www.lbl.gov/)."
32 > *       Alternately, this acknowledgment may appear in the software itself,
33 > *       if and wherever such third-party acknowledgments normally appear.
34 > *
35 > * 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
36 > *       and "The Regents of the University of California" must
37 > *       not be used to endorse or promote products derived from this
38 > *       software without prior written permission. For written
39 > *       permission, please contact [email protected].
40 > *
41 > * 5. Products derived from this software may not be called "Radiance",
42 > *       nor may "Radiance" appear in their name, without prior written
43 > *       permission of Lawrence Berkeley National Laboratory.
44 > *
45 > * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
46 > * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
47 > * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
48 > * DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
49 > * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
50 > * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
51 > * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
52 > * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
53 > * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
54 > * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
55 > * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
56 > * SUCH DAMAGE.
57 > * ====================================================================
58 > *
59 > * This software consists of voluntary contributions made by many
60 > * individuals on behalf of Lawrence Berkeley National Laboratory.   For more
61 > * information on Lawrence Berkeley National Laboratory, please see
62 > * <http://www.lbl.gov/>.
63   */
64  
65   #include  "ray.h"
# Line 50 | Line 102 | static int  lambda();
102  
103   #define  MINCOS         0.997           /* minimum dot product for dispersion */
104  
53 – extern COLOR  cextinction;              /* global coefficient of extinction */
54 – extern COLOR  salbedo;                  /* global scattering albedo */
105  
56 –
106   static double
107   mylog(x)                /* special log for extinction coefficients */
108   double  x;
# Line 73 | Line 122 | register RAY  *r;
122          double  cos1, cos2, nratio;
123          COLOR  ctrans;
124          COLOR  talb;
125 +        int  hastexture;
126          double  refl, trans;
127          FVECT  dnorm;
128          double  d1, d2;
# Line 84 | Line 134 | register RAY  *r;
134  
135          raytexture(r, m->omod);                 /* get modifiers */
136  
137 <        cos1 = raynormal(dnorm, r);             /* cosine of theta1 */
137 >        if (hastexture = DOT(r->pert,r->pert) > FTINY*FTINY)
138 >                cos1 = raynormal(dnorm, r);     /* perturb normal */
139 >        else {
140 >                VCOPY(dnorm, r->ron);
141 >                cos1 = r->rod;
142 >        }
143                                                  /* index of refraction */
144          if (m->otype == MAT_DIELECTRIC)
145                  nratio = m->oargs.farg[3] + m->oargs.farg[4]/MLAMBDA;
# Line 92 | Line 147 | register RAY  *r;
147                  nratio = m->oargs.farg[3] / m->oargs.farg[7];
148          
149          if (cos1 < 0.0) {                       /* inside */
150 +                hastexture = -hastexture;
151                  cos1 = -cos1;
152                  dnorm[0] = -dnorm[0];
153                  dnorm[1] = -dnorm[1];
# Line 150 | Line 206 | register RAY  *r;
206                  refl *= 0.5;
207                  trans = 1.0 - refl;
208  
209 +                trans *= nratio*nratio;         /* solid angle ratio */
210 +
211                  if (rayorigin(&p, r, REFRACTED, trans) == 0) {
212  
213                                                  /* compute refracted ray */
214                          d1 = nratio*cos1 - cos2;
215                          for (i = 0; i < 3; i++)
216                                  p.rdir[i] = nratio*r->rdir[i] + d1*dnorm[i];
217 <
217 >                                                /* accidental reflection? */
218 >                        if (hastexture &&
219 >                                DOT(p.rdir,r->ron)*hastexture >= -FTINY) {
220 >                                d1 *= (double)hastexture;
221 >                                for (i = 0; i < 3; i++) /* ignore texture */
222 >                                        p.rdir[i] = nratio*r->rdir[i] +
223 >                                                        d1*r->ron[i];
224 >                                normalize(p.rdir);      /* not exact */
225 >                        }
226   #ifdef  DISPERSE
227                          if (m->otype != MAT_DIELECTRIC
228                                          || r->rod > 0.0
# Line 184 | Line 250 | register RAY  *r;
250                                          /* compute reflected ray */
251                  for (i = 0; i < 3; i++)
252                          p.rdir[i] = r->rdir[i] + 2.0*cos1*dnorm[i];
253 +                                        /* accidental penetration? */
254 +                if (hastexture && DOT(p.rdir,r->ron)*hastexture <= FTINY)
255 +                        for (i = 0; i < 3; i++)         /* ignore texture */
256 +                                p.rdir[i] = r->rdir[i] + 2.0*r->rod*r->ron[i];
257  
258                  rayvalue(&p);                   /* reflected ray value */
259  

Diff Legend

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