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.19 by greg, Wed Jun 19 15:08:39 1991 UTC vs.
Revision 2.4 by greg, Sat Jan 4 19:53:25 1992 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  REFVNAME       "`FILE_REFCNT"
22 + #define  CALSUF         ".cal"
23 + #define  LCALSUF        4
24 +
25   XF  unitxf = {                  /* identity transform */
26          1.0, 0.0, 0.0, 0.0,
27          0.0, 1.0, 0.0, 0.0,
# 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(), l_erf(), l_erfc();
48 <        extern long  eclock;
49 <        static char  *initfile = "rayinit.cal";
50 <        static long  lastrno = -1;
51 <                                        /* check to see if already set */
52 <        if (m == fobj && r->rno == lastrno)
53 <                return(0);
54 <                                        /* initialize if first call */
55 <        if (initfile != NULL) {
56 <                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 >        if (initfile[0]) {              /* initialize on first call */
58 >                setcontext("");
59                  scompile("Dx=$1;Dy=$2;Dz=$3;", NULL, 0);
60                  scompile("Nx=$4;Ny=$5;Nz=$6;", NULL, 0);
61                  scompile("Px=$7;Py=$8;Pz=$9;", NULL, 0);
# Line 55 | Line 68 | XF  *bx;
68                  funset("erf", 1, ':', l_erf);
69                  funset("erfc", 1, ':', l_erfc);
70                  setnoisefuncs();
71 <                initfile = NULL;
71 >                loadfunc(initfile);
72 >                initfile[0] = '\0';
73          }
74 +        if ((na = m->oargs.nsargs) <= ff)
75 +                goto toofew;
76 +        arg = m->oargs.sarg;
77 +        if ((f = (MFUNC *)calloc(1, sizeof(MFUNC))) == NULL)
78 +                goto memerr;
79 +        i = strlen(arg[ff]);                    /* set up context */
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 */
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]);
89 +                        varset(REFVNAME, '=', 1.0);
90 +                } else                          /* reference_count++ */
91 +                        varset(REFVNAME, '=', varvalue(REFVNAME)+1.0);
92 +        }
93 +        curfunc = NULL;                 /* parse expressions */
94 +        sprintf(sbuf, "%s \"%s\"", ofun[m->otype].funame, m->oname);
95 +        for (i=0, ne=0; ef && i < na; i++, ef>>=1)
96 +                if (ef & 1) {                   /* flagged as an expression? */
97 +                        if (ne >= MAXEXPR)
98 +                                objerror(m, INTERNAL, "too many expressions");
99 +                        initstr(arg[i], sbuf, 0);
100 +                        f->ep[ne++] = getE1();
101 +                        if (nextc != EOF)
102 +                                syntax("unexpected character");
103 +                }
104 +        if (ef)
105 +                goto toofew;
106 +        if (i <= ff)                    /* find transform args */
107 +                i = ff+1;
108 +        while (i < na && arg[i][0] != '-')
109 +                i++;
110 +        if (i == na)                    /* no transform */
111 +                f->f = f->b = &unitxf;
112 +        else {                          /* get transform */
113 +                if ((f->b = (XF *)malloc(sizeof(XF))) == NULL)
114 +                        goto memerr;
115 +                if (invxf(f->b, na-i, arg+i) != na-i)
116 +                        objerror(m, USER, "bad transform");
117 +                if (f->b->sca < 0.0)
118 +                        f->b->sca = -f->b->sca;
119 +                if (dofwd) {                    /* do both transforms */
120 +                        if ((f->f = (XF *)malloc(sizeof(XF))) == NULL)
121 +                                goto memerr;
122 +                        xf(f->f, na-i, arg+i);
123 +                        if (f->f->sca < 0.0)
124 +                                f->f->sca = -f->f->sca;
125 +                }
126 +        }
127 +        m->os = (char *)f;
128 +        return(f);
129 + toofew:
130 +        objerror(m, USER, "too few string arguments");
131 + memerr:
132 +        error(SYSTEM, "out of memory in getfunc");
133 + }
134 +
135 +
136 + freefunc(m)                     /* free memory associated with modifier */
137 + OBJREC  *m;
138 + {
139 +        register MFUNC  *f;
140 +        register int  i;
141 +
142 +        if ((f = (MFUNC *)m->os) == NULL)
143 +                return;
144 +        for (i = 0; f->ep[i] != NULL; i++)
145 +                epfree(f->ep[i]);
146 +        if (f->ctx[0]) {                        /* done with definitions */
147 +                setcontext(f->ctx);
148 +                i = varvalue(REFVNAME)-.5;      /* reference_count-- */
149 +                if (i > 0)
150 +                        varset(REFVNAME, '=', (double)i);
151 +                else
152 +                        dcleanup(2);            /* remove definitions */
153 +                freestr(f->ctx);
154 +        }
155 +        if (f->b != &unitxf)
156 +                free((char *)f->b);
157 +        if (f->f != NULL && f->f != &unitxf)
158 +                free((char *)f->f);
159 +        free((char *)f);
160 +        m->os = NULL;
161 + }
162 +
163 +
164 + setfunc(m, r)                   /* set channels for function call */
165 + OBJREC  *m;
166 + register RAY  *r;
167 + {
168 +        static long  lastrno = -1;
169 +        register MFUNC  *f;
170 +                                        /* get function */
171 +        if ((f = (MFUNC *)m->os) == NULL)
172 +                error(m, CONSISTENCY, "setfunc called before getfunc");
173 +                                        /* set evaluator context */
174 +        setcontext(f->ctx);
175 +                                        /* check to see if matrix set */
176 +        if (m == fobj && r->rno == lastrno)
177 +                return(0);
178          fobj = m;
179          fray = r;
180          if (r->rox != NULL)
181 <                if (bx != &unitxf) {
182 <                        funcxf.sca = r->rox->b.sca * bx->sca;
183 <                        multmat4(funcxf.xfm, r->rox->b.xfm, bx->xfm);
181 >                if (f->b != &unitxf) {
182 >                        funcxf.sca = r->rox->b.sca * f->b->sca;
183 >                        multmat4(funcxf.xfm, r->rox->b.xfm, f->b->xfm);
184                  } else
185                          copystruct(&funcxf, &r->rox->b);
186          else
187 <                copystruct(&funcxf, bx);
187 >                copystruct(&funcxf, f->b);
188          lastrno = r->rno;
189          eclock++;               /* notify expression evaluator */
190          return(1);
191   }
192  
193  
76 – setfunc(m, r)                           /* simplified interface to setmap */
77 – register OBJREC  *m;
78 – RAY  *r;
79 – {
80 –        register XF  *mxf;
81 –
82 –        if ((mxf = (XF *)m->os) == NULL) {
83 –                register int  n;
84 –                register char  **sa;
85 –
86 –                for (n = m->oargs.nsargs, sa = m->oargs.sarg;
87 –                                n > 0 && **sa != '-'; n--, sa++)
88 –                        ;
89 –                if (n == 0)
90 –                        mxf = &unitxf;
91 –                else {
92 –                        mxf = (XF *)malloc(sizeof(XF));
93 –                        if (mxf == NULL)
94 –                                goto memerr;
95 –                        if (invxf(mxf, n, sa) != n)
96 –                                objerror(m, USER, "bad transform");
97 –                        if (mxf->sca < 0.0)
98 –                                mxf->sca = -mxf->sca;
99 –                }
100 –                m->os = (char *)mxf;
101 –        }
102 –        return(setmap(m, r, mxf));
103 – memerr:
104 –        error(SYSTEM, "out of memory in setfunc");
105 – }
106 –
107 –
194   loadfunc(fname)                 /* load definition file */
195   char  *fname;
196   {
# Line 119 | Line 205 | char  *fname;
205   }
206  
207  
208 < double
208 > static double
209   l_arg()                         /* return nth real argument */
210   {
211          extern double  argument();
# Line 138 | Line 224 | l_arg()                                /* return nth real argument */
224   }
225  
226  
227 < double
227 > static double
228   l_erf()                         /* error function */
229   {
230          extern double  erf();
# Line 147 | Line 233 | l_erf()                                /* error function */
233   }
234  
235  
236 < double
236 > static double
237   l_erfc()                        /* cumulative error function */
238   {
239          extern double  erfc();

Diff Legend

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