char *minform = NULL;
workspace_t *programs = NULL;
+int nb_programs = 0;
/* compare codes */
return ERROR_OP;
}
element_t **prog = NULL;
- new = newelement (Prog, 0, 5);
+ new = newelement (Code, 0, 5);
root = new;
prog = &root;
delelement (root);
return ERROR_OP;
}
- element_t *newprog = newelement (Prog, (*prog)->nbops + 1, 5);
+ element_t *newprog = newelement (Code, (*prog)->nbops + 1, 5);
for (i = 0; i < (*prog)->nbops; i++) {
newprog->ops[i] = (*prog)->ops[i];
(*prog)->ops[i] = NULL;
case Not: func = "Not"; break;
case Cond: func = "Condition"; break;
case While: func = "While"; break;
- case Prog: func = "Program"; break;
+ case Code: func = "Code"; break;
case Print: func = "Print"; break;
case Prog: func = "Program"; break;
case Call: func = "Call"; break;
/* program function */
-double program_do (element_t **prog, int nbcalls)
+double execute_code (element_t **prog, int nbcalls)
{
double ret = 0;
int i;
/* program function */
void prog (int id, int nbmems, element_t *root)
-{ }
+{
+ int i, n = -1;
+
+ if (programs == NULL) {
+
+ /* initial memory allocation */
+ programs = (workspace_t *) calloc (1, sizeof (workspace_t));
+ if (programs == NULL) {
+ VERBOSE (ERROR, fprintf (stderr, "can't allocate memory\n"));
+ exit (1);
+ }
+ nb_programs = 1;
+ n = 0;
+
+ } else {
+
+ /* look for existing program */
+ for (i = 0; i < nb_programs; i++) {
+ if ((programs + i)->id == id) {
+ n = i;
+ break;
+ }
+ }
+ if (n == -1) {
-double call (int id, int nbobs, element_t **ops)
+ /* new program */
+ n = nb_programs++;
+ programs = (workspace_t *) realloc (programs, nb_programs * sizeof (workspace_t));
+ } else {
+
+ /* clean old program */
+ if ((programs + n)->storage) {
+ free ((programs + n)->storage);
+ }
+ if ((programs + n)->root) {
+ delelement ((programs + n)->root);
+ }
+ }
+ }
+
+ /* set program */
+ (programs + n)->id = id;
+ (programs + n)->answer = 0;
+ (programs + n)->storage = (double *) calloc (nbmems, sizeof (double));
+ if ((programs + n)->storage == NULL) {
+ VERBOSE (ERROR, fprintf (stderr, "can't allocate memory\n"));
+ exit (1);
+ }
+ (programs + n)->storage_size = nbmems;
+ (programs + n)->root = dupelement (root);
+}
+
+double call (int id, int nbops, element_t **ops)
{
- return 0;
+ workspace_t tmp = {0};
+ int i, n = -1;
+ double ret = 0;
+
+ if (programs) {
+
+ /* look for program */
+ for (i = 0; i < nb_programs; i++) {
+ if ((programs + i)->id == id) {
+ n = i;
+ break;
+ }
+ }
+ if (n == -1) {
+ VERBOSE (WARNING, fprintf (stdout, "error unknown program (%d)\n", id));
+ return 0;
+ }
+
+ /* store context */
+ tmp.answer = answer;
+ tmp.storage = storage;
+ tmp.storage_size = storage_size;
+
+ /* change context */
+ answer = 0;
+ storage = (programs + n)->storage;
+ storage_size = (programs + n)->storage_size;
+ if (nbops > storage_size) {
+ storage = realloc (storage, nbops * sizeof (double));
+ storage_size = nbops;
+ }
+ for (i = 0; i < nbops; i++) {
+ double val = evaluate_element (ops[i], 0);
+ store (i + 1, val);
+ }
+
+ /* evaluate program */
+ element_t *elements = dupelement ((programs + n)->root);
+ ret = evaluate_element (elements, 0);
+
+ /* restore context */
+ answer = tmp.answer;
+ storage = tmp.storage;
+ storage_size = tmp.storage_size;
+ }
+
+ return ret;
}
void list ()
-{ }
+{
+ int i;
+ if (programs != NULL) {
+ fprintf (stdout, "programs:");
+ for (i = 0; i < nb_programs; i++) {
+ fprintf (stdout, " %d", (programs + i)->id);
+ }
+ fprintf (stdout, "\n");
+ }
+}
void edit (int id)
{ }
{
double op0 = 0, op1 = 0;
char nextmask = mask;
+ int i, nb;
if ((root == NULL) || (root == ERROR_OP)) {
VERBOSE (WARNING, fprintf (stdout, "error while evaluating\n"));
case Mem:
case Cond:
case Call:
- case List:
case Edit:
case Del:
if (root->ops[0]) {
case Ans:
case Pi:
case E:
- case Prog:
+ case Code:
+ case List:
break;
case While:
if (root->ops[0] == NULL) {
return 0;
}
case While: return while_do (root->ops[0], root->ops[1]);
- case Prog: return program_do (root->ops, root->nbops);
+ case Code: return execute_code (root->ops, root->nbops);
case Print: return print (op0);
case Prog: prog ((int)op0, (int)op1, root->ops[2]); break;
- case Call: return call ((int)op0, root->nbops + 1, root->ops + 1);
+ 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 Del: del ((int)op0); break;
}
return 0;