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

Diff for /OpenXM_contrib2/asir2018/builtin/gr.c between version 1.1 and 1.2

version 1.1, 2018/09/19 05:45:06 version 1.2, 2018/09/28 08:20:27
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: OpenXM_contrib2/asir2018/builtin/gr.c,v 1.1 2018/09/19 05:45:06 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
Line 388  void dp_gr_main(LIST f,LIST v,Num homo,int modular,int
Line 388  void dp_gr_main(LIST f,LIST v,Num homo,int modular,int
     else      else
       dtop(CO,vv,ps[(long)BDY(x)],(Obj *)&BDY(r));        dtop(CO,vv,ps[(long)BDY(x)],(Obj *)&BDY(r));
     NEXTNODE(ind0,ind);      NEXTNODE(ind0,ind);
     STOQ((long)BDY(x),q); BDY(ind) = q;      STOZ((long)BDY(x),q); BDY(ind) = q;
   }    }
   if ( r0 ) NEXT(r) = 0;    if ( r0 ) NEXT(r) = 0;
   if ( ind0 ) NEXT(ind) = 0;    if ( ind0 ) NEXT(ind) = 0;
Line 448  void dp_interreduce(LIST f,LIST v,int field,struct ord
Line 448  void dp_interreduce(LIST f,LIST v,int field,struct ord
     else      else
       dtop(CO,vv,ps[(long)BDY(x)],(Obj *)&BDY(r));        dtop(CO,vv,ps[(long)BDY(x)],(Obj *)&BDY(r));
     NEXTNODE(ind0,ind);      NEXTNODE(ind0,ind);
     STOQ((long)BDY(x),q); BDY(ind) = q;      STOZ((long)BDY(x),q); BDY(ind) = q;
   }    }
   if ( r0 ) NEXT(r) = 0;    if ( r0 ) NEXT(r) = 0;
   if ( ind0 ) NEXT(ind) = 0;    if ( ind0 ) NEXT(ind) = 0;
Line 1211  void makesubst(VL v,NODE *s)
Line 1211  void makesubst(VL v,NODE *s)
   for ( r0 = 0; v; v = NEXT(v) ) {    for ( r0 = 0; v; v = NEXT(v) ) {
     NEXTNODE(r0,r); BDY(r) = (pointer)v->v;      NEXTNODE(r0,r); BDY(r) = (pointer)v->v;
 #if defined(_PA_RISC1_1)  #if defined(_PA_RISC1_1)
     n = mrand48()&BMASK; UTOQ(n,q);      n = mrand48()&BMASK; UTOZ(n,q);
 #else  #else
     n = random(); UTOQ(n,q);      n = random(); UTOZ(n,q);
 #endif  #endif
     NEXTNODE(r0,r); BDY(r) = (pointer)q;      NEXTNODE(r0,r); BDY(r) = (pointer)q;
   }    }
Line 1226  void printsubst(NODE s)
Line 1226  void printsubst(NODE s)
   fputc('[',asir_out);    fputc('[',asir_out);
   while ( s ) {    while ( s ) {
     printv(CO,(V)BDY(s)); s = NEXT(s);      printv(CO,(V)BDY(s)); s = NEXT(s);
     fprintf(asir_out,"->%d",QTOS((Q)BDY(s)));      fprintf(asir_out,"->%d",ZTOS((Q)BDY(s)));
     if ( NEXT(s) ) {      if ( NEXT(s) ) {
       fputc(',',asir_out); s = NEXT(s);        fputc(',',asir_out); s = NEXT(s);
     } else      } else
Line 1318  void setup_arrays(NODE f,int m,NODE *r)
Line 1318  void setup_arrays(NODE f,int m,NODE *r)
     MKLIST(input,t);      MKLIST(input,t);
   
     if ( OXCheck >= 0 ) {      if ( OXCheck >= 0 ) {
       STOQ(OXCheck,q);        STOZ(OXCheck,q);
       MKSTR(fname,"register_input");        MKSTR(fname,"register_input");
       arg = mknode(3,q,fname,input);        arg = mknode(3,q,fname,input);
       Pox_cmo_rpc(arg,&obj);        Pox_cmo_rpc(arg,&obj);
Line 1446  void reduceall(NODE in,NODE *h)
Line 1446  void reduceall(NODE in,NODE *h)
       NODE node;        NODE node;
       LIST hist;        LIST hist;
   
       STOQ(w[i],q);        STOZ(w[i],q);
       node = mknode(4,ONE,q,ONE,ONE);        node = mknode(4,ONE,q,ONE,ONE);
       MKLIST(hist,node);        MKLIST(hist,node);
       MKNODE(TraceList,hist,0);        MKNODE(TraceList,hist,0);
Line 1566  int newps(DP a,int m,NODE subst)
Line 1566  int newps(DP a,int m,NODE subst)
     }      }
     MKLIST(trace,tr);      MKLIST(trace,tr);
     if ( OXCheck >= 0 ) {      if ( OXCheck >= 0 ) {
       STOQ(OXCheck,q1);        STOZ(OXCheck,q1);
       MKSTR(fname,"check_trace");        MKSTR(fname,"check_trace");
       STOQ(psn,q2);        STOZ(psn,q2);
       arg = mknode(5,q1,fname,a,q2,trace);        arg = mknode(5,q1,fname,a,q2,trace);
       Pox_cmo_rpc(arg,&obj);        Pox_cmo_rpc(arg,&obj);
     } else if ( OXCheck < 0 ) {      } else if ( OXCheck < 0 ) {
       STOQ(psn,q1);        STOZ(psn,q1);
       tn = mknode(2,q1,trace);        tn = mknode(2,q1,trace);
       MKLIST(trace1,tn);        MKLIST(trace1,tn);
       MKNODE(tr,trace1,AllTraceList);        MKNODE(tr,trace1,AllTraceList);
Line 1657  void reducebase_dehomo(NODE f,NODE *g)
Line 1657  void reducebase_dehomo(NODE f,NODE *g)
         LIST hist;          LIST hist;
         NODE node;          NODE node;
   
         STOQ(r[i],q);          STOZ(r[i],q);
         node = mknode(4,NULLP,q,NULLP,NULLP);          node = mknode(4,NULLP,q,NULLP,NULLP);
         MKLIST(hist,node);          MKLIST(hist,node);
         MKNODE(TraceList,hist,0);          MKNODE(TraceList,hist,0);
Line 1789  NODE gb(NODE f,int m,NODE subst)
Line 1789  NODE gb(NODE f,int m,NODE subst)
       } else        } else
         dp_sp(ps[l->dp1],ps[l->dp2],&h);          dp_sp(ps[l->dp1],ps[l->dp2],&h);
       if ( GenTrace ) {        if ( GenTrace ) {
         STOQ(l->dp1,q); ARG1(BDY((LIST)BDY(NEXT(TraceList)))) = q;          STOZ(l->dp1,q); ARG1(BDY((LIST)BDY(NEXT(TraceList)))) = q;
         STOQ(l->dp2,q); ARG1(BDY((LIST)BDY(TraceList))) = q;          STOZ(l->dp2,q); ARG1(BDY((LIST)BDY(TraceList))) = q;
       }        }
       if ( h )        if ( h )
         new_sugar = h->sugar;          new_sugar = h->sugar;
Line 2222  void gbcheck_list(NODE f,int n,VECT *gp,LIST *pp)
Line 2222  void gbcheck_list(NODE f,int n,VECT *gp,LIST *pp)
   
   for ( u0 = 0, l = d; l; l = NEXT(l) ) {    for ( u0 = 0, l = d; l; l = NEXT(l) ) {
     NEXTNODE(u0,u);      NEXTNODE(u0,u);
     STOQ(l->dp1,q1);      STOZ(l->dp1,q1);
     STOQ(l->dp2,q2);      STOZ(l->dp2,q2);
     t = mknode(2,q1,q2);      t = mknode(2,q1,q2);
     MKLIST(pair,t);      MKLIST(pair,t);
     BDY(u) = (pointer)pair;      BDY(u) = (pointer)pair;
Line 2280  void dp_set_flag(Obj name,Obj value)
Line 2280  void dp_set_flag(Obj name,Obj value)
     ratio = (Q)value;      ratio = (Q)value;
     if ( ratio ) {      if ( ratio ) {
       nmq(ratio,&t);        nmq(ratio,&t);
       DP_Multiple = QTOS(t);        DP_Multiple = ZTOS(t);
       dnq(ratio,&t);        dnq(ratio,&t);
       Denominator = QTOS(t);        Denominator = ZTOS(t);
     } else {      } else {
       DP_Multiple = 0;        DP_Multiple = 0;
       Denominator = 1;        Denominator = 1;
Line 2290  void dp_set_flag(Obj name,Obj value)
Line 2290  void dp_set_flag(Obj name,Obj value)
   }    }
   if ( value && OID(value) != O_N )    if ( value && OID(value) != O_N )
     return;      return;
   v = QTOS((Q)value);    v = ZTOS((Q)value);
   if ( !strcmp(n,"NoSugar") )    if ( !strcmp(n,"NoSugar") )
     NoSugar = v;      NoSugar = v;
   else if ( !strcmp(n,"NoCriB") )    else if ( !strcmp(n,"NoCriB") )
Line 2343  void dp_make_flaglist(LIST *list)
Line 2343  void dp_make_flaglist(LIST *list)
   NODE n,n1;    NODE n,n1;
   
 #if 0  #if 0
   STOQ(DP_Multiple,v); MKNODE(n,v,0); MKSTR(name,"DP_Multiple"); MKNODE(n1,name,n); n = n1;    STOZ(DP_Multiple,v); MKNODE(n,v,0); MKSTR(name,"DP_Multiple"); MKNODE(n1,name,n); n = n1;
   STOQ(Denominator,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Denominator"); MKNODE(n1,name,n); n = n1;    STOZ(Denominator,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Denominator"); MKNODE(n1,name,n); n = n1;
 #else  #else
   if ( DP_Multiple ) {    if ( DP_Multiple ) {
     STOQ(DP_Multiple,nm); STOQ(Denominator,dn); divq((Q)nm,(Q)dn,&r);      STOZ(DP_Multiple,nm); STOZ(Denominator,dn); divq((Q)nm,(Q)dn,&r);
   } else    } else
     v = 0;      v = 0;
   MKNODE(n,v,0); MKSTR(name,"Content"); MKNODE(n1,name,n); n = n1;    MKNODE(n,v,0); MKSTR(name,"Content"); MKNODE(n1,name,n); n = n1;
 #endif  #endif
   MKNODE(n1,Dist,n); n = n1; MKSTR(name,"Dist"); MKNODE(n1,name,n); n = n1;    MKNODE(n1,Dist,n); n = n1; MKSTR(name,"Dist"); MKNODE(n1,name,n); n = n1;
   STOQ(Reverse,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Reverse"); MKNODE(n1,name,n); n = n1;    STOZ(Reverse,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Reverse"); MKNODE(n1,name,n); n = n1;
   STOQ(Stat,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Stat"); MKNODE(n1,name,n); n = n1;    STOZ(Stat,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Stat"); MKNODE(n1,name,n); n = n1;
   STOQ(DP_Print,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Print"); MKNODE(n1,name,n); n = n1;    STOZ(DP_Print,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Print"); MKNODE(n1,name,n); n = n1;
   STOQ(DP_PrintShort,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"PrintShort"); MKNODE(n1,name,n); n = n1;    STOZ(DP_PrintShort,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"PrintShort"); MKNODE(n1,name,n); n = n1;
   STOQ(DP_NFStat,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NFStat"); MKNODE(n1,name,n); n = n1;    STOZ(DP_NFStat,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NFStat"); MKNODE(n1,name,n); n = n1;
   STOQ(OXCheck,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"OXCheck"); MKNODE(n1,name,n); n = n1;    STOZ(OXCheck,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"OXCheck"); MKNODE(n1,name,n); n = n1;
   STOQ(GenTrace,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"GenTrace"); MKNODE(n1,name,n); n = n1;    STOZ(GenTrace,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"GenTrace"); MKNODE(n1,name,n); n = n1;
   STOQ(GenSyz,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"GenSyz"); MKNODE(n1,name,n); n = n1;    STOZ(GenSyz,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"GenSyz"); MKNODE(n1,name,n); n = n1;
   STOQ(MaxDeg,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"MaxDeg"); MKNODE(n1,name,n); n = n1;    STOZ(MaxDeg,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"MaxDeg"); MKNODE(n1,name,n); n = n1;
   STOQ(OneZeroHomo,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"OneZeroHomo"); MKNODE(n1,name,n); n = n1;    STOZ(OneZeroHomo,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"OneZeroHomo"); MKNODE(n1,name,n); n = n1;
   STOQ(PtozpRA,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"PtozpRA"); MKNODE(n1,name,n); n = n1;    STOZ(PtozpRA,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"PtozpRA"); MKNODE(n1,name,n); n = n1;
   STOQ(ShowMag,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"ShowMag"); MKNODE(n1,name,n); n = n1;    STOZ(ShowMag,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"ShowMag"); MKNODE(n1,name,n); n = n1;
   STOQ(Top,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Top"); MKNODE(n1,name,n); n = n1;    STOZ(Top,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Top"); MKNODE(n1,name,n); n = n1;
   STOQ(NoGCD,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoGCD"); MKNODE(n1,name,n); n = n1;    STOZ(NoGCD,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoGCD"); MKNODE(n1,name,n); n = n1;
   STOQ(NoRA,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoRA"); MKNODE(n1,name,n); n = n1;    STOZ(NoRA,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoRA"); MKNODE(n1,name,n); n = n1;
   STOQ(NoMC,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoMC"); MKNODE(n1,name,n); n = n1;    STOZ(NoMC,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoMC"); MKNODE(n1,name,n); n = n1;
   STOQ(NoGC,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoGC"); MKNODE(n1,name,n); n = n1;    STOZ(NoGC,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoGC"); MKNODE(n1,name,n); n = n1;
   STOQ(NoCriB,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoCriB"); MKNODE(n1,name,n); n = n1;    STOZ(NoCriB,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoCriB"); MKNODE(n1,name,n); n = n1;
   STOQ(NoSugar,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoSugar"); MKNODE(n1,name,n); n = n1;    STOZ(NoSugar,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoSugar"); MKNODE(n1,name,n); n = n1;
   if ( Demand )    if ( Demand )
     MKSTR(path,Demand);      MKSTR(path,Demand);
   else    else
Line 2497  void _dp_nf(NODE b,DP g,DP *ps,int full,DP *rp)
Line 2497  void _dp_nf(NODE b,DP g,DP *ps,int full,DP *rp)
           Z cq;            Z cq;
           NODE node,node0;            NODE node,node0;
   
           STOQ((long)BDY(l),cq);            STOZ((long)BDY(l),cq);
           node0 = mknode(4,coef,cq,mult,ONE);            node0 = mknode(4,coef,cq,mult,ONE);
           MKLIST(hist,node0);            MKLIST(hist,node0);
           MKNODE(node,hist,TraceList); TraceList = node;            MKNODE(node,hist,TraceList); TraceList = node;
Line 2579  void _dp_nf_z(NODE b,DP g,DP *ps,int full,int multiple
Line 2579  void _dp_nf_z(NODE b,DP g,DP *ps,int full,int multiple
   
         if ( GenTrace ) {          if ( GenTrace ) {
           /* u = cr*rp + (-cred)*shift*red */            /* u = cr*rp + (-cred)*shift*red */
           STOQ((long)BDY(l),cq);            STOZ((long)BDY(l),cq);
           node = mknode(4,cr,cq,NULLP,NULLP);            node = mknode(4,cr,cq,NULLP,NULLP);
           mulz(cred,rc,&rcred);            mulz(cred,rc,&rcred);
           chsgnz(rcred,&mrcred);            chsgnz(rcred,&mrcred);

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>