[BACK]Return to quote.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / parse

Diff for /OpenXM_contrib2/asir2000/parse/quote.c between version 1.1 and 1.4

version 1.1, 2001/08/21 00:47:27 version 1.4, 2001/09/04 05:14:04
Line 1 
Line 1 
   /* $OpenXM: OpenXM_contrib2/asir2000/parse/quote.c,v 1.3 2001/09/04 02:46:15 noro Exp $ */
   
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
   
Line 66  QUOTE *c;
Line 68  QUOTE *c;
   
         fn = mkfnode(3,I_BOP,subfs,0,BDY(a));          fn = mkfnode(3,I_BOP,subfs,0,BDY(a));
         MKQUOTE(*c,fn);          MKQUOTE(*c,fn);
   }
   
   void polytoquote(), dctoquote(), vartoquote();
   void dptoquote(), mptoquote();
   
   void objtoquote(a,c)
   Obj a;
   QUOTE *c;
   {
           QUOTE nm,dn;
           NODE arg,t0,t,t1,t2,t3;
           FNODE fn;
           Obj *b;
           Obj **m;
           int i,j,len,row,col;
           Q q,qrow,qcol;
           FUNC f;
   
           if ( !a ) {
                   MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a));
                   return;
           }
           switch ( OID(a) ) {
                   case O_N:
                   case O_STR:
                           MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a));
                           break;
                   case O_P:
                           polytoquote((P)a,c);
                           break;
                   case O_R:
                           polytoquote(NM((R)a),&nm);
                           polytoquote(DN((R)a),&dn);
                           divquote(CO,nm,dn,c);
                           break;
                   case O_LIST:
                           t0 = 0;
                           for ( arg = BDY((LIST)a); arg; arg = NEXT(arg) ) {
                                   NEXTNODE(t0,t);
                                   objtoquote(BDY(arg),&nm);
                                   BDY(t) = BDY(nm);
                           }
                           if ( t0 )
                                   NEXT(t) = 0;
                           MKQUOTE(*c,mkfnode(1,I_LIST,t0));
                           break;
                   case O_VECT:
                           len = ((VECT)a)->len;
                           b = (Obj *)BDY(((VECT)a));
                           t = 0;
                           for ( i = len-1; i >= 0; i-- ) {
                                   objtoquote(b[i],&nm);
                                   MKNODE(t1,BDY(nm),t);
                                   t = t1;
                           }
                           STOQ(len,q);
                           t = mknode(2,mkfnode(1,I_FORMULA,q),mkfnode(1,I_LIST,t));
                           gen_searchf("vector",&f);
                           MKQUOTE(*c,mkfnode(2,I_FUNC,f,mkfnode(1,I_LIST,t)));
                           break;
                   case O_MAT:
                           row = ((MAT)a)->row;
                           col = ((MAT)a)->row;
                           m = (Obj **)BDY(((MAT)a));
                           t2 = 0;
                           for ( i = row-1; i >= 0; i-- ) {
                                   t = 0;
                                   for ( j = col-1; j >= 0; j-- ) {
                                           objtoquote(m[i][j],&nm);
                                           MKNODE(t1,BDY(nm),t);
                                           t = t1;
                                   }
                                   fn = mkfnode(1,I_LIST,t);
                                   MKNODE(t3,fn,t2);
                                   t2 = t3;
                           }
                           fn = mkfnode(1,I_LIST,t2);
   
                           STOQ(row,qrow);
                           STOQ(col,qcol);
                           t = mknode(3,
                                   mkfnode(1,I_FORMULA,qrow),mkfnode(1,I_FORMULA,qcol),fn);
                           gen_searchf("matrix",&f);
                           MKQUOTE(*c,mkfnode(2,I_FUNC,f,mkfnode(1,I_LIST,t)));
                           break;
                   case O_DP:
                           dptoquote((DP)a,c);
                           break;
                   case O_QUOTE:
                           *c = (QUOTE)a;
                           break;
                   default:
                           error("objtoquote : not implemented");
           }
   }
   
   void polytoquote(a,c)
   P a;
   QUOTE *c;
   {
           DCP dc,t;
           DCP *dca;
           int n,i;
           QUOTE v,r,s,u;
   
           if ( !a || (OID(a) == O_N) ) {
                   MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a));
                   return;
           }
           dc = DC((P)a);
           vartoquote(VR((P)a),&v);
           for ( t = dc, n = 0; t; t = NEXT(t), n++ );
           dca = (DCP *)ALLOCA(n*sizeof(DCP));
           for ( t = dc, i = 0; t; t = NEXT(t), i++ )
                   dca[i] = t;
           dctoquote(dca[n-1],v,&r);
           for ( i = n-2; i >= 0; i-- ) {
                   dctoquote(dca[i],v,&s);
                   addquote(CO,s,r,&u);
                   r = u;
           }
           *c = r;
   }
   
   void dptoquote(a,c)
   DP a;
   QUOTE *c;
   {
           MP t;
           MP *m;
           int i,n,nv;
           QUOTE s,r,u;
   
           if ( !a ) {
                   MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)a));
                   return;
           }
           for ( t = BDY(a), n = 0; t; t = NEXT(t), n++ );
           m = (MP *)ALLOCA(n*sizeof(MP));
           for ( t = BDY(a), i = 0; t; t = NEXT(t), i++ )
                   m[i] = t;
           nv = NV(a);
           mptoquote(m[n-1],nv,&r);
           for ( i = n-2; i >= 0; i-- ) {
                   mptoquote(m[i],nv,&s);
                   addquote(CO,s,r,&u);
                   r = u;
           }
           *c = r;
   }
   
   void dctoquote(dc,v,c)
   DCP dc;
   QUOTE v;
   QUOTE *c;
   {
           QUOTE r,d,s,u;
   
           if ( UNIQ(COEF(dc)) ) {
                   if ( DEG(dc) ) {
                           if ( UNIQ(DEG(dc)) )
                                   *c = v;
                           else {
                                   objtoquote(DEG(dc),&d);
                                   pwrquote(CO,v,d,c);
                           }
                   } else
                           objtoquote(ONE,c);
           } else {
                   objtoquote(COEF(dc),&u);
                   if ( DEG(dc) ) {
                           if ( UNIQ(DEG(dc)) )
                                   s = v;
                           else {
                                   objtoquote(DEG(dc),&d);
                                   pwrquote(CO,v,d,&s);
                           }
                           mulquote(CO,u,s,c);
                   } else
                           *c = u;
           }
   }
   
   void mptoquote(m,n,c)
   MP m;
   int n;
   QUOTE *c;
   {
           QUOTE s,u;
           NODE t,t1;
           FNODE f;
           Q q;
           DL dl;
           int i;
   
           objtoquote(C(m),&s);
           dl = m->dl;
           for ( i = n-1; i >= 0; i-- ) {
                   STOQ(dl->d[i],q);
                   f = mkfnode(1,I_FORMULA,q);
                   MKNODE(t1,f,t);
                   t = t1;
           }
           MKQUOTE(u,mkfnode(1,I_EV,t));
           mulquote(CO,s,u,c);
   }
   
   void vartoquote(v,c)
   V v;
   QUOTE *c;
   {
           P x;
           PF pf;
           PFAD ad;
           QUOTE a,b;
           int i;
           FUNC f;
           NODE t,t1;
   
           if ( NAME(v) ) {
                   MKV(v,x);
                   MKQUOTE(*c,mkfnode(1,I_FORMULA,(pointer)x));
           } else if ( (vid)v->attr == V_PF ) {
                   /* pure function */
                   pf = ((PFINS)v->priv)->pf;
                   ad = ((PFINS)v->priv)->ad;
                   if ( !strcmp(NAME(pf),"pow") ) {
                           /* pow(a,b) = a^b */
                           objtoquote(ad[0].arg,&a); objtoquote(ad[1].arg,&b);
                           pwrquote(CO,a,b,c);
                   } else {
                           for ( i = 0; i < pf->argc; i++ )
                                   if ( ad[i].d )
                                           break;
                           if ( i < pf->argc )
                                   error("vartoquote : not implemented");
                           gen_searchf(NAME(pf),&f);
                           t = 0;
                           for ( i = pf->argc-1; i >= 0; i-- ) {
                                   objtoquote(ad[i].arg,&a);
                                   MKNODE(t1,BDY(a),t);
                                   t = t1;
                           }
                           MKQUOTE(*c,mkfnode(2,I_FUNC,f,mkfnode(1,I_LIST,t)));
                   }
           }
 }  }
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>