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

Comparing ray/src/rt/m_direct.c (file contents):
Revision 1.4 by greg, Tue Aug 6 14:07:06 1991 UTC vs.
Revision 2.8 by greg, Sat Feb 22 02:07:28 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1991 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   * Routines for light-redirecting materials and
6   *   their associated virtual light sources
7   */
8  
9 + /* ====================================================================
10 + * The Radiance Software License, Version 1.0
11 + *
12 + * Copyright (c) 1990 - 2002 The Regents of the University of California,
13 + * through Lawrence Berkeley National Laboratory.   All rights reserved.
14 + *
15 + * Redistribution and use in source and binary forms, with or without
16 + * modification, are permitted provided that the following conditions
17 + * are met:
18 + *
19 + * 1. Redistributions of source code must retain the above copyright
20 + *         notice, this list of conditions and the following disclaimer.
21 + *
22 + * 2. Redistributions in binary form must reproduce the above copyright
23 + *       notice, this list of conditions and the following disclaimer in
24 + *       the documentation and/or other materials provided with the
25 + *       distribution.
26 + *
27 + * 3. The end-user documentation included with the redistribution,
28 + *           if any, must include the following acknowledgment:
29 + *             "This product includes Radiance software
30 + *                 (http://radsite.lbl.gov/)
31 + *                 developed by the Lawrence Berkeley National Laboratory
32 + *               (http://www.lbl.gov/)."
33 + *       Alternately, this acknowledgment may appear in the software itself,
34 + *       if and wherever such third-party acknowledgments normally appear.
35 + *
36 + * 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
37 + *       and "The Regents of the University of California" must
38 + *       not be used to endorse or promote products derived from this
39 + *       software without prior written permission. For written
40 + *       permission, please contact [email protected].
41 + *
42 + * 5. Products derived from this software may not be called "Radiance",
43 + *       nor may "Radiance" appear in their name, without prior written
44 + *       permission of Lawrence Berkeley National Laboratory.
45 + *
46 + * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
47 + * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
48 + * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
49 + * DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
50 + * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
51 + * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
52 + * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
53 + * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
54 + * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
55 + * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
56 + * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
57 + * SUCH DAMAGE.
58 + * ====================================================================
59 + *
60 + * This software consists of voluntary contributions made by many
61 + * individuals on behalf of Lawrence Berkeley National Laboratory.   For more
62 + * information on Lawrence Berkeley National Laboratory, please see
63 + * <http://www.lbl.gov/>.
64 + */
65 +
66   #include  "ray.h"
67  
68   #include  "otypes.h"
69  
70   #include  "source.h"
71  
72 + #include  "func.h"
73 +
74   /*
75   * The arguments for MAT_DIRECT1 are:
76   *
# Line 30 | Line 86 | static char SCCSid[] = "$SunId$ LBL";
86   */
87  
88  
89 < extern double  varvalue();
89 > static int  dir_proj();
90  
35 – int  dir_proj();
91   VSMATERIAL  direct1_vs = {dir_proj, 1};
92   VSMATERIAL  direct2_vs = {dir_proj, 2};
93  
94 + #define getdfunc(m)     ( (m)->otype == MAT_DIRECT1 ? \
95 +                                getfunc(m, 4, 0xf, 1) : \
96 +                                getfunc(m, 8, 0xff, 1) )
97  
98 +
99 + int
100   m_direct(m, r)                  /* shade redirected ray */
101   register OBJREC  *m;
102   register RAY  *r;
103   {
104                                          /* check if source ray */
105          if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
106 <                return;                         /* got the wrong guy */
47 <        dir_check(m);
106 >                return(1);                      /* got the wrong guy */
107                                          /* compute first projection */
108          if (m->otype == MAT_DIRECT1 ||
109                          (r->rsrc < 0 || source[r->rsrc].sa.sv.pn == 0))
# Line 53 | Line 112 | register RAY  *r;
112          if (m->otype == MAT_DIRECT2 &&
113                          (r->rsrc < 0 || source[r->rsrc].sa.sv.pn == 1))
114                  redirect(m, r, 1);
115 +        return(1);
116   }
117  
118  
119 + int
120   redirect(m, r, n)               /* compute n'th ray redirection */
121   OBJREC  *m;
122   RAY  *r;
123   int  n;
124   {
125 <        register char  **sa;
125 >        MFUNC  *mf;
126 >        register EPNODE  **va;
127 >        FVECT  nsdir;
128          RAY  nr;
129          double  coef;
130          register int  j;
131                                          /* set up function */
132 <        setmap(m, r, &((FULLXF *)m->os)->b);
133 <        sa = m->oargs.sarg + 4*n;
132 >        mf = getdfunc(m);
133 >        setfunc(m, r);
134 >                                        /* assign direction variable */
135 >        if (r->rsrc >= 0) {
136 >                register SRCREC  *sp = source + source[r->rsrc].sa.sv.sn;
137 >
138 >                if (sp->sflags & SDISTANT)
139 >                        VCOPY(nsdir, sp->sloc);
140 >                else {
141 >                        for (j = 0; j < 3; j++)
142 >                                nsdir[j] = sp->sloc[j] - r->rop[j];
143 >                        normalize(nsdir);
144 >                }
145 >                multv3(nsdir, nsdir, funcxf.xfm);
146 >                varset("DxA", '=', nsdir[0]/funcxf.sca);
147 >                varset("DyA", '=', nsdir[1]/funcxf.sca);
148 >                varset("DzA", '=', nsdir[2]/funcxf.sca);
149 >        } else {
150 >                varset("DxA", '=', 0.0);
151 >                varset("DyA", '=', 0.0);
152 >                varset("DzA", '=', 0.0);
153 >        }
154                                          /* compute coefficient */
155          errno = 0;
156 <        coef = varvalue(sa[0]);
156 >        va = mf->ep + 4*n;
157 >        coef = evalue(va[0]);
158          if (errno)
159                  goto computerr;
160          if (coef <= FTINY || rayorigin(&nr, r, TRANS, coef) < 0)
161                  return(0);
162 <                                        /* compute direction */
163 <        errno = 0;
164 <        for (j = 0; j < 3; j++)
165 <                nr.rdir[j] = varvalue(sa[j+1]);
166 <        if (errno)
167 <                goto computerr;
168 <        multv3(nr.rdir, nr.rdir, ((FULLXF *)m->os)->f.xfm);
162 >        va++;                           /* compute direction */
163 >        for (j = 0; j < 3; j++) {
164 >                nr.rdir[j] = evalue(va[j]);
165 >                if (errno)
166 >                        goto computerr;
167 >        }
168 >        if (mf->f != &unitxf)
169 >                multv3(nr.rdir, nr.rdir, mf->f->xfm);
170          if (r->rox != NULL)
171                  multv3(nr.rdir, nr.rdir, r->rox->f.xfm);
172          if (normalize(nr.rdir) == 0.0)
# Line 92 | Line 177 | int  n;
177          rayvalue(&nr);
178          scalecolor(nr.rcol, coef);
179          addcolor(r->rcol, nr.rcol);
180 +        if (r->ro != NULL && isflat(r->ro->otype))
181 +                r->rt = r->rot + nr.rt;
182          return(1);
183   computerr:
184          objerror(m, WARNING, "compute error");
# Line 99 | Line 186 | computerr:
186   }
187  
188  
189 + static int
190   dir_proj(pm, o, s, n)           /* compute a director's projection */
191   MAT4  pm;
192   OBJREC  *o;
# Line 106 | Line 194 | SRCREC  *s;
194   int  n;
195   {
196          RAY  tr;
197 <        register OBJREC  *m;
198 <        char  **sa;
197 >        OBJREC  *m;
198 >        MFUNC  *mf;
199 >        EPNODE  **va;
200          FVECT  cent, newdir, nv, h;
201 <        double  olddot, newdot, od;
201 >        double  coef, olddot, newdot, od;
202          register int  i, j;
114 –                                /* get material arguments */
115 –        m = objptr(o->omod);
116 –        dir_check(m);
117 –        sa = m->oargs.sarg + 4*n;
203                                  /* initialize test ray */
204          getmaxdisk(cent, o);
205          if (s->sflags & SDISTANT)
# Line 134 | Line 219 | int  n;
219          olddot = DOT(tr.rdir, nv);
220          if (olddot <= FTINY && olddot >= -FTINY)
221                  return(0);              /* old dir parallels plane */
222 +        tr.rmax = 0.0;
223          rayorigin(&tr, NULL, PRIMARY, 1.0);
224          if (!(*ofun[o->otype].funp)(o, &tr))
225                  return(0);              /* no intersection! */
226                                  /* compute redirection */
227 <        setmap(m, &tr, &((FULLXF *)m->os)->b);
227 >        m = vsmaterial(o);
228 >        mf = getdfunc(m);
229 >        setfunc(m, &tr);
230 >        varset("DxA", '=', 0.0);
231 >        varset("DyA", '=', 0.0);
232 >        varset("DzA", '=', 0.0);
233          errno = 0;
234 <        if (varvalue(sa[0]) <= FTINY)
235 <                return(0);              /* insignificant */
234 >        va = mf->ep + 4*n;
235 >        coef = evalue(va[0]);
236          if (errno)
237                  goto computerr;
238 <        for (i = 0; i < 3; i++)
239 <                newdir[i] = varvalue(sa[i+1]);
240 <        if (errno)
241 <                goto computerr;
242 <        multv3(newdir, newdir, ((FULLXF *)m->os)->f.xfm);
238 >        if (coef <= FTINY)
239 >                return(0);              /* insignificant */
240 >        va++;
241 >        for (i = 0; i < 3; i++) {
242 >                newdir[i] = evalue(va[i]);
243 >                if (errno)
244 >                        goto computerr;
245 >        }
246 >        if (mf->f != &unitxf)
247 >                multv3(newdir, newdir, mf->f->xfm);
248                                          /* normalization unnecessary */
249          newdot = DOT(newdir, nv);
250          if (newdot <= FTINY && newdot >= -FTINY)
# Line 172 | Line 268 | int  n;
268   computerr:
269          objerror(m, WARNING, "projection compute error");
270          return(0);
175 – }
176 –
177 –
178 – static
179 – dir_check(m)                    /* check arguments and load function file */
180 – register OBJREC  *m;
181 – {
182 –        register FULLXF  *mxf;
183 –        register int  ff;
184 –
185 –        ff = m->otype == MAT_DIRECT1 ? 5 : 9;
186 –        if (ff > m->oargs.nsargs)
187 –                objerror(m, USER, "too few arguments");
188 –        if (!vardefined(m->oargs.sarg[0]))
189 –                loadfunc(m->oargs.sarg[ff-1]);
190 –        if (m->os == NULL) {
191 –                mxf = (FULLXF *)malloc(sizeof(FULLXF));
192 –                if (mxf == NULL)
193 –                        error(SYSTEM, "out of memory in dir_check");
194 –                if (fullxf(mxf, m->oargs.nsargs-ff, m->oargs.sarg+ff) !=
195 –                                m->oargs.nsargs-ff)
196 –                        objerror(m, USER, "bad transform");
197 –                if (mxf->f.sca < 0.0)
198 –                        mxf->f.sca = -mxf->f.sca;
199 –                if (mxf->b.sca < 0.0)
200 –                        mxf->b.sca = -mxf->b.sca;
201 –                m->os = (char *)mxf;
202 –        }
271   }

Diff Legend

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