ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cal/calc.c
Revision: 1.1
Committed: Sat Feb 22 02:07:20 2003 UTC (21 years, 2 months ago) by greg
Content type: text/plain
Branch: MAIN
CVS Tags: rad3R5
Log Message:
Changes and check-in for 3.5 release
Includes new source files and modifications not recorded for many years
See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release

File Contents

# User Rev Content
1 greg 1.1 #ifndef lint
2     static const char RCSid[] = "$Id$";
3     #endif
4     /*
5     * calc.c - simple algebraic desk calculator program.
6     *
7     * 4/1/86
8     */
9    
10     #include <stdio.h>
11    
12     #include <setjmp.h>
13    
14     #include <ctype.h>
15    
16     #include "calcomp.h"
17    
18    
19     #define MAXRES 100
20    
21     double result[MAXRES];
22     int nres = 0;
23    
24     #ifndef BSD
25     #define index strchr
26     #endif
27     extern char *index(), *fgets();
28    
29    
30     jmp_buf env;
31     int recover = 0;
32    
33    
34     main(argc, argv)
35     int argc;
36     char *argv[];
37     {
38     char expr[512];
39     FILE *fp;
40     int i;
41     register char *cp;
42    
43     esupport |= E_VARIABLE|E_INCHAN|E_FUNCTION;
44     esupport &= ~(E_REDEFW|E_RCONST|E_OUTCHAN);
45     #ifdef BIGGERLIB
46     biggerlib();
47     #endif
48     varset("PI", ':', 3.14159265358979323846);
49    
50     for (i = 1; i < argc; i++)
51     fcompile(argv[i]);
52    
53     setjmp(env);
54     recover = 1;
55     eclock++;
56    
57     while (fgets(expr, sizeof(expr), stdin) != NULL) {
58     for (cp = expr; *cp && *cp != '\n'; cp++)
59     ;
60     *cp = '\0';
61     switch (expr[0]) {
62     case '\0':
63     continue;
64     case '?':
65     for (cp = expr+1; isspace(*cp); cp++)
66     ;
67     if (*cp)
68     dprint(cp, stdout);
69     else
70     dprint(NULL, stdout);
71     continue;
72     case '>':
73     for (cp = expr+1; isspace(*cp); cp++)
74     ;
75     if (!*cp) {
76     eputs("file name required\n");
77     continue;
78     }
79     if ((fp = fopen(cp, "w")) == NULL) {
80     eputs(cp);
81     eputs(": cannot open\n");
82     continue;
83     }
84     dprint(NULL, fp);
85     fclose(fp);
86     continue;
87     case '<':
88     for (cp = expr+1; isspace(*cp); cp++)
89     ;
90     if (!*cp) {
91     eputs("file name required\n");
92     continue;
93     }
94     fcompile(cp);
95     eclock++;
96     continue;
97     }
98     if ((cp = index(expr, '=')) != NULL ||
99     (cp = index(expr, ':')) != NULL) {
100     if (cp[1])
101     scompile(expr, NULL, 0);
102     else if (*cp == '=') {
103     *cp = '\0';
104     if (!strcmp(expr, "*"))
105     dcleanup(1);
106     else
107     dclear(expr);
108     } else {
109     *cp = '\0';
110     if (!strcmp(expr, "*"))
111     dcleanup(2);
112     else
113     dremove(expr);
114     }
115     eclock++;
116     } else {
117     printf("$%d=%.9g\n", nres+1,
118     result[nres%MAXRES] = eval(expr));
119     nres++;
120     }
121     }
122    
123     recover = 0;
124     quit(0);
125     }
126    
127    
128     double
129     chanvalue(n) /* return channel value */
130     int n;
131     {
132     if (n == 0)
133     n = nres;
134     else if (n > nres || nres-n >= MAXRES) {
135     fprintf(stderr, "$%d: illegal result\n", n);
136     return(0.0);
137     }
138     return(result[(n-1)%MAXRES]);
139     }
140    
141    
142     void
143     eputs(msg)
144     char *msg;
145     {
146     fputs(msg, stderr);
147     }
148    
149    
150     void
151     wputs(msg)
152     char *msg;
153     {
154     eputs(msg);
155     }
156    
157    
158     void
159     quit(code)
160     int code;
161     {
162     if (recover) /* a cavalier approach */
163     longjmp(env, 1);
164     exit(code);
165     }