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 1.22 by greg, Thu Aug 8 15:41:14 1991 UTC vs.
Revision 2.10 by greg, Tue Aug 3 12:04:58 1993 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1991 Regents of the University of California */
1 > /* Copyright (c) 1992 Regents of the University of California */
2  
3   #ifndef lint
4   static char SCCSid[] = "$SunId$ LBL";
# Line 14 | Line 14 | static char SCCSid[] = "$SunId$ LBL";
14  
15   #include  "otypes.h"
16  
17 + #include  "func.h"
18  
19 +
20   #define  INITFILE       "rayinit.cal"
21 < #define  DEFVNAME       "FILE_LOADED`"
21 > #define  CALSUF         ".cal"
22 > #define  LCALSUF        4
23 > char  REFVNAME[] = "`FILE_REFCNT";
24  
25   XF  unitxf = {                  /* identity transform */
26 <        1.0, 0.0, 0.0, 0.0,
27 <        0.0, 1.0, 0.0, 0.0,
28 <        0.0, 0.0, 1.0, 0.0,
29 <        0.0, 0.0, 0.0, 1.0,
26 >        {{1.0, 0.0, 0.0, 0.0},
27 >        {0.0, 1.0, 0.0, 0.0},
28 >        {0.0, 0.0, 1.0, 0.0},
29 >        {0.0, 0.0, 0.0, 1.0}},
30          1.0
31   };
32  
# Line 30 | Line 34 | XF  funcxf;                    /* current transformation */
34   static OBJREC  *fobj = NULL;    /* current function object */
35   static RAY  *fray = NULL;       /* current function ray */
36  
37 + static double  l_erf(), l_erfc(), l_arg();
38  
39 < setmap(m, r, bx)                /* set channels for function call */
39 >
40 > MFUNC *
41 > getfunc(m, ff, ef, dofwd)       /* get function for this modifier */
42   OBJREC  *m;
43 + int  ff;
44 + unsigned  ef;
45 + int  dofwd;
46 + {
47 +        extern EPNODE  *curfunc;
48 +        static char  initfile[] = INITFILE;
49 +        char  sbuf[MAXSTR];
50 +        register char  **arg;
51 +        register MFUNC  *f;
52 +        int  ne, na;
53 +        register int  i;
54 +                                        /* check to see if done already */
55 +        if ((f = (MFUNC *)m->os) != NULL)
56 +                return(f);
57 +        fobj = NULL; fray = NULL;
58 +        if (initfile[0]) {              /* initialize on first call */
59 +                setcontext("");
60 +                scompile("Dx=$1;Dy=$2;Dz=$3;", NULL, 0);
61 +                scompile("Nx=$4;Ny=$5;Nz=$6;", NULL, 0);
62 +                scompile("Px=$7;Py=$8;Pz=$9;", NULL, 0);
63 +                scompile("T=$10;Ts=$25;Rdot=$11;", NULL, 0);
64 +                scompile("S=$12;Tx=$13;Ty=$14;Tz=$15;", NULL, 0);
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 +                funset("arg", 1, '=', l_arg);
69 +                funset("erf", 1, ':', l_erf);
70 +                funset("erfc", 1, ':', l_erfc);
71 +                setnoisefuncs();
72 +                loadfunc(initfile);
73 +                initfile[0] = '\0';
74 +        }
75 +        if ((na = m->oargs.nsargs) <= ff)
76 +                goto toofew;
77 +        arg = m->oargs.sarg;
78 +        if ((f = (MFUNC *)calloc(1, sizeof(MFUNC))) == NULL)
79 +                goto memerr;
80 +        i = strlen(arg[ff]);                    /* set up context */
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 */
85 +                if (i > LCALSUF && !strcmp(sbuf+i-LCALSUF, CALSUF))
86 +                        sbuf[i-LCALSUF] = '\0'; /* remove suffix */
87 +                setcontext(f->ctx = savestr(sbuf));
88 +                if (!vardefined(REFVNAME)) {    /* file loaded? */
89 +                        loadfunc(arg[ff]);
90 +                        varset(REFVNAME, '=', 1.0);
91 +                } else                          /* reference_count++ */
92 +                        varset(REFVNAME, '=', varvalue(REFVNAME)+1.0);
93 +        }
94 +        curfunc = NULL;                 /* parse expressions */
95 +        sprintf(sbuf, "%s \"%s\"", ofun[m->otype].funame, m->oname);
96 +        for (i=0, ne=0; ef && i < na; i++, ef>>=1)
97 +                if (ef & 1) {                   /* flagged as an expression? */
98 +                        if (ne >= MAXEXPR)
99 +                                objerror(m, INTERNAL, "too many expressions");
100 +                        initstr(arg[i], sbuf, 0);
101 +                        f->ep[ne++] = getE1();
102 +                        if (nextc != EOF)
103 +                                syntax("unexpected character");
104 +                }
105 +        if (ef)
106 +                goto toofew;
107 +        if (i <= ff)                    /* find transform args */
108 +                i = ff+1;
109 +        while (i < na && arg[i][0] != '-')
110 +                i++;
111 +        if (i == na)                    /* no transform */
112 +                f->f = f->b = &unitxf;
113 +        else {                          /* get transform */
114 +                if ((f->b = (XF *)malloc(sizeof(XF))) == NULL)
115 +                        goto memerr;
116 +                if (invxf(f->b, na-i, arg+i) != na-i)
117 +                        objerror(m, USER, "bad transform");
118 +                if (f->b->sca < 0.0)
119 +                        f->b->sca = -f->b->sca;
120 +                if (dofwd) {                    /* do both transforms */
121 +                        if ((f->f = (XF *)malloc(sizeof(XF))) == NULL)
122 +                                goto memerr;
123 +                        xf(f->f, na-i, arg+i);
124 +                        if (f->f->sca < 0.0)
125 +                                f->f->sca = -f->f->sca;
126 +                }
127 +        }
128 +        m->os = (char *)f;
129 +        return(f);
130 + toofew:
131 +        objerror(m, USER, "too few string arguments");
132 + memerr:
133 +        error(SYSTEM, "out of memory in getfunc");
134 + }
135 +
136 +
137 + freefunc(m)                     /* free memory associated with modifier */
138 + OBJREC  *m;
139 + {
140 +        register MFUNC  *f;
141 +        register int  i;
142 +
143 +        if ((f = (MFUNC *)m->os) == NULL)
144 +                return;
145 +        for (i = 0; f->ep[i] != NULL; i++)
146 +                epfree(f->ep[i]);
147 +        if (f->ctx[0]) {                        /* done with definitions */
148 +                setcontext(f->ctx);
149 +                i = varvalue(REFVNAME)-.5;      /* reference_count-- */
150 +                if (i > 0)
151 +                        varset(REFVNAME, '=', (double)i);
152 +                else
153 +                        dcleanup(2);            /* remove definitions */
154 +                freestr(f->ctx);
155 +        }
156 +        if (f->b != &unitxf)
157 +                free((char *)f->b);
158 +        if (f->f != NULL && f->f != &unitxf)
159 +                free((char *)f->f);
160 +        free((char *)f);
161 +        m->os = NULL;
162 + }
163 +
164 +
165 + setfunc(m, r)                   /* set channels for function call */
166 + OBJREC  *m;
167   register RAY  *r;
37 XF  *bx;
168   {
169 <        extern long  eclock;
170 <        static long  lastrno = -1;
171 <                                        /* check to see if already set */
169 >        static unsigned long  lastrno = ~0;
170 >        register MFUNC  *f;
171 >                                        /* get function */
172 >        if ((f = (MFUNC *)m->os) == NULL)
173 >                error(m, CONSISTENCY, "setfunc called before getfunc");
174 >                                        /* set evaluator context */
175 >        setcontext(f->ctx);
176 >                                        /* check to see if matrix set */
177          if (m == fobj && r->rno == lastrno)
178                  return(0);
179          fobj = m;
180          fray = r;
181          if (r->rox != NULL)
182 <                if (bx != &unitxf) {
183 <                        funcxf.sca = r->rox->b.sca * bx->sca;
184 <                        multmat4(funcxf.xfm, r->rox->b.xfm, bx->xfm);
182 >                if (f->b != &unitxf) {
183 >                        funcxf.sca = r->rox->b.sca * f->b->sca;
184 >                        multmat4(funcxf.xfm, r->rox->b.xfm, f->b->xfm);
185                  } else
186                          copystruct(&funcxf, &r->rox->b);
187          else
188 <                copystruct(&funcxf, bx);
188 >                copystruct(&funcxf, f->b);
189          lastrno = r->rno;
190          eclock++;               /* notify expression evaluator */
191          return(1);
192   }
193  
194  
60 setfunc(m, r)                           /* simplified interface to setmap */
61 register OBJREC  *m;
62 RAY  *r;
63 {
64        register XF  *mxf;
65
66        if ((mxf = (XF *)m->os) == NULL) {
67                register int  n;
68                register char  **sa;
69
70                for (n = m->oargs.nsargs, sa = m->oargs.sarg;
71                                n > 0 && **sa != '-'; n--, sa++)
72                        ;
73                if (n == 0)
74                        mxf = &unitxf;
75                else {
76                        mxf = (XF *)malloc(sizeof(XF));
77                        if (mxf == NULL)
78                                goto memerr;
79                        if (invxf(mxf, n, sa) != n)
80                                objerror(m, USER, "bad transform");
81                        if (mxf->sca < 0.0)
82                                mxf->sca = -mxf->sca;
83                }
84                m->os = (char *)mxf;
85        }
86        return(setmap(m, r, mxf));
87 memerr:
88        error(SYSTEM, "out of memory in setfunc");
89 }
90
91
195   loadfunc(fname)                 /* load definition file */
196   char  *fname;
197   {
# Line 103 | Line 206 | char  *fname;
206   }
207  
208  
209 < double
209 > static double
210   l_arg()                         /* return nth real argument */
211   {
212          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 122 | Line 228 | l_arg()                                /* return nth real argument */
228   }
229  
230  
231 < double
231 > static double
232   l_erf()                         /* error function */
233   {
234          extern double  erf();
# Line 131 | Line 237 | l_erf()                                /* error function */
237   }
238  
239  
240 < double
240 > static double
241   l_erfc()                        /* cumulative error function */
242   {
243          extern double  erfc();
# Line 140 | Line 246 | l_erfc()                       /* cumulative error function */
246   }
247  
248  
143 funcfile(fname)                 /* set context, load file if necessary */
144 register char  *fname;
145 {
146        extern char  *setcontext();
147        static char  initfile[] = INITFILE;
148
149        if (initfile[0]) {              /* initialize on first call */
150                setcontext("");
151                scompile("Dx=$1;Dy=$2;Dz=$3;", NULL, 0);
152                scompile("Nx=$4;Ny=$5;Nz=$6;", NULL, 0);
153                scompile("Px=$7;Py=$8;Pz=$9;", NULL, 0);
154                scompile("T=$10;Rdot=$11;", NULL, 0);
155                scompile("S=$12;Tx=$13;Ty=$14;Tz=$15;", NULL, 0);
156                scompile("Ix=$16;Iy=$17;Iz=$18;", NULL, 0);
157                scompile("Jx=$19;Jy=$20;Jz=$21;", NULL, 0);
158                scompile("Kx=$22;Ky=$23;Kz=$24;", NULL, 0);
159                funset("arg", 1, '=', l_arg);
160                funset("erf", 1, ':', l_erf);
161                funset("erfc", 1, ':', l_erfc);
162                setnoisefuncs();
163                loadfunc(initfile);
164                initfile[0] = '\0';
165        }
166        if (fname[0] == '.' && fname[1] == '\0')
167                setcontext("");                 /* "." means no file */
168        else {
169                setcontext(fname);
170                if (!vardefined(DEFVNAME)) {
171                        loadfunc(fname);
172                        varset(DEFVNAME, ':', 1.0);
173                }
174        }
175 }
176
177
249   double
250   chanvalue(n)                    /* return channel n to calcomp */
251   register int  n;
# Line 182 | Line 253 | register int  n;
253          double  sum;
254          register RAY  *r;
255  
256 +        if (fray == NULL)
257 +                syntax("ray parameter used in constant expression");
258 +
259          if (--n < 0)
260                  goto badchan;
261  
# Line 206 | Line 280 | register int  n;
280                                  fray->rop[2]*funcxf.xfm[2][n-6] +
281                                               funcxf.xfm[3][n-6] );
282  
283 <        if (n == 9) {                   /* distance */
210 <
283 >        if (n == 9) {                   /* total distance */
284                  sum = fray->rot;
285                  for (r = fray->parent; r != NULL; r = r->parent)
286                          sum += r->rot;
287                  return(sum * funcxf.sca);
288  
289          }
290 +
291          if (n == 10)                    /* dot product (range [-1,1]) */
292                  return( fray->rod <= -1.0 ? -1.0 :
293                          fray->rod >= 1.0 ? 1.0 :
# Line 233 | Line 307 | register int  n;
307  
308          if (n < 24)                     /* k unit vector */
309                  return(funcxf.xfm[2][n-21] / funcxf.sca);
310 +
311 +        if (n == 24) {                  /* single ray (shadow) distance */
312 +                sum = fray->rot;
313 +                for (r = fray->parent; r != NULL && r->crtype&SHADOW;
314 +                                r = r->parent)
315 +                        sum += r->rot;
316 +                return(sum * funcxf.sca);
317 +        }
318   badchan:
319          error(USER, "illegal channel number");
320   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines