[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.9 and 1.13

version 1.9, 2010/01/31 03:25:54 version 1.13, 2018/03/29 01:32:50
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.8 2010/01/28 08:56:26 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/subst.c,v 1.12 2018/03/27 06:29:19 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 54  void Psubst(), Ppsubst(), Psubstf(), Psubst_quote();
Line 54  void Psubst(), Ppsubst(), Psubstf(), Psubst_quote();
 void Psubstr2np();  void Psubstr2np();
   
 struct ftab subst_tab[] = {  struct ftab subst_tab[] = {
         {"subst",Psubst,-99999999},    {"subst",Psubst,-99999999},
         {"substr2np",Psubstr2np,2},    {"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;  extern Obj VOIDobj;
   
   /* substr2np(P,[[v,a],...]) or substr2np(P,[v1,...],[a1,...]) */
   
 void Psubstr2np(NODE arg,Obj *rp)  void Psubstr2np(NODE arg,Obj *rp)
 {  {
         Obj a;    Obj a,b;
         P nm,dn,p;    P nm,dn,p,q;
         R r;    R r;
         VL vl,tvl;    VL vl,tvl;
         int nv,i;    int nv,i,ac,j;
         NODE slist,t,ps;    NODE slist,t,ps,u,vlist;
         P s;    P s;
         P *svect;    P *svect;
         V v;    V v;
         V *vvect;    V *vvect;
   
         a = (Obj)ARG0(arg);    a = (Obj)ARG0(arg);
         if ( !a || NUM(a) ) {    if ( !a || NUM(a) ) {
                 *rp = a;      *rp = a;
                 return;      return;
         }    }
         asir_assert(ARG0(arg),O_R,"substr2np");    asir_assert(ARG0(arg),O_R,"substr2np");
         asir_assert(ARG1(arg),O_LIST,"substr2np");    asir_assert(ARG1(arg),O_LIST,"substr2np");
         get_vars(a,&vl);    get_vars(a,&vl);
         for ( i = 0, tvl = vl; tvl; tvl = NEXT(tvl), i++ );    for ( i = 0, tvl = vl; tvl; tvl = NEXT(tvl), i++ );
         nv = i;    nv = i;
         vvect = (V *)MALLOC(nv*sizeof(V));    vvect = (V *)MALLOC((nv)*sizeof(V));
         for ( i = 0, tvl = vl; tvl; tvl = NEXT(tvl), i++ ) vvect[i] = tvl->v;    for ( i = 0, tvl = vl; tvl; tvl = NEXT(tvl), i++ ) vvect[i] = tvl->v;
         svect = (P *)MALLOC(nv*sizeof(P));    svect = (P *)MALLOC((nv)*sizeof(P));
         slist = BDY((LIST)ARG1(arg));    ac = argc(arg);
         for ( i = 0; i < nv; i++ ) svect[i] = (P)VOIDobj;    if ( ac == 2 ) {
         for ( t = slist; t; t = NEXT(t) ) {      slist = BDY((LIST)ARG1(arg));
                 ps = BDY((LIST)BDY(t)); p = (P)BDY(ps); s = (P)BDY(NEXT(ps));      for ( i = 0; i < nv; i++ ) svect[i] = (P)VOIDobj;
                 asir_assert(p,O_P,"substr2np"); asir_assert(s,O_N,"substr2np");      for ( t = slist; t; t = NEXT(t) ) {
                 v = VR(p);        ps = BDY((LIST)BDY(t)); p = (P)BDY(ps); s = (P)BDY(NEXT(ps));
                 for ( i = 0; i < nv; i++ ) if ( vvect[i] == v ) break;        asir_assert(p,O_P,"substr2np"); asir_assert(s,O_P,"substr2np");
                 svect[i] = s;        v = VR(p);
         }        for ( i = 0; i < nv; i++ ) if ( vvect[i] == v ) break;
         svect[i] = s;
         switch ( OID(a) ) {      }
                 case O_P:    } else if ( ac == 3 ) {
                         substpp(vl,(P)a,vvect,svect,nv,&nm); *rp = (Obj)nm;      asir_assert(ARG2(arg),O_LIST,"substr2np");
                         return;      vlist = BDY((LIST)ARG1(arg));
                 case O_R:      slist = BDY((LIST)ARG2(arg));
                         substpp(vl,(P)NM((R)a),vvect,svect,nv,&nm);      for ( i = 0; i < nv; i++ ) svect[i] = (P)VOIDobj;
                         substpp(vl,(P)DN((R)a),vvect,svect,nv,&dn);        for ( u = vlist, t = slist; u && t; u = NEXT(u), t = NEXT(t) ) {
                         if ( !dn )          v = VR((P)BDY(u));
                                 error("substr2np: division by 0");          for ( i = 0; i < nv; i++ ) if ( vvect[i] == v ) break;
                         else if ( !nm )          svect[i] = (P)BDY(t);
                                 *rp = 0;        }
                         else if ( NUM(dn) ) {      } else
                                 divsp(vl,nm,dn,&p);        error("substr2np : argument mismatch");
                                 *rp = (Obj)p;    for ( i = 0; i < nv; i++ ) {
                         } else {      if ( (int)(vvect[i]->attr) == V_PF ) {
                                 MKRAT(nm,dn,0,r);        MKV(vvect[i],p);
                                 *rp = (Obj)r;        for ( j = 0; j < nv; j++ )
                         }          if ( j != i ) {
                         return;            substr(CO,0,(Obj)p,vvect[j],(Obj)svect[j],&b); p = (P)b;
                 default:          }
                         error("substr2np: invalid argument");        if ( OID(svect[i]) == O_VOID ) svect[i] = p;
         }        else if ( arf_comp(CO,(Obj)p,(Obj)svect[i]) )
           error("substr2np : inconsistent values for substitution");
       }
     }
     switch ( OID(a) ) {
       case O_P:
         substpp(CO,(P)a,vvect,svect,nv,&nm); *rp = (Obj)nm;
         return;
       case O_R:
         substpp(CO,(P)NM((R)a),vvect,svect,nv,&nm);
         substpp(CO,(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(CO,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;
 {  {
         Obj a,b,t;    Obj a,b,t;
         LIST l;    LIST l;
         V v;    V v;
         int row,col,len;    int row,col,len;
         VECT vect;    VECT vect;
         MAT mat;    MAT mat;
         int i,j;    int i,j;
         NODE n0,n,nd;    NODE n0,n,nd;
         struct oNODE arg0;    struct oNODE arg0;
         MP m,mp,mp0;    MP m,mp,mp0;
         DP d;    DP d;
     VL lastvl,vl,tvl,prev,cur;
   
         if ( !arg ) {    if ( !arg ) {
                 *rp = 0; return;      *rp = 0; return;
         }    }
         a = (Obj)ARG0(arg);    a = (Obj)ARG0(arg);
         if ( !a ) {    if ( !a ) {
                 *rp = 0;      *rp = 0;
                 return;      return;
         }    }
         switch ( OID(a) ) {    lastvl = LASTCO;
                 case O_N: case O_P: case O_R:    switch ( OID(a) ) {
                         reductr(CO,(Obj)ARG0(arg),&a);      case O_N: case O_P: case O_R:
                         arg = NEXT(arg);        reductr(CO,(Obj)ARG0(arg),&a);
                         if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )        arg = NEXT(arg);
                                 arg = BDY(l);        if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )
                         while ( arg ) {          arg = BDY(l);
                                 if ( !BDY(arg) || OID((Obj)BDY(arg)) != O_P )        while ( arg ) {
                                         error("subst : invalid argument");          if ( !BDY(arg) || OID((Obj)BDY(arg)) != O_P )
                                 v = VR((P)BDY(arg)); arg = NEXT(arg);            error("subst : invalid argument");
                                 if ( !arg )          v = VR((P)BDY(arg)); arg = NEXT(arg);
                                         error("subst : invalid argument");          if ( !arg )
                                 asir_assert(ARG0(arg),O_R,"subst");            error("subst : invalid argument");
                                 reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);          asir_assert(ARG0(arg),O_R,"subst");
                 /*              b = (Obj)BDY(arg); arg = NEXT(arg); */          reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);
                                 substr(CO,0,a,v,b,&t); a = t;      /*    b = (Obj)BDY(arg); arg = NEXT(arg); */
                         }          substr(CO,0,a,v,b,&t); a = t;
                         *rp = a;        }
                         break;        *rp = a;
                 case O_LIST:        break;
                         n0 = 0;      case O_LIST:
                         for ( nd = BDY((LIST)a); nd; nd = NEXT(nd) ) {        n0 = 0;
                                 NEXTNODE(n0,n);        for ( nd = BDY((LIST)a); nd; nd = NEXT(nd) ) {
                                 arg0.body = (pointer)BDY(nd);          NEXTNODE(n0,n);
                                 arg0.next = NEXT(arg);          arg0.body = (pointer)BDY(nd);
                                 Psubst(&arg0,&b);          arg0.next = NEXT(arg);
                                 BDY(n) = (pointer)b;          Psubst(&arg0,&b);
                         }          BDY(n) = (pointer)b;
                         if ( n0 )        }
                                 NEXT(n) = 0;        if ( n0 )
                         MKLIST(l,n0);          NEXT(n) = 0;
                         *rp = (Obj)l;        MKLIST(l,n0);
                         break;        *rp = (Obj)l;
                 case O_VECT:        break;
                         len = ((VECT)a)->len;      case O_VECT:
                         MKVECT(vect,len);        len = ((VECT)a)->len;
                         for ( i = 0; i < len; i++ ) {        MKVECT(vect,len);
                                 arg0.body = (pointer)BDY((VECT)a)[i];        for ( i = 0; i < len; i++ ) {
                                 arg0.next = NEXT(arg);          arg0.body = (pointer)BDY((VECT)a)[i];
                                 Psubst(&arg0,&b);          arg0.next = NEXT(arg);
                                 BDY(vect)[i] = (pointer)b;          Psubst(&arg0,&b);
                         }          BDY(vect)[i] = (pointer)b;
                         *rp = (Obj)vect;        }
                         break;        *rp = (Obj)vect;
                 case O_MAT:        break;
                         row = ((MAT)a)->row;      case O_MAT:
                         col = ((MAT)a)->col;        row = ((MAT)a)->row;
                         MKMAT(mat,row,col);        col = ((MAT)a)->col;
                         for ( i = 0; i < row; i++ )        MKMAT(mat,row,col);
                                 for ( j = 0; j < col; j++ ) {        for ( i = 0; i < row; i++ )
                                         arg0.body = (pointer)BDY((MAT)a)[i][j];          for ( j = 0; j < col; j++ ) {
                                         arg0.next = NEXT(arg);            arg0.body = (pointer)BDY((MAT)a)[i][j];
                                         Psubst(&arg0,&b);            arg0.next = NEXT(arg);
                                         BDY(mat)[i][j] = (pointer)b;            Psubst(&arg0,&b);
                                 }            BDY(mat)[i][j] = (pointer)b;
                         *rp = (Obj)mat;          }
                         break;        *rp = (Obj)mat;
                 case O_DP:        break;
                         mp0 = 0;      case O_DP:
                         for ( m = BDY((DP)a); m; m = NEXT(m) ) {        mp0 = 0;
                                 arg0.body = (pointer)C(m);        for ( m = BDY((DP)a); m; m = NEXT(m) ) {
                                 arg0.next = NEXT(arg);          arg0.body = (pointer)C(m);
                                 Psubst(&arg0,&b);          arg0.next = NEXT(arg);
                                 if ( b ) {          Psubst(&arg0,&b);
                                         NEXTMP(mp0,mp);          if ( b ) {
                                         C(mp) = (P)b;            NEXTMP(mp0,mp);
                                         mp->dl = m->dl;            C(mp) = (P)b;
                                 }            mp->dl = m->dl;
                         }          }
                         if ( mp0 ) {        }
                                 MKDP(NV((DP)a),mp0,d);        if ( mp0 ) {
                                 d->sugar = ((DP)a)->sugar;          MKDP(NV((DP)a),mp0,d);
                                 *rp = (Obj)d;          d->sugar = ((DP)a)->sugar;
                         } else          *rp = (Obj)d;
                                 *rp = 0;        } else
           *rp = 0;
   
                         break;        break;
                 default:      default:
                         error("subst : invalid argument");        error("subst : invalid argument");
         }    }
     if ( lastvl != LASTCO ) {
       get_vars_recursive(*rp,&vl);
       prev = lastvl; cur = NEXT(prev);
       while ( cur ) {
         v = cur->v;
         for ( tvl = vl; tvl && tvl->v != v; tvl = NEXT(tvl) );
         if ( !tvl ) NEXT(prev) = NEXT(cur);
         else prev = cur;
         cur = NEXT(cur);
       }
       update_LASTCO();
     }
 }  }
   
 FNODE subst_in_fnode();  FNODE subst_in_fnode();
Line 236  void Psubst_quote(arg,rp)
Line 277  void Psubst_quote(arg,rp)
 NODE arg;  NODE arg;
 QUOTE *rp;  QUOTE *rp;
 {  {
         QUOTE h;    QUOTE h;
         FNODE fn;    FNODE fn;
         Obj g;    Obj g;
         LIST l;    LIST l;
         V v;    V v;
   
         if ( !arg ) {    if ( !arg ) {
                 *rp = 0; return;      *rp = 0; return;
         }    }
         asir_assert(ARG0(arg),O_QUOTE,"subst_quote");    asir_assert(ARG0(arg),O_QUOTE,"subst_quote");
         fn = BDY((QUOTE)ARG0(arg)); arg = NEXT(arg);    fn = BDY((QUOTE)ARG0(arg)); arg = NEXT(arg);
         if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )    if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )
                 arg = BDY(l);      arg = BDY(l);
         while ( arg ) {    while ( arg ) {
                 asir_assert(BDY(arg),O_P,"subst_quote");      asir_assert(BDY(arg),O_P,"subst_quote");
                 v = VR((P)BDY(arg)); arg = NEXT(arg);      v = VR((P)BDY(arg)); arg = NEXT(arg);
                 if ( !arg )      if ( !arg )
                         error("subst_quote : invalid argument");        error("subst_quote : invalid argument");
                 g = (Obj)ARG0(arg); arg = NEXT(arg);      g = (Obj)ARG0(arg); arg = NEXT(arg);
                 if ( !g || OID(g) != O_QUOTE )      if ( !g || OID(g) != O_QUOTE )
                         objtoquote(g,&h);        objtoquote(g,&h);
                 else      else
                         h = (QUOTE)g;        h = (QUOTE)g;
                 fn = subst_in_fnode(fn,v,BDY(h));      fn = subst_in_fnode(fn,v,BDY(h));
         }    }
         MKQUOTE(*rp,fn);    MKQUOTE(*rp,fn);
 }  }
   
 void Ppsubst(arg,rp)  void Ppsubst(arg,rp)
 NODE arg;  NODE arg;
 Obj *rp;  Obj *rp;
 {  {
         Obj a,b,t;    Obj a,b,t;
         LIST l;    LIST l;
         V v;    V v;
   
         if ( !arg ) {    if ( !arg ) {
                 *rp = 0; return;      *rp = 0; return;
         }    }
         asir_assert(ARG0(arg),O_R,"psubst");    asir_assert(ARG0(arg),O_R,"psubst");
         reductr(CO,(Obj)ARG0(arg),&a);    reductr(CO,(Obj)ARG0(arg),&a);
 /*      a = (Obj)ARG0(arg); */  /*  a = (Obj)ARG0(arg); */
         arg = NEXT(arg);    arg = NEXT(arg);
         if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )    if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )
                 arg = BDY(l);      arg = BDY(l);
         while ( arg ) {    while ( arg ) {
                 asir_assert(BDY(arg),O_P,"psubst");      asir_assert(BDY(arg),O_P,"psubst");
                 v = VR((P)BDY(arg)); arg = NEXT(arg);      v = VR((P)BDY(arg)); arg = NEXT(arg);
                 if ( !arg )      if ( !arg )
                         error("psubst : invalid argument");        error("psubst : invalid argument");
                 asir_assert(ARG0(arg),O_R,"psubst");      asir_assert(ARG0(arg),O_R,"psubst");
                 reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);      reductr(CO,(Obj)BDY(arg),&b); arg = NEXT(arg);
 /*              b = (Obj)BDY(arg); arg = NEXT(arg); */  /*    b = (Obj)BDY(arg); arg = NEXT(arg); */
                 substr(CO,1,a,v,b,&t); a = t;      substr(CO,1,a,v,b,&t); a = t;
         }    }
         *rp = a;    *rp = a;
 }  }
   
 void Psubstf(arg,rp)  void Psubstf(arg,rp)
 NODE arg;  NODE arg;
 Obj *rp;  Obj *rp;
 {  {
         Obj a,t;    Obj a,t;
         LIST l;    LIST l;
         V v,f;    V v,f;
   
         if ( !arg ) {    if ( !arg ) {
                 *rp = 0; return;      *rp = 0; return;
         }    }
         asir_assert(ARG0(arg),O_R,"substf");    asir_assert(ARG0(arg),O_R,"substf");
         reductr(CO,(Obj)ARG0(arg),&a);    reductr(CO,(Obj)ARG0(arg),&a);
 /*      a = (Obj)ARG0(arg); */  /*  a = (Obj)ARG0(arg); */
         arg = NEXT(arg);    arg = NEXT(arg);
         if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )    if ( arg && (l = (LIST)ARG0(arg)) && OID(l) == O_LIST )
                 arg = BDY(l);      arg = BDY(l);
         while ( arg ) {    while ( arg ) {
                 asir_assert(BDY(arg),O_P,"substf");      asir_assert(BDY(arg),O_P,"substf");
                 v = VR((P)BDY(arg)); arg = NEXT(arg);      v = VR((P)BDY(arg)); arg = NEXT(arg);
                 if ( !arg || (int)v->attr != V_SR )      if ( !arg || (int)v->attr != V_SR )
                         error("substf : invalid argument");        error("substf : invalid argument");
                 f = VR((P)BDY(arg)); arg = NEXT(arg);      f = VR((P)BDY(arg)); arg = NEXT(arg);
                 if ( (int)f->attr != V_SR )      if ( (int)f->attr != V_SR )
                         error("substf : invalid argument\n");        error("substf : invalid argument\n");
                 substfr(CO,a,((FUNC)v->priv)->f.puref,((FUNC)f->priv)->f.puref,&t);      substfr(CO,a,((FUNC)v->priv)->f.puref,((FUNC)f->priv)->f.puref,&t);
                 a = t;      a = t;
         }    }
         *rp = a;    *rp = a;
 }  }
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.13
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>