generate error on mixing terminal functions
[calc.git] / parser.c
CommitLineData
bc97a989 1#include <malloc.h>
f2927108 2#include <math.h>
87621fe1 3#include <stdio.h>
bc97a989
LM
4#include <stdlib.h>
5
6#include "debug.h"
bc97a989
LM
7
8#include "parser.h"
9
5075f6ea
LM
10/* global variables */
11
12double answer = 0;
13
6ba1dd0f
LM
14#define STORAGE_SIZE 10
15double storage[STORAGE_SIZE] = {0};
16
bc97a989
LM
17/* compare codes */
18
19int codecmp (char *ref, char *str)
20{
21 int sig;
22
23 while (*ref != '\0') {
0b9cc9b0
LM
24 if (*ref == '\t') {
25 sig = (*str == '.') ? -1 : ((*str >= '0') && (*str <= '9'));
26 } else {
27 sig = *str - *ref;
28 }
bc97a989
LM
29 if (sig != 0) {
30 return (sig > 0) ? 1 : -1;
31 }
0b9cc9b0
LM
32 str++;
33 ref++;
bc97a989
LM
34 }
35
36 return 0;
37}
38
39/* allocate new element */
40
11cda8d7 41element_t *newelement (func_t function, int nbops, int prio)
bc97a989
LM
42{
43 element_t *new = (element_t *) calloc (1, sizeof (element_t));
44 if (new == NULL) {
04d68907 45 VERBOSE (ERROR, fprintf (stderr, "can't allocate memory\n"));
bc97a989
LM
46 return NULL;
47 }
48 new->func = function;
49 new->nbops = nbops;
11cda8d7 50 new->prio = prio;
bc97a989
LM
51
52 return new;
53}
54
031d7bba
LM
55/* desallocate element */
56
57void delelement (element_t *root)
58{
59 int i;
60 if ((root != NULL) && (root != ERROR_OP)) {
61 for (i = 0; i < root->nbops; i++) {
62 if ((root->ops[i] != NULL) && (root->ops[i] != ERROR_OP)) {
63 delelement (root->ops[i]);
64 }
65 }
66 free (root);
67 }
68}
69
bc97a989
LM
70/* functions */
71
ca3e2a2f 72#define NB_OPERATORS 14
bc97a989
LM
73
74keyword_t operators[NB_OPERATORS] = {
11cda8d7
LM
75 { "+\t", Add, 2, 1, 1},
76 { "-\t", Sub, 2, 1, 1},
77 { "*", Mul, 2, 1, 2},
78 { "/", Div, 2, 1, 2},
c47a9298 79 { "%", Mod, 2, 1, 3},
7fe742c9 80 { "^", Pow, 2, 1, 4},
ca3e2a2f
LM
81 { "==", Equal, 2, 2, -1},
82 { "!=", Diff, 2, 2, -1},
83 { ">=", Ge, 2, 2, -1},
84 { "<=", Le, 2, 2, -1},
85 { ">", Gt, 2, 1, -1},
86 { "<", Lt, 2, 1, -1},
87 { "&", And, 2, 1, -2},
88 { "|", Or, 2, 1, -2}
bc97a989
LM
89};
90
ca3e2a2f 91#define NB_FUNCTIONS 13
bc97a989 92keyword_t functions[NB_FUNCTIONS] = {
c47a9298
LM
93 { "sqrt", Sqr, 1, 4, 5},
94 { "pow", Pow, 2, 3, 5},
95 { "cos", Cos, 1, 3, 5},
96 { "sin", Sin, 1, 3, 5},
471da7c9 97 { "atan", Atan, 1, 4, 5},
c47a9298 98 { "exp", Exp, 1, 3, 5},
89cf0955 99 { "log", Log, 1, 3, 5},
ca3e2a2f 100 { "sto", Store, 2, 3, 9},
6ba1dd0f 101 { "rcl", Recall, 1, 3, 5},
ca3e2a2f
LM
102 { "disp", Disp, 0, 4, 9},
103 { "quit", Quit, 0, 4, 9},
104 { "help", Help, 0, 4, 9},
105 { "!", Not, 1, 1, 5}
471da7c9
LM
106};
107
5075f6ea 108#define NB_CONSTANTS 3
471da7c9 109keyword_t constants[NB_CONSTANTS] = {
5075f6ea
LM
110 { "ans", Ans, 0, 3, 5},
111 { "e", E, 0, 1, 5},
112 { "pi", Pi, 0, 2, 5}
bc97a989
LM
113};
114
d2ff8478
LM
115/* subparser function */
116
117element_t *subparser (element_t **proot, char **pstr, func_t func, int nbops, int prio)
118{
119 element_t *new = newelement (func, nbops, prio);
120 if (new == NULL) {
121 return ERROR_OP;
122 }
123 new->ops[0] = *proot;
124 new->ops[1] = parser (*pstr, pstr, new->prio);
67b13e3d
LM
125 if ((new->ops[1] != NULL) && (new->ops[1] != ERROR_OP) && (new->ops[1]->prio == 9)) {
126 delelement (new->ops[1]);
127 new->ops[1] = ERROR_OP;
128 }
d2ff8478 129 if (new->ops[1] == ERROR_OP) {
031d7bba
LM
130 delelement (new);
131 *proot = NULL;
d2ff8478
LM
132 return ERROR_OP;
133 }
3b4b0bbe 134 *proot = newelement (Val, 1, 5);
031d7bba
LM
135 if (*proot == NULL) {
136 delelement (new);
d2ff8478
LM
137 return ERROR_OP;
138 }
139 (*proot)->ops[0] = new;
140
141 return *proot;
142}
143
bc97a989
LM
144/* parser function */
145
ef37d966 146element_t *parser (char *str, char **next, int prio)
0b489a77 147{
bc97a989 148 element_t *root = NULL;
85b4a72c 149 int i;
bc97a989 150
04d68907 151 VERBOSE (DEBUG, fprintf (stdout, "Starting parsing\n"));
49223129 152
bc97a989
LM
153 /* main loop */
154 while (*str != '\0') {
155 int found = 0;
156 element_t *new = NULL;
04d68907 157 VERBOSE (INFO, fprintf (stdout, "Processing: %s\n", str));
bc97a989
LM
158
159 /* skip spaces and tabs */
160
161 if ((*str == ' ') || (*str == '\t')) {
162 str++;
163 continue;
164 }
165
85b4a72c 166 /* check for open bracket */
0b489a77 167
0b489a77 168 if (*str == '(') {
04d68907 169 VERBOSE (DEBUG, fprintf (stdout, "start processing bracket\n"));
0b489a77 170 if (root) {
85b4a72c
LM
171 do {
172 found = 0;
ef37d966 173 new = parser (str + 1, &str, 0);
67b13e3d
LM
174 if ((new != NULL) && (new != ERROR_OP) && (new->prio == 9)) {
175 delelement (new);
176 new = ERROR_OP;
177 }
85b4a72c 178 if (new == ERROR_OP) {
031d7bba 179 delelement (root);
85b4a72c
LM
180 return ERROR_OP;
181 }
182 for (i = 0; i < root->nbops; i++) {
183 if (root->ops[i] == NULL) {
0b489a77
LM
184 root->ops[i] = new;
185 found = 1;
186 break;
187 }
188 }
85b4a72c 189 if (!found) {
031d7bba
LM
190 delelement (new);
191 delelement (root);
85b4a72c
LM
192 return ERROR_OP;
193 }
194 } while (*str == ',');
195 } else {
3b4b0bbe 196 root = newelement (Val, 1, 5);
85b4a72c 197 if (root == NULL) {
0b489a77
LM
198 return ERROR_OP;
199 }
ef37d966 200 new = parser (str + 1, &str, 0);
67b13e3d
LM
201 if ((new != NULL) && (new != ERROR_OP) && (new->prio == 9)) {
202 delelement (new);
203 new = ERROR_OP;
204 }
85b4a72c 205 if ((new == ERROR_OP) || (*str == ',')) {
031d7bba
LM
206 delelement (new);
207 delelement (root);
efdfb543
LM
208 return ERROR_OP;
209 }
85b4a72c 210 root->ops[0] = new;
0b489a77 211 }
85b4a72c 212 str++;
04d68907 213 VERBOSE (DEBUG, fprintf (stdout, "stop processing bracket\n"));
efdfb543 214 continue;
49223129 215 }
85b4a72c
LM
216
217 /* check for closing bracket or koma */
218
219 if ((*str == ')') || (*str == ',')) {
49223129 220 if (next != NULL) {
85b4a72c 221 *next = str;
49223129
LM
222 }
223 return root;
0b489a77
LM
224 }
225
bc97a989
LM
226 /* look for operators */
227
228 for (i = 0; i < NB_OPERATORS; i++) {
229 keyword_t *operator = operators + i;
230 if (codecmp (operator->keyword, str) == 0) {
04d68907 231 VERBOSE (DEBUG, fprintf (stdout, "start processing operator\n"));
11cda8d7 232 if (root) {
ef37d966 233 if ((prio) && (prio > operator->prio)) {
04d68907 234 VERBOSE (DEBUG, fprintf (stdout, "stop because operator priority\n"));
ef37d966
LM
235 *next = str;
236 return root;
237 }
238 str += operator->offset;
04d68907 239 VERBOSE (INFO, fprintf (stdout, "Oper: %d\n", operator->func));
d2ff8478 240 if (subparser (&root, &str, operator->func, operator->nbops, operator->prio) == ERROR_OP) {
031d7bba 241 delelement (root);
ef37d966 242 return ERROR_OP;
11cda8d7 243 }
3b4b0bbe
LM
244 } else if (*str == '-') {
245 new = newelement (Sig, 1, 9);
246 if (new == NULL) {
247 return ERROR_OP;
248 }
249 root = new;
49223129 250 } else {
0b489a77 251 return ERROR_OP;
bc97a989 252 }
bc97a989 253 found = 1;
04d68907 254 VERBOSE (DEBUG, fprintf (stdout, "stop processing operator\n"));
bc97a989
LM
255 break;
256 }
257 }
258 if (found) {
bc97a989
LM
259 continue;
260 }
49223129 261
bc97a989
LM
262 /* look for functions */
263
264 for (i = 0; i < NB_FUNCTIONS; i++) {
265 keyword_t *function = functions + i;
266 if (codecmp (function->keyword, str) == 0) {
04d68907 267 VERBOSE (DEBUG, fprintf (stdout, "start processing function\n"));
bc97a989 268 if (root == NULL) {
04d68907 269 VERBOSE (INFO, fprintf (stdout, "Func: %d\n", function->func));
11cda8d7 270 new = newelement (function->func, function->nbops, function->prio);
bc97a989 271 if (new == NULL) {
0b489a77 272 return ERROR_OP;
bc97a989
LM
273 }
274 root = new;
031d7bba
LM
275 } else {
276 delelement (root);
0b489a77 277 return ERROR_OP;
bc97a989
LM
278 }
279 str += function->offset;
280 found = 1;
04d68907 281 VERBOSE (DEBUG, fprintf (stdout, "stop processing function\n"));
bc97a989
LM
282 break;
283 }
284 }
285 if (found) {
bc97a989
LM
286 continue;
287 }
288
471da7c9
LM
289 /* look for constant */
290
291 for (i = 0; i < NB_CONSTANTS; i++) {
292 keyword_t *constant = constants + i;
293 if (codecmp (constant->keyword, str) == 0) {
294 VERBOSE (DEBUG, fprintf (stdout, "start processing constant\n"));
295 if (root == NULL) {
5075f6ea 296 VERBOSE (INFO, fprintf (stdout, "Const: %d\n", constant->func));
471da7c9
LM
297 new = newelement (constant->func, constant->nbops, constant->prio);
298 if (new == NULL) {
299 return ERROR_OP;
300 }
301 root = new;
302 } else {
303 delelement (root);
304 return ERROR_OP;
305 }
306 str += constant->offset;
307 found = 1;
308 VERBOSE (DEBUG, fprintf (stdout, "stop processing constant\n"));
309 break;
310 }
311 }
312 if (found) {
313 continue;
314 }
315
bc97a989
LM
316 /* look for number */
317
0b9cc9b0
LM
318 if (((*str >= '0') && (*str <= '9')) ||
319 (*str == '.') || (*str == '+') || (*str == '-')) {
04d68907 320 VERBOSE (DEBUG, fprintf (stdout, "start processing value\n"));
bc97a989 321 char *pt;
fd88e359 322 double value = strtod (str, &pt);
04d68907 323 VERBOSE (INFO, fprintf (stdout, "Value: %f\n", value));
bc97a989 324 if (str != pt) {
bc97a989 325 if (root == NULL) {
3b4b0bbe 326 new = newelement (Val, 1, 5);
c47a9298
LM
327 if (new == NULL) {
328 return ERROR_OP;
329 }
330 new->value = value;
bc97a989 331 root = new;
c47a9298 332 str = pt;
32741902
LM
333 } else if ((*str == '+') || (*str == '-')) {
334 if ((prio) && (prio > 1)) {
335 VERBOSE (DEBUG, fprintf (stdout, "stop because operator priority\n"));
336 *next = str;
337 return root;
338 }
339 if (subparser (&root, &str, Add, 2, 1) == ERROR_OP) {
031d7bba 340 delelement (root);
0b489a77 341 return ERROR_OP;
bc97a989 342 }
0b9cc9b0 343 } else {
031d7bba 344 delelement (root);
0b9cc9b0 345 return ERROR_OP;
bc97a989 346 }
bc97a989
LM
347 found = 1;
348 }
04d68907 349 VERBOSE (DEBUG, fprintf (stdout, "stop processing value\n"));
bc97a989
LM
350 }
351
352 /* error */
353
354 if (!found) {
031d7bba 355 delelement (root);
0b489a77 356 return ERROR_OP;
bc97a989
LM
357 }
358
359 }
360
49223129
LM
361 if (next != NULL) {
362 *next = str;
363 }
031d7bba 364
bc97a989
LM
365 return root;
366}
367
368/* print element tree */
369
370void print_element (element_t *root, int level)
371{
372 char *func = NULL;
373 int i;
374
49223129 375 if ((root == NULL) || (root == ERROR_OP)) {
bc97a989 376 return;
49223129
LM
377 }
378
bc97a989 379 for (i = 0; i < level; i++) {
04d68907 380 fprintf (stdout, " ");
bc97a989
LM
381 }
382
383 switch (root->func) {
384 case Val: func = "Value"; break;
89cf0955 385 case Sig: func = "Sign"; break;
bc97a989
LM
386 case Add: func = "Addition"; break;
387 case Sub: func = "Subtraction"; break;
388 case Mul: func = "Multiplication"; break;
389 case Div: func = "Division"; break;
c47a9298 390 case Mod: func = "Modulo"; break;
bc97a989
LM
391 case Pow: func = "Power"; break;
392 case Sqr: func = "Square Root"; break;
393 case Cos: func = "Cosine"; break;
394 case Sin: func = "Sine"; break;
471da7c9 395 case Atan: func = "Arc Tangent"; break;
bc97a989
LM
396 case Log: func = "Logarithm"; break;
397 case Exp: func = "Exponantial"; break;
6ba1dd0f
LM
398 case Store: func = "Store"; break;
399 case Recall: func = "Recall"; break;
400 case Disp: func = "Display"; break;
471da7c9
LM
401 case Quit: func = "Quit"; break;
402 case Help: func = "Help"; break;
5075f6ea 403 case Ans: func = "Ans"; break;
471da7c9
LM
404 case Pi: func = "Pi"; break;
405 case E: func = "E"; break;
7fe742c9
LM
406 case Equal: func = "Equal"; break;
407 case Diff: func = "Different"; break;
408 case Ge: func = "Greater or equal"; break;
409 case Le: func = "Lesser or equal"; break;
410 case Gt: func = "Greater"; break;
411 case Lt: func = "Lesser"; break;
ca3e2a2f
LM
412 case And: func = "And"; break;
413 case Or: func = "Or"; break;
414 case Not: func = "Not"; break;
bc97a989
LM
415 }
416
04d68907 417 fprintf (stdout, "Function: %s\n", func);
bc97a989 418
85b4a72c 419 if ((root->func == Val) && (root->ops[0] == NULL)) {
bc97a989 420 for (i = 0; i < level; i++) {
04d68907 421 fprintf (stdout, " ");
bc97a989 422 }
04d68907 423 fprintf (stdout, "value: %f\n", root->value);
bc97a989
LM
424 } else {
425 for (i = 0; i < root->nbops; i++) {
426 print_element (root->ops[i], level + 1);
427 }
428 }
429}
430
6ba1dd0f
LM
431/* storage functions */
432
3639df66 433double store (int index, double value)
6ba1dd0f
LM
434{
435 if ((index > 0) && (index <= STORAGE_SIZE)) {
3639df66 436 storage[index - 1] = value;
6ba1dd0f
LM
437 } else {
438 VERBOSE (WARNING, fprintf (stdout, "invalid index (%d) [1, %d]\n", index, STORAGE_SIZE));
439 }
3639df66 440 return value;
6ba1dd0f
LM
441}
442
443double recall (int index)
444{
445 if ((index > 0) && (index <= STORAGE_SIZE)) {
446 return storage[index - 1];
447 } else {
448 VERBOSE (WARNING, fprintf (stdout, "invalid index (%d) [1, %d]\n", index, STORAGE_SIZE));
449 }
450 return 0;
451}
452
453void display (void)
454{
455 int i;
456 fprintf (stdout, "storage:");
457 for (i = 0; i < STORAGE_SIZE; i++) {
458 fprintf (stdout, " %g", storage[i]);
459 }
460 fprintf (stdout, "\n");
461}
462
89cf0955
LM
463/* quit function */
464
465void quit (void)
466{
04d68907 467 fprintf (stdout, "bye\n");
89cf0955
LM
468 exit (0);
469}
470
471/* help message */
472
473void help (void)
474{
04d68907
LM
475 fprintf (stdout, "calc is a simple calculator\n\n");
476 fprintf (stdout, "supported operators:\n");
477 fprintf (stdout, " + - * / %% ^\n\n");
7fe742c9
LM
478 fprintf (stdout, "camparison operators:\n");
479 fprintf (stdout, " == != >= <= > <\n\n");
04d68907
LM
480 fprintf (stdout, "supported functions:\n");
481 fprintf (stdout, " pow sqrt cos sin atan log exp\n\n");
6ba1dd0f
LM
482 fprintf (stdout, "storage functions:\n");
483 fprintf (stdout, " sto rcl\n\n");
04d68907 484 fprintf (stdout, "miscellaneous functions:\n");
471da7c9
LM
485 fprintf (stdout, " quit help\n\n");
486 fprintf (stdout, "supported constants:\n");
487 fprintf (stdout, " e pi\n");
89cf0955 488}
3b4b0bbe 489
f2927108
LM
490/* evaluate element tree */
491
89cf0955
LM
492#define MASK_SUB 0x1
493#define MASK_DIV 0x2
494
3b4b0bbe 495double evaluate_element (element_t *root, char mask)
f2927108
LM
496{
497 double op0 = 0, op1 = 0;
3b4b0bbe 498 char nextmask = mask;
f2927108
LM
499
500 if ((root == NULL) || (root == ERROR_OP)) {
04d68907 501 VERBOSE (WARNING, fprintf (stdout, "error while evaluating\n"));
f2927108
LM
502 return 0;
503 }
504
3b4b0bbe
LM
505 /* mask to manage sub operator sub and div */
506 switch (root->func) {
507 case Add:
508 nextmask &= ~MASK_SUB;
509 nextmask &= ~MASK_DIV;
510 break;
511 case Sub:
512 nextmask |= MASK_SUB;
513 nextmask &= ~MASK_DIV;
514 break;
515 case Mul:
516 nextmask &= ~MASK_DIV;
517 break;
518 case Div:
519 nextmask |= MASK_DIV;
520 break;
521 default:
522 nextmask = mask;
523 }
524
f2927108
LM
525 switch (root->func) {
526 case Val:
3b4b0bbe
LM
527 case Sig:
528 op0 = (root->ops[0]) ? evaluate_element (root->ops[0], nextmask) : root->value;
f2927108
LM
529 break;
530 case Add:
531 case Sub:
532 case Mul:
533 case Div:
c47a9298 534 case Mod:
f2927108 535 case Pow:
3639df66 536 case Store:
7fe742c9
LM
537 case Equal:
538 case Diff:
539 case Ge:
540 case Le:
541 case Gt:
542 case Lt:
ca3e2a2f
LM
543 case And:
544 case Or:
f2927108 545 if (root->ops[1]) {
3b4b0bbe 546 op1 = evaluate_element (root->ops[1], nextmask);
3639df66 547 } else if (root->func != Store) {
04d68907 548 VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[1])\n"));
f2927108
LM
549 return 0;
550 }
551 /* fallthrough */
552 case Sqr:
553 case Cos:
554 case Sin:
471da7c9 555 case Atan:
f2927108
LM
556 case Log:
557 case Exp:
6ba1dd0f 558 case Recall:
ca3e2a2f 559 case Not:
f2927108 560 if (root->ops[0]) {
3b4b0bbe 561 op0 = evaluate_element (root->ops[0], 0);
f2927108 562 } else {
04d68907 563 VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[0])\n"));
f2927108
LM
564 return 0;
565 }
89cf0955 566 break;
6ba1dd0f 567 case Disp:
471da7c9
LM
568 case Quit:
569 case Help:
5075f6ea 570 case Ans:
471da7c9
LM
571 case Pi:
572 case E:
89cf0955 573 break;
f2927108
LM
574 }
575
576 switch (root->func) {
3b4b0bbe
LM
577 case Val: return op0;
578 case Sig: return -op0;
579 case Add: return ((mask & MASK_SUB) == 0) ? op0 + op1 : op0 - op1;
580 case Sub: return ((mask & MASK_SUB) == 0) ? op0 - op1 : op0 + op1;
581 case Mul: return ((mask & MASK_DIV) == 0) ? op0 * op1 : op0 / op1;
582 case Div: return ((mask & MASK_DIV) == 0) ? op0 / op1 : op0 * op1;
c47a9298 583 case Mod: return fmod (op0, op1);
f2927108
LM
584 case Pow: return pow (op0, op1);
585 case Sqr: return sqrt (op0);
586 case Cos: return cos (op0);
587 case Sin: return sin (op0);
471da7c9 588 case Atan: return atan (op0);
f2927108
LM
589 case Log: return log (op0);
590 case Exp: return exp (op0);
3639df66 591 case Store: return store ((int)op0, (op1) ? op1 : answer);
6ba1dd0f
LM
592 case Recall: return recall ((int)op0);
593 case Disp: display (); break;
471da7c9
LM
594 case Quit: quit (); break;
595 case Help: help (); break;
5075f6ea 596 case Ans: return answer;
471da7c9
LM
597 case Pi: return M_PI;
598 case E: return M_E;
7fe742c9
LM
599 case Equal: return op0 == op1;
600 case Diff: return op0 != op1;
601 case Ge: return op0 >= op1;
602 case Le: return op0 <= op1;
603 case Gt: return op0 > op1;
604 case Lt: return op0 < op1;
ca3e2a2f
LM
605 case And: return (op0 != 0) && (op1 != 0);
606 case Or: return (op0 != 0) || (op1 != 0);
607 case Not: return (op0 == 0);
f2927108
LM
608 }
609
610 return 0;
611}
612
bc97a989 613/* vim: set ts=4 sw=4 et: */