ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/Development/ray/src/rt/func.c
(Generate patch)

Comparing ray/src/rt/func.c (file contents):
Revision 1.15 by greg, Fri May 24 13:52:10 1991 UTC vs.
Revision 2.13 by greg, Thu Nov 2 17:38:05 1995 UTC

# Line 1 | Line 1
1 < /* Copyright (c) 1991 Regents of the University of California */
1 > /* Copyright (c) 1995 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  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 27 | 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 < register RAY  *r;
44 < XF  *bx;
43 > int  ff;
44 > unsigned  ef;
45 > int  dofwd;
46   {
47 <        extern double  l_arg();
48 <        extern long  eclock;
49 <        static char  *initfile = "rayinit.cal";
50 <                                        /* initialize if first call */
51 <        if (initfile != NULL) {
52 <                loadfunc(initfile);
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;Rdot=$11;", 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 <                initfile = NULL;
72 >                setprismfuncs();
73 >                loadfunc(initfile);
74 >                initfile[0] = '\0';
75          }
76 +        if ((na = m->oargs.nsargs) <= ff)
77 +                goto toofew;
78 +        arg = m->oargs.sarg;
79 +        if ((f = (MFUNC *)calloc(1, sizeof(MFUNC))) == NULL)
80 +                goto memerr;
81 +        i = strlen(arg[ff]);                    /* set up context */
82 +        if (i == 1 && arg[ff][0] == '.')
83 +                setcontext(f->ctx = "");        /* "." means no file */
84 +        else {
85 +                strcpy(sbuf,arg[ff]);   /* file name is context */
86 +                if (i > LCALSUF && !strcmp(sbuf+i-LCALSUF, CALSUF))
87 +                        sbuf[i-LCALSUF] = '\0'; /* remove suffix */
88 +                setcontext(f->ctx = savestr(sbuf));
89 +                if (!vardefined(REFVNAME)) {    /* file loaded? */
90 +                        loadfunc(arg[ff]);
91 +                        varset(REFVNAME, '=', 1.0);
92 +                } else                          /* reference_count++ */
93 +                        varset(REFVNAME, '=', varvalue(REFVNAME)+1.0);
94 +        }
95 +        curfunc = NULL;                 /* parse expressions */
96 +        sprintf(sbuf, "%s \"%s\"", ofun[m->otype].funame, m->oname);
97 +        for (i=0, ne=0; ef && i < na; i++, ef>>=1)
98 +                if (ef & 1) {                   /* flagged as an expression? */
99 +                        if (ne >= MAXEXPR)
100 +                                objerror(m, INTERNAL, "too many expressions");
101 +                        initstr(arg[i], sbuf, 0);
102 +                        f->ep[ne++] = getE1();
103 +                        if (nextc != EOF)
104 +                                syntax("unexpected character");
105 +                }
106 +        if (ef)
107 +                goto toofew;
108 +        if (i <= ff)                    /* find transform args */
109 +                i = ff+1;
110 +        while (i < na && arg[i][0] != '-')
111 +                i++;
112 +        if (i == na)                    /* no transform */
113 +                f->f = f->b = &unitxf;
114 +        else {                          /* get transform */
115 +                if ((f->b = (XF *)malloc(sizeof(XF))) == NULL)
116 +                        goto memerr;
117 +                if (invxf(f->b, na-i, arg+i) != na-i)
118 +                        objerror(m, USER, "bad transform");
119 +                if (f->b->sca < 0.0)
120 +                        f->b->sca = -f->b->sca;
121 +                if (dofwd) {                    /* do both transforms */
122 +                        if ((f->f = (XF *)malloc(sizeof(XF))) == NULL)
123 +                                goto memerr;
124 +                        xf(f->f, na-i, arg+i);
125 +                        if (f->f->sca < 0.0)
126 +                                f->f->sca = -f->f->sca;
127 +                }
128 +        }
129 +        m->os = (char *)f;
130 +        return(f);
131 + toofew:
132 +        objerror(m, USER, "too few string arguments");
133 + memerr:
134 +        error(SYSTEM, "out of memory in getfunc");
135 + }
136 +
137 +
138 + freefunc(m)                     /* free memory associated with modifier */
139 + OBJREC  *m;
140 + {
141 +        register MFUNC  *f;
142 +        register int  i;
143 +
144 +        if ((f = (MFUNC *)m->os) == NULL)
145 +                return;
146 +        for (i = 0; f->ep[i] != NULL; i++)
147 +                epfree(f->ep[i]);
148 +        if (f->ctx[0]) {                        /* done with definitions */
149 +                setcontext(f->ctx);
150 +                i = varvalue(REFVNAME)-.5;      /* reference_count-- */
151 +                if (i > 0)
152 +                        varset(REFVNAME, '=', (double)i);
153 +                else
154 +                        dcleanup(2);            /* remove definitions */
155 +                freestr(f->ctx);
156 +        }
157 +        if (f->b != &unitxf)
158 +                free((char *)f->b);
159 +        if (f->f != NULL && f->f != &unitxf)
160 +                free((char *)f->f);
161 +        free((char *)f);
162 +        m->os = NULL;
163 + }
164 +
165 +
166 + setfunc(m, r)                   /* set channels for function call */
167 + OBJREC  *m;
168 + register RAY  *r;
169 + {
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");
175 +                                        /* set evaluator context */
176 +        setcontext(f->ctx);
177 +                                        /* check to see if matrix set */
178 +        if (m == fobj && r->rno == lastrno)
179 +                return(0);
180          fobj = m;
181          fray = r;
182          if (r->rox != NULL)
183 <                if (bx != &unitxf) {
184 <                        funcxf.sca = r->rox->b.sca * bx->sca;
185 <                        multmat4(funcxf.xfm, r->rox->b.xfm, bx->xfm);
183 >                if (f->b != &unitxf) {
184 >                        funcxf.sca = r->rox->b.sca * f->b->sca;
185 >                        multmat4(funcxf.xfm, r->rox->b.xfm, f->b->xfm);
186                  } else
187                          copystruct(&funcxf, &r->rox->b);
188          else
189 <                copystruct(&funcxf, bx);
189 >                copystruct(&funcxf, f->b);
190 >        lastrno = r->rno;
191          eclock++;               /* notify expression evaluator */
192 +        return(1);
193   }
194  
195  
68 – setfunc(m, r)                           /* simplified interface to setmap */
69 – register OBJREC  *m;
70 – RAY  *r;
71 – {
72 –        register XF  *mxf;
73 –
74 –        if ((mxf = (XF *)m->os) == NULL) {
75 –                register int  n;
76 –                register char  **sa;
77 –
78 –                for (n = m->oargs.nsargs, sa = m->oargs.sarg;
79 –                                n > 0 && **sa != '-'; n--, sa++)
80 –                        ;
81 –                if (n == 0)
82 –                        mxf = &unitxf;
83 –                else {
84 –                        mxf = (XF *)malloc(sizeof(XF));
85 –                        if (mxf == NULL)
86 –                                goto memerr;
87 –                        if (invxf(mxf, n, sa) != n)
88 –                                objerror(m, USER, "bad transform");
89 –                        if (mxf->sca < 0.0)
90 –                                mxf->sca = -mxf->sca;
91 –                }
92 –                m->os = (char *)mxf;
93 –        }
94 –        setmap(m, r, mxf);
95 –        return;
96 – memerr:
97 –        error(SYSTEM, "out of memory in setfunc");
98 – }
99 –
100 –
196   loadfunc(fname)                 /* load definition file */
197   char  *fname;
198   {
199 <        extern char  *libpath;          /* library search path */
199 >        extern char  *getlibpath();     /* 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 112 | Line 207 | char  *fname;
207   }
208  
209  
210 < double
210 > static double
211   l_arg()                         /* return nth real argument */
212   {
213          extern double  argument();
214          register int  n;
215  
216 +        if (fobj == NULL)
217 +                syntax("arg(n) used in constant expression");
218 +
219          n = argument(1) + .5;           /* round to integer */
220  
221          if (n < 1)
# Line 131 | Line 229 | l_arg()                                /* return nth real argument */
229   }
230  
231  
232 + static double
233 + l_erf()                         /* error function */
234 + {
235 +        extern double  erf();
236 +
237 +        return(erf(argument(1)));
238 + }
239 +
240 +
241 + static double
242 + l_erfc()                        /* cumulative error function */
243 + {
244 +        extern double  erfc();
245 +
246 +        return(erfc(argument(1)));
247 + }
248 +
249 +
250   double
251   chanvalue(n)                    /* return channel n to calcomp */
252   register int  n;
253   {
254 <        double  sum;
255 <        register RAY  *r;
254 >        if (fray == NULL)
255 >                syntax("ray parameter used in constant expression");
256  
257          if (--n < 0)
258                  goto badchan;
# Line 162 | Line 278 | register int  n;
278                                  fray->rop[2]*funcxf.xfm[2][n-6] +
279                                               funcxf.xfm[3][n-6] );
280  
281 <        if (n == 9) {                   /* distance */
281 >        if (n == 9)                     /* total distance */
282 >                return(raydist(fray,PRIMARY) * funcxf.sca);
283  
284 <                sum = fray->rot;
285 <                for (r = fray->parent; r != NULL; r = r->parent)
286 <                        sum += r->rot;
287 <                return(sum * funcxf.sca);
284 >        if (n == 10)                    /* dot product (range [-1,1]) */
285 >                return( fray->rod <= -1.0 ? -1.0 :
286 >                        fray->rod >= 1.0 ? 1.0 :
287 >                        fray->rod );
288  
172 –        }
173 –        if (n == 10)                    /* dot product */
174 –                return(fray->rod);
175 –
289          if (n == 11)                    /* scale */
290                  return(funcxf.sca);
291  
# Line 187 | Line 300 | register int  n;
300  
301          if (n < 24)                     /* k unit vector */
302                  return(funcxf.xfm[2][n-21] / funcxf.sca);
303 +
304 +        if (n == 24)                    /* single ray (shadow) distance */
305 +                return((fray->rot+raydist(fray->parent,SHADOW)) * funcxf.sca);
306   badchan:
307          error(USER, "illegal channel number");
308   }

Diff Legend

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