[BACK]Return to ox_eval.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / ox_gsl

Diff for /OpenXM/src/ox_gsl/ox_eval.c between version 1.3 and 1.4

version 1.3, 2018/04/05 13:02:39 version 1.4, 2018/04/06 10:44:51
Line 1 
Line 1 
 /* $OpenXM: OpenXM/src/ox_gsl/ox_eval.c,v 1.2 2018/04/05 05:53:52 ohara Exp $ */  /* $OpenXM: OpenXM/src/ox_gsl/ox_eval.c,v 1.3 2018/04/05 13:02:39 ohara Exp $ */
   
 #include <stdio.h>  #include <stdio.h>
 #include <stdlib.h>  #include <stdlib.h>
Line 67  typedef struct {
Line 67  typedef struct {
 entry global_dic[512] = {  entry global_dic[512] = {
     {"sin",0,sin,1},      {"sin",0,sin,1},
     {"cos",0,cos,1},      {"cos",0,cos,1},
       {"tan",0,tan,1},
       {"sinh",0,sinh,1},
       {"cosh",0,cosh,1},
       {"tanh",0,tanh,1},
       {"asin",0,asin,1},
       {"acos",0,acos,1},
       {"asinh",0,asinh,1},
       {"acosh",0,acosh,1},
       {"erf",0,erf,1},
     {"exp",0,exp,1},      {"exp",0,exp,1},
       {"exp2",0,exp2,1},
     {"log",0,log,1},      {"log",0,log,1},
     {"negative",0,op_negative,1},      {"log2",0,log2,1},
       {"log10",0,log10,1},
       {"gamma",0,gamma,1},
       {"lgamma",0,lgamma,1},
       {"sqrt",0,sqrt,1},
       {"cbrt",0,cbrt,1},
       {"fabs",0,fabs,1},
       {"j0",0,j0,1},
       {"j1",0,j1,1},
       {"y0",0,y0,1},
       {"y1",0,y1,1},
       {"-",  0,op_negative,1},
     {"+",  0,op_add,2},      {"+",  0,op_add,2},
     {"-",  0,op_sub,2},      {"-",  0,op_sub,2},
     {"*",  0,op_mul,2},      {"*",  0,op_mul,2},
     {"/",  0,op_div,2},      {"/",  0,op_div,2},
     {"^",  0,pow,2},      {"^",  0,pow,2},
     {"e",  M_E, NULL,0},      {"pow",  0,pow,2},
     {"pi", M_PI,NULL,0},      {"@e",  M_E, NULL,0},
       {"@pi", M_PI,NULL,0},
     {NULL,0,NULL,0}      {NULL,0,NULL,0}
 };  };
   
Line 132  void replace2(int n, char *s[], double v[]) 
Line 154  void replace2(int n, char *s[], double v[]) 
     }      }
 }  }
   
 static entry *find_entry(cmo *node, entry *dic)  static entry *find_entry(cmo *node, int len, entry *dic)
 {  {
     char *s;      char *s;
     entry *e;      entry *e;
Line 144  static entry *find_entry(cmo *node, entry *dic)
Line 166  static entry *find_entry(cmo *node, entry *dic)
         return NULL;          return NULL;
     }      }
     for(e=dic; e->name != NULL; e++) {      for(e=dic; e->name != NULL; e++) {
         if(strcmp(e->name,s)==0) {          if(strcmp(e->name,s)==0 && (len<0 || len==e->n_args)) {
             return e;              return e;
         }          }
     }      }
Line 207  int entry_value(entry *e, double *retval)
Line 229  int entry_value(entry *e, double *retval)
 */  */
 static int eval_cmo_tree(cmo_tree* t, double *retval)  static int eval_cmo_tree(cmo_tree* t, double *retval)
 {  {
     entry *e = find_entry((cmo *)t->name,global_dic);      entry *e = find_entry((cmo *)t->name,list_length(t->leaves),global_dic);
         if (FUNCTION_P(e)) {      if (FUNCTION_P(e)) {
                 return entry_function(e,t->leaves,retval);          return entry_function(e,t->leaves,retval);
         }else if (VALUE_P(e)) {      }else if (VALUE_P(e)) {
                 return entry_value(e, retval);          return entry_value(e, retval);
         }      }
     return 0;      return 0;
 }  }
   
 static int eval_cmo_indeterminate(cmo_indeterminate *c, double *retval)  static int eval_cmo_indeterminate(cmo_indeterminate *c, double *retval)
 {  {
     entry *e = find_entry((cmo *)c,local_dic);      entry *e = find_entry((cmo *)c,-1,local_dic);
     if (VALUE_P(e)) {      if (VALUE_P(e)) {
         return entry_value(e,retval);          return entry_value(e,retval);
     }      }
Line 287  static int eval_cmo_recursive_polynomial(cmo_recursive
Line 309  static int eval_cmo_recursive_polynomial(cmo_recursive
                 }                  }
                 vars=(double *)calloc(n,sizeof(double));                  vars=(double *)calloc(n,sizeof(double));
                 for(i=0; i<n; i++) {                  for(i=0; i<n; i++) {
                         e = find_entry(list_nth(c->ringdef,i),local_dic);                          e = find_entry(list_nth(c->ringdef,i),-1,local_dic);
                         if(e == NULL) {                          if(e == NULL) {
                                 free(vars);                                  free(vars);
                                 return 0; /* failed */                                  return 0; /* failed */
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>