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.25 by greg, Thu Jun 7 18:56:06 2012 UTC vs.
Revision 2.34 by greg, Sat May 13 00:26:27 2017 UTC

# Line 11 | Line 11 | static const char      RCSid[] = "$Id$";
11   #include  "paths.h"
12   #include  "otypes.h"
13   #include  "func.h"
14 + #include <ctype.h>
15  
16  
17   #define  INITFILE       "rayinit.cal"
# Line 30 | 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;
71 +        static int      lplen = 0;
72 +        int             len;
73 +        char            vname[RMAXWORD];
74 +        double          value;
75 +        char            *cpd;
76 +                                        /* check if already set */
77 +        if (prms == NULL || !*prms)
78 +                return;
79 +        if (lplen && !strcmp(prms, last_params))
80 +                return;
81 +        len = strlen(prms);             /* record new settings */
82 +        if ((lplen != 0) & (lplen <= len)) {
83 +                free(last_params);
84 +                lplen = 0;
85 +        }
86 +        if (!lplen) {
87 +                lplen = len + 100;
88 +                last_params = (char *)malloc(lplen);
89 +                if (last_params == NULL)
90 +                        error(SYSTEM, "out of memory in set_eparams()");
91 +        }
92 +        strcpy(last_params, prms);
93 +                                        /* assign each variable */
94 +        while (*prms) {
95 +                if (isspace(*prms)) {
96 +                        ++prms; continue;
97 +                }
98 +                if (!isalpha(*prms))
99 +                        goto bad_params;
100 +                cpd = vname;
101 +                while (*prms && (*prms != '=') & !isspace(*prms)) {
102 +                        if (!isid(*prms) | (cpd-vname >= RMAXWORD-1))
103 +                                goto bad_params;
104 +                        *cpd++ = *prms++;
105 +                }
106 +                *cpd = '\0';
107 +                while (isspace(*prms)) prms++;
108 +                if (*prms++ != '=')
109 +                        goto bad_params;
110 +                value = atof(prms);
111 +                if ((prms = fskip(prms)) == NULL)
112 +                        goto bad_params;
113 +                while (isspace(*prms)) prms++;
114 +                prms += (*prms == ',') | (*prms == ';');
115 +                varset(vname, '=', value);
116 +        }
117 +        eclock++;               /* notify expression evaluator */
118 +        return;
119 + bad_params:
120 +        sprintf(errmsg, "bad parameter list '%s'", last_params);
121 +        error(USER, errmsg);
122 + }
123 +
124 +
125   MFUNC *
126   getfunc(        /* get function for this modifier */
127          OBJREC  *m,
# Line 41 | Line 130 | getfunc(       /* get function for this modifier */
130          int  dofwd
131   )
132   {
44 –        static char  initfile[] = INITFILE;
133          char  sbuf[MAXSTR];
134          char  **arg;
135          MFUNC  *f;
# Line 51 | Line 139 | getfunc(       /* get function for this modifier */
139          if ((f = (MFUNC *)m->os) != NULL)
140                  return(f);
141          fobj = NULL; fray = NULL;
142 <        if (initfile[0]) {              /* initialize on first call */
143 <                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 <        }
142 >        if (rayinitcal[0])              /* initialize on first call */
143 >                initfunc();
144          if ((na = m->oargs.nsargs) <= ff)
145                  goto toofew;
146          arg = m->oargs.sarg;
147          if ((f = (MFUNC *)calloc(1, sizeof(MFUNC))) == NULL)
148                  goto memerr;
149          i = strlen(arg[ff]);                    /* set up context */
150 <        if (i == 1 && arg[ff][0] == '.')
150 >        if (i == 1 && arg[ff][0] == '.') {
151                  setcontext(f->ctx = "");        /* "." means no file */
152 <        else {
152 >        } else {
153                  strcpy(sbuf,arg[ff]);           /* file name is context */
154                  if (i > LCALSUF && !strcmp(sbuf+i-LCALSUF, CALSUF))
155                          sbuf[i-LCALSUF] = '\0'; /* remove suffix */
# Line 108 | Line 177 | getfunc(       /* get function for this modifier */
177                  i = ff+1;
178          while (i < na && arg[i][0] != '-')
179                  i++;
180 <        if (i == na)                    /* no transform */
181 <                f->f = f->b = &unitxf;
182 <        else {                          /* get transform */
183 <                if ((f->b = (XF *)malloc(sizeof(XF))) == NULL)
180 >        if (i == na) {                  /* no transform */
181 >                f->fxp = f->bxp = &unitxf;
182 >        } else {                        /* get transform */
183 >                if ((f->bxp = (XF *)malloc(sizeof(XF))) == NULL)
184                          goto memerr;
185 <                if (invxf(f->b, na-i, arg+i) != na-i)
185 >                if (invxf(f->bxp, na-i, arg+i) != na-i)
186                          objerror(m, USER, "bad transform");
187 <                if (f->b->sca < 0.0)
188 <                        f->b->sca = -f->b->sca;
187 >                if (f->bxp->sca < 0.0)
188 >                        f->bxp->sca = -f->bxp->sca;
189                  if (dofwd) {                    /* do both transforms */
190 <                        if ((f->f = (XF *)malloc(sizeof(XF))) == NULL)
190 >                        if ((f->fxp = (XF *)malloc(sizeof(XF))) == NULL)
191                                  goto memerr;
192 <                        xf(f->f, na-i, arg+i);
193 <                        if (f->f->sca < 0.0)
194 <                                f->f->sca = -f->f->sca;
192 >                        xf(f->fxp, na-i, arg+i);
193 >                        if (f->fxp->sca < 0.0)
194 >                                f->fxp->sca = -f->fxp->sca;
195                  }
196          }
197          m->os = (char *)f;
# Line 156 | Line 225 | freefunc(                      /* free memory associated with modifier */
225                          dcleanup(2);            /* remove definitions */
226                  freestr(f->ctx);
227          }
228 <        if (f->b != &unitxf)
229 <                free((void *)f->b);
230 <        if ((f->f != NULL) & (f->f != &unitxf))
231 <                free((void *)f->f);
228 >        if (f->bxp != &unitxf)
229 >                free((void *)f->bxp);
230 >        if ((f->fxp != NULL) & (f->fxp != &unitxf))
231 >                free((void *)f->fxp);
232          free((void *)f);
233          m->os = NULL;
234   }
# Line 167 | Line 236 | freefunc(                      /* free memory associated with modifier */
236  
237   int
238   setfunc(                        /* set channels for function call */
239 <        OBJREC  *m,                     /* can be NULL */
240 <        RAY  *r
239 >        OBJREC  *m,
240 >        RAY     *r
241   )
242   {
243 <        static RNUMBER  lastrno = ~0;
244 <        MFUNC  *f;
243 >        static RNUMBER  lastrno = ~0;
244 >        MFUNC           *f;
245                                          /* get function if any */
246 <        if (m != NULL) {
247 <                if ((f = (MFUNC *)m->os) == NULL)
248 <                        objerror(m, CONSISTENCY, "setfunc called before getfunc");
249 <                setcontext(f->ctx);     /* set evaluator context */
181 <        } else
182 <                setcontext("");
246 >        if ((f = (MFUNC *)m->os) == NULL)
247 >                objerror(m, CONSISTENCY, "setfunc called before getfunc");
248 >                
249 >        setcontext(f->ctx);             /* set evaluator context */
250                                          /* check to see if matrix set */
251          if ((m == fobj) & (r->rno == lastrno))
252                  return(0);
253          fobj = m;
254          fray = r;
255 <        if (r->rox != NULL)
256 <                if (f->b != &unitxf) {
257 <                        funcxf.sca = r->rox->b.sca * f->b->sca;
258 <                        multmat4(funcxf.xfm, r->rox->b.xfm, f->b->xfm);
255 >        if (r->rox != NULL) {
256 >                if (f->bxp != &unitxf) {
257 >                        funcxf.sca = r->rox->b.sca * f->bxp->sca;
258 >                        multmat4(funcxf.xfm, r->rox->b.xfm, f->bxp->xfm);
259                  } else
260                          funcxf = r->rox->b;
261 <        else
262 <                funcxf = *(f->b);
261 >        } else
262 >                funcxf = *f->bxp;
263          lastrno = r->rno;
264          eclock++;               /* notify expression evaluator */
265          return(1);
266   }
267  
268  
269 + int
270 + worldfunc(                      /* special function context sans object */
271 +        char    *ctx,
272 +        RAY     *r
273 + )
274 + {
275 +        static RNUMBER  lastrno = ~0;
276 +
277 +        if (rayinitcal[0])              /* initialize on first call */
278 +                initfunc();
279 +                                        /* set evaluator context */
280 +        setcontext(ctx);
281 +                                        /* check if ray already set */
282 +        if ((fobj == NULL) & (r->rno == lastrno))
283 +                return(0);
284 +        fobj = NULL;
285 +        fray = r;
286 +        funcxf = unitxf;
287 +        lastrno = r->rno;
288 +        eclock++;               /* notify expression evaluator */
289 +        return(1);
290 + }
291 +
292 +
293   void
294   loadfunc(                       /* load definition file */
295          char  *fname
# Line 208 | Line 299 | loadfunc(                      /* load definition file */
299  
300          if ((ffname = getpath(fname, getrlibpath(), R_OK)) == NULL) {
301                  sprintf(errmsg, "cannot find function file \"%s\"", fname);
302 <                error(USER, errmsg);
302 >                error(SYSTEM, errmsg);
303          }
304          fcompile(ffname);
305   }
# Line 220 | Line 311 | l_arg(char *nm)                        /* return nth real argument */
311          int  n;
312  
313          if (fobj == NULL)
314 <                error(USER, "arg(n) called without a context");
314 >                error(INTERNAL, "arg(n) called without a modifier context");
315  
316          n = argument(1) + .5;           /* round to integer */
317  
# Line 274 | Line 365 | chanvalue(                     /* return channel n to calcomp */
365                                  fray->ron[2]*funcxf.xfm[2][n-3] )
366                           / funcxf.sca );
367  
368 <        if (n <= 8)                     /* intersection */
368 >        if (n <= 8) {                   /* intersection point */
369 >                if (fray->rot >= FHUGE)
370 >                        return(0.0);    /* XXX should be runtime error? */
371  
372                  return( fray->rop[0]*funcxf.xfm[0][n-6] +
373                                  fray->rop[1]*funcxf.xfm[1][n-6] +
374                                  fray->rop[2]*funcxf.xfm[2][n-6] +
375                                               funcxf.xfm[3][n-6] );
376 +        }
377  
378          if (n == 9)                     /* total distance */
379                  return(raydist(fray,PRIMARY) * funcxf.sca);

Diff Legend

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