WebSVN – wimsdev – Diff – /trunk/wims/src/Lib/evalue.c

 /* log(-1) does not make sense in real */
 #ifndef NAN
 #define NAN log(-1)
 #endif
-        /* Only two decimal points, less than 1 million.
+    /* Only two decimal points, less than 1 million.
-         * No check of buffer length which should be at least 12.
+     * No check of buffer length which should be at least 12.
-         * returns the end of buffer. */
+     * returns the end of buffer. */
 char *moneyprint(char *p, double s)
+{
     char *p1, *p2, buf[16];
     int t, t1, t2;
     if(s<0) {*p++='-'; s=-s;}
     if(t==0) {*p++='0'; *p=0; return p;}
     t1=t/100; t2=t%100; p1=buf+10;
     for(*p1--=t1%10+'0',t1/=10;t1>0;*p1--=t1%10+'0',t1/=10);
     p2=buf+11;
     if(t2) {
-        *p2++='.';
+      *p2++='.';
-        *p2++=t2/10+'0'; t2%=10;
+      *p2++=t2/10+'0'; t2%=10;
-        if(t2) *p2++=t2+'0';
+      if(t2) *p2++=t2+'0';
+    }
     p1++; *p2=0; memmove(p,p1,p2-p1+1); p+=p2-p1;
     return p;
+}
 /* #define RAND_BUF_SIZE 4096
 static char rand_buf[RAND_BUF_SIZE];
+*/
-*/      /* The trouble here is that httpd does not initialize
+/* The trouble here is that httpd does not initialize
-         * the variable RANDOM.
+     * the variable RANDOM.
-         * So I use time (microseconds) to get a quick solution. */
+     * So I use time (microseconds) to get a quick solution. */
 void init_random(void)
+{
     int r;
     struct timeval t;
 /*    initstate(1,rand_buf,RAND_BUF_SIZE); */
     r=t.tv_usec+t.tv_sec*1000;
     if(r<0) r=-r; if(r==0) r=1;
     srandom(r);
+}
-        /* gives a double random number between 0 and m */
+/* gives a double random number between 0 and m */
 double drand(double m)
+{
     double r;
     r=((double) random()+(double) random()/(double) RAND_MAX);
     return (r/(double) RAND_MAX)*m;
+}
-        /* gives a random integer between 0 and n.
+/* gives a random integer between 0 and n.
-         * n maybe floating, but will be rounded */
+ * n maybe floating, but will be rounded */
 double irand(double n)
+{
     int  end,r;
     if(n==0) return 0;
     if(n>0) end=n; else end=-n;
     r=(double) random()*end/RAND_MAX;
     if(r==n) r--;
     if(n>0) return r; else return -r;
+}
-        /* sign of d */
+/* sign of d */
 double sign(double d)
+{
     if(d==0) return 0;
     if(d<0) return -1;
     else return 1;
+}
-        /* rounding to integer: problem with half-way rounding */
+/* rounding to integer: problem with half-way rounding */
 double myround(double d)
+{
     long int t;
     if(d<0) t=d-0.5; else t=d+0.5;
     return t;
+}
-        /* log of base 2 */
+/* log of base 2 */
 double mylog2(double d)
+{
     return log(d)/log(2);
+}
-        /* sec */
+/* sec */
 double sec(double d)
-{       return 1/cos(d);        }
+{    return 1/cos(d);}
-        /* csc */
+/* csc */
 double csc(double d)
-{       return 1/sin(d);        }
+{    return 1/sin(d);}
-        /* cotangent function */
+/* cotangent function */
 double cotan(double d)
+{
     return 1/tan(d);
+}
-        /* hyperbolic cotangent */
+/* hyperbolic cotangent */
 double cotanh(double d)
+{
     return 1/tanh(d);
+}
-        /* factorial of an integer */
+/* factorial of an integer */
 double factorial(double d)
+{
     int i,n; double t;
     n=d;
     if(n<0 || n!=d) return NAN;
     if(n>1000) return HUGE_VAL;
     t=1; for(i=1;i<=n;i++) t=t*i;
     return t;
+}
-        /* binomial coefficient */
+/* binomial coefficient */
 double binomial(double d1,double d2)
+{
     return factorial(d1)/(factorial(d2)*factorial(d1-d2));
+}
-        /* max and min */
+/* max and min */
 double max(double d1, double d2)
+{
     if(!finite(d1) || !finite(d2)) return NAN;
     if(d1<d2) return d2; else return d1;
+}
+{
     if(!finite(d1) || !finite(d2)) return NAN;
     if(d1<d2) return d1; else return d2;
+}
-        /* gcd and lcm, not really checking errors. */
+/* gcd and lcm, not really checking errors. */
 double gcd(double n1, double n2)
+{
     unsigned long long int l1, l2, ll;
     n1=abs(n1); n2=abs(n2);
     if(!finite(n1) || !finite(n2) || n1<0 || n2<0 ||
        n1>1E18 || n2>1E18) return NAN;
     l1=n1; l2=n2;
     if(l1<l2) {
-        ll=l1;l1=l2;l2=ll;
+      ll=l1;l1=l2;l2=ll;
+    }
     if(l1==0) return NAN;
     while(l2>0) {
-        ll=l2;l2=l1%l2;l1=ll;
+      ll=l2;l2=l1%l2;l1=ll;
+    }
     return l1;
+}
 double lcm(double n1, double n2)
 void set_evalue_pointer(char *p)
+{
     evalue_pt=p;
+}
-        /* get position of name in nametable */
+/* get position of name in nametable */
 int eval_getpos(char *name)
+{
     return search_list(evalname,evalname_no,sizeof(evalname[0]),name);
+}
-        /* set value to name */
+/* set value to name */
 void eval_setval(int pos, double v)
+{
     if(pos>=0 && pos<evalname_no) evalname[pos].val=v;
+}
-        /* get string pointer (after evaluation) */
+/* get string pointer (after evaluation) */
 char *get_evalue_pointer(void)
+{
     return evalue_pt;
+}
+{
     double d,dd;
     int i,k;
     char buf[32];
     if(evalue_error) return NAN;
     d=0;
     while(*evalue_pt=='+') evalue_pt++;
     if(*evalue_pt==0) return 0; /* empty string */
     switch(*evalue_pt) {
-        case '(':
+      case '(':
-        evalue_pt++; d=_evalue(')');goto vld;
+       evalue_pt++; d=_evalue(')');goto vld;
-        case '|':
+      case '|':
-        if(ord=='|') {
+       if(ord=='|') {
-            evalue_pt++; return 0;
+        evalue_pt++; return 0;
-        }
+       }
-        evalue_pt++; d=fabs(_evalue('|'));goto vld;
+       evalue_pt++; d=fabs(_evalue('|'));goto vld;
-        case '-':
+      case '-':
-        evalue_pt++; d=-_evalue(6);goto vld;
+       evalue_pt++; d=-_evalue(6);goto vld;
+    }
-    if((128&*evalue_pt)!=0) {
+    if((128&*evalue_pt)!=0) {/* special character */
-        k=(*evalue_pt)&255; evalue_pt++;
+      k=(*evalue_pt)&255; evalue_pt++;
-        if(k>=130 && k<140) {
+      if(k>=130 && k<140) {
-            i=(k-130)*200; k=(*evalue_pt)&255; evalue_pt++;
+          i=(k-130)*200; k=(*evalue_pt)&255; evalue_pt++;
-            if(k<33 || k>=233) goto badeval;
+          if(k<33 || k>=233) goto badeval;
-            i+=k-33; if(i<0 || i>=evalname_no) goto badeval;
+          i+=k-33; if(i<0 || i>=evalname_no) goto badeval;
-            goto ename;
+          goto ename;
-        }
+      }
-        if(k>=140 && k<150) {
+      if(k>=140 && k<150) {
-            i=(k-140)*200; k=(*evalue_pt)&255; evalue_pt++;
+          i=(k-140)*200; k=(*evalue_pt)&255; evalue_pt++;
-            if(k<33 || k>=233) goto badeval;
+          if(k<33 || k>=233) goto badeval;
-            if(ev_var==NULL || ev_varcnt==NULL) goto badeval;
+          if(ev_var==NULL || ev_varcnt==NULL) goto badeval;
-            i+=k-33; if(i<0 || i>=*ev_varcnt) goto badeval;
+          i+=k-33; if(i<0 || i>=*ev_varcnt) goto badeval;
-            goto vname;
+          goto vname;
-        }
+      }
-        evalue_pt++; goto badeval;
+      evalue_pt++; goto badeval;
+    }
     if(*evalue_pt=='.' || myisdigit(*evalue_pt))
-        {d=strtod(evalue_pt,&evalue_pt);goto binary;}
+      {d=strtod(evalue_pt,&evalue_pt);goto binary;}
     for(i=0;myisalnum(*(evalue_pt+i)) && i<16; i++)
       buf[i]=*(evalue_pt+i);
     buf[i]=0; evalue_pt+=i;
     if(i==0) goto badeval;
     if(ev_varcnt!=NULL && ev_var!=NULL && *ev_varcnt>0)
       for(i=0;i<*ev_varcnt;i++) {
-          if(strcmp(buf,ev_var[i].name)==0) {
+        if(strcmp(buf,ev_var[i].name)==0) {
-              vname: d=ev_var[i].value; goto vld;
+            vname: d=ev_var[i].value; goto vld;
-          }
+        }
       }
     i=search_list(evalname,evalname_no,sizeof(evalname[0]),buf);
-    ename: if(i>=0) switch(evalname[i].type) {
+    ename: if(i>=0) switch(evalname[i].type) { /* evaluation of expressions */
-        case 0: {
+      case 0: {
-            d=evalname[i].val;
+          d=evalname[i].val;
-            if(evalname[i].f1!=NULL) {
+          if(evalname[i].f1!=NULL) {
-                if(d==0) d=NAN; if(d==1) d=HUGE_VAL;
+            if(d==0) d=NAN; if(d==1) d=HUGE_VAL;
-            }
+          }
-            break;
+          break;
-        }
+      }
-        case 1: {
+      case 1: {
-            if(*evalue_pt!='(') return NAN;
+          if(*evalue_pt!='(') return NAN;
-            evalue_pt++;
+          evalue_pt++;
-            d=evalname[i].f1(_evalue(')')); break;
+          d=evalname[i].f1(_evalue(')')); break;/* evaluation of function */
-        }
+      }
-        case 2: {
+      case 2: {
-            double parm1,parm2;
+          double parm1,parm2;
-            if(*evalue_pt!='(') return NAN;
+          if(*evalue_pt!='(') return NAN;
-            evalue_pt++;
+          evalue_pt++;
-            parm1=_evalue(',');parm2=_evalue(')');
+          parm1=_evalue(',');parm2=_evalue(')');
-            d=evalname[i].f2(parm1,parm2); break;
+          d=evalname[i].f2(parm1,parm2); break;
-        }
+      }
-        default: {      /* This is impossible. */
+      default: {      /* This is impossible. */
-            return NAN;
+          return NAN;
-        }
+      }
     }
     else {
-        badeval: evalue_error=-1; return NAN;
+      badeval: evalue_error=-1; return NAN;
+    }
   vld:
     if(evalue_error) return NAN;
-  binary:
+  binary: /*evaluation des expressions */
     if(*evalue_pt=='!') {
-        evalue_pt++; d=factorial(d);
+      evalue_pt++; d=factorial(d);
     }
-    if(*evalue_pt==ord) {evalue_pt++;goto ok;}
+    if(*evalue_pt==ord) {evalue_pt++;goto ok;}/* */
-    if(*evalue_pt==0 ||
+    if(*evalue_pt==0 || /* chaine de caractere finie*/
        (ord<10 && (*evalue_pt==',' || *evalue_pt==';' || *evalue_pt==')'
-                   || *evalue_pt=='|')))
+               || *evalue_pt=='|')))
        goto ok;
     switch(*evalue_pt) {
-                case '+':
+      case '+':
-                        if(ord<=8) break;
+        if(ord<=8) break;
-                        evalue_pt++; d+=_evalue(8);goto vld;
+        evalue_pt++; d+=_evalue(8);goto vld;
-                case '-':
+      case '-':
-                        if(ord<=8) break;
+        if(ord<=8) break;
-                        evalue_pt++; d-=_evalue(8);goto vld;
+        evalue_pt++; d-=_evalue(8);goto vld;
-                case '*':
+      case '*':
-                        if(ord<=6) break;
+        if(ord<=6) break;
-                        evalue_pt++; d*=_evalue(6);goto vld;
+        evalue_pt++; d*=_evalue(6);goto vld;
-                case '/':
+      case '/':
-                        if(ord<=6) break;
+        if(ord<=6) break;
-                        evalue_pt++; dd=_evalue(6);
+        evalue_pt++; dd=_evalue(6);
-                        if(dd==0) {evalue_error=10;return NAN;}
+        if(dd==0) {evalue_error=10;return NAN;}
-                        d/=dd;goto vld;
+            d/=dd;goto vld;
-        case '%': {
+      case '%': {
-            int di, ddi;
+        int di, ddi;
-            if(ord<=6) break;
+        if(ord<=6) break;
-            evalue_pt++; dd=_evalue(6);
+        evalue_pt++; dd=_evalue(6);
-            if(dd==0) {evalue_error=10;return NAN;}
+        if(dd==0) {evalue_error=10;return NAN;}
-            di=d; ddi=dd; d=di%ddi;goto vld;
+          di=d; ddi=dd; d=di%ddi;goto vld;
-        }
+      }
-        case '^': {
+      case '^': {
-            if(ord<5) break;
+       if(ord<5) break;
-            evalue_pt++; d=pow(d,_evalue(5));goto vld;
+       evalue_pt++; d=pow(d,_evalue(5));goto vld;
-        }
+      }
-        default : {
+      default : {
-            /*if(ord<=6) break;
-            d*=_evalue(6);goto vld;*/
-            return NAN;
+          return NAN;
-        }
+      }
+    }
     ok: return d;
 }
-        /* substitute variable names by their environment strings
+/* substitute variable names by their environment strings
-         * The buffer pointed to by p must have enough space
+ * The buffer pointed to by p must have enough space
-         * (defined by MAX_LINELEN). */
+ * (defined by MAX_LINELEN). */
 char *_substit(char *p)
+{
     return p;
+}
 char *(*substitute) (char *p)=_substit;
-        /* evalue a string to double */
+/* evalue a string to double */
 double strevalue(char *p)
+{
     char buf[MAX_LINELEN+1];
     if(p==NULL) return 0;
     mystrncpy(buf,p,sizeof(buf));
     substitute(buf); nospace(buf);
     if(check_parentheses(buf,0)) return NAN;
     set_evalue_error(0);
     set_evalue_pointer(buf);
     return _evalue(10);
+}
-        /* compile an expression for faster evaluation
+/* compile an expression for faster evaluation
-         * returns -1 if cannot be compiled.
+ * returns -1 if cannot be compiled.
-         * else returns the number of compilations. */
+ * else returns the number of compilations.
+ */
 int evalue_compile(char *p)
+{
     char *p1, *p2, *pe, name[256], buf[8];
     int i,k;
     k=0;
     for(p1=p; *p1; p1++) if((128&*p1)!=0) return -1;
     nospace(p);
     for(p1=find_mathvar_start(p); *p1; p1=find_mathvar_start(pe)) {
-        pe=find_mathvar_end(p1);
+      pe=find_mathvar_end(p1);
-        if(!myisalpha(*p1)) continue;
+      if(!myisalpha(*p1)) continue;
-        p2=pe; if(p2-p1>16) continue;
+      p2=pe; if(p2-p1>16) continue;
-        memmove(name,p1,p2-p1); name[p2-p1]=0;
+      memmove(name,p1,p2-p1); name[p2-p1]=0;
-        if(ev_varcnt!=NULL && ev_var!=NULL && *ev_varcnt>0) {
+/* replace the variables by a number
-            for(i=0;i<*ev_varcnt && strcmp(name,ev_var[i].name)!=0;i++);
-            if(i<*ev_varcnt && i<2000) {
+ * at most 2000 variables on two characters :
-                buf[0]=i/200+140; buf[1]=i%200+33; buf[2]=0;
+ * variable: 140 <= integer <150,  number between 33 and 233
-                string_modify(p,p1,p2,"%s",buf);
+ * function: 130 <= integer < 140, number between 33 and 233
-                pe=p1+2; k++; continue;
-            }
-        }
+ */
+    if(ev_varcnt!=NULL && ev_var!=NULL && *ev_varcnt>0) {
-        i=search_list(evalname,evalname_no,sizeof(evalname[0]),name);
+        for(i=0;i<*ev_varcnt && strcmp(name,ev_var[i].name)!=0;i++);
-        if(i>=0 && i<2000) {
+        if(i<*ev_varcnt && i<2000) {
-            buf[0]=i/200+130; buf[1]=i%200+33; buf[2]=0;
+        buf[0]=i/200+140; buf[1]=i%200+33; buf[2]=0;
-            string_modify(p,p1,p2,"%s",buf);
+        string_modify(p,p1,p2,"%s",buf);
-            pe=p1+2; k++; continue;
+        pe=p1+2; k++; continue;
         }
+    }
+    i=search_list(evalname,evalname_no,sizeof(evalname[0]),name);
+    if(i>=0 && i<2000) {
+        buf[0]=i/200+130; buf[1]=i%200+33; buf[2]=0;
+        string_modify(p,p1,p2,"%s",buf);
+        pe=p1+2; k++; continue;
+    }
+  }
-    return k;
+  return k;
+}

Subversion Repositories wimsdev

(root)/trunk/wims/src/Lib/evalue.c – Rev 10 → 7798