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.8 by greg, Mon Mar 8 12:37:23 1993 UTC vs.
Revision 2.20 by greg, Mon Aug 4 22:37:53 2003 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1992 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   *  func.c - interface to calcomp functions.
9 – *
10 – *     4/7/86
6   */
7  
8 + #include "copyright.h"
9 +
10   #include  "ray.h"
11  
12   #include  "otypes.h"
# Line 18 | Line 15 | static char SCCSid[] = "$SunId$ LBL";
15  
16  
17   #define  INITFILE       "rayinit.cal"
21 – #define  REFVNAME       "`FILE_REFCNT"
18   #define  CALSUF         ".cal"
19   #define  LCALSUF        4
20 + char  REFVNAME[] = "`FILE_REFCNT";
21  
22   XF  unitxf = {                  /* identity transform */
23          {{1.0, 0.0, 0.0, 0.0},
# Line 34 | Line 31 | XF  funcxf;                    /* current transformation */
31   static OBJREC  *fobj = NULL;    /* current function object */
32   static RAY  *fray = NULL;       /* current function ray */
33  
34 < static double  l_erf(), l_erfc(), l_arg();
34 > static double  l_erf(char *), l_erfc(char *), l_arg(char *);
35  
36  
37   MFUNC *
38   getfunc(m, ff, ef, dofwd)       /* get function for this modifier */
39   OBJREC  *m;
40   int  ff;
41 < unsigned  ef;
41 > unsigned int  ef;
42   int  dofwd;
43   {
47 –        extern EPNODE  *curfunc;
44          static char  initfile[] = INITFILE;
45          char  sbuf[MAXSTR];
46          register char  **arg;
# Line 56 | Line 52 | int  dofwd;
52                  return(f);
53          fobj = NULL; fray = NULL;
54          if (initfile[0]) {              /* initialize on first call */
55 +                esupport |= E_VARIABLE|E_FUNCTION|E_INCHAN|E_RCONST|E_REDEFW;
56 +                esupport &= ~(E_OUTCHAN);
57                  setcontext("");
58                  scompile("Dx=$1;Dy=$2;Dz=$3;", NULL, 0);
59                  scompile("Nx=$4;Ny=$5;Nz=$6;", NULL, 0);
# Line 65 | Line 63 | int  dofwd;
63                  scompile("Ix=$16;Iy=$17;Iz=$18;", NULL, 0);
64                  scompile("Jx=$19;Jy=$20;Jz=$21;", NULL, 0);
65                  scompile("Kx=$22;Ky=$23;Kz=$24;", NULL, 0);
66 +                scompile("Lu=$26;Lv=$27;", NULL, 0);
67                  funset("arg", 1, '=', l_arg);
68                  funset("erf", 1, ':', l_erf);
69                  funset("erfc", 1, ':', l_erfc);
70                  setnoisefuncs();
71 +                setprismfuncs();
72                  loadfunc(initfile);
73                  initfile[0] = '\0';
74          }
# Line 81 | Line 81 | int  dofwd;
81          if (i == 1 && arg[ff][0] == '.')
82                  setcontext(f->ctx = "");        /* "." means no file */
83          else {
84 <                strcpy(sbuf,arg[ff]);   /* file name is context */
84 >                strcpy(sbuf,arg[ff]);           /* file name is context */
85                  if (i > LCALSUF && !strcmp(sbuf+i-LCALSUF, CALSUF))
86                          sbuf[i-LCALSUF] = '\0'; /* remove suffix */
87                  setcontext(f->ctx = savestr(sbuf));
# Line 134 | Line 134 | memerr:
134   }
135  
136  
137 + void
138   freefunc(m)                     /* free memory associated with modifier */
139   OBJREC  *m;
140   {
# Line 154 | Line 155 | OBJREC  *m;
155                  freestr(f->ctx);
156          }
157          if (f->b != &unitxf)
158 <                free((char *)f->b);
158 >                free((void *)f->b);
159          if (f->f != NULL && f->f != &unitxf)
160 <                free((char *)f->f);
161 <        free((char *)f);
160 >                free((void *)f->f);
161 >        free((void *)f);
162          m->os = NULL;
163   }
164  
165  
166 + int
167   setfunc(m, r)                   /* set channels for function call */
168   OBJREC  *m;
169   register RAY  *r;
170   {
171 <        static long  lastrno = -1;
171 >        static unsigned long  lastrno = ~0;
172          register MFUNC  *f;
173                                          /* get function */
174          if ((f = (MFUNC *)m->os) == NULL)
175 <                error(m, CONSISTENCY, "setfunc called before getfunc");
175 >                objerror(m, CONSISTENCY, "setfunc called before getfunc");
176                                          /* set evaluator context */
177          setcontext(f->ctx);
178                                          /* check to see if matrix set */
# Line 183 | Line 185 | register RAY  *r;
185                          funcxf.sca = r->rox->b.sca * f->b->sca;
186                          multmat4(funcxf.xfm, r->rox->b.xfm, f->b->xfm);
187                  } else
188 <                        copystruct(&funcxf, &r->rox->b);
188 >                        funcxf = r->rox->b;
189          else
190 <                copystruct(&funcxf, f->b);
190 >                funcxf = *(f->b);
191          lastrno = r->rno;
192          eclock++;               /* notify expression evaluator */
193          return(1);
194   }
195  
196  
197 + void
198   loadfunc(fname)                 /* load definition file */
199   char  *fname;
200   {
198 –        extern char  *libpath;          /* library search path */
201          char  *ffname;
202  
203 <        if ((ffname = getpath(fname, libpath, R_OK)) == NULL) {
203 >        if ((ffname = getpath(fname, getrlibpath(), R_OK)) == NULL) {
204                  sprintf(errmsg, "cannot find function file \"%s\"", fname);
205                  error(USER, errmsg);
206          }
# Line 207 | Line 209 | char  *fname;
209  
210  
211   static double
212 < l_arg()                         /* return nth real argument */
212 > l_arg(char *nm)                 /* return nth real argument */
213   {
212 –        extern double  argument();
214          register int  n;
215  
216          if (fobj == NULL)
# Line 229 | Line 230 | l_arg()                                /* return nth real argument */
230  
231  
232   static double
233 < l_erf()                         /* error function */
233 > l_erf(char *nm)                 /* error function */
234   {
235          extern double  erf();
236  
# Line 238 | Line 239 | l_erf()                                /* error function */
239  
240  
241   static double
242 < l_erfc()                        /* cumulative error function */
242 > l_erfc(char *nm)                /* cumulative error function */
243   {
244          extern double  erfc();
245  
# Line 250 | Line 251 | double
251   chanvalue(n)                    /* return channel n to calcomp */
252   register int  n;
253   {
253 –        double  sum;
254 –        register RAY  *r;
255 –
254          if (fray == NULL)
255                  syntax("ray parameter used in constant expression");
256  
257          if (--n < 0)
258                  goto badchan;
259  
260 <        if (n < 3)                      /* ray direction */
260 >        if (n <= 2)                     /* ray direction */
261  
262                  return( (       fray->rdir[0]*funcxf.xfm[0][n] +
263                                  fray->rdir[1]*funcxf.xfm[1][n] +
264                                  fray->rdir[2]*funcxf.xfm[2][n]  )
265                           / funcxf.sca );
266  
267 <        if (n < 6)                      /* surface normal */
267 >        if (n <= 5)                     /* surface normal */
268  
269                  return( (       fray->ron[0]*funcxf.xfm[0][n-3] +
270                                  fray->ron[1]*funcxf.xfm[1][n-3] +
271                                  fray->ron[2]*funcxf.xfm[2][n-3] )
272                           / funcxf.sca );
273  
274 <        if (n < 9)                      /* intersection */
274 >        if (n <= 8)                     /* intersection */
275  
276                  return( fray->rop[0]*funcxf.xfm[0][n-6] +
277                                  fray->rop[1]*funcxf.xfm[1][n-6] +
278                                  fray->rop[2]*funcxf.xfm[2][n-6] +
279                                               funcxf.xfm[3][n-6] );
280  
281 <        if (n == 9) {                   /* total distance */
282 <                sum = fray->rot;
285 <                for (r = fray->parent; r != NULL; r = r->parent)
286 <                        sum += r->rot;
287 <                return(sum * funcxf.sca);
281 >        if (n == 9)                     /* total distance */
282 >                return(raydist(fray,PRIMARY) * funcxf.sca);
283  
289 –        }
290 –
284          if (n == 10)                    /* dot product (range [-1,1]) */
285                  return( fray->rod <= -1.0 ? -1.0 :
286                          fray->rod >= 1.0 ? 1.0 :
# Line 296 | Line 289 | register int  n;
289          if (n == 11)                    /* scale */
290                  return(funcxf.sca);
291  
292 <        if (n < 15)                     /* origin */
292 >        if (n <= 14)                    /* origin */
293                  return(funcxf.xfm[3][n-12]);
294  
295 <        if (n < 18)                     /* i unit vector */
295 >        if (n <= 17)                    /* i unit vector */
296                  return(funcxf.xfm[0][n-15] / funcxf.sca);
297  
298 <        if (n < 21)                     /* j unit vector */
299 <                return(funcxf.xfm[1][n-15] / funcxf.sca);
298 >        if (n <= 20)                    /* j unit vector */
299 >                return(funcxf.xfm[1][n-18] / funcxf.sca);
300  
301 <        if (n < 24)                     /* k unit vector */
301 >        if (n <= 23)                    /* k unit vector */
302                  return(funcxf.xfm[2][n-21] / funcxf.sca);
303  
304 <        if (n == 24) {                  /* single ray (shadow) distance */
305 <                sum = fray->rot;
306 <                for (r = fray->parent; r != NULL && r->crtype&SHADOW;
307 <                                r = r->parent)
308 <                        sum += r->rot;
316 <                return(sum * funcxf.sca);
317 <        }
304 >        if (n == 24)                    /* single ray (shadow) distance */
305 >                return((fray->rot+raydist(fray->parent,SHADOW)) * funcxf.sca);
306 >
307 >        if (n <= 26)                    /* local (u,v) coordinates */
308 >                return(fray->uv[n-25]);
309   badchan:
310          error(USER, "illegal channel number");
311 +        return(0.0);
312   }

Diff Legend

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