9796048af19e5f8d552a682fef0c1afc564fd799
10 /* global variables */
14 #define STORAGE_SIZE 10
15 double storage
[STORAGE_SIZE
] = {0};
19 int codecmp (char *ref
, char *str
)
23 while (*ref
!= '\0') {
25 sig
= (*str
== '.') ? -1 : ((*str
>= '0') && (*str
<= '9'));
30 return (sig
> 0) ? 1 : -1;
39 /* allocate new element */
41 element_t
*newelement (func_t function
, int nbops
, int prio
)
43 element_t
*new = (element_t
*) calloc (1, sizeof (element_t
));
45 VERBOSE (ERROR
, fprintf (stderr
, "can't allocate memory\n"));
49 new->ops
= (element_t
**) calloc (nbops
, sizeof (element_t
*));
50 if (new->ops
== NULL
) {
51 VERBOSE (ERROR
, fprintf (stderr
, "can't allocate memory\n"));
62 /* desallocate element */
64 void delelement (element_t
*root
)
66 if ((root
!= NULL
) && (root
!= ERROR_OP
)) {
68 for (i
= 0; i
< root
->nbops
; i
++) {
69 if ((root
->ops
[i
] != NULL
) && (root
->ops
[i
] != ERROR_OP
)) {
70 delelement (root
->ops
[i
]);
80 /* duplicate element */
82 element_t
*dupelement (element_t
*root
)
84 element_t
*tmp
= NULL
;
87 if ((root
== NULL
) || (root
== ERROR_OP
)) {
90 tmp
= newelement (root
->func
, root
->nbops
, root
->prio
);
91 tmp
->value
= root
->value
;
92 for (i
= 0; i
< root
->nbops
; i
++) {
93 tmp
->ops
[i
] = dupelement (root
->ops
[i
]);
100 #define NB_OPERATORS 14
101 keyword_t operators
[NB_OPERATORS
] = {
102 { "+\t", Add
, 2, 1, 1},
103 { "-\t", Sub
, 2, 1, 1},
104 { "*", Mul
, 2, 1, 2},
105 { "/", Div
, 2, 1, 2},
106 { "%", Mod
, 2, 1, 3},
107 { "^", Pow
, 2, 1, 4},
108 { "==", Equal
, 2, 2, -1},
109 { "!=", Diff
, 2, 2, -1},
110 { ">=", Ge
, 2, 2, -1},
111 { "<=", Le
, 2, 2, -1},
112 { ">", Gt
, 2, 1, -1},
113 { "<", Lt
, 2, 1, -1},
114 { "&", And
, 2, 1, -2},
118 #define NB_FUNCTIONS 17
119 keyword_t functions
[NB_FUNCTIONS
] = {
120 { "sqrt", Sqr
, 1, 4, 5},
121 { "pow", Pow
, 2, 3, 5},
122 { "cos", Cos
, 1, 3, 5},
123 { "sin", Sin
, 1, 3, 5},
124 { "atan", Atan
, 1, 4, 5},
125 { "exp", Exp
, 1, 3, 5},
126 { "log", Log
, 1, 3, 5},
127 { "sto", Store
, 2, 3, 5},
128 { "rcl", Recall
, 1, 3, 5},
129 { "inc", Inc
, 1, 3, 5},
130 { "dec", Dec
, 1, 3, 5},
131 { "disp", Disp
, 0, 4, 9},
132 { "quit", Quit
, 0, 4, 9},
133 { "help", Help
, 0, 4, 9},
134 { "!", Not
, 1, 1, 6},
135 { "cond", Cond
, 3, 4, 5},
136 { "whl", While
, 2, 3, 5}
139 #define NB_CONSTANTS 3
140 keyword_t constants
[NB_CONSTANTS
] = {
141 { "ans", Ans
, 0, 3, 5},
146 /* subparser function */
148 element_t
*subparser (element_t
**proot
, char **pstr
, func_t func
, int nbops
, int prio
)
150 element_t
*new = newelement (func
, nbops
, prio
);
151 new->ops
[0] = *proot
;
152 new->ops
[1] = parser (*pstr
, pstr
, new->prio
);
153 if ((new->ops
[1] == NULL
) || ((new->ops
[1] != ERROR_OP
) && (new->ops
[1]->prio
== 9))) {
154 delelement (new->ops
[1]);
155 new->ops
[1] = ERROR_OP
;
157 if (new->ops
[1] == ERROR_OP
) {
162 *proot
= newelement (Val
, 1, 5);
163 (*proot
)->ops
[0] = new;
168 /* parser function */
170 element_t
*parser (char *str
, char **next
, int prio
)
172 element_t
*root
= NULL
;
175 VERBOSE (DEBUG
, fprintf (stdout
, "Starting parsing\n"));
178 while (*str
!= '\0') {
180 element_t
*new = NULL
;
181 VERBOSE (INFO
, fprintf (stdout
, "Processing: %s\n", str
));
183 /* skip spaces and tabs */
185 if ((*str
== ' ') || (*str
== '\t')) {
190 /* check for open brace */
193 VERBOSE (DEBUG
, fprintf (stdout
, "start processing brace\n"));
198 element_t
**prog
= NULL
;
199 new = newelement (Prog
, 0, 5);
204 new = parser (str
+ 1, &str
, 0);
205 if ((new == NULL
) || ((new != ERROR_OP
) && (new->prio
== 9))) {
209 if (new == ERROR_OP
) {
213 element_t
*newprog
= newelement (Prog
, (*prog
)->nbops
+ 1, 5);
214 for (i
= 0; i
< (*prog
)->nbops
; i
++) {
215 newprog
->ops
[i
] = (*prog
)->ops
[i
];
216 (*prog
)->ops
[i
] = NULL
;
218 newprog
->ops
[(*prog
)->nbops
] = new;
221 } while (*str
== ',');
228 VERBOSE (DEBUG
, fprintf (stdout
, "stop processing brace\n"));
232 /* check for open bracket */
235 VERBOSE (DEBUG
, fprintf (stdout
, "start processing bracket\n"));
239 new = parser (str
+ 1, &str
, 0);
240 if ((new == NULL
) || ((new != ERROR_OP
) && (new->prio
== 9))) {
244 if ((new == NULL
) || (new == ERROR_OP
)) {
248 for (i
= 0; i
< root
->nbops
; i
++) {
249 if (root
->ops
[i
] == NULL
) {
260 } while (*str
== ',');
262 root
= newelement (Val
, 1, 5);
263 new = parser (str
+ 1, &str
, 0);
264 if ((new == NULL
) || ((new != ERROR_OP
) && (new->prio
== 9))) {
268 if ((new == NULL
) || (new == ERROR_OP
) || (*str
== ',')) {
280 VERBOSE (DEBUG
, fprintf (stdout
, "stop processing bracket\n"));
284 /* check for closing bracket, closing brace or koma */
286 if ((*str
== ')') || (*str
== '}') || (*str
== ',')) {
297 /* look for operators */
299 for (i
= 0; i
< NB_OPERATORS
; i
++) {
300 keyword_t
*operator = operators
+ i
;
301 if (codecmp (operator->keyword
, str
) == 0) {
302 VERBOSE (DEBUG
, fprintf (stdout
, "start processing operator\n"));
304 if ((prio
) && (prio
> operator->prio
)) {
305 VERBOSE (DEBUG
, fprintf (stdout
, "stop because operator priority\n"));
309 str
+= operator->offset
;
310 VERBOSE (INFO
, fprintf (stdout
, "Oper: %d\n", operator->func
));
311 if (subparser (&root
, &str
, operator->func
, operator->nbops
, operator->prio
) == ERROR_OP
) {
315 } else if (*str
== '-') {
316 root
= newelement (Sig
, 1, 9);
321 VERBOSE (DEBUG
, fprintf (stdout
, "stop processing operator\n"));
329 /* look for functions */
331 for (i
= 0; i
< NB_FUNCTIONS
; i
++) {
332 keyword_t
*function
= functions
+ i
;
333 if (codecmp (function
->keyword
, str
) == 0) {
334 VERBOSE (DEBUG
, fprintf (stdout
, "start processing function\n"));
336 VERBOSE (INFO
, fprintf (stdout
, "Func: %d\n", function
->func
));
337 root
= newelement (function
->func
, function
->nbops
, function
->prio
);
342 str
+= function
->offset
;
344 VERBOSE (DEBUG
, fprintf (stdout
, "stop processing function\n"));
352 /* look for constant */
354 for (i
= 0; i
< NB_CONSTANTS
; i
++) {
355 keyword_t
*constant
= constants
+ i
;
356 if (codecmp (constant
->keyword
, str
) == 0) {
357 VERBOSE (DEBUG
, fprintf (stdout
, "start processing constant\n"));
359 VERBOSE (INFO
, fprintf (stdout
, "Const: %d\n", constant
->func
));
360 root
= newelement (constant
->func
, constant
->nbops
, constant
->prio
);
365 str
+= constant
->offset
;
367 VERBOSE (DEBUG
, fprintf (stdout
, "stop processing constant\n"));
375 /* look for number */
377 if (((*str
>= '0') && (*str
<= '9')) ||
378 (*str
== '.') || (*str
== '+') || (*str
== '-')) {
379 VERBOSE (DEBUG
, fprintf (stdout
, "start processing value\n"));
381 double value
= strtod (str
, &pt
);
382 VERBOSE (INFO
, fprintf (stdout
, "Value: %f\n", value
));
384 if ((root
== NULL
) || (root
->prio
== 6)) {
385 new = newelement (Val
, 1, 5);
390 for (i
= 0; i
< root
->nbops
; i
++) {
391 if (root
->ops
[i
] == NULL
) {
404 } else if ((*str
== '+') || (*str
== '-')) {
405 if ((prio
) && (prio
> 1)) {
406 VERBOSE (DEBUG
, fprintf (stdout
, "stop because operator priority\n"));
410 if (subparser (&root
, &str
, Add
, 2, 1) == ERROR_OP
) {
420 VERBOSE (DEBUG
, fprintf (stdout
, "stop processing value\n"));
439 /* print element tree */
441 void print_element (element_t
*root
, int level
)
446 if ((root
== NULL
) || (root
== ERROR_OP
)) {
450 for (i
= 0; i
< level
; i
++) {
451 fprintf (stdout
, " ");
454 switch (root
->func
) {
455 case Val
: func
= "Value"; break;
456 case Sig
: func
= "Sign"; break;
457 case Add
: func
= "Addition"; break;
458 case Sub
: func
= "Subtraction"; break;
459 case Mul
: func
= "Multiplication"; break;
460 case Div
: func
= "Division"; break;
461 case Mod
: func
= "Modulo"; break;
462 case Pow
: func
= "Power"; break;
463 case Sqr
: func
= "Square Root"; break;
464 case Cos
: func
= "Cosine"; break;
465 case Sin
: func
= "Sine"; break;
466 case Atan
: func
= "Arc Tangent"; break;
467 case Log
: func
= "Logarithm"; break;
468 case Exp
: func
= "Exponantial"; break;
469 case Store
: func
= "Store"; break;
470 case Recall
: func
= "Recall"; break;
471 case Inc
: func
= "Increase"; break;
472 case Dec
: func
= "Decrease"; break;
473 case Disp
: func
= "Display"; break;
474 case Quit
: func
= "Quit"; break;
475 case Help
: func
= "Help"; break;
476 case Ans
: func
= "Ans"; break;
477 case Pi
: func
= "Pi"; break;
478 case E
: func
= "E"; break;
479 case Equal
: func
= "Equal"; break;
480 case Diff
: func
= "Different"; break;
481 case Ge
: func
= "Greater or equal"; break;
482 case Le
: func
= "Lesser or equal"; break;
483 case Gt
: func
= "Greater"; break;
484 case Lt
: func
= "Lesser"; break;
485 case And
: func
= "And"; break;
486 case Or
: func
= "Or"; break;
487 case Not
: func
= "Not"; break;
488 case Cond
: func
= "Condition"; break;
489 case While
: func
= "While"; break;
490 case Prog
: func
= "Program"; break;
493 fprintf (stdout
, "Function: %s\n", func
);
495 if ((root
->func
== Val
) && (root
->ops
[0] == NULL
)) {
496 for (i
= 0; i
< level
; i
++) {
497 fprintf (stdout
, " ");
499 fprintf (stdout
, "value: %f\n", root
->value
);
501 for (i
= 0; i
< root
->nbops
; i
++) {
502 print_element (root
->ops
[i
], level
+ 1);
507 /* storage functions */
509 double store (int index
, double value
)
511 if ((index
> 0) && (index
<= STORAGE_SIZE
)) {
512 storage
[index
- 1] = value
;
514 VERBOSE (WARNING
, fprintf (stdout
, "invalid index (%d) [1, %d]\n", index
, STORAGE_SIZE
));
519 double recall (int index
)
521 if ((index
> 0) && (index
<= STORAGE_SIZE
)) {
522 return storage
[index
- 1];
524 VERBOSE (WARNING
, fprintf (stdout
, "invalid index (%d) [1, %d]\n", index
, STORAGE_SIZE
));
529 double increase (int index
)
531 if ((index
> 0) && (index
<= STORAGE_SIZE
)) {
532 return storage
[index
- 1]++;
534 VERBOSE (WARNING
, fprintf (stdout
, "invalid index (%d) [1, %d]\n", index
, STORAGE_SIZE
));
539 double decrease (int index
)
541 if ((index
> 0) && (index
<= STORAGE_SIZE
)) {
542 return storage
[index
- 1]--;
544 VERBOSE (WARNING
, fprintf (stdout
, "invalid index (%d) [1, %d]\n", index
, STORAGE_SIZE
));
552 fprintf (stdout
, "storage:");
553 for (i
= 0; i
< STORAGE_SIZE
; i
++) {
554 fprintf (stdout
, " %g", storage
[i
]);
556 fprintf (stdout
, "\n");
559 /* While do function */
561 double while_do (element_t
*cond
, element_t
*action
)
564 element_t
*temp
= NULL
;
566 VERBOSE (DEBUG
, fprintf (stdout
, "starting while loop\n"));
568 VERBOSE (DEBUG
, fprintf (stdout
, "loop...\n"));
570 temp
= dupelement (cond
);
571 if (!evaluate_element (temp
, 0)) {
575 temp
= dupelement (action
);
576 ret
= evaluate_element (temp
, 0);
580 VERBOSE (DEBUG
, fprintf (stdout
, "ending while loop\n"));
585 /* program function */
587 double program_do (element_t
**prog
, int nbcalls
)
591 for (i
= 0; i
< nbcalls
; i
++) {
592 ret
= evaluate_element (prog
[i
], 0);
602 fprintf (stdout
, "bye\n");
610 fprintf (stdout
, "calc is a simple calculator\n\n");
611 fprintf (stdout
, "supported operators:");
612 fprintf (stdout
, " + - * / %% ^\n");
613 fprintf (stdout
, "camparison operators:");
614 fprintf (stdout
, " == != >= <= > <\n");
615 fprintf (stdout
, "logical operators:");
616 fprintf (stdout
, " & | !\n");
617 fprintf (stdout
, "supported functions:");
618 fprintf (stdout
, " pow sqrt cos sin atan log exp\n");
619 fprintf (stdout
, "storage functions:");
620 fprintf (stdout
, " sto rcl inc dec\n");
621 fprintf (stdout
, "conditional functions:");
622 fprintf (stdout
, " cond\n");
623 fprintf (stdout
, "miscellaneous functions:");
624 fprintf (stdout
, " quit help\n");
625 fprintf (stdout
, "supported constants:");
626 fprintf (stdout
, " e pi\n");
629 /* evaluate element tree */
634 double evaluate_element (element_t
*root
, char mask
)
636 double op0
= 0, op1
= 0;
637 char nextmask
= mask
;
639 if ((root
== NULL
) || (root
== ERROR_OP
)) {
640 VERBOSE (WARNING
, fprintf (stdout
, "error while evaluating\n"));
644 /* mask to manage sub operator sub and div */
645 switch (root
->func
) {
647 nextmask
&= ~MASK_SUB
;
648 nextmask
&= ~MASK_DIV
;
651 nextmask
|= MASK_SUB
;
652 nextmask
&= ~MASK_DIV
;
655 nextmask
&= ~MASK_DIV
;
658 nextmask
|= MASK_DIV
;
664 switch (root
->func
) {
667 op0
= (root
->ops
[0]) ? evaluate_element (root
->ops
[0], nextmask
) : root
->value
;
685 op1
= evaluate_element (root
->ops
[1], nextmask
);
686 } else if (root
->func
!= Store
) {
687 VERBOSE (WARNING
, fprintf (stdout
, "error while evaluating (op[1])\n"));
703 op0
= evaluate_element (root
->ops
[0], 0);
705 VERBOSE (WARNING
, fprintf (stdout
, "error while evaluating (op[0])\n"));
718 if (root
->ops
[0] == NULL
) {
719 VERBOSE (WARNING
, fprintf (stdout
, "error while evaluating (op[0])\n"));
725 switch (root
->func
) {
726 case Val
: return op0
;
727 case Sig
: return -op0
;
728 case Add
: return ((mask
& MASK_SUB
) == 0) ? op0
+ op1
: op0
- op1
;
729 case Sub
: return ((mask
& MASK_SUB
) == 0) ? op0
- op1
: op0
+ op1
;
730 case Mul
: return ((mask
& MASK_DIV
) == 0) ? op0
* op1
: op0
/ op1
;
731 case Div
: return ((mask
& MASK_DIV
) == 0) ? op0
/ op1
: op0
* op1
;
732 case Mod
: return fmod (op0
, op1
);
733 case Pow
: return pow (op0
, op1
);
734 case Sqr
: return sqrt (op0
);
735 case Cos
: return cos (op0
);
736 case Sin
: return sin (op0
);
737 case Atan
: return atan (op0
);
738 case Log
: return log (op0
);
739 case Exp
: return exp (op0
);
740 case Store
: return store ((int)op0
, (op1
) ? op1
: answer
);
741 case Recall
: return recall ((int)op0
);
742 case Inc
: return increase ((int)op0
);
743 case Dec
: return decrease ((int)op0
);
744 case Disp
: display (); break;
745 case Quit
: quit (); break;
746 case Help
: help (); break;
747 case Ans
: return answer
;
748 case Pi
: return M_PI
;
750 case Equal
: return op0
== op1
;
751 case Diff
: return op0
!= op1
;
752 case Ge
: return op0
>= op1
;
753 case Le
: return op0
<= op1
;
754 case Gt
: return op0
> op1
;
755 case Lt
: return op0
< op1
;
756 case And
: return (op0
!= 0) && (op1
!= 0);
757 case Or
: return (op0
!= 0) || (op1
!= 0);
758 case Not
: return (op0
== 0);
760 if ((op0
) && (root
->ops
[1])) {
761 return evaluate_element (root
->ops
[1], 0);
762 } else if ((!op0
) && (root
->ops
[2])) {
763 return evaluate_element (root
->ops
[2], 0);
767 case While
: return while_do (root
->ops
[0], root
->ops
[1]);
768 case Prog
: return program_do (root
->ops
, root
->nbops
);
774 /* vim: set ts=4 sw=4 et: */