ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/common/calexpr.c
(Generate patch)

Comparing ray/src/common/calexpr.c (file contents):
Revision 1.16 by greg, Wed Aug 14 16:26:31 1991 UTC vs.
Revision 2.9 by greg, Tue Oct 6 12:29:36 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 24 | Line 24 | static char SCCSid[] = "$SunId$ LBL";
24  
25   #include  <errno.h>
26  
27 + #include  <math.h>
28 +
29   #include  "calcomp.h"
30  
31 < #define  MAXLINE        256             /* maximum line length */
31 > #define  MAXLINE        256             /* maximum line length */
32  
33 < #define  newnode()      (EPNODE *)ecalloc(1, sizeof(EPNODE))
33 > #define  newnode()      (EPNODE *)ecalloc(1, sizeof(EPNODE))
34  
35 < #define  isdecimal(c)   (isdigit(c) || (c) == '.')
35 > #define  isdecimal(c)   (isdigit(c) || (c) == '.')
36  
37 < extern double  atof(), pow();
37 > #ifndef atof
38 > extern double  atof();
39 > #endif
40   extern char  *fgets(), *savestr();
41   extern char  *emalloc(), *ecalloc();
42   extern EPNODE  *curfunc;
# Line 40 | Line 44 | extern double  efunc(), evariable();
44   static double  euminus(), echannel(), eargument(), enumber();
45   static double  eadd(), esubtr(), emult(), edivi(), epow();
46   static double  ebotch();
43 – extern int  errno;
47  
48   int  nextc;                             /* lookahead character */
49  
50 < double  (*eoper[])() = {                /* expression operations */
50 > double  (*eoper[])() = {                /* expression operations */
51          ebotch,
52 < #ifdef  VARIABLE
52 > #ifdef  VARIABLE
53          evariable,
54   #else
55          ebotch,
56   #endif
57          enumber,
58          euminus,
59 < #ifdef  INCHAN
59 > #ifdef  INCHAN
60          echannel,
61   #else
62          ebotch,
63   #endif
64 < #ifdef  FUNCTION
64 > #ifdef  FUNCTION
65          efunc,
66          eargument,
67   #else
# Line 123 | Line 126 | char  *expr;
126  
127  
128   epfree(epar)                    /* free a parse tree */
129 < register EPNODE  *epar;
129 > register EPNODE  *epar;
130   {
131      register EPNODE  *ep;
132  
# Line 133 | Line 136 | register EPNODE  *epar;
136          case VAR:
137              varfree(epar->v.ln);
138              break;
136 – #endif
139              
140          case SYM:
141              freestr(epar->v.name);
142              break;
143 + #endif
144  
145          case NUM:
146          case CHAN:
# Line 156 | Line 159 | register EPNODE  *epar;
159   }
160  
161                                  /* the following used to be a switch */
162 < #ifdef  FUNCTION
162 > #ifdef  FUNCTION
163   static double
164   eargument(ep)
165 < EPNODE  *ep;
165 > EPNODE  *ep;
166   {
167      return(argument(ep->v.chan));
168   }
# Line 167 | Line 170 | EPNODE  *ep;
170  
171   static double
172   enumber(ep)
173 < EPNODE  *ep;
173 > EPNODE  *ep;
174   {
175      return(ep->v.num);
176   }
177  
178   static double
179   euminus(ep)
180 < EPNODE  *ep;
180 > EPNODE  *ep;
181   {
182      register EPNODE  *ep1 = ep->v.kid;
183  
184      return(-evalue(ep1));
185   }
186  
187 < #ifdef  INCHAN
187 > #ifdef  INCHAN
188   static double
189   echannel(ep)
190 < EPNODE  *ep;
190 > EPNODE  *ep;
191   {
192      return(chanvalue(ep->v.chan));
193   }
# Line 192 | Line 195 | EPNODE  *ep;
195  
196   static double
197   eadd(ep)
198 < EPNODE  *ep;
198 > EPNODE  *ep;
199   {
200      register EPNODE  *ep1 = ep->v.kid;
201  
# Line 201 | Line 204 | EPNODE  *ep;
204  
205   static double
206   esubtr(ep)
207 < EPNODE  *ep;
207 > EPNODE  *ep;
208   {
209      register EPNODE  *ep1 = ep->v.kid;
210  
# Line 210 | Line 213 | EPNODE  *ep;
213  
214   static double
215   emult(ep)
216 < EPNODE  *ep;
216 > EPNODE  *ep;
217   {
218      register EPNODE  *ep1 = ep->v.kid;
219  
# Line 219 | Line 222 | EPNODE  *ep;
222  
223   static double
224   edivi(ep)
225 < EPNODE  *ep;
225 > EPNODE  *ep;
226   {
227      register EPNODE  *ep1 = ep->v.kid;
228      double  d;
# Line 235 | Line 238 | EPNODE  *ep;
238  
239   static double
240   epow(ep)
241 < EPNODE  *ep;
241 > EPNODE  *ep;
242   {
243      register EPNODE  *ep1 = ep->v.kid;
244      double  d;
245 <    int  lasterrno;
245 >    int  lasterrno;
246  
247      lasterrno = errno;
248      errno = 0;
249      d = pow(evalue(ep1), evalue(ep1->sibling));
250 < #ifdef  IEEE
250 > #ifdef  IEEE
251      if (!finite(d))
252          errno = EDOM;
253   #endif
# Line 258 | Line 261 | EPNODE  *ep;
261  
262   static double
263   ebotch(ep)
264 < EPNODE  *ep;
264 > EPNODE  *ep;
265   {
266      eputs("Bad expression!\n");
267      quit(1);
# Line 267 | Line 270 | EPNODE  *ep;
270  
271   EPNODE *
272   ekid(ep, n)                     /* return pointer to a node's nth kid */
273 < register EPNODE  *ep;
273 > register EPNODE  *ep;
274   register int  n;
275   {
276  
# Line 281 | Line 284 | register int  n;
284  
285   int
286   nekids(ep)                      /* return # of kids for node ep */
287 < register EPNODE  *ep;
287 > register EPNODE  *ep;
288   {
289      register int  n = 0;
290  
# Line 297 | Line 300 | FILE  *fp;
300   char  *fn;
301   int  ln;
302   {
303 <    static char  inpbuf[MAXLINE];
303 >    static char  inpbuf[MAXLINE];
304  
305      infp = fp;
306      linbuf = inpbuf;
# Line 369 | Line 372 | scan()                         /* scan next character, return literal next
372  
373  
374   char *
375 < ltoa(l)                         /* convert long to ascii */
375 > long2ascii(l)                         /* convert long to ascii */
376   long  l;
377   {
378 <    static char  buf[16];
378 >    static char  buf[16];
379      register char  *cp;
380 <    int  neg = 0;
380 >    int  neg = 0;
381  
382      if (l == 0)
383          return("0");
# Line 403 | Line 406 | char  *err;
406          if (infile != NULL) eputs(infile);
407          if (lineno != 0) {
408              eputs(infile != NULL ? ", line " : "line ");
409 <            eputs(ltoa((long)lineno));
409 >            eputs(long2ascii((long)lineno));
410          }
411 <        eputs(": syntax error:\n");
411 >        eputs(":\n");
412      }
413      eputs(linbuf);
414      if (linbuf[strlen(linbuf)-1] != '\n')
# Line 420 | Line 423 | char  *err;
423  
424  
425   addekid(ep, ekid)                       /* add a child to ep */
426 < register EPNODE  *ep;
427 < EPNODE  *ekid;
426 > register EPNODE  *ep;
427 > EPNODE  *ekid;
428   {
429      if (ep->v.kid == NULL)
430          ep->v.kid = ekid;
# Line 434 | Line 437 | EPNODE  *ekid;
437   }
438  
439  
440 + #if  defined(VARIABLE) || defined(FUNCTION)
441   char *
442   getname()                       /* scan an identifier */
443   {
444 <    static char  str[MAXWORD+1];
444 >    static char  str[MAXWORD+1];
445      register int  i, lnext;
446  
447      lnext = nextc;
# Line 449 | Line 453 | getname()                      /* scan an identifier */
453  
454      return(str);
455   }
456 + #endif
457  
458  
459   int
# Line 479 | Line 484 | getnum()                       /* scan a positive float */
484          lnext = scan();
485      }
486      if (lnext == '.' && i < MAXWORD) {
487 <        str[i++] = lnext;
488 <        lnext = scan();
487 >        str[i++] = lnext;
488 >        lnext = scan();
489          while (isdigit(lnext) && i < MAXWORD) {
490              str[i++] = lnext;
491              lnext = scan();
492          }
493      }
494      if ((lnext == 'e' || lnext == 'E') && i < MAXWORD) {
495 <        str[i++] = lnext;
496 <        lnext = scan();
495 >        str[i++] = lnext;
496 >        lnext = scan();
497          if ((lnext == '-' || lnext == '+') && i < MAXWORD) {
498              str[i++] = lnext;
499              lnext = scan();
# Line 506 | Line 511 | getnum()                       /* scan a positive float */
511  
512   EPNODE *
513   getE1()                         /* E1 -> E1 ADDOP E2 */
514 <                                /*       E2 */
514 >                                /*       E2 */
515   {
516      register EPNODE  *ep1, *ep2;
517  
# Line 517 | Line 522 | getE1()                                /* E1 -> E1 ADDOP E2 */
522          scan();
523          addekid(ep2, ep1);
524          addekid(ep2, getE2());
525 < #ifdef  RCONST
525 > #ifdef  RCONST
526          if (ep1->type == NUM && ep1->sibling->type == NUM)
527                  ep2 = rconst(ep2);
528   #endif
# Line 529 | Line 534 | getE1()                                /* E1 -> E1 ADDOP E2 */
534  
535   EPNODE *
536   getE2()                         /* E2 -> E2 MULOP E3 */
537 <                                /*       E3 */
537 >                                /*       E3 */
538   {
539      register EPNODE  *ep1, *ep2;
540  
# Line 540 | Line 545 | getE2()                                /* E2 -> E2 MULOP E3 */
545          scan();
546          addekid(ep2, ep1);
547          addekid(ep2, getE3());
548 < #ifdef  RCONST
548 > #ifdef  RCONST
549          if (ep1->type == NUM && ep1->sibling->type == NUM)
550                  ep2 = rconst(ep2);
551   #endif
# Line 552 | Line 557 | getE2()                                /* E2 -> E2 MULOP E3 */
557  
558   EPNODE *
559   getE3()                         /* E3 -> E4 ^ E3 */
560 <                                /*       E4 */
560 >                                /*       E4 */
561   {
562      register EPNODE  *ep1, *ep2;
563  
# Line 563 | Line 568 | getE3()                                /* E3 -> E4 ^ E3 */
568          scan();
569          addekid(ep2, ep1);
570          addekid(ep2, getE3());
571 < #ifdef  RCONST
571 > #ifdef  RCONST
572          if (ep1->type == NUM && ep1->sibling->type == NUM)
573                  ep2 = rconst(ep2);
574   #endif
# Line 575 | Line 580 | getE3()                                /* E3 -> E4 ^ E3 */
580  
581   EPNODE *
582   getE4()                         /* E4 -> ADDOP E5 */
583 <                                /*       E5 */
583 >                                /*       E5 */
584   {
585      register EPNODE  *ep1, *ep2;
586  
# Line 603 | Line 608 | getE4()                                /* E4 -> ADDOP E5 */
608  
609   EPNODE *
610   getE5()                         /* E5 -> (E1) */
611 <                                /*       VAR */
612 <                                /*       NUM */
613 <                                /*       $N */
614 <                                /*       FUNC(E1,..) */
615 <                                /*       ARG */
611 >                                /*       VAR */
612 >                                /*       NUM */
613 >                                /*       $N */
614 >                                /*       FUNC(E1,..) */
615 >                                /*       ARG */
616   {
617 <    int  i;
617 >    int  i;
618      char  *nam;
619      register EPNODE  *ep1, *ep2;
620  
# Line 622 | Line 627 | getE5()                                /* E5 -> (E1) */
627          return(ep1);
628      }
629  
630 < #ifdef  INCHAN
630 > #ifdef  INCHAN
631      if (nextc == '$') {
632          scan();
633          ep1 = newnode();
# Line 639 | Line 644 | getE5()                                /* E5 -> (E1) */
644          ep1 = NULL;
645          if (curfunc != NULL)
646              for (i = 1, ep2 = curfunc->v.kid->sibling;
647 <                                ep2 != NULL; i++, ep2 = ep2->sibling)
647 >                                ep2 != NULL; i++, ep2 = ep2->sibling)
648                  if (!strcmp(ep2->v.name, nam)) {
649                      ep1 = newnode();
650                      ep1->type = ARG;
# Line 653 | Line 658 | getE5()                                /* E5 -> (E1) */
658              ep1->type = VAR;
659              ep1->v.ln = varinsert(nam);
660          }
661 < #ifdef  FUNCTION
661 > #ifdef  FUNCTION
662          if (nextc == '(') {
663              ep2 = newnode();
664              ep2->type = FUNC;
# Line 667 | Line 672 | getE5()                                /* E5 -> (E1) */
672                  syntax("')' expected");
673              scan();
674          }
675 < #ifndef  VARIABLE
675 > #ifndef  VARIABLE
676          else
677              syntax("'(' expected");
678   #endif
679   #endif
680 < #ifdef  RCONST
680 > #ifdef  RCONST
681          if (isconstvar(ep1))
682              ep1 = rconst(ep1);
683   #endif
# Line 690 | Line 695 | getE5()                                /* E5 -> (E1) */
695   }
696  
697  
698 < #ifdef  RCONST
698 > #ifdef  RCONST
699   EPNODE *
700   rconst(epar)                    /* reduce a constant expression */
701 < register EPNODE  *epar;
701 > register EPNODE  *epar;
702   {
703      register EPNODE  *ep;
704  
# Line 702 | Line 707 | register EPNODE  *epar;
707      errno = 0;
708      ep->v.num = evalue(epar);
709      if (errno)
710 <        syntax("bad constant expression");
710 >        syntax("bad constant expression");
711      epfree(epar);
712  
713      return(ep);
# Line 710 | Line 715 | register EPNODE  *epar;
715  
716  
717   isconstvar(ep)                  /* is ep linked to a constant expression? */
718 < register EPNODE  *ep;
718 > register EPNODE  *ep;
719   {
720 < #ifdef  VARIABLE
720 > #ifdef  VARIABLE
721      register EPNODE  *ep1;
722 < #ifdef  FUNCTION
722 > #ifdef  FUNCTION
723  
724      if (ep->type == FUNC) {
725          if (!isconstfun(ep->v.kid))
# Line 730 | Line 735 | register EPNODE  *ep;
735      ep1 = ep->v.ln->def;
736      if (ep1 == NULL || ep1->type != ':')
737          return(0);
738 < #ifdef  FUNCTION
738 > #ifdef  FUNCTION
739      if (ep1->v.kid->type != SYM)
740          return(0);
741   #endif
# Line 743 | Line 748 | register EPNODE  *ep;
748  
749   #if  defined(FUNCTION) && defined(VARIABLE)
750   isconstfun(ep)                  /* is ep linked to a constant function? */
751 < register EPNODE  *ep;
751 > register EPNODE  *ep;
752   {
753      register EPNODE  *dp;
754      register LIBR  *lp;
755  
756      if (ep->type != VAR)
757          return(0);
758 <    dp = ep->v.ln->def;
759 <    if (dp != NULL && dp->type != ':')
760 <        return(0);
761 <    if ((dp == NULL || dp->v.kid->type != FUNC)
762 <            && ((lp = liblookup(ep->v.ln->name)) == NULL
758 <                    || lp->atyp != ':'))
759 <        return(0);
760 <    return(1);
758 >    if ((dp = ep->v.ln->def) != NULL && dp->v.kid->type == FUNC)
759 >        return(dp->type == ':');
760 >    if ((lp = ep->v.ln->lib) != NULL)
761 >        return(lp->atyp == ':');
762 >    return(0);
763   }
764   #endif
765   #endif

Diff Legend

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