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

Diff for /OpenXM_contrib2/asir2018/engine/nd.c between version 1.42 and 1.45

version 1.42, 2020/12/03 07:58:54 version 1.45, 2021/01/11 08:37:44
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.41 2020/11/26 03:55:23 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.44 2020/12/15 07:40:09 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
   void print_siglist(NODE l);
   
 int Nnd_add,Nf4_red,NcriB,NcriMF,Ncri2,Npairs;  int Nnd_add,Nf4_red,NcriB,NcriMF,Ncri2,Npairs;
 struct oEGT eg_search,f4_symb,f4_conv,f4_elim1,f4_elim2;  struct oEGT eg_search,f4_symb,f4_conv,f4_elim1,f4_elim2;
   
Line 67  static NODE nd_tracelist;
Line 69  static NODE nd_tracelist;
 static NODE nd_alltracelist;  static NODE nd_alltracelist;
 static int nd_gentrace,nd_gensyz,nd_nora,nd_newelim,nd_intersect,nd_lf,nd_norb;  static int nd_gentrace,nd_gensyz,nd_nora,nd_newelim,nd_intersect,nd_lf,nd_norb;
 static int nd_f4_td,nd_sba_f4step,nd_sba_pot,nd_sba_largelcm,nd_sba_dontsort,nd_sba_redundant_check;  static int nd_f4_td,nd_sba_f4step,nd_sba_pot,nd_sba_largelcm,nd_sba_dontsort,nd_sba_redundant_check;
 static int nd_top;  static int nd_top,nd_sba_syz,nd_sba_fixord,nd_sba_grevlexgb;
 static int *nd_gbblock;  static int *nd_gbblock;
 static NODE nd_nzlist,nd_check_splist;  static NODE nd_nzlist,nd_check_splist;
 static int nd_splist;  static int nd_splist;
Line 98  void dltondl(int n,DL dl,UINT *r);
Line 100  void dltondl(int n,DL dl,UINT *r);
 DP ndvtodp(int mod,NDV p);  DP ndvtodp(int mod,NDV p);
 DP ndtodp(int mod,ND p);  DP ndtodp(int mod,ND p);
 DPM ndvtodpm(int mod,NDV p);  DPM ndvtodpm(int mod,NDV p);
   NDV dptondv(int mod,DP p);
 NDV dpmtondv(int mod,DPM p);  NDV dpmtondv(int mod,DPM p);
   int dp_getdeg(DP p);
 int dpm_getdeg(DPM p,int *rank);  int dpm_getdeg(DPM p,int *rank);
 void dpm_ptozp(DPM p,Z *cont,DPM *r);  void dpm_ptozp(DPM p,Z *cont,DPM *r);
 int compdmm(int nv,DMM a,DMM b);  int compdmm(int nv,DMM a,DMM b);
   DPM sigtodpm(SIG s);
   SIG dup_sig(SIG sig);
   
 void Pdp_set_weight(NODE,VECT *);  void Pdp_set_weight(NODE,VECT *);
 void Pox_cmo_rpc(NODE,Obj *);  void Pox_cmo_rpc(NODE,Obj *);
Line 794  int _eqdl(int n,DL d1,DL d2);
Line 800  int _eqdl(int n,DL d1,DL d2);
   
 int ndl_module_schreyer_compare(UINT *m1,UINT *m2)  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
 {  {
   int pos1,pos2,t,j;    int pos1,pos2,t,j,retpot;
   DMM *in;    DMM *in;
   DMMstack s;    DMMstack s;
   static DL d1=0;    static DL d1=0;
   static DL d2=0;    static DL d2=0;
   static int dlen=0;    static int dlen=0;
     extern int ReversePOT;
   
     if ( ReversePOT ) retpot = -1;
     else retpot = 1;
   pos1 = MPOS(m1); pos2 = MPOS(m2);    pos1 = MPOS(m1); pos2 = MPOS(m2);
   if ( pos1 == pos2 ) return (*ndl_base_compare_function)(m1,m2);    if ( pos1 == pos2 ) return (*ndl_base_compare_function)(m1,m2);
   if ( nd_nvar > dlen ) {    if ( nd_nvar > dlen ) {
Line 817  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
Line 826  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
     _addtodl(nd_nvar,in[pos1]->dl,d1);      _addtodl(nd_nvar,in[pos1]->dl,d1);
     _addtodl(nd_nvar,in[pos2]->dl,d2);      _addtodl(nd_nvar,in[pos2]->dl,d2);
     if ( in[pos1]->pos == in[pos2]->pos && _eqdl(nd_nvar,d1,d2)) {      if ( in[pos1]->pos == in[pos2]->pos && _eqdl(nd_nvar,d1,d2)) {
       if ( pos1 < pos2 ) return 1;        if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return 0;        else return 0;
     }      }
     pos1 = in[pos1]->pos;      pos1 = in[pos1]->pos;
Line 831  LAST:
Line 840  LAST:
     case 0:      case 0:
       t = (*dl_base_compare_function)(nd_nvar,d1,d2);        t = (*dl_base_compare_function)(nd_nvar,d1,d2);
       if ( t ) return t;        if ( t ) return t;
       else if ( pos1 < pos2 ) return 1;        else if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return 0;        else return 0;
       break;        break;
     case 1:      case 1:
       if ( pos1 < pos2 ) return 1;        if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return (*dl_base_compare_function)(nd_nvar,d1,d2);        else return (*dl_base_compare_function)(nd_nvar,d1,d2);
       break;        break;
     case 2:      case 2:
       if ( d1->td > d2->td  ) return 1;        if ( d1->td > d2->td  ) return 1;
       else if ( d1->td < d2->td ) return -1;        else if ( d1->td < d2->td ) return -1;
       else if ( pos1 < pos2 ) return 1;        else if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return (*dl_base_compare_function)(nd_nvar,d1,d2);        else return (*dl_base_compare_function)(nd_nvar,d1,d2);
       break;        break;
     default:      default:
Line 1238  void print_sig(SIG s)
Line 1247  void print_sig(SIG s)
   fprintf(asir_out,">>*e%d",s->pos);    fprintf(asir_out,">>*e%d",s->pos);
 }  }
   
   void print_siglist(NODE l)
   {
     for ( ; l; l = NEXT(l) )
      print_sig((SIG)l->body);
   }
   
   
 // assuming increasing order wrt signature  // assuming increasing order wrt signature
   
 INLINE int ndl_find_reducer_s(UINT *dg,SIG sig)  INLINE int ndl_find_reducer_s(UINT *dg,SIG sig)
Line 2589  NODE insert_sig(NODE l,SIG s)
Line 2605  NODE insert_sig(NODE l,SIG s)
       else if ( _dl_redble(sig,DL(t),nd_nvar) )        else if ( _dl_redble(sig,DL(t),nd_nvar) )
         // remove p          // remove p
         prev->next = p->next;          prev->next = p->next;
         else
           prev = p;
     } else      } else
       prev = p;        prev = p;
   }    }
Line 2629  int small_lcm(ND_pairs l)
Line 2647  int small_lcm(ND_pairs l)
   int i;    int i;
   NODE t;    NODE t;
   static DL lcm,mul,quo;    static DL lcm,mul,quo;
   static int nvar;    static int nvar = 0;
   
   if ( nd_sba_largelcm ) return 0;    if ( nd_sba_largelcm ) return 0;
   if ( nvar < nd_nvar ) {    if ( nvar < nd_nvar ) {
Line 2672  ND_pairs find_smallest_lcm(ND_pairs l)
Line 2690  ND_pairs find_smallest_lcm(ND_pairs l)
   ND_pairs r;    ND_pairs r;
   struct oSIG sig1;    struct oSIG sig1;
   static DL mul,quo,minlm;    static DL mul,quo,minlm;
   static int nvar;    static int nvar = 0;
   
   if ( nvar < nd_nvar ) {    if ( nvar < nd_nvar ) {
     nvar = nd_nvar;      nvar = nd_nvar;
Line 2747  struct oEGT eg_create,eg_newpairs,eg_merge;
Line 2765  struct oEGT eg_create,eg_newpairs,eg_merge;
   
 NODE conv_ilist_s(int demand,int trace,int **indp);  NODE conv_ilist_s(int demand,int trace,int **indp);
   
 NODE nd_sba_buch(int m,int ishomo,int **indp)  // S(fj*ei-fi*ej)
   
   void _subdl(int,DL,DL,DL);
   
   SIG trivial_sig(int i,int j)
 {  {
     static DL lcm;
     static struct oSIG sigi,sigj;
     static int nvar = 0;
     SIG sig;
   
     if ( nvar != nd_nvar ) {
       nvar = nd_nvar; NEWDL(lcm,nvar); NEWDL(sigi.dl,nvar); NEWDL(sigj.dl,nvar);
     }
     if ( nd_sba_grevlexgb != 0 ) {
       lcm_of_DL(nd_nvar,nd_sba_hm[i],nd_sba_hm[j],lcm);
       sigi.pos = i; _subdl(nd_nvar,lcm,nd_sba_hm[i],sigi.dl);
       sigj.pos = j; _subdl(nd_nvar,lcm,nd_sba_hm[j],sigj.dl);
       if ( comp_sig(&sigi,&sigj) > 0 ) sig = dup_sig(&sigi);
       else sig = dup_sig(&sigj);
     } else {
       sigi.pos = i; _copydl(nd_nvar,nd_sba_hm[j],sigi.dl);
       sigj.pos = j; _copydl(nd_nvar,nd_sba_hm[i],sigj.dl);
       if ( comp_sig(&sigi,&sigj) > 0 ) sig = dup_sig(&sigi);
       else sig = dup_sig(&sigj);
     }
     return sig;
   }
   
   NODE nd_sba_buch(int m,int ishomo,int **indp,NODE *syzp)
   {
   int i,j,nh,sugar,stat,pos;    int i,j,nh,sugar,stat,pos;
   NODE r,t,g;    NODE r,t,g;
   ND_pairs d;    ND_pairs d;
Line 2761  NODE nd_sba_buch(int m,int ishomo,int **indp)
Line 2808  NODE nd_sba_buch(int m,int ishomo,int **indp)
   LIST list;    LIST list;
   SIG sig;    SIG sig;
   NODE *syzlist;    NODE *syzlist;
     int ngen;
   int Nnominimal,Nredundant;    int Nnominimal,Nredundant;
   DL lcm,quo,mul;    DL lcm,quo,mul;
   struct oEGT eg1,eg2,eg_update,eg_remove,eg_large,eg_nf,eg_nfzero;    struct oEGT eg1,eg2,eg_update,eg_remove,eg_large,eg_nf,eg_nfzero;
Line 2773  init_eg(&eg_remove);
Line 2821  init_eg(&eg_remove);
   Nnominimal = 0;    Nnominimal = 0;
   Nredundant = 0;    Nredundant = 0;
   d = 0;    d = 0;
     ngen = nd_psn;
   for ( i = 0; i < nd_psn; i++ )    for ( i = 0; i < nd_psn; i++ )
     for ( j = i+1; j < nd_psn; j++ ) {      for ( j = i+1; j < nd_psn; j++ ) {
       NEWSIG(sig); sig->pos = j;        sig = trivial_sig(i,j);
       _copydl(nd_nvar,nd_sba_hm[i],sig->dl);  
       syzlist[sig->pos] = insert_sig(syzlist[sig->pos],sig);        syzlist[sig->pos] = insert_sig(syzlist[sig->pos],sig);
     }      }
   for ( i = 0; i < nd_psn; i++ ) {    for ( i = 0; i < nd_psn; i++ ) {
Line 2907  get_eg(&eg2); add_eg(&eg_remove,&eg1,&eg2);
Line 2955  get_eg(&eg2); add_eg(&eg_remove,&eg1,&eg2);
    print_eg("nfzero",&eg_nfzero);     print_eg("nfzero",&eg_nfzero);
    printf("\n");     printf("\n");
  }   }
    if ( nd_sba_syz ) {
      NODE hsyz,tsyz,prev;
   
      hsyz = 0;
      for ( i = 0; i < ngen; i++ ) {
        tsyz = syzlist[i];
        for ( prev = 0; tsyz != 0; prev = tsyz, tsyz = NEXT(tsyz))
          BDY(tsyz) = (pointer)sigtodpm((SIG)BDY(tsyz));
        if ( prev != 0 ) {
          prev->next = hsyz; hsyz = syzlist[i];
        }
      }
      *syzp = hsyz;
    } else *syzp = 0;
  return g;   return g;
 }  }
   
Line 3398  int comp_sig(SIG s1,SIG s2)
Line 3460  int comp_sig(SIG s1,SIG s2)
     else return (*cmpdl)(nd_nvar,s1->dl,s2->dl);      else return (*cmpdl)(nd_nvar,s1->dl,s2->dl);
   } else {    } else {
     static DL m1,m2;      static DL m1,m2;
     static int nvar;      static int nvar = 0;
     int ret;      int ret;
   
     if ( nvar != nd_nvar ) {      if ( nvar != nd_nvar ) {
Line 3406  int comp_sig(SIG s1,SIG s2)
Line 3468  int comp_sig(SIG s1,SIG s2)
     }      }
     _adddl(nd_nvar,s1->dl,nd_sba_hm[s1->pos],m1);      _adddl(nd_nvar,s1->dl,nd_sba_hm[s1->pos],m1);
     _adddl(nd_nvar,s2->dl,nd_sba_hm[s2->pos],m2);      _adddl(nd_nvar,s2->dl,nd_sba_hm[s2->pos],m2);
     ret = (*cmpdl)(nd_nvar,m1,m2);      if ( nd_sba_fixord )
         ret = cmpdl_revgradlex(nd_nvar,m1,m2);
       else
         ret = (*cmpdl)(nd_nvar,m1,m2);
     if ( ret != 0 ) return ret;      if ( ret != 0 ) return ret;
     else if ( s1->pos > s2->pos ) return 1;      else if ( s1->pos > s2->pos ) return 1;
     else if ( s1->pos < s2->pos ) return -1;      else if ( s1->pos < s2->pos ) return -1;
Line 3569  ND_pairs nd_newpairs_s(int t, NODE *syz)
Line 3634  ND_pairs nd_newpairs_s(int t, NODE *syz)
   int ts,ret,i;    int ts,ret,i;
   ND_pairs r,r0,_sp,sp;    ND_pairs r,r0,_sp,sp;
   SIG spsig,tsig;    SIG spsig,tsig;
   static int nvar;    static int nvar = 0;
   static SIG _sig1,_sig2;    static SIG _sig1,_sig2;
   struct oEGT eg1,eg2,eg3,eg4;    struct oEGT eg1,eg2,eg3,eg4;
   
Line 3968  int ndv_newps(int m,NDV a,NDV aq)
Line 4033  int ndv_newps(int m,NDV a,NDV aq)
     return nd_psn++;      return nd_psn++;
 }  }
   
   // find LM wrt grevlex
   void ndv_lm_fixord(NDV p,DL d)
   {
     NMV m;
     DL tmp;
     int len,i,ret;
   
     NEWDL(tmp,nd_nvar);
     m = BDY(p); len = LEN(p);
     _ndltodl(DL(m),d); printdl(d); printf("->");
     for ( i = 1, NMV_ADV(m); i < len; i++, NMV_ADV(m) ) {
       _ndltodl(DL(m),tmp);
       ret = cmpdl_revgradlex(nd_nvar,tmp,d);
       if ( ret > 0 ) _copydl(nd_nvar,tmp,d);
     }
     printdl(d); printf("\n");
   }
   
 /* nd_tracelist = [[0,index,div],...,[nd_psn-1,index,div]] */  /* nd_tracelist = [[0,index,div],...,[nd_psn-1,index,div]] */
 /* return 1 if success, 0 if failure (HC(a mod p)) */  /* return 1 if success, 0 if failure (HC(a mod p)) */
   
Line 4086  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
Line 4169  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
         if ( nd_demand ) nd_ps_trace_sym[i]->sig = sig;          if ( nd_demand ) nd_ps_trace_sym[i]->sig = sig;
       }        }
       NEWDL(nd_sba_hm[i],nd_nvar);        NEWDL(nd_sba_hm[i],nd_nvar);
       _ndltodl(DL(nd_psh[i]),nd_sba_hm[i]);        if ( nd_sba_fixord )
           ndv_lm_fixord(nd_ps[i],nd_sba_hm[i]);
         else
           _ndltodl(DL(nd_psh[i]),nd_sba_hm[i]);
     }      }
     nd_sba_pos = (NODE *)MALLOC(nd_psn*sizeof(NODE));      nd_sba_pos = (NODE *)MALLOC(nd_psn*sizeof(NODE));
     for ( i = 0; i < nd_psn; i++ ) {      for ( i = 0; i < nd_psn; i++ ) {
Line 4285  void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int 
Line 4371  void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int 
                   }                    }
                 }                  }
             } else {              } else {
                 e = getdeg(tv->v,(P)BDY(t));                  if ( OID(BDY(t)) == O_DP ) {
                 max = MAX(e,max);                    e = dp_getdeg((DP)BDY(t));
                     max = MAX(e,max);
                   } else {
                     e = getdeg(tv->v,(P)BDY(t));
                     max = MAX(e,max);
                   }
             }              }
         }          }
     nd_setup_parameters(nvar,nd_nzlist?0:max);      nd_setup_parameters(nvar,nd_nzlist?0:max);
Line 4307  void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int 
Line 4398  void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int 
           b = (pointer)pltondv(CO,vv,zpl);            b = (pointer)pltondv(CO,vv,zpl);
         }          }
       } else {        } else {
         if ( !m && !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);          if ( OID(BDY(t)) == O_DP ) {
         else zp = (P)BDY(t);            DP zdp;
         b = (pointer)ptondv(CO,vv,zp);  
             if ( !m && !nd_gentrace ) dp_ptozp((DP)BDY(t),&zdp);
             else zdp = (DP)BDY(t);
             b = (pointer)dptondv(m,zdp);
           } else {
             if ( !m && !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);
             else zp = (P)BDY(t);
             b = (pointer)ptondv(CO,vv,zp);
           }
       }        }
       if ( ishomo )        if ( ishomo )
         ishomo = ishomo && ndv_ishomo(b);          ishomo = ishomo && ndv_ishomo(b);
Line 4439  NODE nd_sba_f4(int m,int **indp);
Line 4538  NODE nd_sba_f4(int m,int **indp);
 void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int f4,struct order_spec *ord,LIST *rp)  void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int f4,struct order_spec *ord,LIST *rp)
 {  {
   VL tv,fv,vv,vc,av;    VL tv,fv,vv,vc,av;
   NODE fd,fd0,r,r0,t,x,s,xx;    NODE fd,fd0,r,r0,t,x,s,xx,nd,syz;
   int e,max,nvar,i;    int e,max,nvar,i;
   NDV b;    NDV b;
   int ishomo,nalg,wmax,len;    int ishomo,nalg,wmax,len;
Line 4477  void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int
Line 4576  void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int
   initd(ord);    initd(ord);
   for ( t = BDY(f), max = 1; t; t = NEXT(t) ) {    for ( t = BDY(f), max = 1; t; t = NEXT(t) ) {
     for ( tv = vv; tv; tv = NEXT(tv) ) {      for ( tv = vv; tv; tv = NEXT(tv) ) {
       e = getdeg(tv->v,(P)BDY(t));        if ( OID(BDY(t)) == O_DP ) {
       max = MAX(e,max);          e = dp_getdeg((DP)BDY(t));
           max = MAX(e,max);
         } else {
           e = getdeg(tv->v,(P)BDY(t));
           max = MAX(e,max);
         }
     }      }
   }    }
   nd_setup_parameters(nvar,max);    nd_setup_parameters(nvar,max);
   obpe = nd_bpe; oadv = nmv_adv; oepos = nd_epos; ompos = nd_mpos;    obpe = nd_bpe; oadv = nmv_adv; oepos = nd_epos; ompos = nd_mpos;
   ishomo = 1;    ishomo = 1;
   for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {    for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {
     if ( !m ) ptozp((P)BDY(t),1,&dmy,&zp);      if ( OID(BDY(t)) == O_DP ) {
     else zp = (P)BDY(t);        DP zdp;
     b = (pointer)ptondv(CO,vv,zp);  
         if ( !m ) dp_ptozp((DP)BDY(t),&zdp);
         else zdp = (DP)BDY(t);
         b = (pointer)dptondv(m,zdp);
       } else {
         if ( !m ) ptozp((P)BDY(t),1,&dmy,&zp);
         else zp = (P)BDY(t);
         b = (pointer)ptondv(CO,vv,zp);
       }
     if ( ishomo )      if ( ishomo )
       ishomo = ishomo && ndv_ishomo(b);        ishomo = ishomo && ndv_ishomo(b);
     if ( m ) ndv_mod(m,b);      if ( m ) ndv_mod(m,b);
Line 4511  void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int
Line 4623  void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int
   }    }
   
   ndv_setup(m,0,fd0,nd_sba_dontsort,0,1);    ndv_setup(m,0,fd0,nd_sba_dontsort,0,1);
   x = f4 ? nd_sba_f4(m,&perm) : nd_sba_buch(m,ishomo || homo,&perm);    x = f4 ? nd_sba_f4(m,&perm) : nd_sba_buch(m,ishomo || homo,&perm,&syz);
   if ( !x ) {    if ( !x ) {
     *rp = 0; return;      *rp = 0; return;
   }    }
Line 4534  void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int
Line 4646  void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int
     else BDY(r) = ndvtop(m,CO,vv,BDY(t));      else BDY(r) = ndvtop(m,CO,vv,BDY(t));
   }    }
   if ( r0 ) NEXT(r) = 0;    if ( r0 ) NEXT(r) = 0;
   MKLIST(*rp,r0);    if ( nd_sba_syz ) {
       LIST gb,hsyz;
       NODE nd;
   
       MKLIST(gb,r0);
       MKLIST(hsyz,syz);
       nd = mknode(2,gb,hsyz);
       MKLIST(*rp,nd);
     } else
       MKLIST(*rp,r0);
   get_eg(&eg1); init_eg(&egconv); add_eg(&egconv,&eg0,&eg1);    get_eg(&eg1); init_eg(&egconv); add_eg(&egconv,&eg0,&eg1);
   print_eg("conv",&egconv); fprintf(asir_out,"\n");    print_eg("conv",&egconv); fprintf(asir_out,"\n");
 }  }
Line 4868  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 4989  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
                 }                  }
               }                }
             } else {              } else {
                 if ( OID(BDY(t)) == O_DP ) {
                   e = dp_getdeg((DP)BDY(t));
                   max = MAX(e,max);
                 } else {
                 e = getdeg(tv->v,(P)BDY(t));                  e = getdeg(tv->v,(P)BDY(t));
                 max = MAX(e,max);                  max = MAX(e,max);
                 }
             }              }
         }          }
     nd_setup_parameters(nvar,max);      nd_setup_parameters(nvar,max);
Line 4890  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5016  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
             c = (pointer)pltondv(CO,vv,zpl);              c = (pointer)pltondv(CO,vv,zpl);
           }            }
         } else {          } else {
           if ( !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);            if ( OID(BDY(t)) == O_DP ) {
           else zp = (P)BDY(t);              DP zdp;
           c = (pointer)ptondv(CO,vv,zp);  
               if ( !nd_gentrace ) dp_ptozp((DP)BDY(t),&zdp);
               else zdp = (DP)BDY(t);
               c = (pointer)dptondv(m,zdp);
             } else {
               if ( !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);
               else zp = (P)BDY(t);
               c = (pointer)ptondv(CO,vv,zp);
             }
         }          }
         if ( ishomo )          if ( ishomo )
             ishomo = ishomo && ndv_ishomo(c);              ishomo = ishomo && ndv_ishomo(c);
Line 5628  int nd_get_exporigin(struct order_spec *ord)
Line 5762  int nd_get_exporigin(struct order_spec *ord)
 void nd_setup_parameters(int nvar,int max) {  void nd_setup_parameters(int nvar,int max) {
     int i,j,n,elen,ord_o,ord_l,l,s,wpd;      int i,j,n,elen,ord_o,ord_l,l,s,wpd;
     struct order_pair *op;      struct order_pair *op;
       extern int CNVars;
   
     nd_nvar = nvar;      nd_nvar = nvar;
       CNVars = nvar;
     if ( max ) {      if ( max ) {
         /* XXX */          /* XXX */
         if ( do_weyl ) nd_bpe = 32;          if ( do_weyl ) nd_bpe = 32;
Line 6759  int nmv_comp(NMV a,NMV b)
Line 6895  int nmv_comp(NMV a,NMV b)
   return -t;    return -t;
 }  }
   
   NDV dptondv(int mod,DP p)
   {
     NDV d;
     NMV m,m0;
     MP t;
     MP *a;
     int i,len,n;
   
     if ( !p ) return 0;
     for ( t = BDY(p), len = 0; t; t = NEXT(t), len++ );
     n = p->nv;
     if ( mod > 0 || mod == -1 )
       m0 = m = (NMV)MALLOC_ATOMIC_IGNORE_OFF_PAGE(len*nmv_adv);
     else
       m0 = m = MALLOC(len*nmv_adv);
     for ( i = 0, t = BDY(p); i < len; i++, NMV_ADV(m), t = NEXT(t) ) {
       dltondl(n,t->dl,DL(m));
       TD(DL(m)) = ndl_weight(DL(m));
       CZ(m) = (Z)t->c;
     }
     qsort(m0,len,nmv_adv,(int (*)(const void *,const void *))nmv_comp);
     MKNDV(NV(p),m0,len,d);
     SG(d) = SG(p);
     return d;
   }
   
 NDV dpmtondv(int mod,DPM p)  NDV dpmtondv(int mod,DPM p)
 {  {
   NDV d;    NDV d;
Line 6835  DP ndvtodp(int mod,NDV p)
Line 6997  DP ndvtodp(int mod,NDV p)
     return d;      return d;
 }  }
   
   DPM sigtodpm(SIG s)
   {
     DMM m;
     DPM d;
   
     NEWDMM(m);
     m->c = (Obj)ONE;
     m->dl = s->dl;
     m->pos = s->pos+1;
     m->next = 0;
     MKDPM(nd_nvar,m,d);
     SG(d) = s->dl->td;
     return d;
   }
   
 DPM ndvtodpm(int mod,NDV p)  DPM ndvtodpm(int mod,NDV p)
 {  {
   DMM m,m0;    DMM m,m0;
Line 9920  void parse_nd_option(NODE opt)
Line 10097  void parse_nd_option(NODE opt)
   nd_sugarweight = 0; nd_f4red =0; nd_rank0 = 0;    nd_sugarweight = 0; nd_f4red =0; nd_rank0 = 0;
   nd_f4_td = 0; nd_sba_f4step = 2; nd_sba_pot = 0; nd_sba_largelcm = 0;    nd_f4_td = 0; nd_sba_f4step = 2; nd_sba_pot = 0; nd_sba_largelcm = 0;
   nd_sba_dontsort = 0; nd_top = 0; nd_sba_redundant_check = 0;    nd_sba_dontsort = 0; nd_top = 0; nd_sba_redundant_check = 0;
     nd_sba_syz = 0; nd_sba_fixord = 0; nd_sba_grevlexgb = 0;
   
   for ( t = opt; t; t = NEXT(t) ) {    for ( t = opt; t; t = NEXT(t) ) {
     p = BDY((LIST)BDY(t));      p = BDY((LIST)BDY(t));
Line 9978  void parse_nd_option(NODE opt)
Line 10156  void parse_nd_option(NODE opt)
     } else if ( !strcmp(key,"sba_f4step") ) {      } else if ( !strcmp(key,"sba_f4step") ) {
       nd_sba_f4step = value?ZTOS((Q)value):0;        nd_sba_f4step = value?ZTOS((Q)value):0;
     } else if ( !strcmp(key,"sba_pot") ) {      } else if ( !strcmp(key,"sba_pot") ) {
       nd_sba_pot = value?1:0;        nd_sba_pot = ZTOS((Q)value);
     } else if ( !strcmp(key,"sba_largelcm") ) {      } else if ( !strcmp(key,"sba_largelcm") ) {
       nd_sba_largelcm = value?1:0;        nd_sba_largelcm = value?1:0;
     } else if ( !strcmp(key,"sba_dontsort") ) {      } else if ( !strcmp(key,"sba_dontsort") ) {
       nd_sba_dontsort = value?1:0;        nd_sba_dontsort = value?1:0;
       } else if ( !strcmp(key,"sba_syz") ) {
         nd_sba_syz = value?1:0;
       } else if ( !strcmp(key,"sba_fixord") ) {
         nd_sba_fixord = value?1:0;
       } else if ( !strcmp(key,"sba_grevlexgb") ) {
         nd_sba_grevlexgb = value?1:0;
     } else if ( !strcmp(key,"sba_redundant_check") ) {      } else if ( !strcmp(key,"sba_redundant_check") ) {
       nd_sba_redundant_check = value?1:0;        nd_sba_redundant_check = value?1:0;
     } else if ( !strcmp(key,"top") ) {      } else if ( !strcmp(key,"top") ) {
       nd_top = value?1:0;        nd_top = value?1:0;
     }      }
   }    }
     if ( nd_sba_syz ) nd_sba_dontsort = 1;
 }  }
   
 ND mdptond(DP d);  ND mdptond(DP d);
Legend:
Removed from v.1.42  
changed lines
  Added in v.1.45
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>