remove unnecessary code
[calc.git] / parser.c
CommitLineData
bc97a989 1#include <malloc.h>
f2927108 2#include <math.h>
87621fe1 3#include <stdio.h>
bc97a989 4#include <stdlib.h>
b9c1c40d 5#include <string.h>
bc97a989
LM
6
7#include "debug.h"
bc97a989
LM
8
9#include "parser.h"
10
5075f6ea
LM
11/* global variables */
12
13double answer = 0;
14
f438579a
LM
15#define DEFAULT_STORAGE_SIZE 10
16int storage_size = -1;
6ba1dd0f 17
f438579a
LM
18double *storage = NULL;
19
20#define DEFAULT_FORMAT "=> %.6g\n"
2a5ec9d1 21char *format = NULL;
9b759869 22char *minform = NULL;
715580ff 23
c7cbb833 24workspace_t *programs = NULL;
05aabb97 25int nb_programs = 0;
c7cbb833 26
bc97a989
LM
27/* compare codes */
28
29int codecmp (char *ref, char *str)
30{
31 int sig;
32
33 while (*ref != '\0') {
0b9cc9b0
LM
34 if (*ref == '\t') {
35 sig = (*str == '.') ? -1 : ((*str >= '0') && (*str <= '9'));
36 } else {
37 sig = *str - *ref;
38 }
bc97a989
LM
39 if (sig != 0) {
40 return (sig > 0) ? 1 : -1;
41 }
0b9cc9b0
LM
42 str++;
43 ref++;
bc97a989
LM
44 }
45
46 return 0;
47}
48
6f14d1cc 49/* calloc or die function */
bc97a989 50
6f14d1cc 51void *callocordie (size_t count, size_t size)
bc97a989 52{
d59a8079
LM
53 if (count * size == 0) {
54 return NULL;
55 }
6f14d1cc 56 void *new = calloc (count, size);
bc97a989 57 if (new == NULL) {
04d68907 58 VERBOSE (ERROR, fprintf (stderr, "can't allocate memory\n"));
d32ea0b0 59 exit (1);
bc97a989 60 }
6f14d1cc
LM
61 return new;
62}
63
64/* allocate new element */
65
66element_t *newelement (func_t function, int nbops, int prio)
67{
68 element_t *new = (element_t *) callocordie (1, sizeof (element_t));
62d56da6 69 if (nbops) {
6f14d1cc 70 new->ops = (element_t **) callocordie (nbops, sizeof (element_t *));
0c95a3d3 71 }
bc97a989
LM
72 new->func = function;
73 new->nbops = nbops;
11cda8d7 74 new->prio = prio;
bc97a989
LM
75
76 return new;
77}
78
031d7bba
LM
79/* desallocate element */
80
81void delelement (element_t *root)
82{
031d7bba 83 if ((root != NULL) && (root != ERROR_OP)) {
0c95a3d3 84 int i;
031d7bba
LM
85 for (i = 0; i < root->nbops; i++) {
86 if ((root->ops[i] != NULL) && (root->ops[i] != ERROR_OP)) {
87 delelement (root->ops[i]);
88 }
89 }
0c95a3d3
LM
90 if (root->nbops) {
91 free (root->ops);
92 }
031d7bba
LM
93 free (root);
94 }
95}
96
2a688642
LM
97/* duplicate element */
98
99element_t *dupelement (element_t *root)
100{
101 element_t *tmp = NULL;
102 int i;
103
104 if ((root == NULL) || (root == ERROR_OP)) {
105 return root;
106 }
107 tmp = newelement (root->func, root->nbops, root->prio);
2a688642
LM
108 tmp->value = root->value;
109 for (i = 0; i < root->nbops; i++) {
110 tmp->ops[i] = dupelement (root->ops[i]);
2a688642
LM
111 }
112 return tmp;
113}
114
bc97a989
LM
115/* functions */
116
ca3e2a2f 117#define NB_OPERATORS 14
bc97a989 118keyword_t operators[NB_OPERATORS] = {
11cda8d7
LM
119 { "+\t", Add, 2, 1, 1},
120 { "-\t", Sub, 2, 1, 1},
121 { "*", Mul, 2, 1, 2},
122 { "/", Div, 2, 1, 2},
c47a9298 123 { "%", Mod, 2, 1, 3},
7fe742c9 124 { "^", Pow, 2, 1, 4},
ca3e2a2f
LM
125 { "==", Equal, 2, 2, -1},
126 { "!=", Diff, 2, 2, -1},
127 { ">=", Ge, 2, 2, -1},
128 { "<=", Le, 2, 2, -1},
129 { ">", Gt, 2, 1, -1},
130 { "<", Lt, 2, 1, -1},
131 { "&", And, 2, 1, -2},
132 { "|", Or, 2, 1, -2}
bc97a989
LM
133};
134
c7cbb833 135#define NB_FUNCTIONS 31
bc97a989 136keyword_t functions[NB_FUNCTIONS] = {
c47a9298
LM
137 { "sqrt", Sqr, 1, 4, 5},
138 { "pow", Pow, 2, 3, 5},
139 { "cos", Cos, 1, 3, 5},
140 { "sin", Sin, 1, 3, 5},
9e52bbc5
LM
141 { "tan", Tan, 1, 3, 5},
142 { "acos", Acos, 1, 4, 5},
143 { "asin", Asin, 1, 4, 5},
471da7c9 144 { "atan", Atan, 1, 4, 5},
89cf0955 145 { "log", Log, 1, 3, 5},
9e52bbc5
LM
146 { "exp", Exp, 1, 3, 5},
147 { "abs", Abs, 1, 3, 5},
148 { "floor", Floor, 1, 5, 5},
149 { "ceil", Ceil, 1, 4, 5},
2a688642 150 { "sto", Store, 2, 3, 5},
6ba1dd0f 151 { "rcl", Recall, 1, 3, 5},
a8cf32ba
LM
152 { "inc", Inc, 1, 3, 5},
153 { "dec", Dec, 1, 3, 5},
ca3e2a2f 154 { "disp", Disp, 0, 4, 9},
f438579a 155 { "mem", Mem, 1, 3, 9},
143f91ae 156 { "clr", Clear, 0, 3, 9},
ca3e2a2f
LM
157 { "quit", Quit, 0, 4, 9},
158 { "help", Help, 0, 4, 9},
94b4e517 159 { "!", Not, 1, 1, 6},
2a688642 160 { "cond", Cond, 3, 4, 5},
715580ff 161 { "while", While, 2, 5, 5},
c7cbb833
LM
162 { "print", Print, 1, 5, 5},
163 { "prog", Prog, 3, 4, 9},
164 { "call", Call, 10, 4, 5},
165 { "ls", List, 0, 2, 9},
166 { "edit", Edit, 1, 4, 9},
167 { "del", Del, 1, 3, 9}
471da7c9
LM
168};
169
5075f6ea 170#define NB_CONSTANTS 3
471da7c9 171keyword_t constants[NB_CONSTANTS] = {
5075f6ea
LM
172 { "ans", Ans, 0, 3, 5},
173 { "e", E, 0, 1, 5},
174 { "pi", Pi, 0, 2, 5}
bc97a989
LM
175};
176
db660e60
LM
177#define NB_SYMBOLS 4
178char *symbols[NB_SYMBOLS] = {
179 "(", ")", "{", "}"
180};
181
d2ff8478
LM
182/* subparser function */
183
184element_t *subparser (element_t **proot, char **pstr, func_t func, int nbops, int prio)
185{
186 element_t *new = newelement (func, nbops, prio);
d2ff8478
LM
187 new->ops[0] = *proot;
188 new->ops[1] = parser (*pstr, pstr, new->prio);
686c4bc2 189 if ((new->ops[1] == NULL) || ((new->ops[1] != ERROR_OP) && (new->ops[1]->prio == 9))) {
67b13e3d
LM
190 delelement (new->ops[1]);
191 new->ops[1] = ERROR_OP;
192 }
d2ff8478 193 if (new->ops[1] == ERROR_OP) {
031d7bba
LM
194 delelement (new);
195 *proot = NULL;
d2ff8478
LM
196 return ERROR_OP;
197 }
3b4b0bbe 198 *proot = newelement (Val, 1, 5);
d2ff8478
LM
199 (*proot)->ops[0] = new;
200
201 return *proot;
202}
203
bc97a989
LM
204/* parser function */
205
ef37d966 206element_t *parser (char *str, char **next, int prio)
0b489a77 207{
bc97a989 208 element_t *root = NULL;
779282bb 209 char *string = str;
85b4a72c 210 int i;
bc97a989 211
04d68907 212 VERBOSE (DEBUG, fprintf (stdout, "Starting parsing\n"));
49223129 213
bc97a989
LM
214 /* main loop */
215 while (*str != '\0') {
216 int found = 0;
217 element_t *new = NULL;
04d68907 218 VERBOSE (INFO, fprintf (stdout, "Processing: %s\n", str));
bc97a989 219
2d0cd54c
LM
220 /* end without printing */
221
222 if (*str == ';') {
223 if (root) {
224 root->hidden = 1;
225 }
226 break;
227 }
228
bc97a989
LM
229 /* skip spaces and tabs */
230
231 if ((*str == ' ') || (*str == '\t')) {
232 str++;
233 continue;
234 }
235
124da7fd
LM
236 /* check for open brace */
237
238 if (*str == '{') {
239 VERBOSE (DEBUG, fprintf (stdout, "start processing brace\n"));
62d56da6 240 if (root != NULL) {
124da7fd
LM
241 delelement (root);
242 return ERROR_OP;
243 }
62d56da6 244 element_t **prog = NULL;
05aabb97 245 new = newelement (Code, 0, 5);
62d56da6
LM
246 root = new;
247 prog = &root;
124da7fd
LM
248
249 do {
124da7fd
LM
250 new = parser (str + 1, &str, 0);
251 if ((new == NULL) || ((new != ERROR_OP) && (new->prio == 9))) {
252 delelement (new);
253 new = ERROR_OP;
254 }
7929d368 255 if (new == ERROR_OP) {
124da7fd 256 delelement (root);
62d56da6 257 return ERROR_OP;
124da7fd 258 }
05aabb97 259 element_t *newprog = newelement (Code, (*prog)->nbops + 1, 5);
62d56da6
LM
260 for (i = 0; i < (*prog)->nbops; i++) {
261 newprog->ops[i] = (*prog)->ops[i];
262 (*prog)->ops[i] = NULL;
124da7fd 263 }
62d56da6
LM
264 newprog->ops[(*prog)->nbops] = new;
265 delelement (*prog);
266 (*prog) = newprog;
124da7fd
LM
267 } while (*str == ',');
268
269 if (*str != '}') {
270 delelement (root);
271 return ERROR_OP;
272 }
273 str++;
274 VERBOSE (DEBUG, fprintf (stdout, "stop processing brace\n"));
275 continue;
276 }
277
85b4a72c 278 /* check for open bracket */
0b489a77 279
0b489a77 280 if (*str == '(') {
04d68907 281 VERBOSE (DEBUG, fprintf (stdout, "start processing bracket\n"));
0b489a77 282 if (root) {
85b4a72c
LM
283 do {
284 found = 0;
ef37d966 285 new = parser (str + 1, &str, 0);
686c4bc2 286 if ((new == NULL) || ((new != ERROR_OP) && (new->prio == 9))) {
67b13e3d
LM
287 delelement (new);
288 new = ERROR_OP;
289 }
686c4bc2 290 if ((new == NULL) || (new == ERROR_OP)) {
031d7bba 291 delelement (root);
85b4a72c
LM
292 return ERROR_OP;
293 }
294 for (i = 0; i < root->nbops; i++) {
295 if (root->ops[i] == NULL) {
0b489a77
LM
296 root->ops[i] = new;
297 found = 1;
298 break;
299 }
300 }
85b4a72c 301 if (!found) {
031d7bba
LM
302 delelement (new);
303 delelement (root);
85b4a72c
LM
304 return ERROR_OP;
305 }
306 } while (*str == ',');
307 } else {
3b4b0bbe 308 root = newelement (Val, 1, 5);
ef37d966 309 new = parser (str + 1, &str, 0);
686c4bc2 310 if ((new == NULL) || ((new != ERROR_OP) && (new->prio == 9))) {
67b13e3d
LM
311 delelement (new);
312 new = ERROR_OP;
313 }
686c4bc2 314 if ((new == NULL) || (new == ERROR_OP) || (*str == ',')) {
031d7bba
LM
315 delelement (new);
316 delelement (root);
efdfb543
LM
317 return ERROR_OP;
318 }
85b4a72c 319 root->ops[0] = new;
0b489a77 320 }
3c0db5bc
LM
321 if (*str != ')') {
322 delelement (root);
323 return ERROR_OP;
324 }
85b4a72c 325 str++;
04d68907 326 VERBOSE (DEBUG, fprintf (stdout, "stop processing bracket\n"));
efdfb543 327 continue;
49223129 328 }
85b4a72c 329
124da7fd 330 /* check for closing bracket, closing brace or koma */
85b4a72c 331
124da7fd 332 if ((*str == ')') || (*str == '}') || (*str == ',')) {
45a631a8
LM
333 if (prio == -9) {
334 delelement (root);
335 return ERROR_OP;
336 }
49223129 337 if (next != NULL) {
85b4a72c 338 *next = str;
49223129
LM
339 }
340 return root;
0b489a77
LM
341 }
342
bc97a989
LM
343 /* look for operators */
344
345 for (i = 0; i < NB_OPERATORS; i++) {
346 keyword_t *operator = operators + i;
347 if (codecmp (operator->keyword, str) == 0) {
04d68907 348 VERBOSE (DEBUG, fprintf (stdout, "start processing operator\n"));
11cda8d7 349 if (root) {
ef37d966 350 if ((prio) && (prio > operator->prio)) {
04d68907 351 VERBOSE (DEBUG, fprintf (stdout, "stop because operator priority\n"));
ef37d966
LM
352 *next = str;
353 return root;
354 }
355 str += operator->offset;
04d68907 356 VERBOSE (INFO, fprintf (stdout, "Oper: %d\n", operator->func));
d2ff8478 357 if (subparser (&root, &str, operator->func, operator->nbops, operator->prio) == ERROR_OP) {
031d7bba 358 delelement (root);
ef37d966 359 return ERROR_OP;
11cda8d7 360 }
3b4b0bbe 361 } else if (*str == '-') {
d32ea0b0 362 root = newelement (Sig, 1, 9);
49223129 363 } else {
0b489a77 364 return ERROR_OP;
bc97a989 365 }
bc97a989 366 found = 1;
04d68907 367 VERBOSE (DEBUG, fprintf (stdout, "stop processing operator\n"));
bc97a989
LM
368 break;
369 }
370 }
371 if (found) {
bc97a989
LM
372 continue;
373 }
49223129 374
bc97a989
LM
375 /* look for functions */
376
377 for (i = 0; i < NB_FUNCTIONS; i++) {
378 keyword_t *function = functions + i;
379 if (codecmp (function->keyword, str) == 0) {
04d68907 380 VERBOSE (DEBUG, fprintf (stdout, "start processing function\n"));
bc97a989 381 if (root == NULL) {
04d68907 382 VERBOSE (INFO, fprintf (stdout, "Func: %d\n", function->func));
d32ea0b0 383 root = newelement (function->func, function->nbops, function->prio);
031d7bba
LM
384 } else {
385 delelement (root);
0b489a77 386 return ERROR_OP;
bc97a989
LM
387 }
388 str += function->offset;
389 found = 1;
04d68907 390 VERBOSE (DEBUG, fprintf (stdout, "stop processing function\n"));
bc97a989
LM
391 break;
392 }
393 }
394 if (found) {
bc97a989
LM
395 continue;
396 }
397
471da7c9
LM
398 /* look for constant */
399
400 for (i = 0; i < NB_CONSTANTS; i++) {
401 keyword_t *constant = constants + i;
402 if (codecmp (constant->keyword, str) == 0) {
403 VERBOSE (DEBUG, fprintf (stdout, "start processing constant\n"));
404 if (root == NULL) {
5075f6ea 405 VERBOSE (INFO, fprintf (stdout, "Const: %d\n", constant->func));
d32ea0b0 406 root = newelement (constant->func, constant->nbops, constant->prio);
471da7c9
LM
407 } else {
408 delelement (root);
409 return ERROR_OP;
410 }
411 str += constant->offset;
412 found = 1;
413 VERBOSE (DEBUG, fprintf (stdout, "stop processing constant\n"));
414 break;
415 }
416 }
417 if (found) {
418 continue;
419 }
420
bc97a989
LM
421 /* look for number */
422
0b9cc9b0
LM
423 if (((*str >= '0') && (*str <= '9')) ||
424 (*str == '.') || (*str == '+') || (*str == '-')) {
04d68907 425 VERBOSE (DEBUG, fprintf (stdout, "start processing value\n"));
bc97a989 426 char *pt;
fd88e359 427 double value = strtod (str, &pt);
04d68907 428 VERBOSE (INFO, fprintf (stdout, "Value: %f\n", value));
bc97a989 429 if (str != pt) {
ce6627f2 430 if ((root == NULL) || (root->prio == 6)) {
3b4b0bbe 431 new = newelement (Val, 1, 5);
c47a9298 432 new->value = value;
ce6627f2
LM
433 if (root == NULL) {
434 root = new;
435 } else {
436 for (i = 0; i < root->nbops; i++) {
437 if (root->ops[i] == NULL) {
438 root->ops[i] = new;
439 found = 1;
440 break;
441 }
442 }
443 if (!found) {
444 delelement (new);
445 delelement (root);
446 return ERROR_OP;
447 }
448 }
c47a9298 449 str = pt;
32741902
LM
450 } else if ((*str == '+') || (*str == '-')) {
451 if ((prio) && (prio > 1)) {
452 VERBOSE (DEBUG, fprintf (stdout, "stop because operator priority\n"));
453 *next = str;
454 return root;
455 }
456 if (subparser (&root, &str, Add, 2, 1) == ERROR_OP) {
031d7bba 457 delelement (root);
0b489a77 458 return ERROR_OP;
bc97a989 459 }
0b9cc9b0 460 } else {
031d7bba 461 delelement (root);
0b9cc9b0 462 return ERROR_OP;
bc97a989 463 }
bc97a989
LM
464 found = 1;
465 }
04d68907 466 VERBOSE (DEBUG, fprintf (stdout, "stop processing value\n"));
bc97a989
LM
467 }
468
469 /* error */
470
471 if (!found) {
031d7bba 472 delelement (root);
0b489a77 473 return ERROR_OP;
bc97a989
LM
474 }
475
476 }
477
49223129
LM
478 if (next != NULL) {
479 *next = str;
480 }
031d7bba 481
779282bb 482 /* save string */
0a8a1de3
ML
483 if (root != NULL) {
484 root->string = string;
485 }
779282bb 486
bc97a989
LM
487 return root;
488}
489
490/* print element tree */
491
492void print_element (element_t *root, int level)
493{
494 char *func = NULL;
495 int i;
496
49223129 497 if ((root == NULL) || (root == ERROR_OP)) {
bc97a989 498 return;
49223129
LM
499 }
500
bc97a989 501 for (i = 0; i < level; i++) {
04d68907 502 fprintf (stdout, " ");
bc97a989
LM
503 }
504
505 switch (root->func) {
506 case Val: func = "Value"; break;
89cf0955 507 case Sig: func = "Sign"; break;
bc97a989
LM
508 case Add: func = "Addition"; break;
509 case Sub: func = "Subtraction"; break;
510 case Mul: func = "Multiplication"; break;
511 case Div: func = "Division"; break;
c47a9298 512 case Mod: func = "Modulo"; break;
bc97a989
LM
513 case Pow: func = "Power"; break;
514 case Sqr: func = "Square Root"; break;
515 case Cos: func = "Cosine"; break;
516 case Sin: func = "Sine"; break;
9e52bbc5
LM
517 case Tan: func = "Tangent"; break;
518 case Acos: func = "Arc Cosine"; break;
519 case Asin: func = "Arc Sine"; break;
471da7c9 520 case Atan: func = "Arc Tangent"; break;
bc97a989
LM
521 case Log: func = "Logarithm"; break;
522 case Exp: func = "Exponantial"; break;
9e52bbc5
LM
523 case Abs: func = "Absolute value"; break;
524 case Ceil: func = "Ceil value"; break;
525 case Floor: func = "Floor value"; break;
6ba1dd0f
LM
526 case Store: func = "Store"; break;
527 case Recall: func = "Recall"; break;
a8cf32ba
LM
528 case Inc: func = "Increase"; break;
529 case Dec: func = "Decrease"; break;
6ba1dd0f 530 case Disp: func = "Display"; break;
f438579a 531 case Mem: func = "Memory"; break;
143f91ae 532 case Clear: func = "Clear"; break;
471da7c9
LM
533 case Quit: func = "Quit"; break;
534 case Help: func = "Help"; break;
5075f6ea 535 case Ans: func = "Ans"; break;
471da7c9
LM
536 case Pi: func = "Pi"; break;
537 case E: func = "E"; break;
7fe742c9
LM
538 case Equal: func = "Equal"; break;
539 case Diff: func = "Different"; break;
540 case Ge: func = "Greater or equal"; break;
541 case Le: func = "Lesser or equal"; break;
542 case Gt: func = "Greater"; break;
543 case Lt: func = "Lesser"; break;
ca3e2a2f
LM
544 case And: func = "And"; break;
545 case Or: func = "Or"; break;
546 case Not: func = "Not"; break;
94b4e517 547 case Cond: func = "Condition"; break;
2a688642 548 case While: func = "While"; break;
05aabb97 549 case Code: func = "Code"; break;
715580ff 550 case Print: func = "Print"; break;
c7cbb833
LM
551 case Prog: func = "Program"; break;
552 case Call: func = "Call"; break;
553 case List: func = "List"; break;
554 case Edit: func = "Edit"; break;
555 case Del: func = "Del"; break;
bc97a989
LM
556 }
557
04d68907 558 fprintf (stdout, "Function: %s\n", func);
bc97a989 559
85b4a72c 560 if ((root->func == Val) && (root->ops[0] == NULL)) {
bc97a989 561 for (i = 0; i < level; i++) {
04d68907 562 fprintf (stdout, " ");
bc97a989 563 }
04d68907 564 fprintf (stdout, "value: %f\n", root->value);
bc97a989
LM
565 } else {
566 for (i = 0; i < root->nbops; i++) {
567 print_element (root->ops[i], level + 1);
568 }
569 }
570}
571
6ba1dd0f
LM
572/* storage functions */
573
f438579a
LM
574void memory (int nb)
575{
576 int i, l;
577 double *tmp = NULL;
578 if (nb != storage_size) {
579 l = (nb < storage_size) ? nb : storage_size;
6f14d1cc 580 tmp = (double *) callocordie (nb, sizeof (double));
f438579a
LM
581 for (i = 0; i < l; i++) {
582 tmp[i] = storage[i];
583 }
584 if (storage != NULL) {
585 free (storage);
586 }
587 storage = tmp;
588 storage_size = nb;
589 }
590}
591
3639df66 592double store (int index, double value)
6ba1dd0f 593{
f438579a
LM
594 if (storage_size == -1) {
595 memory (DEFAULT_STORAGE_SIZE);
596 }
597 if ((index > 0) && (index <= storage_size)) {
3639df66 598 storage[index - 1] = value;
6ba1dd0f 599 } else {
9b759869 600 VERBOSE (WARNING, fprintf (stdout, "invalid index (%d) [%d, %d]\n", index, (storage_size) ? 1 : 0, storage_size));
6ba1dd0f 601 }
3639df66 602 return value;
6ba1dd0f
LM
603}
604
605double recall (int index)
606{
f438579a
LM
607 if (storage_size == -1) {
608 memory (DEFAULT_STORAGE_SIZE);
609 }
610 if ((index > 0) && (index <= storage_size)) {
6ba1dd0f
LM
611 return storage[index - 1];
612 } else {
9b759869 613 VERBOSE (WARNING, fprintf (stdout, "invalid index (%d) [%d, %d]\n", index, (storage_size) ? 1 : 0, storage_size));
6ba1dd0f
LM
614 }
615 return 0;
616}
617
a8cf32ba
LM
618double increase (int index)
619{
f438579a
LM
620 if (storage_size == -1) {
621 memory (DEFAULT_STORAGE_SIZE);
622 }
623 if ((index > 0) && (index <= storage_size)) {
a8cf32ba
LM
624 return storage[index - 1]++;
625 } else {
9b759869 626 VERBOSE (WARNING, fprintf (stdout, "invalid index (%d) [%d, %d]\n", index, (storage_size) ? 1 : 0, storage_size));
a8cf32ba
LM
627 }
628 return 0;
629}
630
631double decrease (int index)
632{
f438579a
LM
633 if (storage_size == -1) {
634 memory (DEFAULT_STORAGE_SIZE);
635 }
636 if ((index > 0) && (index <= storage_size)) {
a8cf32ba
LM
637 return storage[index - 1]--;
638 } else {
9b759869 639 VERBOSE (WARNING, fprintf (stdout, "invalid index (%d) [%d, %d]\n", index, (storage_size) ? 1 : 0, storage_size));
a8cf32ba
LM
640 }
641 return 0;
642}
643
6ba1dd0f
LM
644void display (void)
645{
646 int i;
f438579a
LM
647 if (storage_size == -1) {
648 memory (DEFAULT_STORAGE_SIZE);
649 }
6ba1dd0f 650 fprintf (stdout, "storage:");
f438579a 651 for (i = 0; i < storage_size; i++) {
9b759869
LM
652 fprintf (stdout, " ");
653 fprintf (stdout, minform, storage[i]);
6ba1dd0f
LM
654 }
655 fprintf (stdout, "\n");
656}
657
143f91ae
LM
658void clear ()
659{
660 int i;
661 for (i = 0; i < storage_size; i++) {
662 storage[i] = 0;
663 }
664}
665
2a688642
LM
666/* While do function */
667
668double while_do (element_t *cond, element_t *action)
669{
670 double ret = 0;
671 element_t *temp = NULL;
672
673 VERBOSE (DEBUG, fprintf (stdout, "starting while loop\n"));
2a688642
LM
674 while (1) {
675 VERBOSE (DEBUG, fprintf (stdout, "loop...\n"));
676
677 temp = dupelement (cond);
678 if (!evaluate_element (temp, 0)) {
679 break;
680 }
681 if (action) {
682 temp = dupelement (action);
683 ret = evaluate_element (temp, 0);
684 }
685 }
686
687 VERBOSE (DEBUG, fprintf (stdout, "ending while loop\n"));
688
689 return ret;
690}
691
124da7fd
LM
692/* program function */
693
05aabb97 694double execute_code (element_t **prog, int nbcalls)
124da7fd
LM
695{
696 double ret = 0;
697 int i;
698 for (i = 0; i < nbcalls; i++) {
699 ret = evaluate_element (prog[i], 0);
700 prog[i] = NULL;
701 }
702 return ret;
703}
704
2a5ec9d1
LM
705/* print function */
706
707void set_format (char *prompt, int precision)
708{
709 char buffer[128] = {0};
2a5ec9d1 710 free_format ();
9b759869 711 sprintf (buffer, "%s%%.%dg\n", prompt, precision);
2a5ec9d1 712 format = strdup (buffer);
9b759869
LM
713 sprintf (buffer, "%%.%dg", precision);
714 minform = strdup (buffer);
2a5ec9d1
LM
715}
716
717void free_format ()
718{
719 if (format) {
720 free (format);
721 format = NULL;
722 }
9b759869
LM
723 if (minform) {
724 free (minform);
725 minform = NULL;
726 }
2a5ec9d1
LM
727}
728
729double print (double value)
730{
f438579a 731 fprintf (stdout, format ? format : DEFAULT_FORMAT, value);
2a5ec9d1
LM
732 fflush (stdout);
733 return value;
734}
735
89cf0955
LM
736/* quit function */
737
738void quit (void)
739{
04d68907 740 fprintf (stdout, "bye\n");
89cf0955
LM
741 exit (0);
742}
743
c7cbb833
LM
744/* program function */
745
746void prog (int id, int nbmems, element_t *root)
05aabb97
LM
747{
748 int i, n = -1;
749
750 if (programs == NULL) {
751
752 /* initial memory allocation */
6f14d1cc 753 programs = (workspace_t *) callocordie (1, sizeof (workspace_t));
05aabb97
LM
754 nb_programs = 1;
755 n = 0;
756
757 } else {
758
759 /* look for existing program */
760 for (i = 0; i < nb_programs; i++) {
761 if ((programs + i)->id == id) {
762 n = i;
763 break;
764 }
765 }
766 if (n == -1) {
767
768 /* new program */
769 n = nb_programs++;
6f14d1cc
LM
770 workspace_t *tmp = (workspace_t *) callocordie (nb_programs, sizeof (workspace_t));
771 memcpy (tmp, programs, (nb_programs - 1) * sizeof (workspace_t));
772 free (programs);
773 programs = tmp;
05aabb97
LM
774 } else {
775
776 /* clean old program */
777 if ((programs + n)->storage) {
778 free ((programs + n)->storage);
779 }
780 if ((programs + n)->root) {
781 delelement ((programs + n)->root);
782 }
2e34899a
LM
783 if ((programs + n)->string) {
784 free ((programs + n)->string);
0a8a1de3 785 (programs + n)->string = NULL;
2e34899a 786 }
05aabb97
LM
787 }
788 }
789
790 /* set program */
791 (programs + n)->id = id;
792 (programs + n)->answer = 0;
6f14d1cc 793 (programs + n)->storage = (double *) callocordie (nbmems, sizeof (double));
05aabb97 794 (programs + n)->storage_size = nbmems;
e7d5c144 795 (programs + n)->root = dupelement (root);
05aabb97 796}
c7cbb833 797
e7d5c144 798double call (int id, int nbops, element_t **ops)
c7cbb833 799{
05aabb97
LM
800 workspace_t tmp = {0};
801 int i, n = -1;
802 double ret = 0;
803
804 if (programs) {
e7d5c144
LM
805
806 /* look for program */
05aabb97 807 for (i = 0; i < nb_programs; i++) {
e7d5c144 808 if ((programs + i)->id == id) {
05aabb97
LM
809 n = i;
810 break;
811 }
812 }
d59a8079
LM
813 }
814 if (n == -1) {
815 VERBOSE (WARNING, fprintf (stdout, "error unknown program (%d)\n", id));
816 return 0;
817 }
05aabb97 818
d59a8079
LM
819 /* store context */
820 tmp.answer = answer;
821 tmp.storage = storage;
822 tmp.storage_size = storage_size;
05aabb97 823
d59a8079
LM
824 /* change context */
825 answer = 0;
826 storage = (programs + n)->storage;
827 storage_size = (programs + n)->storage_size;
828 if (nbops > storage_size) {
829 double *tmp = (double *) callocordie (nbops, sizeof (double));
830 memcpy (tmp, storage, storage_size * sizeof (double));
831 free (storage);
832 storage = tmp;
833 storage_size = nbops;
05aabb97 834 }
d59a8079
LM
835 for (i = 0; i < nbops; i++) {
836 double val = evaluate_element (ops[i], 0);
837 store (i + 1, val);
838 }
839
840 /* evaluate program */
841 element_t *elements = dupelement ((programs + n)->root);
842 ret = evaluate_element (elements, 0);
843
844 /* restore context */
845 answer = tmp.answer;
846 storage = tmp.storage;
847 storage_size = tmp.storage_size;
05aabb97
LM
848
849 return ret;
c7cbb833
LM
850}
851
852void list ()
05aabb97
LM
853{
854 int i;
d59a8079
LM
855 fprintf (stdout, "programs:");
856 for (i = 0; i < nb_programs; i++) {
857 fprintf (stdout, " %d", (programs + i)->id);
05aabb97 858 }
d59a8079 859 fprintf (stdout, "\n");
05aabb97 860}
c7cbb833
LM
861
862void edit (int id)
2e34899a
LM
863{
864 int i, n = -1;
865
866 if (programs) {
867
868 /* look for program */
869 for (i = 0; i < nb_programs; i++) {
870 if ((programs + i)->id == id) {
871 n = i;
872 break;
873 }
874 }
875 }
876 if (n == -1) {
877 VERBOSE (WARNING, fprintf (stdout, "error unknown program (%d)\n", id));
878 return;
879 }
880
881 /* set string program */
779282bb 882 fprintf (stdout, "edit: %s\n", (programs + n)->string);
2e34899a
LM
883}
884
779282bb 885void savestring (int id, char *string)
2e34899a
LM
886{
887 int i, n = -1;
888
889 if (programs) {
890
891 /* look for program */
892 for (i = 0; i < nb_programs; i++) {
893 if ((programs + i)->id == id) {
894 n = i;
895 break;
896 }
897 }
898 }
2e34899a 899
1b068216
LM
900 /* unnecesary code */
901 //if (n == -1) {
902 // VERBOSE (WARNING, fprintf (stdout, "error unknown program (%d)\n", id));
903 // return;
904 //}
905 //if ((programs + n)->string) {
906 // free ((programs + n)->string);
907 //}
908
2e34899a
LM
909 (programs + n)->string = strdup (string);
910}
c7cbb833
LM
911
912void del (int id)
d59a8079
LM
913{
914 int i, j, n = -1;
915
916 if (programs) {
917
918 /* look for program */
919 for (i = 0; i < nb_programs; i++) {
920 if ((programs + i)->id == id) {
921 n = i;
922 break;
923 }
924 }
925 }
926 if (n == -1) {
927 VERBOSE (WARNING, fprintf (stdout, "error unknown program (%d)\n", id));
928 return;
929 }
930
931 /* clean program */
932 if ((programs + n)->storage) {
933 free ((programs + n)->storage);
934 }
935 if ((programs + n)->root) {
936 delelement ((programs + n)->root);
937 }
2e34899a
LM
938 if ((programs + n)->string) {
939 free ((programs + n)->string);
940 }
d59a8079
LM
941
942 /* remove entry */
943 workspace_t *tmp = (workspace_t *) callocordie (nb_programs - 1, sizeof (workspace_t));
944 for (i = 0, j = 0; i < nb_programs; i++) {
945 if (i != n) {
946 memcpy (tmp + j, programs + i, sizeof (workspace_t));
947 j++;
948 }
949 }
950 free (programs);
951 programs = tmp;
952 nb_programs--;
953}
c7cbb833 954
89cf0955
LM
955/* help message */
956
957void help (void)
958{
04d68907 959 fprintf (stdout, "calc is a simple calculator\n\n");
94b4e517
LM
960 fprintf (stdout, "supported operators:");
961 fprintf (stdout, " + - * / %% ^\n");
962 fprintf (stdout, "camparison operators:");
963 fprintf (stdout, " == != >= <= > <\n");
964 fprintf (stdout, "logical operators:");
965 fprintf (stdout, " & | !\n");
9e52bbc5 966 fprintf (stdout, "mathematic functions:");
b9c1c40d 967 fprintf (stdout, " pow sqrt log exp\n");
9049369f 968 fprintf (stdout, "trig. func.:");
b9c1c40d 969 fprintf (stdout, " cos sin tan acos asin atan\n");
94b4e517 970 fprintf (stdout, "supported functions:");
9e52bbc5 971 fprintf (stdout, " abs ceil floor\n");
94b4e517 972 fprintf (stdout, "storage functions:");
f438579a 973 fprintf (stdout, " mem sto rcl inc dec disp\n");
9049369f 974 fprintf (stdout, "prog. functions:");
2d0cd54c 975 fprintf (stdout, " cond while print {} ;\n");
9049369f 976 fprintf (stdout, "misc. functions:");
94b4e517
LM
977 fprintf (stdout, " quit help\n");
978 fprintf (stdout, "supported constants:");
471da7c9 979 fprintf (stdout, " e pi\n");
89cf0955 980}
3b4b0bbe 981
f2927108
LM
982/* evaluate element tree */
983
89cf0955
LM
984#define MASK_SUB 0x1
985#define MASK_DIV 0x2
986
3b4b0bbe 987double evaluate_element (element_t *root, char mask)
f2927108
LM
988{
989 double op0 = 0, op1 = 0;
3b4b0bbe 990 char nextmask = mask;
e7d5c144 991 int i, nb;
f2927108
LM
992
993 if ((root == NULL) || (root == ERROR_OP)) {
04d68907 994 VERBOSE (WARNING, fprintf (stdout, "error while evaluating\n"));
f2927108
LM
995 return 0;
996 }
997
3b4b0bbe
LM
998 /* mask to manage sub operator sub and div */
999 switch (root->func) {
1000 case Add:
1001 nextmask &= ~MASK_SUB;
1002 nextmask &= ~MASK_DIV;
1003 break;
1004 case Sub:
1005 nextmask |= MASK_SUB;
1006 nextmask &= ~MASK_DIV;
1007 break;
1008 case Mul:
1009 nextmask &= ~MASK_DIV;
1010 break;
1011 case Div:
1012 nextmask |= MASK_DIV;
1013 break;
1014 default:
1015 nextmask = mask;
1016 }
1017
f2927108
LM
1018 switch (root->func) {
1019 case Val:
3b4b0bbe
LM
1020 case Sig:
1021 op0 = (root->ops[0]) ? evaluate_element (root->ops[0], nextmask) : root->value;
f2927108
LM
1022 break;
1023 case Add:
1024 case Sub:
1025 case Mul:
1026 case Div:
c47a9298 1027 case Mod:
f2927108 1028 case Pow:
3639df66 1029 case Store:
7fe742c9
LM
1030 case Equal:
1031 case Diff:
1032 case Ge:
1033 case Le:
1034 case Gt:
1035 case Lt:
ca3e2a2f
LM
1036 case And:
1037 case Or:
c7cbb833 1038 case Prog:
f2927108 1039 if (root->ops[1]) {
3b4b0bbe 1040 op1 = evaluate_element (root->ops[1], nextmask);
3639df66 1041 } else if (root->func != Store) {
04d68907 1042 VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[1])\n"));
f2927108
LM
1043 return 0;
1044 }
1045 /* fallthrough */
1046 case Sqr:
1047 case Cos:
1048 case Sin:
9e52bbc5
LM
1049 case Tan:
1050 case Acos:
1051 case Asin:
471da7c9 1052 case Atan:
f2927108
LM
1053 case Log:
1054 case Exp:
9e52bbc5
LM
1055 case Abs:
1056 case Ceil:
1057 case Floor:
6ba1dd0f 1058 case Recall:
a8cf32ba
LM
1059 case Inc:
1060 case Dec:
ca3e2a2f 1061 case Not:
f438579a 1062 case Mem:
94b4e517 1063 case Cond:
c7cbb833 1064 case Call:
c7cbb833
LM
1065 case Edit:
1066 case Del:
f2927108 1067 if (root->ops[0]) {
3b4b0bbe 1068 op0 = evaluate_element (root->ops[0], 0);
f2927108 1069 } else {
04d68907 1070 VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[0])\n"));
f2927108
LM
1071 return 0;
1072 }
89cf0955 1073 break;
6ba1dd0f 1074 case Disp:
143f91ae 1075 case Clear:
471da7c9
LM
1076 case Quit:
1077 case Help:
5075f6ea 1078 case Ans:
471da7c9
LM
1079 case Pi:
1080 case E:
05aabb97 1081 case Code:
e7d5c144 1082 case List:
09d87cae 1083 break;
2a688642 1084 case While:
09d87cae
LM
1085 if (root->ops[0] == NULL) {
1086 VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[0])\n"));
1087 return 0;
1088 }
89cf0955 1089 break;
715580ff
LM
1090 case Print:
1091 op0 = (root->ops[0]) ? evaluate_element (root->ops[0], 0) : answer;
1092 break;
f2927108
LM
1093 }
1094
1095 switch (root->func) {
3b4b0bbe
LM
1096 case Val: return op0;
1097 case Sig: return -op0;
1098 case Add: return ((mask & MASK_SUB) == 0) ? op0 + op1 : op0 - op1;
1099 case Sub: return ((mask & MASK_SUB) == 0) ? op0 - op1 : op0 + op1;
1100 case Mul: return ((mask & MASK_DIV) == 0) ? op0 * op1 : op0 / op1;
1101 case Div: return ((mask & MASK_DIV) == 0) ? op0 / op1 : op0 * op1;
c47a9298 1102 case Mod: return fmod (op0, op1);
f2927108
LM
1103 case Pow: return pow (op0, op1);
1104 case Sqr: return sqrt (op0);
1105 case Cos: return cos (op0);
1106 case Sin: return sin (op0);
9e52bbc5
LM
1107 case Tan: return tan (op0);
1108 case Acos: return acos (op0);
1109 case Asin: return asin (op0);
471da7c9 1110 case Atan: return atan (op0);
f2927108
LM
1111 case Log: return log (op0);
1112 case Exp: return exp (op0);
9e52bbc5
LM
1113 case Abs: return fabs (op0);
1114 case Ceil: return ceil (op0);
1115 case Floor: return floor (op0);
3639df66 1116 case Store: return store ((int)op0, (op1) ? op1 : answer);
6ba1dd0f 1117 case Recall: return recall ((int)op0);
a8cf32ba
LM
1118 case Inc: return increase ((int)op0);
1119 case Dec: return decrease ((int)op0);
6ba1dd0f 1120 case Disp: display (); break;
f438579a 1121 case Mem: memory ((int)op0); break;
143f91ae 1122 case Clear: clear (); break;
471da7c9
LM
1123 case Quit: quit (); break;
1124 case Help: help (); break;
5075f6ea 1125 case Ans: return answer;
471da7c9
LM
1126 case Pi: return M_PI;
1127 case E: return M_E;
7fe742c9
LM
1128 case Equal: return op0 == op1;
1129 case Diff: return op0 != op1;
1130 case Ge: return op0 >= op1;
1131 case Le: return op0 <= op1;
1132 case Gt: return op0 > op1;
1133 case Lt: return op0 < op1;
ca3e2a2f
LM
1134 case And: return (op0 != 0) && (op1 != 0);
1135 case Or: return (op0 != 0) || (op1 != 0);
1136 case Not: return (op0 == 0);
94b4e517
LM
1137 case Cond:
1138 if ((op0) && (root->ops[1])) {
1139 return evaluate_element (root->ops[1], 0);
1140 } else if ((!op0) && (root->ops[2])) {
1141 return evaluate_element (root->ops[2], 0);
1142 } else {
1143 return 0;
1144 }
2a688642 1145 case While: return while_do (root->ops[0], root->ops[1]);
05aabb97 1146 case Code: return execute_code (root->ops, root->nbops);
2a5ec9d1 1147 case Print: return print (op0);
779282bb
LM
1148 case Prog:
1149 prog ((int)op0, (int)op1, root->ops[2]);
1150 savestring ((int)op0, root->string);
1151 break;
e7d5c144
LM
1152 case Call:
1153 for (i = 1, nb =0; i < root->nbops; i++) {
1154 if (root->ops[i]) {
1155 nb++;
1156 }
1157 }
1158 return call ((int)op0, nb, root->ops + 1);
c7cbb833
LM
1159 case List: list (); break;
1160 case Edit: edit ((int)op0); break;
05aabb97 1161 case Del: del ((int)op0); break;
f2927108
LM
1162 }
1163
1164 return 0;
1165}
1166
b9c1c40d
LM
1167char **generate_completion_list ()
1168{
db660e60 1169 int i, j, l = 0;
6f14d1cc 1170 char **list = (char **) callocordie (NB_OPERATORS + NB_FUNCTIONS + NB_CONSTANTS + NB_SYMBOLS + 1, sizeof (char *));
b9c1c40d 1171
db660e60
LM
1172 for (i = 0; i < NB_OPERATORS; i++) {
1173 list[l] = strdup ((operators + i)->keyword);
1174 for (j = 0; j < (int)strlen (list[l]); j++) {
1175 if (list[i][j] == '\t') {
1176 list[i][j] = '\0';
1177 }
1178 }
d8f4df0f
LM
1179 if (list[l] != NULL) {
1180 l++;
db660e60 1181 }
db660e60
LM
1182 }
1183
b9c1c40d
LM
1184 for (i = 0; i < NB_FUNCTIONS; i++) {
1185 list[l] = strdup ((functions + i)->keyword);
d8f4df0f
LM
1186 if (list[l] != NULL) {
1187 l++;
b9c1c40d 1188 }
b9c1c40d
LM
1189 }
1190
1191 for (i = 0; i < NB_CONSTANTS; i++) {
1192 list[l] = strdup ((constants + i)->keyword);
d8f4df0f
LM
1193 if (list[l] != NULL) {
1194 l++;
b9c1c40d 1195 }
b9c1c40d
LM
1196 }
1197
db660e60
LM
1198 for (i = 0; i < NB_SYMBOLS; i++) {
1199 list[l] = strdup (symbols[i]);
d8f4df0f
LM
1200 if (list[l] != NULL) {
1201 l++;
db660e60 1202 }
db660e60
LM
1203 }
1204
b9c1c40d
LM
1205 return (list);
1206}
1207
1208void free_completion_list (char **list)
1209{
1210 int i;
1211
1212 if (list) {
1213 for (i = 0; i < NB_FUNCTIONS + NB_CONSTANTS; i++) {
1214 if (list[i] != NULL) {
1215 free (list[i]);
1216 }
1217 }
1218 free (list);
1219 }
1220}
1221
1222
bc97a989 1223/* vim: set ts=4 sw=4 et: */