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