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.20 by greg, Mon Aug 4 22:37:53 2003 UTC vs.
Revision 2.36 by greg, Thu Mar 12 17:19:18 2020 UTC

# Line 8 | Line 8 | static const char      RCSid[] = "$Id$";
8   #include "copyright.h"
9  
10   #include  "ray.h"
11 <
11 > #include  "paths.h"
12   #include  "otypes.h"
13 –
13   #include  "func.h"
14 + #include <ctype.h>
15  
16  
17   #define  INITFILE       "rayinit.cal"
# Line 31 | 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 char  rayinitcal[] = INITFILE;
35 +
36   static double  l_erf(char *), l_erfc(char *), l_arg(char *);
37  
38  
39 + void
40 + initfunc(void)  /* initialize function evaluation */
41 + {
42 +        if (!rayinitcal[0])     /* already done? */
43 +                return;
44 +        esupport |= E_VARIABLE|E_FUNCTION|E_INCHAN|E_RCONST|E_REDEFW;
45 +        esupport &= ~(E_OUTCHAN);
46 +        setcontext("");
47 +        scompile("Dx=$1;Dy=$2;Dz=$3;", NULL, 0);
48 +        scompile("Nx=$4;Ny=$5;Nz=$6;", NULL, 0);
49 +        scompile("Px=$7;Py=$8;Pz=$9;", NULL, 0);
50 +        scompile("T=$10;Ts=$25;Rdot=$11;", NULL, 0);
51 +        scompile("S=$12;Tx=$13;Ty=$14;Tz=$15;", NULL, 0);
52 +        scompile("Ix=$16;Iy=$17;Iz=$18;", NULL, 0);
53 +        scompile("Jx=$19;Jy=$20;Jz=$21;", NULL, 0);
54 +        scompile("Kx=$22;Ky=$23;Kz=$24;", NULL, 0);
55 +        scompile("Lu=$26;Lv=$27;", NULL, 0);
56 +        funset("arg", 1, '=', l_arg);
57 +        funset("erf", 1, ':', l_erf);
58 +        funset("erfc", 1, ':', l_erfc);
59 +        setnoisefuncs();
60 +        setprismfuncs();
61 +        loadfunc(rayinitcal);
62 +        rayinitcal[0] = '\0';
63 + }
64 +
65 +
66 + /* Set parameters for current evaluation */
67 + void
68 + set_eparams(char *prms)
69 + {
70 +        static char     *last_params = NULL;
71 +        char            vname[RMAXWORD];
72 +        double          value;
73 +        char            *cpd;
74 +                                        /* check if already set */
75 +        if (prms == NULL || !*prms)
76 +                return;
77 +        if (prms == last_params || (last_params != NULL &&
78 +                                        !strcmp(prms, last_params)))
79 +                return;
80 +        last_params = prms;             /* XXX assumes static string */
81 +                                        /* assign each variable */
82 +        while (*prms) {
83 +                if (isspace(*prms)) {
84 +                        ++prms; continue;
85 +                }
86 +                if (!isalpha(*prms))
87 +                        goto bad_params;
88 +                cpd = vname;
89 +                while (*prms && (*prms != '=') & !isspace(*prms)) {
90 +                        if (!isid(*prms) | (cpd-vname >= RMAXWORD-1))
91 +                                goto bad_params;
92 +                        *cpd++ = *prms++;
93 +                }
94 +                *cpd = '\0';
95 +                while (isspace(*prms)) prms++;
96 +                if (*prms++ != '=')
97 +                        goto bad_params;
98 +                value = atof(prms);
99 +                if ((prms = fskip(prms)) == NULL)
100 +                        goto bad_params;
101 +                while (isspace(*prms)) prms++;
102 +                prms += (*prms == ',') | (*prms == ';');
103 +                varset(vname, '=', value);
104 +        }
105 +        eclock++;               /* notify expression evaluator */
106 +        return;
107 + bad_params:
108 +        sprintf(errmsg, "bad parameter list '%s'", last_params);
109 +        error(USER, errmsg);
110 + }
111 +
112 +
113   MFUNC *
114 < getfunc(m, ff, ef, dofwd)       /* get function for this modifier */
115 < OBJREC  *m;
116 < int  ff;
117 < unsigned int  ef;
118 < int  dofwd;
114 > getfunc(        /* get function for this modifier */
115 >        OBJREC  *m,
116 >        int  ff,
117 >        unsigned int  ef,
118 >        int  dofwd
119 > )
120   {
44 –        static char  initfile[] = INITFILE;
121          char  sbuf[MAXSTR];
122 <        register char  **arg;
123 <        register MFUNC  *f;
122 >        char  **arg;
123 >        MFUNC  *f;
124          int  ne, na;
125 <        register int  i;
125 >        int  i;
126                                          /* check to see if done already */
127          if ((f = (MFUNC *)m->os) != NULL)
128                  return(f);
129          fobj = NULL; fray = NULL;
130 <        if (initfile[0]) {              /* initialize on first call */
131 <                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;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);
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 <        }
130 >        if (rayinitcal[0])              /* initialize on first call */
131 >                initfunc();
132          if ((na = m->oargs.nsargs) <= ff)
133                  goto toofew;
134          arg = m->oargs.sarg;
135          if ((f = (MFUNC *)calloc(1, sizeof(MFUNC))) == NULL)
136                  goto memerr;
137          i = strlen(arg[ff]);                    /* set up context */
138 <        if (i == 1 && arg[ff][0] == '.')
138 >        if (i == 1 && arg[ff][0] == '.') {
139                  setcontext(f->ctx = "");        /* "." means no file */
140 <        else {
140 >        } else {
141                  strcpy(sbuf,arg[ff]);           /* file name is context */
142                  if (i > LCALSUF && !strcmp(sbuf+i-LCALSUF, CALSUF))
143                          sbuf[i-LCALSUF] = '\0'; /* remove suffix */
# Line 108 | Line 165 | int  dofwd;
165                  i = ff+1;
166          while (i < na && arg[i][0] != '-')
167                  i++;
168 <        if (i == na)                    /* no transform */
169 <                f->f = f->b = &unitxf;
170 <        else {                          /* get transform */
171 <                if ((f->b = (XF *)malloc(sizeof(XF))) == NULL)
168 >        if (i == na) {                  /* no transform */
169 >                f->fxp = f->bxp = &unitxf;
170 >        } else {                        /* get transform */
171 >                if ((f->bxp = (XF *)malloc(sizeof(XF))) == NULL)
172                          goto memerr;
173 <                if (invxf(f->b, na-i, arg+i) != na-i)
173 >                if (invxf(f->bxp, na-i, arg+i) != na-i)
174                          objerror(m, USER, "bad transform");
175 <                if (f->b->sca < 0.0)
176 <                        f->b->sca = -f->b->sca;
175 >                if (f->bxp->sca < 0.0)
176 >                        f->bxp->sca = -f->bxp->sca;
177                  if (dofwd) {                    /* do both transforms */
178 <                        if ((f->f = (XF *)malloc(sizeof(XF))) == NULL)
178 >                        if ((f->fxp = (XF *)malloc(sizeof(XF))) == NULL)
179                                  goto memerr;
180 <                        xf(f->f, na-i, arg+i);
181 <                        if (f->f->sca < 0.0)
182 <                                f->f->sca = -f->f->sca;
180 >                        xf(f->fxp, na-i, arg+i);
181 >                        if (f->fxp->sca < 0.0)
182 >                                f->fxp->sca = -f->fxp->sca;
183                  }
184          }
185          m->os = (char *)f;
# Line 131 | Line 188 | toofew:
188          objerror(m, USER, "too few string arguments");
189   memerr:
190          error(SYSTEM, "out of memory in getfunc");
191 +        return NULL; /* pro forma return */
192   }
193  
194  
195   void
196 < freefunc(m)                     /* free memory associated with modifier */
197 < OBJREC  *m;
196 > freefunc(                       /* free memory associated with modifier */
197 >        OBJREC  *m
198 > )
199   {
200 <        register MFUNC  *f;
201 <        register int  i;
200 >        MFUNC  *f;
201 >        int  i;
202  
203          if ((f = (MFUNC *)m->os) == NULL)
204                  return;
# Line 154 | Line 213 | OBJREC  *m;
213                          dcleanup(2);            /* remove definitions */
214                  freestr(f->ctx);
215          }
216 <        if (f->b != &unitxf)
217 <                free((void *)f->b);
218 <        if (f->f != NULL && f->f != &unitxf)
219 <                free((void *)f->f);
216 >        if (f->bxp != &unitxf)
217 >                free((void *)f->bxp);
218 >        if ((f->fxp != NULL) & (f->fxp != &unitxf))
219 >                free((void *)f->fxp);
220          free((void *)f);
221          m->os = NULL;
222   }
223  
224  
225   int
226 < setfunc(m, r)                   /* set channels for function call */
227 < OBJREC  *m;
228 < register RAY  *r;
226 > setfunc(                        /* set channels for function call */
227 >        OBJREC  *m,
228 >        RAY     *r
229 > )
230   {
231 <        static unsigned long  lastrno = ~0;
232 <        register MFUNC  *f;
233 <                                        /* get function */
231 >        static RNUMBER  lastrno = ~0;
232 >        MFUNC           *f;
233 >                                        /* get function if any */
234          if ((f = (MFUNC *)m->os) == NULL)
235                  objerror(m, CONSISTENCY, "setfunc called before getfunc");
236 <                                        /* set evaluator context */
237 <        setcontext(f->ctx);
236 >                
237 >        setcontext(f->ctx);             /* set evaluator context */
238                                          /* check to see if matrix set */
239 <        if (m == fobj && r->rno == lastrno)
239 >        if ((m == fobj) & (r->rno == lastrno))
240                  return(0);
241          fobj = m;
242          fray = r;
243 <        if (r->rox != NULL)
244 <                if (f->b != &unitxf) {
245 <                        funcxf.sca = r->rox->b.sca * f->b->sca;
246 <                        multmat4(funcxf.xfm, r->rox->b.xfm, f->b->xfm);
243 >        if (r->rox != NULL) {
244 >                if (f->bxp != &unitxf) {
245 >                        funcxf.sca = r->rox->b.sca * f->bxp->sca;
246 >                        multmat4(funcxf.xfm, r->rox->b.xfm, f->bxp->xfm);
247                  } else
248                          funcxf = r->rox->b;
249 <        else
250 <                funcxf = *(f->b);
249 >        } else
250 >                funcxf = *f->bxp;
251          lastrno = r->rno;
252          eclock++;               /* notify expression evaluator */
253          return(1);
254   }
255  
256  
257 + int
258 + worldfunc(                      /* special function context sans object */
259 +        char    *ctx,
260 +        RAY     *r
261 + )
262 + {
263 +        static RNUMBER  lastrno = ~0;
264 +
265 +        if (rayinitcal[0])              /* initialize on first call */
266 +                initfunc();
267 +                                        /* set evaluator context */
268 +        setcontext(ctx);
269 +                                        /* check if ray already set */
270 +        if ((fobj == NULL) & (r->rno == lastrno))
271 +                return(0);
272 +        fobj = NULL;
273 +        fray = r;
274 +        funcxf = unitxf;
275 +        lastrno = r->rno;
276 +        eclock++;               /* notify expression evaluator */
277 +        return(1);
278 + }
279 +
280 +
281   void
282 < loadfunc(fname)                 /* load definition file */
283 < char  *fname;
282 > loadfunc(                       /* load definition file */
283 >        char  *fname
284 > )
285   {
286          char  *ffname;
287  
288          if ((ffname = getpath(fname, getrlibpath(), R_OK)) == NULL) {
289                  sprintf(errmsg, "cannot find function file \"%s\"", fname);
290 <                error(USER, errmsg);
290 >                error(SYSTEM, errmsg);
291          }
292          fcompile(ffname);
293   }
# Line 211 | Line 296 | char  *fname;
296   static double
297   l_arg(char *nm)                 /* return nth real argument */
298   {
299 <        register int  n;
299 >        int  n;
300  
301          if (fobj == NULL)
302 <                syntax("arg(n) used in constant expression");
302 >                error(INTERNAL, "arg(n) called without a modifier context");
303  
304          n = argument(1) + .5;           /* round to integer */
305  
# Line 232 | Line 317 | l_arg(char *nm)                        /* return nth real argument */
317   static double
318   l_erf(char *nm)                 /* error function */
319   {
235 –        extern double  erf();
236 –
320          return(erf(argument(1)));
321   }
322  
# Line 241 | Line 324 | l_erf(char *nm)                        /* error function */
324   static double
325   l_erfc(char *nm)                /* cumulative error function */
326   {
244 –        extern double  erfc();
245 –
327          return(erfc(argument(1)));
328   }
329  
330  
331   double
332 < chanvalue(n)                    /* return channel n to calcomp */
333 < register int  n;
332 > chanvalue(                      /* return channel n to calcomp */
333 >        int  n
334 > )
335   {
336          if (fray == NULL)
337                  syntax("ray parameter used in constant expression");
# Line 271 | Line 353 | register int  n;
353                                  fray->ron[2]*funcxf.xfm[2][n-3] )
354                           / funcxf.sca );
355  
356 <        if (n <= 8)                     /* intersection */
356 >        if (n <= 8) {                   /* intersection point */
357 >                if (fray->rot >= FHUGE*.99)
358 >                        return(0.0);    /* XXX should be runtime error? */
359  
360                  return( fray->rop[0]*funcxf.xfm[0][n-6] +
361                                  fray->rop[1]*funcxf.xfm[1][n-6] +
362                                  fray->rop[2]*funcxf.xfm[2][n-6] +
363                                               funcxf.xfm[3][n-6] );
364 +        }
365  
366          if (n == 9)                     /* total distance */
367                  return(raydist(fray,PRIMARY) * funcxf.sca);

Diff Legend

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