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

Diff for /OpenXM_contrib2/asir2000/builtin/subst.c between version 1.4 and 1.10

version 1.4, 2001/09/05 04:43:58 version 1.10, 2015/09/24 04:43:13
Line 45 
Line 45 
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,   * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.   * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *   *
  * $OpenXM: OpenXM_contrib2/asir2000/builtin/subst.c,v 1.3 2000/08/22 05:04:00 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/subst.c,v 1.9 2010/01/31 03:25:54 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
   
 void Psubst(), Ppsubst(), Psubstf(), Psubst_quote();  void Psubst(), Ppsubst(), Psubstf(), Psubst_quote();
   void Psubstr2np();
   
 struct ftab subst_tab[] = {  struct ftab subst_tab[] = {
         {"subst",Psubst,-99999999},          {"subst",Psubst,-99999999},
           {"substr2np",Psubstr2np,-3},
         {"subst_quote",Psubst_quote,-99999999},          {"subst_quote",Psubst_quote,-99999999},
         {"psubst",Ppsubst,-99999999},          {"psubst",Ppsubst,-99999999},
         {"substf",Psubstf,-99999999},          {"substf",Psubstf,-99999999},
         {0,0,0},          {0,0,0},
 };  };
   
   extern Obj VOIDobj;
   
   /* substr2np(P,[[v,a],...]) or substr2np(P,[v1,...],[a1,...]) */
   
   void Psubstr2np(NODE arg,Obj *rp)
   {
           Obj a;
           P nm,dn,p;
           R r;
           VL vl,tvl;
           int nv,i,ac;
           NODE slist,t,ps,u,vlist;
           P s;
           P *svect;
           V v;
           V *vvect;
   
           a = (Obj)ARG0(arg);
           if ( !a || NUM(a) ) {
                   *rp = a;
                   return;
           }
           asir_assert(ARG0(arg),O_R,"substr2np");
           asir_assert(ARG1(arg),O_LIST,"substr2np");
           get_vars(a,&vl);
           for ( i = 0, tvl = vl; tvl; tvl = NEXT(tvl), i++ );
           nv = i;
           vvect = (V *)MALLOC(nv*sizeof(V));
           for ( i = 0, tvl = vl; tvl; tvl = NEXT(tvl), i++ ) vvect[i] = tvl->v;
           svect = (P *)MALLOC(nv*sizeof(P));
       ac = argc(arg);
       if ( ac == 2 ) {
             slist = BDY((LIST)ARG1(arg));
             for ( i = 0; i < nv; i++ ) svect[i] = (P)VOIDobj;
             for ( t = slist; t; t = NEXT(t) ) {
                     ps = BDY((LIST)BDY(t)); p = (P)BDY(ps); s = (P)BDY(NEXT(ps));
                     asir_assert(p,O_P,"substr2np"); asir_assert(s,O_N,"substr2np");
                     v = VR(p);
                     for ( i = 0; i < nv; i++ ) if ( vvect[i] == v ) break;
                     svect[i] = s;
             }
       } else if ( ac == 3 ) {
             asir_assert(ARG2(arg),O_LIST,"substr2np");
         vlist = BDY((LIST)ARG1(arg));
         slist = BDY((LIST)ARG2(arg));
             for ( i = 0; i < nv; i++ ) svect[i] = (P)VOIDobj;
         for ( u = vlist, t = slist; u && t; u = NEXT(u), t = NEXT(t) ) {
           v = VR((P)BDY(u));
           for ( i = 0; i < nv; i++ ) if ( vvect[i] == v ) break;
           svect[i] = (P)BDY(t);
         }
       } else
         error("substr2np : argument mismatch");
   
           switch ( OID(a) ) {
                   case O_P:
                           substpp(vl,(P)a,vvect,svect,nv,&nm); *rp = (Obj)nm;
                           return;
                   case O_R:
                           substpp(vl,(P)NM((R)a),vvect,svect,nv,&nm);
                           substpp(vl,(P)DN((R)a),vvect,svect,nv,&dn);
                           if ( !dn )
                                   error("substr2np: division by 0");
                           else if ( !nm )
                                   *rp = 0;
                           else if ( NUM(dn) ) {
                                   divsp(vl,nm,dn,&p);
                                   *rp = (Obj)p;
                           } else {
                                   MKRAT(nm,dn,0,r);
                                   *rp = (Obj)r;
                           }
                           return;
                   default:
                           error("substr2np: invalid argument");
           }
   }
   
 void Psubst(arg,rp)  void Psubst(arg,rp)
 NODE arg;  NODE arg;
 Obj *rp;  Obj *rp;
Line 67  Obj *rp;
Line 147  Obj *rp;
         Obj a,b,t;          Obj a,b,t;
         LIST l;          LIST l;
         V v;          V v;
           int row,col,len;
           VECT vect;
           MAT mat;
           int i,j;
           NODE n0,n,nd;
           struct oNODE arg0;
           MP m,mp,mp0;
           DP d;
   
         if ( !arg ) {          if ( !arg ) {
                 *rp = 0; return;                  *rp = 0; return;
         }          }
         asir_assert(ARG0(arg),O_R,"subst");          a = (Obj)ARG0(arg);
         reductr(CO,(Obj)ARG0(arg),&a);          if ( !a ) {
 /*      a = (Obj)ARG0(arg); */                  *rp = 0;
         arg = NEXT(arg);                  return;
         if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )          }
                 arg = BDY(l);          switch ( OID(a) ) {
         while ( arg ) {                  case O_N: case O_P: case O_R:
                 asir_assert(BDY(arg),O_P,"subst");                          reductr(CO,(Obj)ARG0(arg),&a);
                 v = VR((P)BDY(arg)); arg = NEXT(arg);                          arg = NEXT(arg);
                 if ( !arg )                          if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )
                                   arg = BDY(l);
                           while ( arg ) {
                                   if ( !BDY(arg) || OID((Obj)BDY(arg)) != O_P )
                                           error("subst : invalid argument");
                                   v = VR((P)BDY(arg)); arg = NEXT(arg);
                                   if ( !arg )
                                           error("subst : invalid argument");
                                   asir_assert(ARG0(arg),O_R,"subst");
                                   reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);
                   /*              b = (Obj)BDY(arg); arg = NEXT(arg); */
                                   substr(CO,0,a,v,b,&t); a = t;
                           }
                           *rp = a;
                           break;
                   case O_LIST:
                           n0 = 0;
                           for ( nd = BDY((LIST)a); nd; nd = NEXT(nd) ) {
                                   NEXTNODE(n0,n);
                                   arg0.body = (pointer)BDY(nd);
                                   arg0.next = NEXT(arg);
                                   Psubst(&arg0,&b);
                                   BDY(n) = (pointer)b;
                           }
                           if ( n0 )
                                   NEXT(n) = 0;
                           MKLIST(l,n0);
                           *rp = (Obj)l;
                           break;
                   case O_VECT:
                           len = ((VECT)a)->len;
                           MKVECT(vect,len);
                           for ( i = 0; i < len; i++ ) {
                                   arg0.body = (pointer)BDY((VECT)a)[i];
                                   arg0.next = NEXT(arg);
                                   Psubst(&arg0,&b);
                                   BDY(vect)[i] = (pointer)b;
                           }
                           *rp = (Obj)vect;
                           break;
                   case O_MAT:
                           row = ((MAT)a)->row;
                           col = ((MAT)a)->col;
                           MKMAT(mat,row,col);
                           for ( i = 0; i < row; i++ )
                                   for ( j = 0; j < col; j++ ) {
                                           arg0.body = (pointer)BDY((MAT)a)[i][j];
                                           arg0.next = NEXT(arg);
                                           Psubst(&arg0,&b);
                                           BDY(mat)[i][j] = (pointer)b;
                                   }
                           *rp = (Obj)mat;
                           break;
                   case O_DP:
                           mp0 = 0;
                           for ( m = BDY((DP)a); m; m = NEXT(m) ) {
                                   arg0.body = (pointer)C(m);
                                   arg0.next = NEXT(arg);
                                   Psubst(&arg0,&b);
                                   if ( b ) {
                                           NEXTMP(mp0,mp);
                                           C(mp) = (P)b;
                                           mp->dl = m->dl;
                                   }
                           }
                           if ( mp0 ) {
                                   MKDP(NV((DP)a),mp0,d);
                                   d->sugar = ((DP)a)->sugar;
                                   *rp = (Obj)d;
                           } else
                                   *rp = 0;
   
                           break;
                   default:
                         error("subst : invalid argument");                          error("subst : invalid argument");
                 asir_assert(ARG0(arg),O_R,"subst");  
                 reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);  
 /*              b = (Obj)BDY(arg); arg = NEXT(arg); */  
                 substr(CO,0,a,v,b,&t); a = t;  
         }          }
         *rp = a;  
 }  }
   
 FNODE subst_in_fnode();  FNODE subst_in_fnode();
Line 96  void Psubst_quote(arg,rp)
Line 252  void Psubst_quote(arg,rp)
 NODE arg;  NODE arg;
 QUOTE *rp;  QUOTE *rp;
 {  {
         QUOTE a,h;          QUOTE h;
         FNODE fn;          FNODE fn;
         Obj g;          Obj g;
         LIST l;          LIST l;
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.10
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>