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 2.2 by greg, Mon Nov 25 09:51:06 1991 UTC vs.
Revision 2.11 by schorsch, Sun Jul 27 22:12:03 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 + #include "copyright.h"
10 +
11   #include  "ray.h"
12  
13   #include  "otypes.h"
# Line 32 | Line 31 | static char SCCSid[] = "$SunId$ LBL";
31   */
32  
33  
34 < int  dir_proj();
34 > static int  dir_proj();
35 >
36   VSMATERIAL  direct1_vs = {dir_proj, 1};
37   VSMATERIAL  direct2_vs = {dir_proj, 2};
38  
# Line 41 | Line 41 | VSMATERIAL  direct2_vs = {dir_proj, 2};
41                                  getfunc(m, 8, 0xff, 1) )
42  
43  
44 + int
45   m_direct(m, r)                  /* shade redirected ray */
46   register OBJREC  *m;
47   register RAY  *r;
48   {
49                                          /* check if source ray */
50          if (r->rsrc >= 0 && source[r->rsrc].so != r->ro)
51 <                return;                         /* got the wrong guy */
51 >                return(1);                      /* got the wrong guy */
52                                          /* compute first projection */
53          if (m->otype == MAT_DIRECT1 ||
54                          (r->rsrc < 0 || source[r->rsrc].sa.sv.pn == 0))
# Line 56 | Line 57 | register RAY  *r;
57          if (m->otype == MAT_DIRECT2 &&
58                          (r->rsrc < 0 || source[r->rsrc].sa.sv.pn == 1))
59                  redirect(m, r, 1);
60 +        return(1);
61   }
62  
63  
64 + int
65   redirect(m, r, n)               /* compute n'th ray redirection */
66   OBJREC  *m;
67   RAY  *r;
# Line 66 | Line 69 | int  n;
69   {
70          MFUNC  *mf;
71          register EPNODE  **va;
72 +        FVECT  nsdir;
73          RAY  nr;
74          double  coef;
75          register int  j;
76                                          /* set up function */
77          mf = getdfunc(m);
78          setfunc(m, r);
79 +                                        /* assign direction variable */
80 +        if (r->rsrc >= 0) {
81 +                register SRCREC  *sp = source + source[r->rsrc].sa.sv.sn;
82 +
83 +                if (sp->sflags & SDISTANT)
84 +                        VCOPY(nsdir, sp->sloc);
85 +                else {
86 +                        for (j = 0; j < 3; j++)
87 +                                nsdir[j] = sp->sloc[j] - r->rop[j];
88 +                        normalize(nsdir);
89 +                }
90 +                multv3(nsdir, nsdir, funcxf.xfm);
91 +                varset("DxA", '=', nsdir[0]/funcxf.sca);
92 +                varset("DyA", '=', nsdir[1]/funcxf.sca);
93 +                varset("DzA", '=', nsdir[2]/funcxf.sca);
94 +        } else {
95 +                varset("DxA", '=', 0.0);
96 +                varset("DyA", '=', 0.0);
97 +                varset("DzA", '=', 0.0);
98 +        }
99                                          /* compute coefficient */
100          errno = 0;
101          va = mf->ep + 4*n;
102          coef = evalue(va[0]);
103 <        if (errno)
103 >        if (errno == EDOM || errno == ERANGE)
104                  goto computerr;
105          if (coef <= FTINY || rayorigin(&nr, r, TRANS, coef) < 0)
106                  return(0);
107          va++;                           /* compute direction */
108          for (j = 0; j < 3; j++) {
109                  nr.rdir[j] = evalue(va[j]);
110 <                if (errno)
110 >                if (errno == EDOM || errno == ERANGE)
111                          goto computerr;
112          }
113          if (mf->f != &unitxf)
# Line 98 | Line 122 | int  n;
122          rayvalue(&nr);
123          scalecolor(nr.rcol, coef);
124          addcolor(r->rcol, nr.rcol);
125 +        if (r->ro != NULL && isflat(r->ro->otype))
126 +                r->rt = r->rot + nr.rt;
127          return(1);
128   computerr:
129          objerror(m, WARNING, "compute error");
# Line 105 | Line 131 | computerr:
131   }
132  
133  
134 + static int
135   dir_proj(pm, o, s, n)           /* compute a director's projection */
136   MAT4  pm;
137   OBJREC  *o;
# Line 137 | Line 164 | int  n;
164          olddot = DOT(tr.rdir, nv);
165          if (olddot <= FTINY && olddot >= -FTINY)
166                  return(0);              /* old dir parallels plane */
167 +        tr.rmax = 0.0;
168          rayorigin(&tr, NULL, PRIMARY, 1.0);
169          if (!(*ofun[o->otype].funp)(o, &tr))
170                  return(0);              /* no intersection! */
# Line 144 | Line 172 | int  n;
172          m = vsmaterial(o);
173          mf = getdfunc(m);
174          setfunc(m, &tr);
175 +        varset("DxA", '=', 0.0);
176 +        varset("DyA", '=', 0.0);
177 +        varset("DzA", '=', 0.0);
178          errno = 0;
179          va = mf->ep + 4*n;
180          coef = evalue(va[0]);
181 <        if (errno)
181 >        if (errno == EDOM || errno == ERANGE)
182                  goto computerr;
183          if (coef <= FTINY)
184                  return(0);              /* insignificant */
185          va++;
186          for (i = 0; i < 3; i++) {
187                  newdir[i] = evalue(va[i]);
188 <                if (errno)
188 >                if (errno == EDOM || errno == ERANGE)
189                          goto computerr;
190          }
191          if (mf->f != &unitxf)
# Line 172 | Line 203 | int  n;
203                          pm[i][j] += nv[i]*h[j];
204                  pm[3][j] = -od*h[j];
205          }
206 <        if (newdot > 0.0 ^ olddot > 0.0)        /* add mirroring */
206 >        if ((newdot > 0.0) ^ (olddot > 0.0))    /* add mirroring */
207                  for (j = 0; j < 3; j++) {
208                          for (i = 0; i < 3; i++)
209                                  pm[i][j] -= 2.*nv[i]*nv[j];

Diff Legend

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