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

Comparing ray/src/rt/func.c (file contents):
Revision 2.16 by greg, Tue Mar 11 17:08:55 2003 UTC vs.
Revision 2.21 by greg, Thu Feb 12 18:55:50 2004 UTC

# Line 9 | Line 9 | static const char      RCSid[] = "$Id$";
9  
10   #include  "ray.h"
11  
12 + #include  "paths.h"
13 +
14   #include  "otypes.h"
15  
16   #include  "func.h"
# Line 31 | Line 33 | XF  funcxf;                    /* current transformation */
33   static OBJREC  *fobj = NULL;    /* current function object */
34   static RAY  *fray = NULL;       /* current function ray */
35  
36 < static double  l_erf(), l_erfc(), l_arg();
36 > static double  l_erf(char *), l_erfc(char *), l_arg(char *);
37  
38  
39   MFUNC *
# Line 63 | Line 65 | int  dofwd;
65                  scompile("Ix=$16;Iy=$17;Iz=$18;", NULL, 0);
66                  scompile("Jx=$19;Jy=$20;Jz=$21;", NULL, 0);
67                  scompile("Kx=$22;Ky=$23;Kz=$24;", NULL, 0);
68 <                scompile("Lu=$25;Lv=$26;", NULL, 0);
68 >                scompile("Lu=$26;Lv=$27;", NULL, 0);
69                  funset("arg", 1, '=', l_arg);
70                  funset("erf", 1, ':', l_erf);
71                  funset("erfc", 1, ':', l_erfc);
# Line 185 | Line 187 | register RAY  *r;
187                          funcxf.sca = r->rox->b.sca * f->b->sca;
188                          multmat4(funcxf.xfm, r->rox->b.xfm, f->b->xfm);
189                  } else
190 <                        copystruct(&funcxf, &r->rox->b);
190 >                        funcxf = r->rox->b;
191          else
192 <                copystruct(&funcxf, f->b);
192 >                funcxf = *(f->b);
193          lastrno = r->rno;
194          eclock++;               /* notify expression evaluator */
195          return(1);
# Line 200 | Line 202 | char  *fname;
202   {
203          char  *ffname;
204  
205 <        if ((ffname = getpath(fname, getlibpath(), R_OK)) == NULL) {
205 >        if ((ffname = getpath(fname, getrlibpath(), R_OK)) == NULL) {
206                  sprintf(errmsg, "cannot find function file \"%s\"", fname);
207                  error(USER, errmsg);
208          }
# Line 209 | Line 211 | char  *fname;
211  
212  
213   static double
214 < l_arg()                         /* return nth real argument */
214 > l_arg(char *nm)                 /* return nth real argument */
215   {
216          register int  n;
217  
# Line 230 | Line 232 | l_arg()                                /* return nth real argument */
232  
233  
234   static double
235 < l_erf()                         /* error function */
235 > l_erf(char *nm)                 /* error function */
236   {
237          extern double  erf();
238  
# Line 239 | Line 241 | l_erf()                                /* error function */
241  
242  
243   static double
244 < l_erfc()                        /* cumulative error function */
244 > l_erfc(char *nm)                /* cumulative error function */
245   {
246          extern double  erfc();
247  
# Line 257 | Line 259 | register int  n;
259          if (--n < 0)
260                  goto badchan;
261  
262 <        if (n < 3)                      /* ray direction */
262 >        if (n <= 2)                     /* ray direction */
263  
264                  return( (       fray->rdir[0]*funcxf.xfm[0][n] +
265                                  fray->rdir[1]*funcxf.xfm[1][n] +
266                                  fray->rdir[2]*funcxf.xfm[2][n]  )
267                           / funcxf.sca );
268  
269 <        if (n < 6)                      /* surface normal */
269 >        if (n <= 5)                     /* surface normal */
270  
271                  return( (       fray->ron[0]*funcxf.xfm[0][n-3] +
272                                  fray->ron[1]*funcxf.xfm[1][n-3] +
273                                  fray->ron[2]*funcxf.xfm[2][n-3] )
274                           / funcxf.sca );
275  
276 <        if (n < 9)                      /* intersection */
276 >        if (n <= 8)                     /* intersection */
277  
278                  return( fray->rop[0]*funcxf.xfm[0][n-6] +
279                                  fray->rop[1]*funcxf.xfm[1][n-6] +
# Line 289 | Line 291 | register int  n;
291          if (n == 11)                    /* scale */
292                  return(funcxf.sca);
293  
294 <        if (n < 15)                     /* origin */
294 >        if (n <= 14)                    /* origin */
295                  return(funcxf.xfm[3][n-12]);
296  
297 <        if (n < 18)                     /* i unit vector */
297 >        if (n <= 17)                    /* i unit vector */
298                  return(funcxf.xfm[0][n-15] / funcxf.sca);
299  
300 <        if (n < 21)                     /* j unit vector */
301 <                return(funcxf.xfm[1][n-15] / funcxf.sca);
300 >        if (n <= 20)                    /* j unit vector */
301 >                return(funcxf.xfm[1][n-18] / funcxf.sca);
302  
303 <        if (n < 24)                     /* k unit vector */
303 >        if (n <= 23)                    /* k unit vector */
304                  return(funcxf.xfm[2][n-21] / funcxf.sca);
305  
306 <        if (n < 25)                     /* single ray (shadow) distance */
306 >        if (n == 24)                    /* single ray (shadow) distance */
307                  return((fray->rot+raydist(fray->parent,SHADOW)) * funcxf.sca);
308  
309 <        if (n < 27)                     /* local (u,v) coordinates */
310 <                return(fray->uv[n-26]);
309 >        if (n <= 26)                    /* local (u,v) coordinates */
310 >                return(fray->uv[n-25]);
311   badchan:
312          error(USER, "illegal channel number");
313          return(0.0);

Diff Legend

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