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 2.7 by greg, Sat Sep 26 08:13:32 1992 UTC vs.
Revision 2.38 by greg, Sat Dec 7 16:38:08 2019 UTC

# Line 1 | Line 1
1 – /* Copyright (c) 1992 Regents of the University of California */
2 –
1   #ifndef lint
2 < static char SCCSid[] = "$SunId$ LBL";
2 > static const char       RCSid[] = "$Id$";
3   #endif
6 –
4   /*
5   *  Compute data values using expression parser
6   *
# Line 16 | Line 13 | static char SCCSid[] = "$SunId$ LBL";
13   *  1/29/87  Made variables conditional (VARIABLE)
14   *
15   *  5/19/88  Added constant subexpression elimination (RCONST)
16 + *
17 + *  2/19/03     Eliminated conditional compiles in favor of esupport extern.
18   */
19  
20 < #include  <stdio.h>
20 > #include "copyright.h"
21  
22 + #include  <stdio.h>
23 + #include  <string.h>
24   #include  <ctype.h>
24 –
25   #include  <errno.h>
26 + #include  <math.h>
27 + #include  <stdlib.h>
28  
29 + #include  "rtmisc.h"
30 + #include  "rtio.h"
31 + #include  "rterror.h"
32   #include  "calcomp.h"
33  
34   #define  MAXLINE        256             /* maximum line length */
# Line 32 | Line 37 | static char SCCSid[] = "$SunId$ LBL";
37  
38   #define  isdecimal(c)   (isdigit(c) || (c) == '.')
39  
40 < #ifndef atof
41 < extern double  atof();
42 < #endif
43 < extern double  pow();
44 < extern char  *fgets(), *savestr();
45 < extern char  *emalloc(), *ecalloc();
41 < extern EPNODE  *curfunc;
42 < extern double  efunc(), evariable();
43 < static double  euminus(), echannel(), eargument(), enumber();
44 < static double  eadd(), esubtr(), emult(), edivi(), epow();
45 < static double  ebotch();
40 > static double  euminus(EPNODE *), eargument(EPNODE *), enumber(EPNODE *);
41 > static double  echannel(EPNODE *);
42 > static double  eadd(EPNODE *), esubtr(EPNODE *),
43 >               emult(EPNODE *), edivi(EPNODE *),
44 >               epow(EPNODE *);
45 > static double  ebotch(EPNODE *);
46  
47 + unsigned int  esupport =                /* what to support */
48 +                E_VARIABLE | E_FUNCTION ;
49 +
50 + int  eofc = 0;                          /* optional end-of-file character */
51   int  nextc;                             /* lookahead character */
52  
53 < double  (*eoper[])() = {                /* expression operations */
53 > double  (*eoper[])(EPNODE *) = {        /* expression operations */
54          ebotch,
51 – #ifdef  VARIABLE
55          evariable,
53 – #else
54 –        ebotch,
55 – #endif
56          enumber,
57          euminus,
58 – #ifdef  INCHAN
58          echannel,
60 – #else
61 –        ebotch,
62 – #endif
63 – #ifdef  FUNCTION
59          efunc,
60          eargument,
66 – #else
61          ebotch,
62          ebotch,
69 – #endif
70 –        ebotch,
71 –        ebotch,
63          0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
64          0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
65          emult,
# Line 94 | Line 85 | static int  linepos;                   /* position in buffer */
85  
86  
87   EPNODE *
88 < eparse(expr)                    /* parse an expression string */
89 < char  *expr;
88 > eparse(                 /* parse an expression string */
89 >    char  *expr
90 > )
91   {
92      EPNODE  *ep;
93  
94      initstr(expr, NULL, 0);
103 – #if  defined(VARIABLE) && defined(FUNCTION)
95      curfunc = NULL;
105 – #endif
96      ep = getE1();
97      if (nextc != EOF)
98          syntax("unexpected character");
# Line 111 | Line 101 | char  *expr;
101  
102  
103   double
104 < eval(expr)                      /* evaluate an expression string */
105 < char  *expr;
104 > eval(                   /* evaluate an expression string */
105 >    char  *expr
106 > )
107   {
108 <    register EPNODE  *ep;
108 >    int  prev_support = esupport;
109 >    EPNODE  *ep;
110      double  rval;
111  
112 +    esupport &= ~E_RCONST;      /* don't bother reducing constant expr */
113      ep = eparse(expr);
114 +    esupport = prev_support;    /* as you were */
115      rval = evalue(ep);
116      epfree(ep);
117      return(rval);
118   }
119  
120  
121 < epfree(epar)                    /* free a parse tree */
122 < register EPNODE  *epar;
121 > int
122 > epcmp(                  /* compare two expressions for equivalence */
123 >    EPNODE  *ep1,
124 >    EPNODE  *ep2
125 > )
126   {
127 <    register EPNODE  *ep;
127 >        double  d;
128  
129 +        if (ep1->type != ep2->type)
130 +                return(1);
131 +
132 +        switch (ep1->type) {
133 +
134 +        case VAR:
135 +                return(ep1->v.ln != ep2->v.ln);
136 +
137 +        case NUM:
138 +                if (ep2->v.num == 0)
139 +                        return(ep1->v.num != 0);
140 +                d = ep1->v.num / ep2->v.num;
141 +                return((d > 1.000000000001) | (d < 0.999999999999));
142 +
143 +        case CHAN:
144 +        case ARG:
145 +                return(ep1->v.chan != ep2->v.chan);
146 +
147 +        case '=':
148 +        case ':':
149 +                return(epcmp(ep1->v.kid->sibling, ep2->v.kid->sibling));
150 +
151 +        case CLKT:
152 +        case SYM:                       /* should never get this one */
153 +                return(0);
154 +
155 +        default:
156 +                ep1 = ep1->v.kid;
157 +                ep2 = ep2->v.kid;
158 +                while (ep1 != NULL) {
159 +                        if (ep2 == NULL)
160 +                                return(1);
161 +                        if (epcmp(ep1, ep2))
162 +                                return(1);
163 +                        ep1 = ep1->sibling;
164 +                        ep2 = ep2->sibling;
165 +                }
166 +                return(ep2 != NULL);
167 +        }
168 + }
169 +
170 +
171 + void
172 + epfree(                 /* free a parse tree */
173 +    EPNODE       *epar
174 + )
175 + {
176 +    EPNODE  *ep;
177 +
178      switch (epar->type) {
179  
134 – #if  defined(VARIABLE) || defined(FUNCTION)
180          case VAR:
181              varfree(epar->v.ln);
182              break;
# Line 139 | Line 184 | register EPNODE         *epar;
184          case SYM:
185              freestr(epar->v.name);
186              break;
142 – #endif
187  
188          case NUM:
189          case CHAN:
190          case ARG:
191 <        case TICK:
191 >        case CLKT:
192              break;
193  
194          default:
195 <            for (ep = epar->v.kid; ep != NULL; ep = ep->sibling)
195 >            while ((ep = epar->v.kid) != NULL) {
196 >                epar->v.kid = ep->sibling;
197                  epfree(ep);
198 +            }
199              break;
200  
201      }
# Line 158 | Line 204 | register EPNODE         *epar;
204   }
205  
206                                  /* the following used to be a switch */
161 – #ifdef  FUNCTION
207   static double
208 < eargument(ep)
209 < EPNODE  *ep;
208 > eargument(
209 >    EPNODE      *ep
210 > )
211   {
212      return(argument(ep->v.chan));
213   }
168 – #endif
214  
215   static double
216 < enumber(ep)
217 < EPNODE  *ep;
216 > enumber(
217 >    EPNODE      *ep
218 > )
219   {
220      return(ep->v.num);
221   }
222  
223   static double
224 < euminus(ep)
225 < EPNODE  *ep;
224 > euminus(
225 >    EPNODE      *ep
226 > )
227   {
228 <    register EPNODE  *ep1 = ep->v.kid;
228 >    EPNODE  *ep1 = ep->v.kid;
229  
230      return(-evalue(ep1));
231   }
232  
186 – #ifdef  INCHAN
233   static double
234 < echannel(ep)
235 < EPNODE  *ep;
234 > echannel(
235 >    EPNODE      *ep
236 > )
237   {
238      return(chanvalue(ep->v.chan));
239   }
193 – #endif
240  
241   static double
242 < eadd(ep)
243 < EPNODE  *ep;
242 > eadd(
243 >    EPNODE      *ep
244 > )
245   {
246 <    register EPNODE  *ep1 = ep->v.kid;
246 >    EPNODE  *ep1 = ep->v.kid;
247  
248      return(evalue(ep1) + evalue(ep1->sibling));
249   }
250  
251   static double
252 < esubtr(ep)
253 < EPNODE  *ep;
252 > esubtr(
253 >    EPNODE      *ep
254 > )
255   {
256 <    register EPNODE  *ep1 = ep->v.kid;
256 >    EPNODE  *ep1 = ep->v.kid;
257  
258      return(evalue(ep1) - evalue(ep1->sibling));
259   }
260  
261   static double
262 < emult(ep)
263 < EPNODE  *ep;
262 > emult(
263 >    EPNODE      *ep
264 > )
265   {
266 <    register EPNODE  *ep1 = ep->v.kid;
266 >    EPNODE  *ep1 = ep->v.kid;
267  
268      return(evalue(ep1) * evalue(ep1->sibling));
269   }
270  
271   static double
272 < edivi(ep)
273 < EPNODE  *ep;
272 > edivi(
273 >    EPNODE      *ep
274 > )
275   {
276 <    register EPNODE  *ep1 = ep->v.kid;
276 >    EPNODE  *ep1 = ep->v.kid;
277      double  d;
278  
279      d = evalue(ep1->sibling);
# Line 236 | Line 286 | EPNODE *ep;
286   }
287  
288   static double
289 < epow(ep)
290 < EPNODE  *ep;
289 > epow(
290 >    EPNODE      *ep
291 > )
292   {
293 <    register EPNODE  *ep1 = ep->v.kid;
293 >    EPNODE  *ep1 = ep->v.kid;
294      double  d;
295      int  lasterrno;
296  
297      lasterrno = errno;
298      errno = 0;
299      d = pow(evalue(ep1), evalue(ep1->sibling));
300 < #ifdef  IEEE
301 <    if (!finite(d))
302 <        errno = EDOM;
300 > #ifdef  isnan
301 >    if (errno == 0) {
302 >        if (isnan(d))
303 >            errno = EDOM;
304 >        else if (isinf(d))
305 >            errno = ERANGE;
306 >    }
307   #endif
308 <    if (errno) {
308 >    if (errno == EDOM || errno == ERANGE) {
309          wputs("Illegal power\n");
310          return(0.0);
311      }
# Line 259 | Line 314 | EPNODE *ep;
314   }
315  
316   static double
317 < ebotch(ep)
318 < EPNODE  *ep;
317 > ebotch(
318 >    EPNODE      *ep
319 > )
320   {
321      eputs("Bad expression!\n");
322      quit(1);
323 +        return 0.0; /* pro forma return */
324   }
325  
326  
327   EPNODE *
328 < ekid(ep, n)                     /* return pointer to a node's nth kid */
329 < register EPNODE  *ep;
330 < register int  n;
328 > ekid(                   /* return pointer to a node's nth kid */
329 >    EPNODE       *ep,
330 >    int  n
331 > )
332   {
333  
334      for (ep = ep->v.kid; ep != NULL; ep = ep->sibling)
# Line 282 | Line 340 | register int  n;
340  
341  
342   int
343 < nekids(ep)                      /* return # of kids for node ep */
344 < register EPNODE  *ep;
343 > nekids(                 /* return # of kids for node ep */
344 >    EPNODE       *ep
345 > )
346   {
347 <    register int  n = 0;
347 >    int  n = 0;
348  
349      for (ep = ep->v.kid; ep != NULL; ep = ep->sibling)
350          n++;
# Line 294 | Line 353 | register EPNODE         *ep;
353   }
354  
355  
356 < initfile(fp, fn, ln)            /* prepare input file */
357 < FILE  *fp;
358 < char  *fn;
359 < int  ln;
356 > void
357 > initfile(               /* prepare input file */
358 >    FILE  *fp,
359 >    char  *fn,
360 >    int  ln
361 > )
362   {
363      static char  inpbuf[MAXLINE];
364  
# Line 311 | Line 372 | int  ln;
372   }
373  
374  
375 < initstr(s, fn, ln)              /* prepare input string */
376 < char  *s;
377 < char  *fn;
378 < int  ln;
375 > void
376 > initstr(                /* prepare input string */
377 >    char  *s,
378 >    char  *fn,
379 >    int  ln
380 > )
381   {
382      infp = NULL;
383      infile = fn;
# Line 325 | Line 388 | int  ln;
388   }
389  
390  
391 < getscanpos(fnp, lnp, spp, fpp)  /* return current scan position */
392 < char  **fnp;
393 < int  *lnp;
394 < char  **spp;
395 < FILE  **fpp;
391 > void
392 > getscanpos(     /* return current scan position */
393 >    char  **fnp,
394 >    int  *lnp,
395 >    char  **spp,
396 >    FILE  **fpp
397 > )
398   {
399      if (fnp != NULL) *fnp = infile;
400      if (lnp != NULL) *lnp = lineno;
# Line 339 | Line 404 | FILE  **fpp;
404  
405  
406   int
407 < scan()                          /* scan next character, return literal next */
407 > scan(void)              /* scan next character, return literal next */
408   {
409 <    register int  lnext = 0;
409 >    int  lnext = 0;
410  
411      do {
412          if (linbuf[linepos] == '\0')
# Line 356 | Line 421 | scan()                         /* scan next character, return literal next
421              nextc = linbuf[linepos++];
422          if (!lnext)
423                  lnext = nextc;
424 +        if (nextc == eofc) {
425 +                nextc = EOF;
426 +                break;
427 +        }
428          if (nextc == '{') {
429              scan();
430              while (nextc != '}')
# Line 371 | Line 440 | scan()                         /* scan next character, return literal next
440  
441  
442   char *
443 < ltoa(l)                         /* convert long to ascii */
444 < long  l;
443 > long2ascii(                           /* convert long to ascii */
444 >    long  l
445 > )
446   {
447      static char  buf[16];
448 <    register char  *cp;
448 >    char  *cp;
449      int  neg = 0;
450  
451      if (l == 0)
# Line 396 | Line 466 | long  l;
466   }
467  
468  
469 < syntax(err)                     /* report syntax error and quit */
470 < char  *err;
469 > void
470 > syntax(                 /* report syntax error and quit */
471 >    char  *err
472 > )
473   {
474 <    register int  i;
474 >    int  i;
475  
476      if (infile != NULL || lineno != 0) {
477          if (infile != NULL) eputs(infile);
478          if (lineno != 0) {
479              eputs(infile != NULL ? ", line " : "line ");
480 <            eputs(ltoa((long)lineno));
480 >            eputs(long2ascii((long)lineno));
481          }
482          eputs(":\n");
483      }
# Line 421 | Line 493 | char  *err;
493   }
494  
495  
496 < addekid(ep, ekid)                       /* add a child to ep */
497 < register EPNODE  *ep;
498 < EPNODE  *ekid;
496 > void
497 > addekid(                        /* add a child to ep */
498 >    EPNODE       *ep,
499 >    EPNODE      *ekid
500 > )
501   {
502      if (ep->v.kid == NULL)
503          ep->v.kid = ekid;
# Line 436 | Line 510 | EPNODE *ekid;
510   }
511  
512  
439 – #if  defined(VARIABLE) || defined(FUNCTION)
513   char *
514 < getname()                       /* scan an identifier */
514 > getname(void)                   /* scan an identifier */
515   {
516 <    static char  str[MAXWORD+1];
517 <    register int  i, lnext;
516 >    static char  str[RMAXWORD+1];
517 >    int  i, lnext;
518  
519      lnext = nextc;
520 <    for (i = 0; i < MAXWORD && isid(lnext); i++, lnext = scan())
520 >    for (i = 0; i < RMAXWORD && isid(lnext); i++, lnext = scan())
521          str[i] = lnext;
522      str[i] = '\0';
523      while (isid(lnext))         /* skip rest of name */
# Line 452 | Line 525 | getname()                      /* scan an identifier */
525  
526      return(str);
527   }
455 – #endif
528  
529  
530   int
531 < getinum()                       /* scan a positive integer */
531 > getinum(void)                   /* scan a positive integer */
532   {
533 <    register int  n, lnext;
533 >    int  n, lnext;
534  
535      n = 0;
536      lnext = nextc;
# Line 471 | Line 543 | getinum()                      /* scan a positive integer */
543  
544  
545   double
546 < getnum()                        /* scan a positive float */
546 > getnum(void)                    /* scan a positive float */
547   {
548 <    register int  i, lnext;
549 <    char  str[MAXWORD+1];
548 >    int  i, lnext;
549 >    char  str[RMAXWORD+1];
550  
551      i = 0;
552      lnext = nextc;
553 <    while (isdigit(lnext) && i < MAXWORD) {
553 >    while (isdigit(lnext) && i < RMAXWORD) {
554          str[i++] = lnext;
555          lnext = scan();
556      }
557 <    if (lnext == '.' && i < MAXWORD) {
557 >    if (lnext == '.' && i < RMAXWORD) {
558          str[i++] = lnext;
559          lnext = scan();
560 <        while (isdigit(lnext) && i < MAXWORD) {
560 >        if (i == 1 && !isdigit(lnext))
561 >            syntax("badly formed number");
562 >        while (isdigit(lnext) && i < RMAXWORD) {
563              str[i++] = lnext;
564              lnext = scan();
565          }
566      }
567 <    if ((lnext == 'e' || lnext == 'E') && i < MAXWORD) {
567 >    if ((lnext == 'e') | (lnext == 'E') && i < RMAXWORD) {
568          str[i++] = lnext;
569          lnext = scan();
570 <        if ((lnext == '-' || lnext == '+') && i < MAXWORD) {
570 >        if ((lnext == '-') | (lnext == '+') && i < RMAXWORD) {
571              str[i++] = lnext;
572              lnext = scan();
573          }
574 <        while (isdigit(lnext) && i < MAXWORD) {
574 >        if (!isdigit(lnext))
575 >            syntax("missing exponent");
576 >        while (isdigit(lnext) && i < RMAXWORD) {
577              str[i++] = lnext;
578              lnext = scan();
579          }
# Line 509 | Line 585 | getnum()                       /* scan a positive float */
585  
586  
587   EPNODE *
588 < getE1()                         /* E1 -> E1 ADDOP E2 */
588 > getE1(void)                     /* E1 -> E1 ADDOP E2 */
589                                  /*       E2 */
590   {
591 <    register EPNODE  *ep1, *ep2;
591 >    EPNODE  *ep1, *ep2;
592  
593      ep1 = getE2();
594      while (nextc == '+' || nextc == '-') {
# Line 521 | Line 597 | getE1()                                /* E1 -> E1 ADDOP E2 */
597          scan();
598          addekid(ep2, ep1);
599          addekid(ep2, getE2());
600 < #ifdef  RCONST
601 <        if (ep1->type == NUM && ep1->sibling->type == NUM)
600 >        if (esupport&E_RCONST &&
601 >                        ep1->type == NUM && ep1->sibling->type == NUM)
602                  ep2 = rconst(ep2);
527 – #endif
603          ep1 = ep2;
604      }
605      return(ep1);
# Line 532 | Line 607 | getE1()                                /* E1 -> E1 ADDOP E2 */
607  
608  
609   EPNODE *
610 < getE2()                         /* E2 -> E2 MULOP E3 */
610 > getE2(void)                     /* E2 -> E2 MULOP E3 */
611                                  /*       E3 */
612   {
613 <    register EPNODE  *ep1, *ep2;
613 >    EPNODE  *ep1, *ep2;
614  
615      ep1 = getE3();
616      while (nextc == '*' || nextc == '/') {
# Line 544 | Line 619 | getE2()                                /* E2 -> E2 MULOP E3 */
619          scan();
620          addekid(ep2, ep1);
621          addekid(ep2, getE3());
622 < #ifdef  RCONST
623 <        if (ep1->type == NUM && ep1->sibling->type == NUM)
624 <                ep2 = rconst(ep2);
625 < #endif
622 >        if (esupport&E_RCONST) {
623 >                EPNODE  *ep3 = ep1->sibling;
624 >                if (ep1->type == NUM && ep3->type == NUM) {
625 >                        ep2 = rconst(ep2);
626 >                } else if (ep3->type == NUM) {
627 >                        if (ep2->type == '/') {
628 >                                if (ep3->v.num == 0)
629 >                                        syntax("divide by zero constant");
630 >                                ep2->type = '*';        /* for speed */
631 >                                ep3->v.num = 1./ep3->v.num;
632 >                        } else if (ep3->v.num == 0) {
633 >                                ep1->sibling = NULL;    /* (E2 * 0) */
634 >                                epfree(ep2);
635 >                                ep2 = ep3;
636 >                        }
637 >                } else if (ep1->type == NUM && ep1->v.num == 0) {
638 >                        epfree(ep3);            /* (0 * E3) or (0 / E3) */
639 >                        ep1->sibling = NULL;
640 >                        efree((char *)ep2);
641 >                        ep2 = ep1;
642 >                }
643 >        }
644          ep1 = ep2;
645      }
646      return(ep1);
# Line 555 | Line 648 | getE2()                                /* E2 -> E2 MULOP E3 */
648  
649  
650   EPNODE *
651 < getE3()                         /* E3 -> E4 ^ E3 */
651 > getE3(void)                     /* E3 -> E4 ^ E3 */
652                                  /*       E4 */
653   {
654 <    register EPNODE  *ep1, *ep2;
654 >        EPNODE  *ep1, *ep2;
655  
656 <    ep1 = getE4();
657 <    if (nextc == '^') {
656 >        ep1 = getE4();
657 >        if (nextc != '^')
658 >                return(ep1);
659          ep2 = newnode();
660          ep2->type = nextc;
661          scan();
662          addekid(ep2, ep1);
663          addekid(ep2, getE3());
664 < #ifdef  RCONST
665 <        if (ep1->type == NUM && ep1->sibling->type == NUM)
666 <                ep2 = rconst(ep2);
667 < #endif
664 >        if (esupport&E_RCONST) {
665 >                EPNODE  *ep3 = ep1->sibling;
666 >                if (ep1->type == NUM && ep3->type == NUM) {
667 >                        ep2 = rconst(ep2);
668 >                } else if (ep1->type == NUM && ep1->v.num == 0) {
669 >                        epfree(ep3);            /* (0 ^ E3) */
670 >                        ep1->sibling = NULL;
671 >                        efree((char *)ep2);
672 >                        ep2 = ep1;
673 >                } else if ((ep3->type == NUM && ep3->v.num == 0) ||
674 >                                (ep1->type == NUM && ep1->v.num == 1)) {
675 >                        epfree(ep2);            /* (E4 ^ 0) or (1 ^ E3) */
676 >                        ep2 = newnode();
677 >                        ep2->type = NUM;
678 >                        ep2->v.num = 1;
679 >                }
680 >        }
681          return(ep2);
575 –    }
576 –    return(ep1);
682   }
683  
684  
685   EPNODE *
686 < getE4()                         /* E4 -> ADDOP E5 */
686 > getE4(void)                     /* E4 -> ADDOP E5 */
687                                  /*       E5 */
688   {
689 <    register EPNODE  *ep1, *ep2;
689 >    EPNODE  *ep1, *ep2;
690  
691      if (nextc == '-') {
692          scan();
# Line 591 | Line 696 | getE4()                                /* E4 -> ADDOP E5 */
696                  return(ep2);
697          }
698          if (ep2->type == UMINUS) {      /* don't generate -(-E5) */
699 +            ep1 = ep2->v.kid;
700              efree((char *)ep2);
701 <            return(ep2->v.kid);
701 >            return(ep1);
702          }
703          ep1 = newnode();
704          ep1->type = UMINUS;
# Line 606 | Line 712 | getE4()                                /* E4 -> ADDOP E5 */
712  
713  
714   EPNODE *
715 < getE5()                         /* E5 -> (E1) */
715 > getE5(void)                     /* E5 -> (E1) */
716                                  /*       VAR */
717                                  /*       NUM */
718                                  /*       $N */
719                                  /*       FUNC(E1,..) */
720                                  /*       ARG */
721   {
722 <    int  i;
723 <    char  *nam;
724 <    register EPNODE  *ep1, *ep2;
722 >        int      i;
723 >        char  *nam;
724 >        EPNODE  *ep1, *ep2;
725  
726 <    if (nextc == '(') {
727 <        scan();
728 <        ep1 = getE1();
729 <        if (nextc != ')')
730 <            syntax("')' expected");
731 <        scan();
732 <        return(ep1);
733 <    }
726 >        if (nextc == '(') {
727 >                scan();
728 >                ep1 = getE1();
729 >                if (nextc != ')')
730 >                        syntax("')' expected");
731 >                scan();
732 >                return(ep1);
733 >        }
734  
735 < #ifdef  INCHAN
736 <    if (nextc == '$') {
737 <        scan();
738 <        ep1 = newnode();
739 <        ep1->type = CHAN;
740 <        ep1->v.chan = getinum();
741 <        return(ep1);
636 <    }
637 < #endif
735 >        if (esupport&E_INCHAN && nextc == '$') {
736 >                scan();
737 >                ep1 = newnode();
738 >                ep1->type = CHAN;
739 >                ep1->v.chan = getinum();
740 >                return(ep1);
741 >        }
742  
743 < #if  defined(VARIABLE) || defined(FUNCTION)
744 <    if (isalpha(nextc) || nextc == CNTXMARK) {
745 <        nam = getname();
746 < #if  defined(VARIABLE) && defined(FUNCTION)
747 <        ep1 = NULL;
748 <        if (curfunc != NULL)
749 <            for (i = 1, ep2 = curfunc->v.kid->sibling;
750 <                                ep2 != NULL; i++, ep2 = ep2->sibling)
751 <                if (!strcmp(ep2->v.name, nam)) {
752 <                    ep1 = newnode();
753 <                    ep1->type = ARG;
754 <                    ep1->v.chan = i;
755 <                    break;
743 >        if (esupport&(E_VARIABLE|E_FUNCTION) &&
744 >                        (isalpha(nextc) || nextc == CNTXMARK)) {
745 >                nam = getname();
746 >                ep1 = NULL;
747 >                if ((esupport&(E_VARIABLE|E_FUNCTION)) == (E_VARIABLE|E_FUNCTION)
748 >                                && curfunc != NULL)
749 >                        for (i = 1, ep2 = curfunc->v.kid->sibling;
750 >                                        ep2 != NULL; i++, ep2 = ep2->sibling)
751 >                                if (!strcmp(ep2->v.name, nam)) {
752 >                                        ep1 = newnode();
753 >                                        ep1->type = ARG;
754 >                                        ep1->v.chan = i;
755 >                                        break;
756 >                                }
757 >                if (ep1 == NULL) {
758 >                        ep1 = newnode();
759 >                        ep1->type = VAR;
760 >                        ep1->v.ln = varinsert(nam);
761                  }
762 <        if (ep1 == NULL)
763 < #endif
764 <        {
765 <            ep1 = newnode();
766 <            ep1->type = VAR;
767 <            ep1->v.ln = varinsert(nam);
762 >                if (esupport&E_FUNCTION && nextc == '(') {
763 >                        ep2 = newnode();
764 >                        ep2->type = FUNC;
765 >                        addekid(ep2, ep1);
766 >                        ep1 = ep2;
767 >                        do {
768 >                                scan();
769 >                                addekid(ep1, getE1());
770 >                        } while (nextc == ',');
771 >                        if (nextc != ')')
772 >                                syntax("')' expected");
773 >                        scan();
774 >                } else if (!(esupport&E_VARIABLE))
775 >                        syntax("'(' expected");
776 >                if (esupport&E_RCONST && isconstvar(ep1))
777 >                        ep1 = rconst(ep1);
778 >                return(ep1);
779          }
660 – #ifdef  FUNCTION
661 –        if (nextc == '(') {
662 –            ep2 = newnode();
663 –            ep2->type = FUNC;
664 –            addekid(ep2, ep1);
665 –            ep1 = ep2;
666 –            do {
667 –                scan();
668 –                addekid(ep1, getE1());
669 –            } while (nextc == ',');
670 –            if (nextc != ')')
671 –                syntax("')' expected");
672 –            scan();
673 –        }
674 – #ifndef  VARIABLE
675 –        else
676 –            syntax("'(' expected");
677 – #endif
678 – #endif
679 – #ifdef  RCONST
680 –        if (isconstvar(ep1))
681 –            ep1 = rconst(ep1);
682 – #endif
683 –        return(ep1);
684 –    }
685 – #endif
780  
781 <    if (isdecimal(nextc)) {
782 <        ep1 = newnode();
783 <        ep1->type = NUM;
784 <        ep1->v.num = getnum();
785 <        return(ep1);
786 <    }
787 <    syntax("unexpected character");
781 >        if (isdecimal(nextc)) {
782 >                ep1 = newnode();
783 >                ep1->type = NUM;
784 >                ep1->v.num = getnum();
785 >                return(ep1);
786 >        }
787 >        syntax("unexpected character");
788 >        return NULL; /* pro forma return */
789   }
790  
791  
697 – #ifdef  RCONST
792   EPNODE *
793 < rconst(epar)                    /* reduce a constant expression */
794 < register EPNODE  *epar;
793 > rconst(                 /* reduce a constant expression */
794 >    EPNODE       *epar
795 > )
796   {
797 <    register EPNODE  *ep;
797 >    EPNODE  *ep;
798  
799      ep = newnode();
800      ep->type = NUM;
801      errno = 0;
802      ep->v.num = evalue(epar);
803 <    if (errno)
803 >    if (errno == EDOM || errno == ERANGE)
804          syntax("bad constant expression");
805      epfree(epar);
806  
# Line 713 | Line 808 | register EPNODE         *epar;
808   }
809  
810  
811 < isconstvar(ep)                  /* is ep linked to a constant expression? */
812 < register EPNODE  *ep;
811 > int
812 > isconstvar(                     /* is ep linked to a constant expression? */
813 >    EPNODE       *ep
814 > )
815   {
816 < #ifdef  VARIABLE
720 <    register EPNODE  *ep1;
721 < #ifdef  FUNCTION
816 >    EPNODE  *ep1;
817  
818 <    if (ep->type == FUNC) {
818 >    if (esupport&E_FUNCTION && ep->type == FUNC) {
819          if (!isconstfun(ep->v.kid))
820                  return(0);
821          for (ep1 = ep->v.kid->sibling; ep1 != NULL; ep1 = ep1->sibling)
# Line 728 | Line 823 | register EPNODE         *ep;
823                  return(0);
824          return(1);
825      }
731 – #endif
826      if (ep->type != VAR)
827          return(0);
828      ep1 = ep->v.ln->def;
829      if (ep1 == NULL || ep1->type != ':')
830          return(0);
831 < #ifdef  FUNCTION
738 <    if (ep1->v.kid->type != SYM)
831 >    if (esupport&E_FUNCTION && ep1->v.kid->type != SYM)
832          return(0);
740 – #endif
833      return(1);
742 – #else
743 –    return(ep->type == FUNC);
744 – #endif
834   }
835  
836  
837 < #if  defined(FUNCTION) && defined(VARIABLE)
838 < isconstfun(ep)                  /* is ep linked to a constant function? */
839 < register EPNODE  *ep;
837 > int
838 > isconstfun(                     /* is ep linked to a constant function? */
839 >    EPNODE       *ep
840 > )
841   {
842 <    register EPNODE  *dp;
843 <    register LIBR  *lp;
842 >    EPNODE  *dp;
843 >    LIBR  *lp;
844  
845      if (ep->type != VAR)
846          return(0);
847 <    if ((dp = ep->v.ln->def) != NULL && dp->v.kid->type == FUNC)
848 <        return(dp->type == ':');
847 >    if ((dp = ep->v.ln->def) != NULL) {
848 >        if (dp->v.kid->type == FUNC)
849 >            return(dp->type == ':');
850 >        else
851 >            return(0);          /* don't identify masked library functions */
852 >    }
853      if ((lp = ep->v.ln->lib) != NULL)
854          return(lp->atyp == ':');
855      return(0);
856   }
763 – #endif
764 – #endif

Diff Legend

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