-#include <malloc.h>
#include <math.h>
+#include <stdio.h>
#include <stdlib.h>
+#include <string.h>
+#include "alloc.h"
+#include "argument.h"
#include "debug.h"
-#include "fdprintf.h"
+#include "element.h"
+#include "format.h"
+#include "program.h"
+#include "stack.h"
+#include "storage.h"
#include "parser.h"
+/* external definition */
+
+extern void history ();
+
+/* global variables */
+
+double answer = 0;
+
+double anglefactor = 1;
+
/* compare codes */
int codecmp (char *ref, char *str)
return 0;
}
-/* allocate new element */
+/* subparser function */
-element_t *newelement (func_t function, int nbops, int prio)
+element_t *subparser (element_t **proot, char **pstr, func_t func, int nbops, int prio)
{
- element_t *new = (element_t *) calloc (1, sizeof (element_t));
- if (new == NULL) {
- VERBOSE (ERROR, fdprintf (stdfderr, "can't allocate memory\n"));
- return NULL;
+ element_t *new = newelement (func, nbops, prio);
+ new->ops[0] = *proot;
+ new->ops[1] = parser (*pstr, pstr, new->prio);
+ if ((new->ops[1] == NULL) || ((new->ops[1] != ERROR_OP) && (new->ops[1]->prio == 9))) {
+ delelement (new->ops[1]);
+ new->ops[1] = ERROR_OP;
+ }
+ if (new->ops[1] == ERROR_OP) {
+ delelement (new);
+ *proot = NULL;
+ return ERROR_OP;
}
- new->func = function;
- new->nbops = nbops;
- new->prio = prio;
+ *proot = newelement (Val, 1, 5);
+ (*proot)->ops[0] = new;
- return new;
+ return *proot;
}
-/* functions */
-
-#define NB_OPERATORS 5
-
-keyword_t operators[NB_OPERATORS] = {
- { "+\t", Add, 2, 1, 1},
- { "-\t", Sub, 2, 1, 1},
- { "*", Mul, 2, 1, 2},
- { "/", Div, 2, 1, 2},
- { "^", Pow, 2, 1, 3}
-};
-
-#define NB_FUNCTIONS 7
-keyword_t functions[NB_FUNCTIONS] = {
- { "sqrt", Sqr, 1, 4, 4},
- { "pow", Pow, 2, 3, 4},
- { "cos", Cos, 1, 3, 4},
- { "sin", Sin, 1, 3, 4},
- { "atan", Atn, 1, 4, 4},
- { "exp", Exp, 1, 3, 4},
- { "log", Log, 1, 3, 4}
-};
-
/* parser function */
element_t *parser (char *str, char **next, int prio)
{
element_t *root = NULL;
+ char *string = str;
int i;
- VERBOSE (DEBUG, PRINTOUT ("Starting parsing\n"));
+ VERBOSE (DEBUG, fprintf (stdout, "Starting parsing\n"));
/* main loop */
while (*str != '\0') {
int found = 0;
element_t *new = NULL;
- VERBOSE (INFO, PRINTOUT ("Processing: %s\n", str));
+ VERBOSE (INFO, fprintf (stdout, "Processing: %s\n", str));
+
+ /* end without printing */
+
+ if (*str == ';') {
+ if (root) {
+ root->hidden = 1;
+ }
+ break;
+ }
/* skip spaces and tabs */
continue;
}
+ /* check for open brace */
+
+ if (*str == '{') {
+ VERBOSE (DEBUG, fprintf (stdout, "start processing brace\n"));
+ if (root != NULL) {
+ delelement (root);
+ return ERROR_OP;
+ }
+ root = newelement (Code, 0, 5);
+
+ do {
+ new = parser (str + 1, &str, 0);
+ if ((new == NULL) || ((new != ERROR_OP) && (new->prio == 9))) {
+ delelement (new);
+ new = ERROR_OP;
+ }
+ if (new == ERROR_OP) {
+ delelement (root);
+ return ERROR_OP;
+ }
+ element_t *newcode = newelement (Code, root->nbops + 1, 5);
+ for (i = 0; i < root->nbops; i++) {
+ newcode->ops[i] = root->ops[i];
+ root->ops[i] = NULL;
+ }
+ newcode->ops[root->nbops] = new;
+ delelement (root);
+ root = newcode;
+ } while (*str == ',');
+
+ if (*str != '}') {
+ delelement (root);
+ return ERROR_OP;
+ }
+ str++;
+ VERBOSE (DEBUG, fprintf (stdout, "stop processing brace\n"));
+ continue;
+ }
+
/* check for open bracket */
if (*str == '(') {
- VERBOSE (DEBUG, PRINTOUT ("start processing bracket\n"));
+ VERBOSE (DEBUG, fprintf (stdout, "start processing bracket\n"));
if (root) {
do {
found = 0;
new = parser (str + 1, &str, 0);
- if (new == ERROR_OP) {
+ if ((new == NULL) || ((new != ERROR_OP) && (new->prio == 9))) {
+ delelement (new);
+ new = ERROR_OP;
+ }
+ if ((new == NULL) || (new == ERROR_OP)) {
+ delelement (root);
return ERROR_OP;
}
for (i = 0; i < root->nbops; i++) {
}
}
if (!found) {
+ delelement (new);
+ delelement (root);
return ERROR_OP;
}
} while (*str == ',');
} else {
- root = newelement (Val, 1, 4);
- if (root == NULL) {
- return ERROR_OP;
- }
+ root = newelement (Val, 1, 5);
new = parser (str + 1, &str, 0);
- if ((new == ERROR_OP) || (*str == ',')) {
+ if ((new == NULL) || ((new != ERROR_OP) && (new->prio == 9))) {
+ delelement (new);
+ new = ERROR_OP;
+ }
+ if ((new == NULL) || (new == ERROR_OP) || (*str == ',')) {
+ delelement (new);
+ delelement (root);
return ERROR_OP;
}
root->ops[0] = new;
}
+ if (*str != ')') {
+ delelement (root);
+ return ERROR_OP;
+ }
str++;
- VERBOSE (DEBUG, PRINTOUT ("stop processing bracket\n"));
+ VERBOSE (DEBUG, fprintf (stdout, "stop processing bracket\n"));
continue;
}
- /* check for closing bracket or koma */
+ /* check for closing bracket, closing brace or koma */
- if ((*str == ')') || (*str == ',')) {
+ if ((*str == ')') || (*str == '}') || (*str == ',')) {
+ if (prio == -9) {
+ delelement (root);
+ return ERROR_OP;
+ }
if (next != NULL) {
*next = str;
}
for (i = 0; i < NB_OPERATORS; i++) {
keyword_t *operator = operators + i;
if (codecmp (operator->keyword, str) == 0) {
- VERBOSE (DEBUG, PRINTOUT ("start processing operator\n"));
- if (root) {
+ VERBOSE (DEBUG, fprintf (stdout, "start processing operator\n"));
+ if ((root) && (root->prio == 9)) {
+ VERBOSE (DEBUG, fprintf (stdout, "terminal function (%d)\n", root->func));
+ delelement (root);
+ return ERROR_OP;
+ } else if (root) {
if ((prio) && (prio > operator->prio)) {
- VERBOSE (DEBUG, PRINTOUT ("stop processing operator because operator priority\n"));
+ VERBOSE (DEBUG, fprintf (stdout, "stop because operator priority\n"));
*next = str;
return root;
}
str += operator->offset;
- VERBOSE (INFO, PRINTOUT ("Oper: %d\n", operator->func));
- new = newelement (operator->func, operator->nbops, operator->prio);
- if (new == NULL) {
- return ERROR_OP;
- }
- new->ops[0] = root;
- new->ops[1] = parser (str, &str, new->prio);
- if (new->ops[1] == ERROR_OP) {
- return ERROR_OP;
- }
- root = newelement (Val, 1, 4);
- if (root == ERROR_OP) {
+ VERBOSE (INFO, fprintf (stdout, "Oper: %d\n", operator->func));
+ if (subparser (&root, &str, operator->func, operator->nbops, operator->prio) == ERROR_OP) {
+ delelement (root);
return ERROR_OP;
}
- root->ops[0] = new;
+ } else if (*str == '-') {
+ root = newelement (Sig, 1, 6);
} else {
return ERROR_OP;
}
found = 1;
- VERBOSE (DEBUG, PRINTOUT ("stop processing operator\n"));
+ VERBOSE (DEBUG, fprintf (stdout, "stop processing operator\n"));
break;
}
}
for (i = 0; i < NB_FUNCTIONS; i++) {
keyword_t *function = functions + i;
if (codecmp (function->keyword, str) == 0) {
- VERBOSE (DEBUG, PRINTOUT ("start processing function\n"));
+ VERBOSE (DEBUG, fprintf (stdout, "start processing function\n"));
if (root == NULL) {
- VERBOSE (INFO, PRINTOUT ("Func: %d\n", function->func));
- new = newelement (function->func, function->nbops, function->prio);
- if (new == NULL) {
- return ERROR_OP;
- }
- root = new;
- } else {
+ VERBOSE (INFO, fprintf (stdout, "Func: %d\n", function->func));
+ root = newelement (function->func, function->nbops, function->prio);
+ } else {
+ delelement (root);
return ERROR_OP;
}
str += function->offset;
found = 1;
- VERBOSE (DEBUG, PRINTOUT ("stop processing function\n"));
+ VERBOSE (DEBUG, fprintf (stdout, "stop processing function\n"));
break;
}
}
continue;
}
- /* look for number */
+ /* look for constant */
- if (((*str >= '0') && (*str <= '9')) ||
- (*str == '.') || (*str == '+') || (*str == '-')) {
- VERBOSE (DEBUG, PRINTOUT ("start processing value\n"));
- char *pt;
- float value = strtof (str, &pt);
- VERBOSE (INFO, PRINTOUT ("Value: %f\n", value));
- if (str != pt) {
- new = newelement (Val, 1, 4);
- if (new == NULL) {
+ for (i = 0; i < NB_CONSTANTS; i++) {
+ keyword_t *constant = constants + i;
+ if (codecmp (constant->keyword, str) == 0) {
+ VERBOSE (DEBUG, fprintf (stdout, "start processing constant\n"));
+ if (root == NULL) {
+ VERBOSE (INFO, fprintf (stdout, "Const: %d\n", constant->func));
+ root = newelement (constant->func, constant->nbops, constant->prio);
+ } else {
+ delelement (root);
return ERROR_OP;
}
+ str += constant->offset;
+ found = 1;
+ VERBOSE (DEBUG, fprintf (stdout, "stop processing constant\n"));
+ break;
+ }
+ }
+ if (found) {
+ continue;
+ }
+
+ /* look for number */
+
+ VERBOSE (DEBUG, fprintf (stdout, "start processing value\n"));
+ char *pt;
+ double value = (get_ibase () == 10) ? strtod (str, &pt) : strtoul (str, &pt, get_ibase ());
+ VERBOSE (INFO, fprintf (stdout, "Value: %f\n", value));
+ if (str != pt) {
+ if ((root == NULL) || (root->prio == 6)) {
+ new = newelement (Val, 1, 5);
new->value = value;
if (root == NULL) {
root = new;
- } else if (root->func == Val) {
- if ((*str == '+') || (*str == '-')) {
- element_t *add = newelement (Add, 2, 1);
- if (add == NULL) {
- return ERROR_OP;
+ } else {
+ for (i = 0; i < root->nbops; i++) {
+ if (root->ops[i] == NULL) {
+ root->ops[i] = new;
+ found = 1;
+ break;
}
- add->ops[0] = root;
- add->ops[1] = new;
- root = add;
- } else {
+ }
+ if (!found) {
+ delelement (new);
+ delelement (root);
return ERROR_OP;
}
- } else {
- return ERROR_OP;
}
str = pt;
- found = 1;
+ } else if ((*str == '+') || (*str == '-')) {
+ if ((prio) && (prio > 1)) {
+ VERBOSE (DEBUG, fprintf (stdout, "stop because operator priority\n"));
+ *next = str;
+ return root;
+ }
+ if (subparser (&root, &str, Add, 2, 1) == ERROR_OP) {
+ delelement (root);
+ return ERROR_OP;
+ }
+ } else {
+ delelement (root);
+ return ERROR_OP;
}
- VERBOSE (DEBUG, PRINTOUT ("stop processing value\n"));
+ found = 1;
}
+ VERBOSE (DEBUG, fprintf (stdout, "stop processing value\n"));
/* error */
if (!found) {
+ delelement (root);
return ERROR_OP;
}
if (next != NULL) {
*next = str;
}
+
+ /* save string */
+ if (root != NULL) {
+ root->string = string;
+ }
+
return root;
}
}
for (i = 0; i < level; i++) {
- PRINTOUT (" ");
+ printf (" ");
}
switch (root->func) {
case Val: func = "Value"; break;
+ case Sig: func = "Sign"; break;
case Add: func = "Addition"; break;
case Sub: func = "Subtraction"; break;
case Mul: func = "Multiplication"; break;
case Div: func = "Division"; break;
+ case Mod: func = "Modulo"; break;
case Pow: func = "Power"; break;
case Sqr: func = "Square Root"; break;
case Cos: func = "Cosine"; break;
case Sin: func = "Sine"; break;
- case Atn: func = "Arc Tangent"; break;
- case Log: func = "Logarithm"; break;
+ case Tan: func = "Tangent"; break;
+ case Acos: func = "Arc Cosine"; break;
+ case Asin: func = "Arc Sine"; break;
+ case Atan: func = "Arc Tangent"; break;
+ case Ln: func = "Logarithm (natural)"; break;
+ case Log: func = "Logarithm (10 base)"; break;
case Exp: func = "Exponantial"; break;
+ case Erfc: func = "Complementary Error Function"; break;
+ case Erf: func = "Error Function"; break;
+ case Abs: func = "Absolute value"; break;
+ case Ceil: func = "Ceil value"; break;
+ case Floor: func = "Floor value"; break;
+ case Store: func = "Store"; break;
+ case Recall: func = "Recall"; break;
+ case Inc: func = "Increase"; break;
+ case Dec: func = "Decrease"; break;
+ case Disp: func = "Display"; break;
+ case Memory: func = "Memory"; break;
+ case Clear: func = "Clear"; break;
+ case Quit: func = "Quit"; break;
+ case Help: func = "Help"; break;
+ case History: func = "History"; break;
+ case Ans: func = "Ans"; break;
+ case Pi: func = "Pi"; break;
+ case E: func = "E"; break;
+ case Equal: func = "Equal"; break;
+ case Diff: func = "Different"; break;
+ case Ge: func = "Greater or equal"; break;
+ case Le: func = "Lesser or equal"; break;
+ case Gt: func = "Greater"; break;
+ case Lt: func = "Lesser"; break;
+ case And: func = "And"; break;
+ case Or: func = "Or"; break;
+ case Not: func = "Not"; break;
+ case Cond: func = "Condition"; break;
+ case While: func = "While"; break;
+ case Code: func = "Code"; break;
+ case Print: func = "Print"; break;
+ case Prog: func = "Program"; break;
+ case Arg: func = "Argument"; break;
+ case Call: func = "Call"; break;
+ case List: func = "List"; break;
+ case Edit: func = "Edit"; break;
+ case Del: func = "Del"; break;
+ case Get: func = "Get"; break;
+ case Length: func = "Length"; break;
+ case Pop: func = "Pop"; break;
+ case Push: func = "Push"; break;
+ case Put: func = "Put"; break;
+ case Set: func = "Set"; break;
+ case Show: func = "Show"; break;
+ case Max: func = "Maximum"; break;
+ case Mean: func = "Mean"; break;
+ case Median: func = "Median"; break;
+ case Min: func = "Minimum"; break;
+ case Order: func = "Order"; break;
+ case Prod: func = "Product"; break;
+ case Sum: func = "Sum"; break;
+ case Variance: func = "Variance"; break;
+ case Precision: func = "Precision"; break;
+ case Base: func = "Base"; break;
+ case Deg: func = "Degree"; break;
+ case Grad: func = "Gradian"; break;
+ case Rad: func = "Radian"; break;
}
- PRINTOUT ("Function: %s\n", func);
+ printf ("Function: %s\n", func);
if ((root->func == Val) && (root->ops[0] == NULL)) {
for (i = 0; i < level; i++) {
- PRINTOUT (" ");
+ printf (" ");
}
- PRINTOUT ("value: %f\n", root->value);
+ printf ("value: %f\n", root->value);
} else {
for (i = 0; i < root->nbops; i++) {
print_element (root->ops[i], level + 1);
}
}
+/* While do function */
+
+double while_do (element_t *cond, element_t *action)
+{
+ double ret = 0;
+ element_t *temp = NULL;
+
+ VERBOSE (DEBUG, fprintf (stdout, "starting while loop\n"));
+ while (1) {
+ VERBOSE (DEBUG, fprintf (stdout, "loop...\n"));
+
+ temp = dupelement (cond);
+ double test = evaluate_element (temp, 0);
+ delelement (temp);
+ if (!test) {
+ break;
+ }
+ if (action) {
+ temp = dupelement (action);
+ ret = evaluate_element (temp, 0);
+ delelement (temp);
+ }
+ }
+
+ VERBOSE (DEBUG, fprintf (stdout, "ending while loop\n"));
+
+ return ret;
+}
+
+/* program function */
+
+double execute_code (element_t **prog, int nbcalls)
+{
+ double ret = 0;
+ int i;
+ for (i = 0; i < nbcalls; i++) {
+ ret = evaluate_element (prog[i], 0);
+ }
+ return ret;
+}
+
+/* quit function */
+
+void quit (void)
+{
+ printf ("bye\n");
+ exit (0);
+}
+
+/* help message */
+
+void help (void)
+{
+ printf ("calc is a simple calculator\n\n");
+ printf ("arithmetic op.:");
+ printf (" + - * / %% ^\n");
+ printf ("comparison op.:");
+ printf (" == != >= <= > <\n");
+ printf ("logical op.:");
+ printf (" & | !\n");
+ printf ("mathematic func.:");
+ printf (" exp ln log pow sqrt\n");
+ printf ("trigonometric func.:");
+ printf (" acos asin atan cos sin tan\n");
+ printf ("error functions:");
+ printf (" erf erfc\n");
+ printf ("miscellaneous func.:");
+ printf (" abs ceil floor\n");
+ printf ("storage func.:");
+ printf (" clear dec disp inc mem rcl sto\n");
+ printf ("control flow prim.:");
+ printf (" cond print while {} ;\n");
+ printf ("program management:");
+ printf (" arg call del edit ls prog\n");
+ printf ("stack management:");
+ printf (" get len pop push put set show\n");
+ printf ("stack func.:");
+ printf (" max mean med min ord prod sum var\n");
+ printf ("control management:");
+ printf (" base format help hist quit\n");
+ printf ("angle management:");
+ printf (" deg grad rad\n");
+ printf ("constants:");
+ printf (" ans e pi\n");
+}
+
+/* format function */
+
+int format (int precision)
+{
+ if (precision > 0) {
+ set_precision (precision);
+ set_format ();
+ } else if (precision != -1) {
+ VERBOSE (WARNING, fprintf (stdout, "error incorrect precision (%d)\n", precision));
+ return 0;
+ }
+ return get_precision ();
+}
+
+/* base function */
+
+void base (int in, int out)
+{
+ if ((in > 0) && (in < 37)) {
+ set_base (in, in);
+ if ((out > 0) && (out < 37)) {
+ set_base (in, out);
+ } else if (out != - 1) {
+ VERBOSE (WARNING, fprintf (stdout, "error incorrect output base (%d)\n", out));
+ }
+ } else if (in != -1 ) {
+ VERBOSE (WARNING, fprintf (stdout, "error incorrect input base (%d)\n", in));
+ } else {
+ printf ("base (I/O): %s\n", show_base ());
+ }
+}
+
/* evaluate element tree */
-double evaluate_element (element_t *root)
+#define MASK_SUB 0x1
+#define MASK_DIV 0x2
+
+double evaluate_element (element_t *root, char mask)
{
double op0 = 0, op1 = 0;
+ char nextmask = mask;
+ int i, nb;
if ((root == NULL) || (root == ERROR_OP)) {
- VERBOSE (WARNING, PRINTOUT ("error while evaluating\n"));
+ VERBOSE (WARNING, fprintf (stdout, "error while evaluating\n"));
return 0;
}
+ /* mask to manage sub operator sub and div */
+ switch (root->func) {
+ case Add:
+ nextmask &= ~MASK_SUB;
+ nextmask &= ~MASK_DIV;
+ break;
+ case Sub:
+ nextmask |= MASK_SUB;
+ nextmask &= ~MASK_DIV;
+ break;
+ case Mul:
+ nextmask &= ~MASK_DIV;
+ break;
+ case Div:
+ nextmask |= MASK_DIV;
+ break;
+ default:
+ nextmask = mask;
+ }
+
switch (root->func) {
case Val:
- if (root->ops[0]) {
- return evaluate_element (root->ops[0]);
- } else {
- return root->value;
- }
+ case Sig:
+ op0 = (root->ops[0]) ? evaluate_element (root->ops[0], nextmask) : root->value;
break;
case Add:
case Sub:
case Mul:
case Div:
+ case Mod:
case Pow:
+ case Equal:
+ case Diff:
+ case Ge:
+ case Le:
+ case Gt:
+ case Lt:
+ case And:
+ case Or:
if (root->ops[1]) {
- op1 = evaluate_element (root->ops[1]);
+ op1 = evaluate_element (root->ops[1], nextmask);
} else {
- VERBOSE (WARNING, PRINTOUT ("error while evaluating (op[1])\n"));
+ VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[1])\n"));
return 0;
}
/* fallthrough */
case Sqr:
case Cos:
case Sin:
- case Atn:
+ case Tan:
+ case Acos:
+ case Asin:
+ case Atan:
+ case Ln:
case Log:
case Exp:
+ case Erfc:
+ case Erf:
+ case Abs:
+ case Ceil:
+ case Floor:
+ case Recall:
+ case Inc:
+ case Dec:
+ case Not:
+ case Cond:
+ case Prog:
+ case Arg:
+ case Call:
+ case Edit:
+ case Del:
+ case Get:
if (root->ops[0]) {
- op0 = evaluate_element (root->ops[0]);
+ op0 = evaluate_element (root->ops[0], 0);
} else {
- VERBOSE (WARNING, PRINTOUT ("error while evaluating (op[0])\n"));
+ VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[0])\n"));
return 0;
}
+ break;
+ case Disp:
+ case Clear:
+ case Quit:
+ case Help:
+ case History:
+ case Ans:
+ case Pi:
+ case E:
+ case Code:
+ case List:
+ case Length:
+ case Pop:
+ case Set:
+ case Show:
+ case Median:
+ case Order:
+ case Prod:
+ case Sum:
+ case Deg:
+ case Grad:
+ case Rad:
+ break;
+ case While:
+ if (root->ops[0] == NULL) {
+ VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[0])\n"));
+ return 0;
+ }
+ break;
+ case Memory:
+ if (root->ops[0] == NULL) {
+ op0 = -1;
+ } else {
+ op0 = (int)evaluate_element (root->ops[0], 0);
+ if (op0 < 0) {
+ VERBOSE (WARNING, fprintf (stdout, "error incorrect memory size (%d)\n", (int)op0));
+ return 0;
+ }
+ }
+ break;
+ case Push:
+ case Print:
+ op0 = (root->ops[0]) ? evaluate_element (root->ops[0], 0) : answer;
+ break;
+ case Store:
+ case Put:
+ if (root->ops[0]) {
+ op0 = evaluate_element (root->ops[0], 0);
+ } else {
+ VERBOSE (WARNING, fprintf (stdout, "error while evaluating (op[0])\n"));
+ return 0;
+ }
+ op1 = (root->ops[1]) ? evaluate_element (root->ops[1], 0) : answer;
+ break;
+ case Max:
+ case Mean:
+ case Min:
+ case Variance:
+ if (root->ops[0]) {
+ op0 = evaluate_element (root->ops[0], 0);
+ op1 = (root->ops[1]) ? evaluate_element (root->ops[1], 0) : answer;
+ }
+ break;
+ case Precision:
+ op0 = (root->ops[0]) ? evaluate_element (root->ops[0], 0) : -1;
+ break;
+ case Base:
+ op0 = (root->ops[0]) ? evaluate_element (root->ops[0], 0) : -1;
+ op1 = (root->ops[1]) ? evaluate_element (root->ops[1], 0) : -1;
+ break;
}
switch (root->func) {
- case Add: return op0 + op1;
- case Sub: return op0 - op1;
- case Mul: return op0 * op1;
- case Div: return op0 / op1;
+ case Val: return op0;
+ case Sig: return -op0;
+ case Add: return ((mask & MASK_SUB) == 0) ? op0 + op1 : op0 - op1;
+ case Sub: return ((mask & MASK_SUB) == 0) ? op0 - op1 : op0 + op1;
+ case Mul: return ((mask & MASK_DIV) == 0) ? op0 * op1 : op0 / op1;
+ case Div: return ((mask & MASK_DIV) == 0) ? op0 / op1 : op0 * op1;
+ case Mod: return fmod (op0, op1);
case Pow: return pow (op0, op1);
case Sqr: return sqrt (op0);
- case Cos: return cos (op0);
- case Sin: return sin (op0);
- case Atn: return atan (op0);
- case Log: return log (op0);
+ case Cos: return cos (op0 / anglefactor);
+ case Sin: return sin (op0 / anglefactor);
+ case Tan: return tan (op0 / anglefactor);
+ case Acos: return acos (op0) * anglefactor;
+ case Asin: return asin (op0) * anglefactor;
+ case Atan: return atan (op0) * anglefactor;
+ case Ln: return log (op0);
+ case Log: return log10 (op0);
case Exp: return exp (op0);
- default: break;
+ case Erfc: return erfc (op0);
+ case Erf: return erf (op0);
+ case Abs: return fabs (op0);
+ case Ceil: return ceil (op0);
+ case Floor: return floor (op0);
+ case Store: return store ((int)op0, op1);
+ case Recall: return recall ((int)op0);
+ case Inc: return increase ((int)op0);
+ case Dec: return decrease ((int)op0);
+ case Disp: display (); break;
+ case Memory: return memory ((root->ops[0]) ? (int)op0 : -1);
+ case Clear: clear (); break;
+ case Quit: quit (); break;
+ case Help: help (); break;
+ case History: history (); break;
+ case Ans: return answer;
+ case Pi: return M_PI;
+ case E: return M_E;
+ case Equal: return op0 == op1;
+ case Diff: return op0 != op1;
+ case Ge: return op0 >= op1;
+ case Le: return op0 <= op1;
+ case Gt: return op0 > op1;
+ case Lt: return op0 < op1;
+ case And: return (op0 != 0) && (op1 != 0);
+ case Or: return (op0 != 0) || (op1 != 0);
+ case Not: return (op0 == 0);
+ case Cond:
+ if ((op0) && (root->ops[1])) {
+ return evaluate_element (root->ops[1], 0);
+ } else if ((!op0) && (root->ops[2])) {
+ return evaluate_element (root->ops[2], 0);
+ } else {
+ return 0;
+ }
+ case While: return while_do (root->ops[0], root->ops[1]);
+ case Code: return execute_code (root->ops, root->nbops);
+ case Print: return print (op0);
+ case Prog:
+ prog ((int)op0, root->ops[1]);
+ savestring ((int)op0, root->string);
+ break;
+ case Arg: return arg ((int)op0);
+ case Call:
+ for (i = 1, nb = 0; i < root->nbops; i++) {
+ if (root->ops[i]) {
+ nb++;
+ }
+ }
+ return call ((int)op0, nb, root->ops + 1);
+ case List: list (); break;
+ case Edit: edit ((int)op0); break;
+ case Del: del ((int)op0); break;
+ case Get: return get ((int)op0);
+ case Length: return length ();
+ case Pop: return pop ();
+ case Push: return push (op0);
+ case Put: return put ((int)op0, op1);
+ case Set:
+ for (i = 0, nb =0; i < root->nbops; i++) {
+ if (root->ops[i]) {
+ nb++;
+ }
+ }
+ return set (nb, root->ops);
+ case Show: show (); break;
+ case Max:
+ if (root->ops[0]) {
+ return op0 > op1 ? op0 : op1;
+ }
+ return max ();
+ case Mean:
+ if (root->ops[0]) {
+ return (op0 + op1) / 2;
+ }
+ return mean ();
+ case Median: return median ();
+ case Min:
+ if (root->ops[0]) {
+ return op0 < op1 ? op0 : op1;
+ }
+ return min ();
+ case Order: order (); break;
+ case Prod: return prod ();
+ case Sum: return sum ();
+ case Variance:
+ if (root->ops[0]) {
+ double m = (op0 + op1) / 2;
+ op0 -= m;
+ op1 -= m;
+ return op0 * op0 + op1 * op1;
+ }
+ return variance ();
+ case Precision:
+ return format ((int)op0);
+ case Base:
+ base ((int)op0, (int)op1);
+ break;
+ case Deg:
+ anglefactor = 180 / M_PI;
+ break;
+ case Grad:
+ anglefactor = 200 / M_PI;
+ break;
+ case Rad:
+ anglefactor = 1;
+ break;
}
return 0;