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.4 by greg, Sat Jan 4 19:53:25 1992 UTC vs.
Revision 2.15 by greg, Tue Feb 25 02:47:22 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,
24 <        0.0, 1.0, 0.0, 0.0,
25 <        0.0, 0.0, 1.0, 0.0,
26 <        0.0, 0.0, 0.0, 1.0,
23 >        {{1.0, 0.0, 0.0, 0.0},
24 >        {0.0, 1.0, 0.0, 0.0},
25 >        {0.0, 0.0, 1.0, 0.0},
26 >        {0.0, 0.0, 0.0, 1.0}},
27          1.0
28   };
29  
# Line 41 | Line 38 | 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 54 | Line 50 | int  dofwd;
50                                          /* check to see if done already */
51          if ((f = (MFUNC *)m->os) != NULL)
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);
60                  scompile("Px=$7;Py=$8;Pz=$9;", NULL, 0);
61 <                scompile("T=$10;Rdot=$11;", NULL, 0);
61 >                scompile("T=$10;Ts=$25;Rdot=$11;", NULL, 0);
62                  scompile("S=$12;Tx=$13;Ty=$14;Tz=$15;", NULL, 0);
63                  scompile("Ix=$16;Iy=$17;Iz=$18;", NULL, 0);
64                  scompile("Jx=$19;Jy=$20;Jz=$21;", NULL, 0);
# Line 68 | Line 67 | int  dofwd;
67                  funset("erf", 1, ':', l_erf);
68                  funset("erfc", 1, ':', l_erfc);
69                  setnoisefuncs();
70 +                setprismfuncs();
71                  loadfunc(initfile);
72                  initfile[0] = '\0';
73          }
# Line 80 | Line 80 | int  dofwd;
80          if (i == 1 && arg[ff][0] == '.')
81                  setcontext(f->ctx = "");        /* "." means no file */
82          else {
83 <                strcpy(sbuf,m->oargs.sarg[ff]); /* file name is context */
83 >                strcpy(sbuf,arg[ff]);           /* file name is context */
84                  if (i > LCALSUF && !strcmp(sbuf+i-LCALSUF, CALSUF))
85                          sbuf[i-LCALSUF] = '\0'; /* remove suffix */
86                  setcontext(f->ctx = savestr(sbuf));
87                  if (!vardefined(REFVNAME)) {    /* file loaded? */
88 <                        loadfunc(m->oargs.sarg[ff]);
88 >                        loadfunc(arg[ff]);
89                          varset(REFVNAME, '=', 1.0);
90                  } else                          /* reference_count++ */
91                          varset(REFVNAME, '=', varvalue(REFVNAME)+1.0);
# Line 133 | Line 133 | memerr:
133   }
134  
135  
136 + void
137   freefunc(m)                     /* free memory associated with modifier */
138   OBJREC  *m;
139   {
# Line 153 | Line 154 | OBJREC  *m;
154                  freestr(f->ctx);
155          }
156          if (f->b != &unitxf)
157 <                free((char *)f->b);
157 >                free((void *)f->b);
158          if (f->f != NULL && f->f != &unitxf)
159 <                free((char *)f->f);
160 <        free((char *)f);
159 >                free((void *)f->f);
160 >        free((void *)f);
161          m->os = NULL;
162   }
163  
164  
165 + int
166   setfunc(m, r)                   /* set channels for function call */
167   OBJREC  *m;
168   register RAY  *r;
169   {
170 <        static long  lastrno = -1;
170 >        static unsigned long  lastrno = ~0;
171          register MFUNC  *f;
172                                          /* get function */
173          if ((f = (MFUNC *)m->os) == NULL)
174 <                error(m, CONSISTENCY, "setfunc called before getfunc");
174 >                objerror(m, CONSISTENCY, "setfunc called before getfunc");
175                                          /* set evaluator context */
176          setcontext(f->ctx);
177                                          /* check to see if matrix set */
# Line 191 | Line 193 | register RAY  *r;
193   }
194  
195  
196 + void
197   loadfunc(fname)                 /* load definition file */
198   char  *fname;
199   {
197 –        extern char  *libpath;          /* library search path */
200          char  *ffname;
201  
202 <        if ((ffname = getpath(fname, libpath, R_OK)) == NULL) {
202 >        if ((ffname = getpath(fname, getlibpath(), R_OK)) == NULL) {
203                  sprintf(errmsg, "cannot find function file \"%s\"", fname);
204                  error(USER, errmsg);
205          }
# Line 208 | Line 210 | char  *fname;
210   static double
211   l_arg()                         /* return nth real argument */
212   {
211 –        extern double  argument();
213          register int  n;
214  
215 +        if (fobj == NULL)
216 +                syntax("arg(n) used in constant expression");
217 +
218          n = argument(1) + .5;           /* round to integer */
219  
220          if (n < 1)
# Line 246 | Line 250 | double
250   chanvalue(n)                    /* return channel n to calcomp */
251   register int  n;
252   {
253 <        double  sum;
254 <        register RAY  *r;
253 >        if (fray == NULL)
254 >                syntax("ray parameter used in constant expression");
255  
256          if (--n < 0)
257                  goto badchan;
# Line 273 | Line 277 | register int  n;
277                                  fray->rop[2]*funcxf.xfm[2][n-6] +
278                                               funcxf.xfm[3][n-6] );
279  
280 <        if (n == 9) {                   /* distance */
280 >        if (n == 9)                     /* total distance */
281 >                return(raydist(fray,PRIMARY) * funcxf.sca);
282  
278 –                sum = fray->rot;
279 –                for (r = fray->parent; r != NULL; r = r->parent)
280 –                        sum += r->rot;
281 –                return(sum * funcxf.sca);
282 –
283 –        }
283          if (n == 10)                    /* dot product (range [-1,1]) */
284                  return( fray->rod <= -1.0 ? -1.0 :
285                          fray->rod >= 1.0 ? 1.0 :
# Line 300 | Line 299 | register int  n;
299  
300          if (n < 24)                     /* k unit vector */
301                  return(funcxf.xfm[2][n-21] / funcxf.sca);
302 +
303 +        if (n == 24)                    /* single ray (shadow) distance */
304 +                return((fray->rot+raydist(fray->parent,SHADOW)) * funcxf.sca);
305   badchan:
306          error(USER, "illegal channel number");
307   }

Diff Legend

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