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

Diff for /OpenXM_contrib2/asir2000/builtin/gr.c between version 1.13 and 1.74

version 1.13, 2000/12/05 08:29:43 version 1.74, 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/gr.c,v 1.12 2000/12/05 06:59:15 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/builtin/gr.c,v 1.73 2017/08/31 02:36:20 noro Exp $
 */  */
 #include "ca.h"  #include "ca.h"
 #include "parse.h"  #include "parse.h"
 #include "base.h"  #include "base.h"
 #include "ox.h"  #include "ox.h"
   
 #define ITOS(p) (((unsigned int)(p))&0x7fffffff)  #if defined(__GNUC__)
 #define STOI(i) ((P)((unsigned int)(i)|0x80000000))  #define INLINE static inline
   #elif defined(VISUAL) || defined(__MINGW32__)
   #define INLINE __inline
   #else
   #define INLINE
   #endif
   
 #define NEXTVL(r,c) \  #define HMAG(p) (p_mag((P)BDY(p)->c))
 if(!(r)){NEWVL(r);(c)=(r);}else{NEWVL(NEXT(c));(c)=NEXT(c);}  
   
 #define HMAG(p) (p_mag(BDY(p)->c))  #define NEWDP_pairs ((DP_pairs)MALLOC(sizeof(struct dp_pairs)))
   
 struct dp_pairs {  static DP_pairs collect_pairs_of_hdlcm( DP_pairs d1, DP_pairs *prest );
         int dp1, dp2;  double get_rtime();
         DL lcm;  
         int sugar;  
         struct dp_pairs *next;  
 };  
   
 typedef struct dp_pairs *DP_pairs;  struct oEGT eg_nf,eg_nfm;
   struct oEGT eg_znfm,eg_pz,eg_np,eg_ra,eg_mc,eg_gc;
   int TP,N_BP,NMP,NFP,NDP,ZR,NZR;
   
 #define NEWDPP(a) ((a)=(DP_pairs)MALLOC(sizeof(struct dp_pairs)))  
   
 #define NEXTDPP(r,c) \  
 if(!(r)){NEWDPP(r);(c)=(r);}else{NEWDPP(NEXT(c));(c)=NEXT(c);}  
   
 struct oEGT eg_up,eg_sp,eg_spm,eg_nf,eg_nfm;  
 struct oEGT eg_znfm,eg_pz,eg_np,eg_mp,eg_ra,eg_mc,eg_gc;  
 int TP,NBP,NMP,NFP,NDP,ZR,NZR;  
   
 #define NEWDP_pairs ((DP_pairs)MALLOC(sizeof(struct dp_pairs)))  
   
 extern int (*cmpdl)();  extern int (*cmpdl)();
 extern int do_weyl;  extern int do_weyl;
   
 extern DP_Print;  extern DP_Print;
   
 void print_stat(void);  
 void init_stat(void);  
 int dp_load_t(int,DP *);  
 void dp_load(int,DP *);  
 void dp_save(int,Obj,char *);  
 void dp_make_flaglist(LIST *);  
 void dp_set_flag(Obj,Obj);  
 int membercheck(NODE,NODE);  
 int gbcheck(NODE);  
 int dl_redble(DL,DL);  
 NODE remove_reducibles(NODE,int);  
 NODE updbase(NODE,int);  
 DP_pairs criterion_F(DP_pairs);  
 int criterion_2(int,int);  
 static DP_pairs collect_pairs_of_hdlcm(DP_pairs,DP_pairs *);  
 DP_pairs criterion_M(DP_pairs);  
 DP_pairs criterion_B(DP_pairs,int);  
 DP_pairs newpairs(NODE,int);  
 DP_pairs updpairs(DP_pairs,NODE,int);  
 void _dp_nf_ptozp(NODE,DP,DP *,int,int,DP *);  
 void _dp_nf_ptozp_mpi(NODE,DP,DP *,int,int,DP *);  
 void _dp_nf(NODE,DP,DP *,int,DP *);  
 NODE gb_mod(NODE,int);  
 NODE gbd(NODE,int,NODE,NODE);  
 NODE gb(NODE,int,NODE);  
 NODE gb_f4(NODE);  
 NODE gb_f4_mod(NODE,int);  
 DP_pairs minp(DP_pairs, DP_pairs *);  
 void minsugar(DP_pairs,DP_pairs *,DP_pairs *);  
 NODE append_one(NODE,int);  
 void reducebase_dehomo(NODE,NODE *);  
 int newps_mod(DP,int);  
 int newps_nosave(DP,int,NODE);  
 int newps(DP,int,NODE);  
 void reduceall_mod(NODE,int,NODE *);  
 void reduceall(NODE,NODE *);  
 NODE NODE_sortbi(NODE,int);  
 NODE NODE_sortbi_insert(int, NODE,int);  
 NODE NODE_sortb(NODE,int);  
 NODE NODE_sortb_insert(DP,NODE,int);  
 void prim_part(DP,int,DP *);  
 void setup_arrays(NODE,int,NODE *);  
 int validhc(P,int,NODE);  
 void vlminus(VL,VL,VL *);  
 void printsubst(NODE);  
 void makesubst(VL,NODE *);  
 void pltovl(LIST,VL *);  
 void printdl(DL);  
 int DPPlength(DP_pairs);  
 void dp_gr_mod_main(LIST,LIST,Num,int,struct order_spec *,LIST *);  
 void dp_gr_main(LIST,LIST,Num,int,struct order_spec *,LIST *);  
 void dp_f4_main(LIST,LIST,struct order_spec *,LIST *);  
 void dp_f4_mod_main(LIST,LIST,int,struct order_spec *,LIST *);  
 double get_rtime();  
 void _dpmod_to_vect(DP,DL *,int *);  
 void dp_to_vect(DP,DL *,Q *);  
 NODE dp_dllist(DP f);  
 NODE symb_merge(NODE,NODE,int),_symb_merge(NODE,NODE,int);  
 extern int dp_nelim;  extern int dp_nelim;
 extern int dp_fcoeffs;  extern int dp_fcoeffs;
 static DP *ps,*psm;  static DP *ps,*psm;
Line 151  static P *psc;
Line 84  static P *psc;
   
 static int *pss;  static int *pss;
 static int psn,pslen;  static int psn,pslen;
 static int NVars,CNVars,PCoeffs;  static int NVars,CNVars;
 static VL VC;  static VL VC;
 static int NoSugar = 0;  
   int PCoeffs;
   int DP_Print = 0;
   int DP_PrintShort = 0;
   int DP_Multiple = 0;
   int DP_NFStat = 0;
   LIST Dist = 0;
   int NoGCD = 0;
   int GenTrace = 0;
   int GenSyz = 0;
   int OXCheck = -1;
   int OneZeroHomo = 0;
   int MaxDeg = 0;
   
   int NoSugar = 0;
 static int NoCriB = 0;  static int NoCriB = 0;
 static int NoGC = 0;  static int NoGC = 0;
 static int NoMC = 0;  static int NoMC = 0;
 static int NoRA = 0;  static int NoRA = 0;
 int DP_Print = 0;  
 static int DP_PrintShort = 0;  
 static int ShowMag = 0;  static int ShowMag = 0;
 static int Stat = 0;  static int Stat = 0;
 static int Multiple = 0;  int Denominator = 1;
 static int Denominator = 1;  int Top = 0;
 static int Top = 0;  int Reverse = 0;
 static int Reverse = 0;  
 static int Max_mag = 0;  static int Max_mag = 0;
 static char *Demand = 0;  static int Max_coef = 0;
   char *Demand = 0;
 static int PtozpRA = 0;  static int PtozpRA = 0;
 LIST Dist = 0;  
 int NoGCD = 0;  
 int GenTrace = 0;  
 int OXCheck = -1;  
 int doing_f4;  int doing_f4;
 NODE TraceList;  NODE TraceList;
   NODE AllTraceList;
   
 int eqdl(nv,dl1,dl2)  void Pox_cmo_rpc(NODE,Obj *);
 int nv;  void Pox_rpc(NODE,Obj *);
 DL dl1,dl2;  void Pox_pop_local(NODE,Obj *);
   
   INLINE int eqdl(int nv,DL dl1,DL dl2)
 {  {
         int i;    int i;
         int *b1,*b2;    int *p1,*p2;
   
         if ( dl1->td != dl2->td )    if ( dl1->td != dl2->td )
                 return 0;      return 0;
         for ( i = 0, b1 = dl1->d, b2 = dl2->d; i < nv; i++ )    i = nv-1;
                 if ( b1[i] != b2[i] )    p1 = dl1->d;
                         break;    p2 = dl2->d;
         if ( i == nv )    while ( i >= 7 ) {
                 return 1;      if ( *p1++ != *p2++ ) return 0;
         else      if ( *p1++ != *p2++ ) return 0;
                 return 0;      if ( *p1++ != *p2++ ) return 0;
       if ( *p1++ != *p2++ ) return 0;
       if ( *p1++ != *p2++ ) return 0;
       if ( *p1++ != *p2++ ) return 0;
       if ( *p1++ != *p2++ ) return 0;
       if ( *p1++ != *p2++ ) return 0;
       i -= 8;
     }
     switch ( i ) {
       case 6:
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         return 1;
       case 5:
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         return 1;
       case 4:
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         return 1;
       case 3:
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         return 1;
       case 2:
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         return 1;
       case 1:
         if ( *p1++ != *p2++ ) return 0;
         if ( *p1++ != *p2++ ) return 0;
         return 1;
       case 0:
         if ( *p1++ != *p2++ ) return 0;
         return 1;
       default:
         return 1;
     }
 }  }
   
 /* b[] should be cleared */  /* b[] should be cleared */
   
 void _dpmod_to_vect(f,at,b)  void _dpmod_to_vect(DP f,DL *at,int *b)
 DP f;  
 DL *at;  
 int *b;  
 {  {
         int i,nv;    int i,nv;
         MP m;    MP m;
   
         nv = f->nv;    nv = f->nv;
         for ( m = BDY(f), i = 0; m; m = NEXT(m), i++ ) {    for ( m = BDY(f), i = 0; m; m = NEXT(m), i++ ) {
                 for ( ; !eqdl(nv,m->dl,at[i]); i++ );      for ( ; !eqdl(nv,m->dl,at[i]); i++ );
                 b[i] = ITOS(m->c);      b[i] = ITOS(m->c);
         }    }
 }  }
   
 void dp_to_vect(f,at,b)  /* [t,findex] -> tf -> compressed vector */
 DP f;  
 DL *at;  void _tf_to_vect_compress(NODE tf,DL *at,CDP *b)
 Q *b;  
 {  {
         int i,nv;    int i,j,k,nv,len;
         MP m;    DL t,s,d1;
     DP f;
     MP m;
     CDP r;
   
         nv = f->nv;    t = (DL)BDY(tf);
         for ( m = BDY(f), i = 0; m; m = NEXT(m), i++ ) {    f = ps[(int)BDY(NEXT(tf))];
                 for ( ; !eqdl(nv,m->dl,at[i]); i++ );  
                 b[i] =(Q)m->c;    nv = f->nv;
         }    for ( m = BDY(f), len = 0; m; m = NEXT(m), len++ );
     r = (CDP)MALLOC(sizeof(struct oCDP));
      r->len = len;
     r->psindex = (int)BDY(NEXT(tf));
     r->body = (unsigned int *)MALLOC_ATOMIC(sizeof(unsigned int)*len);
   
     NEWDL_NOINIT(s,nv);
     for ( m = BDY(f), i = j = 0; m; m = NEXT(m), j++ ) {
       d1 = m->dl;
       s->td = t->td+d1->td;
       for ( k = 0; k < nv; k++ )
         s->d[k] = t->d[k]+d1->d[k];
       for ( ; !eqdl(nv,s,at[i]); i++ );
       r->body[j] = i;
     }
     *b = r;
 }  }
   
 NODE dp_dllist(f)  void dp_to_vect(DP f,DL *at,Q *b)
 DP f;  
 {  {
         MP m;    int i,nv;
         NODE mp,mp0;    MP m;
   
         if ( !f )    nv = f->nv;
                 return 0;    for ( m = BDY(f), i = 0; m; m = NEXT(m), i++ ) {
         mp0 = 0;      for ( ; !eqdl(nv,m->dl,at[i]); i++ );
         for ( m = BDY(f); m; m = NEXT(m) ) {      b[i] =(Q)m->c;
                 NEXTNODE(mp0,mp); BDY(mp) = (pointer)m->dl;    }
         }  
         NEXT(mp) = 0;  
         return mp0;  
 }  }
   
 void pdl(f)  NODE dp_dllist(DP f)
 NODE f;  
 {  {
         while ( f ) {    MP m;
                 printdl(BDY(f)); f = NEXT(f);    NODE mp,mp0;
         }  
         fflush(stdout);    if ( !f )
         printf("\n");      return 0;
     mp0 = 0;
     for ( m = BDY(f); m; m = NEXT(m) ) {
       NEXTNODE(mp0,mp); BDY(mp) = (pointer)m->dl;
     }
     NEXT(mp) = 0;
     return mp0;
 }  }
   
 void dp_gr_main(f,v,homo,modular,ord,rp)  NODE mul_dllist(DL d,DP f)
 LIST f,v;  
 Num homo;  
 int modular;  
 struct order_spec *ord;  
 LIST *rp;  
 {  {
         int i,mindex,m,nochk;    MP m;
         struct order_spec ord1;    NODE mp,mp0;
         VL fv,vv,vc;    DL t,d1;
         NODE fd,fd0,fi,fi0,r,r0,t,subst,x,s,xx;    int i,nv;
   
         mindex = 0; nochk = 0; dp_fcoeffs = 0;    if ( !f )
         get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);      return 0;
         NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;    nv = NV(f);
         CNVars = homo ? NVars+1 : NVars;    mp0 = 0;
         if ( ord->id && NVars != ord->nv )    for ( m = BDY(f); m; m = NEXT(m) ) {
                 error("dp_gr_main : invalid order specification");      NEXTNODE(mp0,mp);
         initd(ord);      NEWDL_NOINIT(t,nv);
         if ( homo ) {      d1 = m->dl;
                 homogenize_order(ord,NVars,&ord1);      t->td = d->td+d1->td;
                 for ( fd0 = fi0 = 0, t = BDY(f); t; t = NEXT(t) ) {      for ( i = 0; i < nv; i++ )
                         NEXTNODE(fd0,fd); NEXTNODE(fi0,fi);        t->d[i] = d->d[i]+d1->d[i];
                         ptod(CO,vv,(P)BDY(t),(DP *)&BDY(fi)); dp_homo((DP)BDY(fi),(DP *)&BDY(fd));      BDY(mp) = (pointer)t;
                 }    }
                 if ( fd0 ) NEXT(fd) = 0;    NEXT(mp) = 0;
                 if ( fi0 ) NEXT(fi) = 0;    return mp0;
                 initd(&ord1);  
         } else {  
                 for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {  
                         NEXTNODE(fd0,fd); ptod(CO,vv,(P)BDY(t),(DP *)&BDY(fd));  
                 }  
                 if ( fd0 ) NEXT(fd) = 0;  
                 fi0 = fd0;  
         }  
         if ( modular < 0 ) {  
                 modular = -modular; nochk = 1;  
         }  
         if ( modular )  
                 m = modular > 1 ? modular : lprime[mindex];  
         else  
                 m = 0;  
         makesubst(vc,&subst);  
         setup_arrays(fd0,0,&s);  
         init_stat();  
         while ( 1 ) {  
                 if ( homo ) {  
                         initd(&ord1); CNVars = NVars+1;  
                 }  
                 if ( DP_Print && modular ) {  
                         fprintf(asir_out,"mod= %d, eval = ",m); printsubst(subst);  
                 }  
                 x = gb(s,m,subst);  
                 if ( x ) {  
                         if ( homo ) {  
                                 reducebase_dehomo(x,&xx); x = xx;  
                                 initd(ord); CNVars = NVars;  
                         }  
                         reduceall(x,&xx); x = xx;  
                         if ( modular ) {  
                                 if ( nochk || (membercheck(fi0,x) && gbcheck(x)) )  
                                         break;  
                         } else  
                                 break;  
                 }  
                 if ( modular )  
                         if ( modular > 1 ) {  
                                 *rp = 0; return;  
                         } else  
                                 m = lprime[++mindex];  
                 makesubst(vc,&subst);  
                 psn = length(s);  
                 for ( i = psn; i < pslen; i++ ) {  
                         pss[i] = 0; psh[i] = 0; psc[i] = 0; ps[i] = 0;  
                 }  
         }  
         for ( r0 = 0; x; x = NEXT(x) ) {  
                 NEXTNODE(r0,r); dp_load((int)BDY(x),&ps[(int)BDY(x)]);  
                 dtop(CO,vv,ps[(int)BDY(x)],(P *)&BDY(r));  
         }  
         if ( r0 ) NEXT(r) = 0;  
         MKLIST(*rp,r0);  
         print_stat();  
         if ( ShowMag )  
                 fprintf(asir_out,"\nMax_mag=%d\n",Max_mag);  
 }  }
   
 void dp_gr_mod_main(f,v,homo,m,ord,rp)  void pdl(NODE f)
 LIST f,v;  
 Num homo;  
 int m;  
 struct order_spec *ord;  
 LIST *rp;  
 {  {
         struct order_spec ord1;    while ( f ) {
         VL fv,vv,vc;      printdl(BDY(f)); f = NEXT(f);
         NODE fd,fd0,r,r0,t,x,s,xx;    }
         DP a,b,c;    fflush(stdout);
     printf("\n");
   }
   
         get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);  void dp_gr_main(LIST f,LIST v,Num homo,int modular,int field,struct order_spec *ord,LIST *rp)
         NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;  {
         CNVars = homo ? NVars+1 : NVars;    int i,mindex,m,nochk;
         if ( ord->id && NVars != ord->nv )    struct order_spec *ord1;
                 error("dp_gr_mod_main : invalid order specification");    Q q;
         initd(ord);    VL fv,vv,vc;
         if ( homo ) {    NODE fd,fd0,fi,fi0,r,r0,t,subst,x,s,xx;
                 for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {    NODE ind,ind0;
                         ptod(CO,vv,(P)BDY(t),&a); dp_homo(a,&b);    LIST trace,gbindex;
                         if ( PCoeffs )    int input_is_dp = 0;
                                 dp_mod(b,m,0,&c);  
                         else    mindex = 0; nochk = 0; dp_fcoeffs = field;
                                 _dp_mod(b,m,(NODE)0,&c);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);
                         if ( c ) {    NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;
                                 NEXTNODE(fd0,fd); BDY(fd) = (pointer)c;    CNVars = homo ? NVars+1 : NVars;
                         }    if ( ord->id && NVars != ord->nv )
                 }      error("dp_gr_main : invalid order specification");
                 homogenize_order(ord,NVars,&ord1); initd(&ord1);    initd(ord);
         } else {    if ( homo ) {
                 for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {      homogenize_order(ord,NVars,&ord1);
                         ptod(CO,vv,(P)BDY(t),&b);      for ( fd0 = fi0 = 0, t = BDY(f); t; t = NEXT(t) ) {
                         if ( PCoeffs )        NEXTNODE(fd0,fd); NEXTNODE(fi0,fi);
                                 dp_mod(b,m,0,&c);        if ( BDY(t) && OID(BDY(t)) == O_DP ) {
                         else          dp_sort((DP)BDY(t),(DP *)&BDY(fi)); input_is_dp = 1;
                                 _dp_mod(b,m,0,&c);        } else
                         if ( c ) {          ptod(CO,vv,(P)BDY(t),(DP *)&BDY(fi));
                                 NEXTNODE(fd0,fd); BDY(fd) = (pointer)c;        dp_homo((DP)BDY(fi),(DP *)&BDY(fd));
                         }      }
                 }      if ( fd0 ) NEXT(fd) = 0;
         }      if ( fi0 ) NEXT(fi) = 0;
         if ( fd0 ) NEXT(fd) = 0;      initd(ord1);
         setup_arrays(fd0,m,&s);    } else {
         init_stat();      for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {
         if ( homo ) {        NEXTNODE(fd0,fd);
                 initd(&ord1); CNVars = NVars+1;        if ( BDY(t) && OID(BDY(t)) == O_DP ) {
         }          dp_sort((DP)BDY(t),(DP *)&BDY(fd)); input_is_dp = 1;
         x = gb_mod(s,m);        } else
         if ( homo ) {          ptod(CO,vv,(P)BDY(t),(DP *)&BDY(fd));
                 reducebase_dehomo(x,&xx); x = xx;      }
                 initd(ord); CNVars = NVars;      if ( fd0 ) NEXT(fd) = 0;
         }      fi0 = fd0;
         reduceall_mod(x,m,&xx); x = xx;    }
         if ( PCoeffs )    if ( modular < 0 ) {
                 for ( r0 = 0; x; x = NEXT(x) ) {      modular = -modular; nochk = 1;
                         NEXTNODE(r0,r); mdtop(CO,m,vv,ps[(int)BDY(x)],(P *)&BDY(r));    }
                 }    if ( modular )
         else      m = modular > 1 ? modular : get_lprime(mindex);
                 for ( r0 = 0; x; x = NEXT(x) ) {    else
                         NEXTNODE(r0,r); _dtop_mod(CO,vv,ps[(int)BDY(x)],(P *)&BDY(r));      m = 0;
                 }    makesubst(vc,&subst);
         print_stat();    setup_arrays(fd0,0,&s);
         if ( r0 ) NEXT(r) = 0;    init_stat();
         MKLIST(*rp,r0);    while ( 1 ) {
       if ( homo ) {
         initd(ord1); CNVars = NVars+1;
       }
       if ( DP_Print && modular ) {
         fprintf(asir_out,"mod= %d, eval = ",m); printsubst(subst);
       }
       x = gb(s,m,subst);
       if ( x ) {
         if ( homo ) {
           reducebase_dehomo(x,&xx); x = xx;
           initd(ord); CNVars = NVars;
         }
         reduceall(x,&xx); x = xx;
         if ( modular ) {
           if ( nochk || (membercheck(fi0,x) && gbcheck(x)) )
             break;
         } else
           break;
       }
       if ( modular )
         if ( modular > 1 ) {
           *rp = 0; return;
         } else
           m = get_lprime(++mindex);
       makesubst(vc,&subst);
       psn = length(s);
       for ( i = psn; i < pslen; i++ ) {
         pss[i] = 0; psh[i] = 0; psc[i] = 0; ps[i] = 0;
       }
     }
     for ( r0 = 0, ind0 = 0; x; x = NEXT(x) ) {
       NEXTNODE(r0,r); dp_load((int)BDY(x),&ps[(int)BDY(x)]);
       if ( input_is_dp )
         BDY(r) = (pointer)ps[(int)BDY(x)];
       else
         dtop(CO,vv,ps[(int)BDY(x)],(Obj *)&BDY(r));
       NEXTNODE(ind0,ind);
       STOQ((int)BDY(x),q); BDY(ind) = q;
     }
     if ( r0 ) NEXT(r) = 0;
     if ( ind0 ) NEXT(ind) = 0;
     MKLIST(*rp,r0);
     MKLIST(gbindex,ind0);
   
     if ( GenTrace && OXCheck < 0 ) {
   
       x = AllTraceList;
       for ( r = 0; x; x = NEXT(x) ) {
         MKNODE(r0,BDY(x),r); r = r0;
       }
       MKLIST(trace,r);
       r0 = mknode(3,*rp,gbindex,trace);
       MKLIST(*rp,r0);
     }
     print_stat();
     if ( ShowMag )
       fprintf(asir_out,"\nMax_mag=%d, Max_coef=%d\n",Max_mag, Max_coef);
 }  }
   
 void dp_f4_main(f,v,ord,rp)  void dp_interreduce(LIST f,LIST v,int field,struct order_spec *ord,LIST *rp)
 LIST f,v;  
 struct order_spec *ord;  
 LIST *rp;  
 {  {
         int i,mindex,m,nochk;    int i,mindex,m,nochk;
         struct order_spec ord1;    struct order_spec *ord1;
         VL fv,vv,vc;    Q q;
         NODE fd,fd0,fi,fi0,r,r0,t,subst,x,s,xx;    VL fv,vv,vc;
     NODE fd,fd0,fi,fi0,r,r0,t,subst,x,s,xx;
     NODE ind,ind0;
     LIST trace,gbindex;
     int input_is_dp = 0;
   
         dp_fcoeffs = 0;    mindex = 0; nochk = 0; dp_fcoeffs = field;
         get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);
         NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;    NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;
         CNVars = NVars;    CNVars = NVars;
         if ( ord->id && NVars != ord->nv )    if ( ord->id && NVars != ord->nv )
                 error("dp_f4_main : invalid order specification");      error("dp_interreduce : invalid order specification");
         initd(ord);    initd(ord);
         for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {    for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {
                 NEXTNODE(fd0,fd); ptod(CO,vv,(P)BDY(t),(DP *)&BDY(fd));      NEXTNODE(fd0,fd);
         }      if ( BDY(t) && OID(BDY(t)) == O_DP ) {
         if ( fd0 ) NEXT(fd) = 0;        dp_sort((DP)BDY(t),(DP *)&BDY(fd)); input_is_dp = 1;
         setup_arrays(fd0,0,&s);      } else
         x = gb_f4(s);        ptod(CO,vv,(P)BDY(t),(DP *)&BDY(fd));
         reduceall(x,&xx); x = xx;    }
         for ( r0 = 0; x; x = NEXT(x) ) {    if ( fd0 ) NEXT(fd) = 0;
                 NEXTNODE(r0,r); dp_load((int)BDY(x),&ps[(int)BDY(x)]);    fi0 = fd0;
                 dtop(CO,vv,ps[(int)BDY(x)],(P *)&BDY(r));    setup_arrays(fd0,0,&s);
         }    init_stat();
         if ( r0 ) NEXT(r) = 0;    x = s;
         MKLIST(*rp,r0);    reduceall(x,&xx); x = xx;
     for ( r0 = 0, ind0 = 0; x; x = NEXT(x) ) {
       NEXTNODE(r0,r); dp_load((int)BDY(x),&ps[(int)BDY(x)]);
       if ( input_is_dp )
         BDY(r) = (pointer)ps[(int)BDY(x)];
       else
         dtop(CO,vv,ps[(int)BDY(x)],(Obj *)&BDY(r));
       NEXTNODE(ind0,ind);
       STOQ((int)BDY(x),q); BDY(ind) = q;
     }
     if ( r0 ) NEXT(r) = 0;
     if ( ind0 ) NEXT(ind) = 0;
     MKLIST(*rp,r0);
     MKLIST(gbindex,ind0);
 }  }
   
 void dp_f4_mod_main(f,v,m,ord,rp)  void dp_gr_mod_main(LIST f,LIST v,Num homo,int m,struct order_spec *ord,LIST *rp)
 LIST f,v;  
 int m;  
 struct order_spec *ord;  
 LIST *rp;  
 {  {
         int i;    struct order_spec *ord1;
         struct order_spec ord1;    VL fv,vv,vc;
         VL fv,vv,vc;    NODE fd,fd0,r,r0,t,x,s,xx;
         DP b,c,c1;    DP a,b,c;
         NODE fd,fd0,fi,fi0,r,r0,t,subst,x,s,xx;    extern struct oEGT eg_red_mod;
     int input_is_dp = 0;
   
         dp_fcoeffs = 0;    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);
         get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);    NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;
         NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;    CNVars = homo ? NVars+1 : NVars;
         CNVars = NVars;    if ( ord->id && NVars != ord->nv )
         if ( ord->id && NVars != ord->nv )      error("dp_gr_mod_main : invalid order specification");
                 error("dp_f4_mod_main : invalid order specification");    initd(ord);
         initd(ord);    if ( homo ) {
         for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {      for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {
                 ptod(CO,vv,(P)BDY(t),&b);        if ( BDY(t) && OID(BDY(t)) == O_DP ) {
                 _dp_mod(b,m,0,&c);          dp_sort((DP)BDY(t),&a); input_is_dp = 1;
                 _dp_monic(c,m,&c1);        } else
                 if ( c ) {          ptod(CO,vv,(P)BDY(t),&a);
                         NEXTNODE(fd0,fd); BDY(fd) = (pointer)c1;        dp_homo(a,&b);
                 }        if ( PCoeffs )
         }          dp_mod(b,m,0,&c);
         if ( fd0 ) NEXT(fd) = 0;        else
         setup_arrays(fd0,m,&s);          _dp_mod(b,m,(NODE)0,&c);
         x = gb_f4_mod(s,m);        if ( c ) {
         reduceall_mod(x,m,&xx); x = xx;          NEXTNODE(fd0,fd); BDY(fd) = (pointer)c;
         for ( r0 = 0; x; x = NEXT(x) ) {        }
                 NEXTNODE(r0,r); _dtop_mod(CO,vv,ps[(int)BDY(x)],(P *)&BDY(r));      }
         }      homogenize_order(ord,NVars,&ord1); initd(ord1);
         if ( r0 ) NEXT(r) = 0;    } else {
         MKLIST(*rp,r0);      for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {
         if ( BDY(t) && OID(BDY(t)) == O_DP ) {
           dp_sort((DP)BDY(t),&b); input_is_dp = 1;
         } else
           ptod(CO,vv,(P)BDY(t),&b);
         if ( PCoeffs )
           dp_mod(b,m,0,&c);
         else
           _dp_mod(b,m,0,&c);
         if ( c ) {
           NEXTNODE(fd0,fd); BDY(fd) = (pointer)c;
         }
       }
     }
     if ( fd0 ) NEXT(fd) = 0;
     setup_arrays(fd0,m,&s);
     init_stat();
     if ( homo ) {
       initd(ord1); CNVars = NVars+1;
     }
   /* init_eg(&eg_red_mod); */
     x = gb_mod(s,m);
   /* print_eg("Red_mod",&eg_red_mod); */
     if ( homo ) {
       reducebase_dehomo(x,&xx); x = xx;
       initd(ord); CNVars = NVars;
     }
     reduceall_mod(x,m,&xx); x = xx;
     if ( PCoeffs )
       for ( r0 = 0; x; x = NEXT(x) ) {
         NEXTNODE(r0,r);
         if ( input_is_dp )
           mdtodp(ps[(int)BDY(x)],(DP *)&BDY(r));
         else
           mdtop(CO,m,vv,ps[(int)BDY(x)],(P *)&BDY(r));
       }
     else
       for ( r0 = 0; x; x = NEXT(x) ) {
         NEXTNODE(r0,r);
         if ( input_is_dp )
           _mdtodp(ps[(int)BDY(x)],(DP *)&BDY(r));
         else
           _dtop_mod(CO,vv,ps[(int)BDY(x)],(P *)&BDY(r));
       }
     print_stat();
     if ( r0 ) NEXT(r) = 0;
     MKLIST(*rp,r0);
 }  }
   
 NODE gb_f4(f)  void dp_f4_main(LIST f,LIST v,struct order_spec *ord,LIST *rp)
 NODE f;  
 {  {
         int i,j,k,nh,row,col,nv;    int homogen;
         NODE r,g,gall;    VL fv,vv,vc;
         NODE s,s0;    NODE fd,fd0,r,r0,t,x,s,xx;
         DP_pairs d,dm,dr,t;    int input_is_dp = 0;
         DP h,nf,nf1,f1,f2,f21,f21r,sp,sp1,sd,sdm,tdp;  
         MP mp,mp0;  
         NODE blist,bt,nt;  
         DL *ht,*at;  
         MAT mat,nm;  
         int *colstat;  
         int *rind,*cind;  
         int rank,nred;  
         Q dn;  
         struct oEGT tmp0,tmp1,tmp2,eg_split_symb,eg_split_elim;  
         extern struct oEGT eg_mod,eg_elim,eg_chrem,eg_gschk,eg_intrat,eg_symb;  
   
         init_eg(&eg_mod); init_eg(&eg_elim); init_eg(&eg_chrem);    dp_fcoeffs = 0;
         init_eg(&eg_gschk); init_eg(&eg_intrat); init_eg(&eg_symb);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);
     NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;
     CNVars = NVars;
     if ( ord->id && NVars != ord->nv )
       error("dp_f4_main : invalid order specification");
     initd(ord);
     for ( fd0 = 0, t = BDY(f), homogen = 1; t; t = NEXT(t) ) {
       NEXTNODE(fd0,fd);
       if ( BDY(t) && OID(BDY(t)) == O_DP ) {
         dp_sort((DP)BDY(t),(DP *)&BDY(fd)); input_is_dp = 1;
       } else
         ptod(CO,vv,(P)BDY(t),(DP *)&BDY(fd));
       if ( homogen )
         homogen = dp_homogeneous(BDY(fd));
     }
     if ( fd0 ) NEXT(fd) = 0;
     setup_arrays(fd0,0,&s);
     x = gb_f4(s);
     if ( !homogen ) {
       reduceall(x,&xx); x = xx;
     }
     for ( r0 = 0; x; x = NEXT(x) ) {
       NEXTNODE(r0,r); dp_load((int)BDY(x),&ps[(int)BDY(x)]);
       if ( input_is_dp )
         BDY(r) = (pointer)ps[(int)BDY(x)];
       else
         dtop(CO,vv,ps[(int)BDY(x)],(Obj *)&BDY(r));
     }
     if ( r0 ) NEXT(r) = 0;
     MKLIST(*rp,r0);
   }
   
         doing_f4 = 1;  void dp_f4_mod_main(LIST f,LIST v,int m,struct order_spec *ord,LIST *rp)
         for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {  {
                 i = (int)BDY(r);    int homogen;
                 d = updpairs(d,g,i);    VL fv,vv,vc;
                 g = updbase(g,i);    DP b,c,c1;
                 gall = append_one(gall,i);    NODE fd,fd0,r,r0,t,x,s,xx;
         }    int input_is_dp = 0;
         if ( gall )  
                 nv = ((DP)ps[(int)BDY(gall)])->nv;  
         while ( d ) {  
                 get_eg(&tmp0);  
                 minsugar(d,&dm,&dr); d = dr;  
                 if ( DP_Print )  
                         fprintf(asir_out,"sugar=%d\n",dm->sugar);  
                 blist = 0; s0 = 0;  
                 /* asph : sum of all head terms of spoly */  
                 for ( t = dm; t; t = NEXT(t) ) {  
                         dp_sp(ps[t->dp1],ps[t->dp2],&sp);  
                         if ( sp ) {  
                                 MKNODE(bt,sp,blist); blist = bt;  
                                 s0 = symb_merge(s0,dp_dllist(sp),nv);  
                         }  
                 }  
                 /* s0 : all the terms appeared in symbolic redunction */  
                 for ( s = s0, nred = 0; s; s = NEXT(s) ) {  
                         for ( r = gall; r; r = NEXT(r) )  
                                 if ( _dl_redble(BDY(ps[(int)BDY(r)])->dl,BDY(s),nv) )  
                                         break;  
                         if ( r ) {  
                                 dltod(BDY(s),nv,&tdp);  
                                 dp_subd(tdp,ps[(int)BDY(r)],&sd);  
                                 muld(CO,sd,ps[(int)BDY(r)],&f2);  
                                 MKNODE(bt,f2,blist); blist = bt;  
                                 s = symb_merge(s,dp_dllist(f2),nv);  
                                 nred++;  
                         }  
                 }  
   
                 /* the first nred polys in blist are reducers */    dp_fcoeffs = 0;
                 /* row = the number of all the polys */    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&vc);
                 for ( r = blist, row = 0; r; r = NEXT(r), row++ );    NVars = length((NODE)vv); PCoeffs = vc ? 1 : 0; VC = vc;
                 ht = (DL *)MALLOC(nred*sizeof(DL));    CNVars = NVars;
                 for ( r = blist, i = 0; i < nred; r = NEXT(r), i++ )    if ( ord->id && NVars != ord->nv )
                         ht[i] = BDY((DP)BDY(r))->dl;      error("dp_f4_mod_main : invalid order specification");
                 for ( s = s0, col = 0; s; s = NEXT(s), col++ );    initd(ord);
                 at = (DL *)MALLOC(col*sizeof(DL));    for ( fd0 = 0, t = BDY(f), homogen = 1; t; t = NEXT(t) ) {
                 for ( s = s0, i = 0; i < col; s = NEXT(s), i++ )      if ( BDY(t) && OID(BDY(t)) == O_DP ) {
                         at[i] = (DL)BDY(s);        dp_sort((DP)BDY(t),&b); input_is_dp = 1;
                 MKMAT(mat,row,col);      } else
                 for ( i = 0, r = blist; i < row; r = NEXT(r), i++ )        ptod(CO,vv,(P)BDY(t),&b);
                         dp_to_vect(BDY(r),at,(Q *)mat->body[i]);      if ( homogen )
                 get_eg(&tmp1); add_eg(&eg_symb,&tmp0,&tmp1);        homogen = dp_homogeneous(b);
                 init_eg(&eg_split_symb); add_eg(&eg_split_symb,&tmp0,&tmp1);      _dp_mod(b,m,0,&c);
                 if ( DP_Print ) {      _dp_monic(c,m,&c1);
                         print_eg("Symb",&eg_split_symb);      if ( c ) {
                         fprintf(asir_out,"mat : %d x %d",row,col);        NEXTNODE(fd0,fd); BDY(fd) = (pointer)c1;
                         fflush(asir_out);      }
                 }    }
     if ( fd0 ) NEXT(fd) = 0;
     setup_arrays(fd0,m,&s);
     init_stat();
     if ( do_weyl )
       x = gb_f4_mod_old(s,m);
     else
       x = gb_f4_mod(s,m);
     if ( !homogen ) {
       reduceall_mod(x,m,&xx); x = xx;
     }
     for ( r0 = 0; x; x = NEXT(x) ) {
       NEXTNODE(r0,r);
       if ( input_is_dp )
         _mdtodp(ps[(int)BDY(x)],(DP *)&BDY(r));
       else
         _dtop_mod(CO,vv,ps[(int)BDY(x)],(P *)&BDY(r));
     }
     if ( r0 ) NEXT(r) = 0;
     MKLIST(*rp,r0);
     print_stat();
   }
   
   NODE gb_f4(NODE f)
   {
     int i,k,nh,row,col,nv;
     NODE r,g,gall;
     NODE s,s0;
     DP_pairs d,dm,dr,t;
     DP nf,nf1,f2,sp,sd,tdp;
     MP mp,mp0;
     NODE blist,bt;
     DL *ht,*at;
     MAT mat,nm;
     int *rind,*cind;
     int rank,nred;
     Q dn;
     struct oEGT tmp0,tmp1,eg_split_symb;
     extern struct oEGT eg_mod,eg_elim,eg_chrem,eg_gschk,eg_intrat,eg_symb;
   
     init_eg(&eg_mod); init_eg(&eg_elim); init_eg(&eg_chrem);
     init_eg(&eg_gschk); init_eg(&eg_intrat); init_eg(&eg_symb);
   
     doing_f4 = 1;
     for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {
       i = (int)BDY(r);
       d = updpairs(d,g,i);
       g = updbase(g,i);
       gall = append_one(gall,i);
     }
     if ( gall )
       nv = ((DP)ps[(int)BDY(gall)])->nv;
     while ( d ) {
       get_eg(&tmp0);
       minsugar(d,&dm,&dr); d = dr;
       if ( DP_Print )
         fprintf(asir_out,"sugar=%d\n",dm->sugar);
       blist = 0; s0 = 0;
       /* asph : sum of all head terms of spoly */
       for ( t = dm; t; t = NEXT(t) ) {
         dp_sp(ps[t->dp1],ps[t->dp2],&sp);
         if ( sp ) {
           MKNODE(bt,sp,blist); blist = bt;
           s0 = symb_merge(s0,dp_dllist(sp),nv);
         }
       }
       /* s0 : all the terms appeared in symbolic redunction */
       for ( s = s0, nred = 0; s; s = NEXT(s) ) {
         for ( r = gall;  r; r = NEXT(r) )
           if ( _dl_redble(BDY(ps[(int)BDY(r)])->dl,BDY(s),nv) )
             break;
         if ( r ) {
           dltod(BDY(s),nv,&tdp);
           dp_subd(tdp,ps[(int)BDY(r)],&sd);
           muld(CO,sd,ps[(int)BDY(r)],&f2);
           MKNODE(bt,f2,blist); blist = bt;
           s = symb_merge(s,dp_dllist(f2),nv);
           nred++;
         }
       }
   
       /* the first nred polys in blist are reducers */
       /* row = the number of all the polys */
       for ( r = blist, row = 0; r; r = NEXT(r), row++ );
       ht = (DL *)MALLOC(nred*sizeof(DL));
       for ( r = blist, i = 0; i < nred; r = NEXT(r), i++ )
         ht[i] = BDY((DP)BDY(r))->dl;
       for ( s = s0, col = 0; s; s = NEXT(s), col++ );
       at = (DL *)MALLOC(col*sizeof(DL));
       for ( s = s0, i = 0; i < col; s = NEXT(s), i++ )
         at[i] = (DL)BDY(s);
       MKMAT(mat,row,col);
       for ( i = 0, r = blist; i < row; r = NEXT(r), i++ )
         dp_to_vect(BDY(r),at,(Q *)mat->body[i]);
       get_eg(&tmp1); add_eg(&eg_symb,&tmp0,&tmp1);
       init_eg(&eg_split_symb); add_eg(&eg_split_symb,&tmp0,&tmp1);
       if ( DP_Print ) {
         print_eg("Symb",&eg_split_symb);
         fprintf(asir_out,"mat : %d x %d",row,col);
         fflush(asir_out);
       }
 #if 0  #if 0
                 rank = generic_gauss_elim_hensel(mat,&nm,&dn,&rind,&cind);      rank = generic_gauss_elim_hensel(mat,&nm,&dn,&rind,&cind);
 #else  #else
                 rank = generic_gauss_elim(mat,&nm,&dn,&rind,&cind);      rank = generic_gauss_elim(mat,&nm,&dn,&rind,&cind);
 #endif  #endif
                 if ( DP_Print )      if ( DP_Print )
                         fprintf(asir_out,"done rank = %d\n",rank,row,col);        fprintf(asir_out,"done rank = %d\n",rank,row,col);
                 for ( i = 0; i < rank; i++ ) {      for ( i = 0; i < rank; i++ ) {
                         for ( k = 0; k < nred; k++ )        for ( k = 0; k < nred; k++ )
                                 if ( !cmpdl(nv,at[rind[i]],ht[k]) )          if ( !cmpdl(nv,at[rind[i]],ht[k]) )
                                         break;            break;
                         if ( k == nred ) {        if ( k == nred ) {
                                 /* this is a new base */          /* this is a new base */
                                 mp0 = 0;          mp0 = 0;
                                 NEXTMP(mp0,mp); mp->dl = at[rind[i]]; mp->c = (P)dn;          NEXTMP(mp0,mp); mp->dl = at[rind[i]]; mp->c = (Obj)dn;
                                 for ( k = 0; k < col-rank; k++ )          for ( k = 0; k < col-rank; k++ )
                                         if ( nm->body[i][k] ) {            if ( nm->body[i][k] ) {
                                                 NEXTMP(mp0,mp); mp->dl = at[cind[k]];              NEXTMP(mp0,mp); mp->dl = at[cind[k]];
                                                 mp->c = (P)nm->body[i][k];              mp->c = (Obj)nm->body[i][k];
                                         }            }
                                 NEXT(mp) = 0;          NEXT(mp) = 0;
                                 MKDP(nv,mp0,nf); nf->sugar = dm->sugar;          MKDP(nv,mp0,nf); nf->sugar = dm->sugar;
                                 dp_ptozp(nf,&nf1);          dp_ptozp(nf,&nf1);
                                 nh = newps(nf1,0,0);          nh = newps(nf1,0,0);
                                 d = updpairs(d,g,nh);          d = updpairs(d,g,nh);
                                 g = updbase(g,nh);          g = updbase(g,nh);
                                 gall = append_one(gall,nh);          gall = append_one(gall,nh);
                         }        }
                 }      }
         }    }
         if ( DP_Print ) {    if ( DP_Print ) {
                 print_eg("Symb",&eg_symb);      print_eg("Symb",&eg_symb);
                 print_eg("Mod",&eg_mod); print_eg("GaussElim",&eg_elim);      print_eg("Mod",&eg_mod); print_eg("GaussElim",&eg_elim);
                 print_eg("ChRem",&eg_chrem); print_eg("IntToRat",&eg_intrat);      print_eg("ChRem",&eg_chrem); print_eg("IntToRat",&eg_intrat);
                 print_eg("Check",&eg_gschk);      print_eg("Check",&eg_gschk);
         }    }
         return g;    return g;
 }  }
   
 /* initial bases are monic */  /* initial bases are monic */
   
 NODE gb_f4_mod(f,m)  unsigned int **psca;
 NODE f;  GeoBucket create_bucket();
 int m;  DL remove_head_bucket(GeoBucket,int);
   
   NODE gb_f4_mod(NODE f,int m)
 {  {
         int i,j,k,nh,row,col,nv;    int i,j,k,nh,row,col,nv;
         NODE r,g,gall;    NODE r,g,gall;
         NODE s,s0;    NODE s,s0;
         DP_pairs d,dm,dr,t;    DP_pairs d,dm,dr,t;
         DP h,nf,f1,f2,f21,f21r,sp,sp1,sd,sdm,tdp;    DP nf,sp,sd,tdp;
         MP mp,mp0;    MP mp,mp0;
         NODE blist,bt,nt;    NODE blist,bt,bt1,dt;
         DL *ht,*at,*st;    DL *at,*st;
         int **spmat,**redmat;    int **spmat;
         int *colstat,*w;    CDP *redmat;
         int rank,nred,nsp,nonzero,spcol;    int *colstat,*w,*w1;
         int *indred,*isred,*ri;    int rank,nred,nsp,nsp0,nonzero,spcol;
         struct oEGT tmp0,tmp1,tmp2,eg_split_symb,eg_split_elim1,eg_split_elim2;    int *indred,*isred;
         extern struct oEGT eg_symb,eg_elim1,eg_elim2;    CDP ri;
     int pscalen;
     GeoBucket bucket;
     DL head;
     struct oEGT tmp0,tmp1,eg_split_symb,eg_split_conv,eg_split_elim1,eg_split_elim2;
     extern struct oEGT eg_symb,eg_conv,eg_elim1,eg_elim2;
   
         init_eg(&eg_symb); init_eg(&eg_elim1); init_eg(&eg_elim2);    /* initialize coeffcient array list of ps[] */
         for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {    pscalen = pslen;
                 i = (int)BDY(r);    psca = (unsigned int **)MALLOC(pscalen*sizeof(unsigned int *));
                 d = updpairs(d,g,i);  
                 g = updbase(g,i);  
                 gall = append_one(gall,i);  
         }  
         if ( gall )  
                 nv = ((DP)ps[(int)BDY(gall)])->nv;  
         while ( d ) {  
                 get_eg(&tmp0);  
                 minsugar(d,&dm,&dr); d = dr;  
                 if ( DP_Print )  
                         fprintf(asir_out,"sugar=%d\n",dm->sugar);  
                 blist = 0; s0 = 0;  
                 /* asph : sum of all head terms of spoly */  
                 for ( t = dm; t; t = NEXT(t) ) {  
                         _dp_sp_mod(ps[t->dp1],ps[t->dp2],m,&sp);  
                         if ( sp ) {  
                                 MKNODE(bt,sp,blist); blist = bt;  
                                 s0 = symb_merge(s0,dp_dllist(sp),nv);  
                         }  
                 }  
                 /* s0 : all the terms appeared in symbolic redunction */  
                 for ( s = s0, nred = 0; s; s = NEXT(s) ) {  
                         for ( r = gall; r; r = NEXT(r) )  
                                 if ( _dl_redble(BDY(ps[(int)BDY(r)])->dl,BDY(s),nv) )  
                                         break;  
                         if ( r ) {  
                                 dltod(BDY(s),nv,&tdp);  
                                 dp_subd(tdp,ps[(int)BDY(r)],&sd);  
                                 _dp_mod(sd,m,0,&sdm);  
                                 mulmd_dup(m,sdm,ps[(int)BDY(r)],&f2);  
                                 MKNODE(bt,f2,blist); blist = bt;  
                                 s = symb_merge(s,dp_dllist(f2),nv);  
                                 nred++;  
                         }  
                 }  
   
                 get_eg(&tmp1); add_eg(&eg_symb,&tmp0,&tmp1);  
                 init_eg(&eg_split_symb); add_eg(&eg_split_symb,&tmp0,&tmp1);  
   
                 /* the first nred polys in blist are reducers */    init_eg(&eg_symb); init_eg(&eg_conv); init_eg(&eg_elim1); init_eg(&eg_elim2);
                 /* row = the number of all the polys */    for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {
                 for ( r = blist, row = 0; r; r = NEXT(r), row++ );      i = (int)BDY(r);
       d = updpairs(d,g,i);
       g = updbase(g,i);
       gall = append_one(gall,i);
       dptoca(ps[i],&psca[i]);
     }
     if ( gall )
       nv = ((DP)ps[(int)BDY(gall)])->nv;
     while ( d ) {
       get_eg(&tmp0);
       minsugar(d,&dm,&dr); d = dr;
       if ( DP_Print )
         fprintf(asir_out,"sugar=%d\n",dm->sugar);
       blist = 0;
       bucket = create_bucket();
       /* asph : sum of all head terms of spoly */
       for ( t = dm; t; t = NEXT(t) ) {
         _dp_sp_mod(ps[t->dp1],ps[t->dp2],m,&sp);
   /*      fprintf(stderr,"splen=%d-",dp_nt(sp)); */
         if ( sp ) {
           MKNODE(bt,sp,blist); blist = bt;
           add_bucket(bucket,dp_dllist(sp),nv);
   /*        fprintf(stderr,"%d-",length(s0)); */
         }
       }
   #if 0
       if ( DP_Print )
         fprintf(asir_out,"initial spmat : %d x %d ",length(blist),length(s0));
   #endif
       /* s0 : all the terms appeared in symbolic reduction */
       nred = 0;
       s0 = 0;
       while ( 1 ) {
         head = remove_head_bucket(bucket,nv);
         if ( !head ) break;
         else {
           NEXTNODE(s0,s);
           BDY(s) = (pointer)head;
         }
         for ( r = gall;  r; r = NEXT(r) )
           if ( _dl_redble(BDY(ps[(int)BDY(r)])->dl,head,nv) )
             break;
         if ( r ) {
           dltod(head,nv,&tdp);
           dp_subd(tdp,ps[(int)BDY(r)],&sd);
           dt = mul_dllist(BDY(sd)->dl,ps[(int)BDY(r)]);
           add_bucket(bucket,NEXT(dt),nv);
   /*        fprintf(stderr,"[%d]",length(dt)); */
           /* list of [t,f] */
           bt1 = mknode(2,BDY(sd)->dl,BDY(r));
           MKNODE(bt,bt1,blist); blist = bt;
   /*        fprintf(stderr,"%d-",length(s0));  */
           nred++;
         }
       }
       if ( s0 ) NEXT(s) = 0;
   /*    fprintf(stderr,"\n"); */
       get_eg(&tmp1); add_eg(&eg_symb,&tmp0,&tmp1);
       init_eg(&eg_split_symb); add_eg(&eg_split_symb,&tmp0,&tmp1);
   
                 /* head terms of reducers */      if ( DP_Print )
                 ht = (DL *)MALLOC(nred*sizeof(DL));        fprintf(asir_out,"number of reducers : %d\n",nred);
                 for ( r = blist, i = 0; i < nred; r = NEXT(r), i++ )  
                         ht[i] = BDY((DP)BDY(r))->dl;  
   
                 /* col = number of all terms */      get_eg(&tmp0);
                 for ( s = s0, col = 0; s; s = NEXT(s), col++ );      /* the first nred polys in blist are reducers */
       /* row = the number of all the polys */
       for ( r = blist, row = 0; r; r = NEXT(r), row++ );
   
                 /* head terms of all terms */      /* col = number of all terms */
                 at = (DL *)MALLOC(col*sizeof(DL));      for ( s = s0, col = 0; s; s = NEXT(s), col++ );
                 for ( s = s0, i = 0; i < col; s = NEXT(s), i++ )  
                         at[i] = (DL)BDY(s);  
   
                 /* store coefficients separately in spmat and redmat */      /* head terms of all terms */
                 nsp = row-nred;      at = (DL *)MALLOC(col*sizeof(DL));
       for ( s = s0, i = 0; i < col; s = NEXT(s), i++ )
         at[i] = (DL)BDY(s);
   
                 /* reducer matrix */      /* store coefficients separately in spmat and redmat */
                 redmat = (int **)almat(nred,col);      nsp = row-nred;
                 for ( i = 0, r = blist; i < nred; r = NEXT(r), i++ )  
                         _dpmod_to_vect(BDY(r),at,redmat[i]);  
                 /* XXX */  
 /*              reduce_reducers_mod(redmat,nred,col,m); */  
                 /* register the position of the head term */  
                 indred = (int *)MALLOC(nred*sizeof(int));  
                 bzero(indred,nred*sizeof(int));  
                 isred = (int *)MALLOC(col*sizeof(int));  
                 bzero(isred,col*sizeof(int));  
                 for ( i = 0; i < nred; i++ ) {  
                         ri = redmat[i];  
                         for ( j = 0; j < col && !ri[j]; j++ );  
                         indred[i] = j;  
                         isred[j] = 1;  
                 }  
   
                 spcol = col-nred;      /* reducer matrix */
                 /* head terms not in ht */      /* indred : register the position of the head term */
                 st = (DL *)MALLOC(spcol*sizeof(DL));      redmat = (CDP *)MALLOC(nred*sizeof(CDP));
                 for ( j = 0, k = 0; j < col; j++ )      for ( i = 0, r = blist; i < nred; r = NEXT(r), i++ )
                         if ( !isred[j] )        _tf_to_vect_compress(BDY(r),at,&redmat[i]);
                                 st[k++] = at[j];  
   
                 /* spoly matrix; stored in reduced form; terms in ht[] are omitted */      /* register the position of the head term */
                 spmat = almat(nsp,spcol);      indred = (int *)MALLOC_ATOMIC(nred*sizeof(int));
                 w = (int *)MALLOC(col*sizeof(int));      bzero(indred,nred*sizeof(int));
                 for ( ; i < row; r = NEXT(r), i++ ) {      isred = (int *)MALLOC_ATOMIC(col*sizeof(int));
                         bzero(w,col*sizeof(int));      bzero(isred,col*sizeof(int));
                         _dpmod_to_vect(BDY(r),at,w);      for ( i = 0; i < nred; i++ ) {
                         reduce_sp_by_red_mod(w,redmat,indred,nred,col,m);        ri = redmat[i];
                         for ( j = 0, k = 0; j < col; j++ )        indred[i] = ri->body[0];
                                 if ( !isred[j] )        isred[indred[i]] = 1;
                                         spmat[i-nred][k++] = w[j];      }
                 }  
   
                 get_eg(&tmp0); add_eg(&eg_elim1,&tmp1,&tmp0);      spcol = col-nred;
                 init_eg(&eg_split_elim1); add_eg(&eg_split_elim1,&tmp1,&tmp0);      /* head terms not in ht */
       st = (DL *)MALLOC(spcol*sizeof(DL));
       for ( j = 0, k = 0; j < col; j++ )
         if ( !isred[j] )
           st[k++] = at[j];
       get_eg(&tmp1); add_eg(&eg_conv,&tmp0,&tmp1);
       init_eg(&eg_split_conv); add_eg(&eg_split_conv,&tmp0,&tmp1);
   
                 colstat = (int *)MALLOC_ATOMIC(spcol*sizeof(int));      get_eg(&tmp1);
                 for ( i = 0, nonzero=0; i < nsp; i++ )      /* spoly matrix; stored in reduced form; terms in ht[] are omitted */
                         for ( j = 0; j < spcol; j++ )      spmat = (int **)MALLOC(nsp*sizeof(int *));
                                 if ( spmat[i][j] )      w = (int *)MALLOC_ATOMIC(col*sizeof(int));
                                         nonzero++;  
                 if ( DP_Print )  
                         fprintf(asir_out,"spmat : %d x %d (nonzero=%f%%)...",  
                                 nsp,spcol,((double)nonzero*100)/(nsp*spcol));  
                 rank = generic_gauss_elim_mod(spmat,nsp,spcol,m,colstat);  
   
                 get_eg(&tmp1); add_eg(&eg_elim2,&tmp0,&tmp1);      /* skip reducers in blist */
                 init_eg(&eg_split_elim2); add_eg(&eg_split_elim2,&tmp0,&tmp1);      for ( i = 0, r = blist; i < nred; r = NEXT(r), i++ );
       for ( i = 0; r; r = NEXT(r) ) {
         bzero(w,col*sizeof(int));
         _dpmod_to_vect(BDY(r),at,w);
         reduce_sp_by_red_mod_compress(w,redmat,indred,nred,col,m);
         for ( j = 0; j < col; j++ )
           if ( w[j] )
             break;
         if ( j < col ) {
           w1 = (int *)MALLOC_ATOMIC(spcol*sizeof(int));
           for ( j = 0, k = 0; j < col; j++ )
             if ( !isred[j] )
               w1[k++] = w[j];
           spmat[i] = w1;
           i++;
         }
       }
       /* update nsp */
       nsp0 = nsp;
       nsp = i;
   
                 if ( DP_Print ) {      /* XXX free redmat explicitly */
                         fprintf(asir_out,"done rank = %d\n",rank,row,col);      for ( k = 0; k < nred; k++ ) {
                         print_eg("Symb",&eg_split_symb);        GCFREE(BDY(redmat[k]));
                         print_eg("Elim1",&eg_split_elim1);        GCFREE(redmat[k]);
                         print_eg("Elim2",&eg_split_elim2);      }
                         fprintf(asir_out,"\n");  
                 }      get_eg(&tmp0); add_eg(&eg_elim1,&tmp1,&tmp0);
                 for ( j = 0, i = 0; j < spcol; j++ )      init_eg(&eg_split_elim1); add_eg(&eg_split_elim1,&tmp1,&tmp0);
                         if ( colstat[j] ) {  
                                 mp0 = 0;      colstat = (int *)MALLOC_ATOMIC(spcol*sizeof(int));
                                 NEXTMP(mp0,mp); mp->dl = st[j]; mp->c = STOI(1);      bzero(colstat,spcol*sizeof(int));
                                 for ( k = j+1; k < spcol; k++ )      for ( i = 0, nonzero=0; i < nsp; i++ )
                                         if ( !colstat[k] && spmat[i][k] ) {        for ( j = 0; j < spcol; j++ )
                                                 NEXTMP(mp0,mp); mp->dl = st[k];          if ( spmat[i][j] )
                                                 mp->c = STOI(spmat[i][k]);            nonzero++;
                                 }      if ( DP_Print && nsp )
                                 NEXT(mp) = 0;        fprintf(asir_out,"spmat : %d x %d (nonzero=%f%%)...",
                                 MKDP(nv,mp0,nf); nf->sugar = dm->sugar;          nsp,spcol,((double)nonzero*100)/(nsp*spcol));
                                 nh = newps_mod(nf,m);      if ( nsp )
                                 d = updpairs(d,g,nh);        rank = generic_gauss_elim_mod(spmat,nsp,spcol,m,colstat);
                                 g = updbase(g,nh);      else
                                 gall = append_one(gall,nh);        rank = 0;
                                 i++;      get_eg(&tmp1); add_eg(&eg_elim2,&tmp0,&tmp1);
                         }      init_eg(&eg_split_elim2); add_eg(&eg_split_elim2,&tmp0,&tmp1);
         }  
         if ( DP_Print ) {      if ( DP_Print ) {
                 print_eg("Symb",&eg_symb);        fprintf(asir_out,"done rank = %d\n",rank,row,col);
                 print_eg("Elim1",&eg_elim1);        print_eg("Symb",&eg_split_symb);
                 print_eg("Elim2",&eg_elim2);        print_eg("Conv",&eg_split_conv);
                 fflush(asir_out);        print_eg("Elim1",&eg_split_elim1);
         }        print_eg("Elim2",&eg_split_elim2);
         return g;        fprintf(asir_out,"\n");
       }
   
       NZR += rank;
       ZR += nsp0-rank;
   
       if ( !rank )
         continue;
   
       for ( j = 0, i = 0; j < spcol; j++ )
         if ( colstat[j] ) {
           mp0 = 0;
           NEXTMP(mp0,mp); mp->dl = st[j]; mp->c = (Obj)STOI(1);
           for ( k = j+1; k < spcol; k++ )
             if ( !colstat[k] && spmat[i][k] ) {
               NEXTMP(mp0,mp); mp->dl = st[k];
               mp->c = (Obj)STOI(spmat[i][k]);
           }
           NEXT(mp) = 0;
           MKDP(nv,mp0,nf); nf->sugar = dm->sugar;
           nh = newps_mod(nf,m);
           if ( nh == pscalen ) {
             psca = (unsigned int **)
               REALLOC(psca,2*pscalen*sizeof(unsigned int *));
             pscalen *= 2;
           }
           dptoca(ps[nh],&psca[nh]);
           d = updpairs(d,g,nh);
           g = updbase(g,nh);
           gall = append_one(gall,nh);
           i++;
         }
   
       /* XXX free spmat[] explicitly */
       for ( j = 0; j < nsp; j++ ) {
         GCFREE(spmat[j]);
       }
     }
     if ( DP_Print ) {
       print_eg("Symb",&eg_symb);
       print_eg("Conv",&eg_conv);
       print_eg("Elim1",&eg_elim1);
       print_eg("Elim2",&eg_elim2);
       fflush(asir_out);
     }
     return g;
 }  }
   
 int DPPlength(n)  NODE gb_f4_mod_old(NODE f,int m)
 DP_pairs n;  
 {  {
         int i;    int i,j,k,nh,row,col,nv;
     NODE r,g,gall;
     NODE s,s0;
     DP_pairs d,dm,dr,t;
     DP nf,f2,sp,sd,sdm,tdp;
     MP mp,mp0;
     NODE blist,bt;
     DL *ht,*at,*st;
     int **spmat,**redmat;
     int *colstat,*w;
     int rank,nred,nsp,nonzero,spcol;
     int *indred,*isred,*ri;
     struct oEGT tmp0,tmp1,eg_split_symb,eg_split_elim1,eg_split_elim2;
     extern struct oEGT eg_symb,eg_elim1,eg_elim2;
   
         for ( i = 0; n; n = NEXT(n), i++ );    init_eg(&eg_symb); init_eg(&eg_elim1); init_eg(&eg_elim2);
         return i;    for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {
       i = (int)BDY(r);
       d = updpairs(d,g,i);
       g = updbase(g,i);
       gall = append_one(gall,i);
     }
     if ( gall )
       nv = ((DP)ps[(int)BDY(gall)])->nv;
     while ( d ) {
       get_eg(&tmp0);
       minsugar(d,&dm,&dr); d = dr;
       if ( DP_Print )
         fprintf(asir_out,"sugar=%d\n",dm->sugar);
       blist = 0; s0 = 0;
       /* asph : sum of all head terms of spoly */
       for ( t = dm; t; t = NEXT(t) ) {
         _dp_sp_mod(ps[t->dp1],ps[t->dp2],m,&sp);
         if ( sp ) {
           MKNODE(bt,sp,blist); blist = bt;
           s0 = symb_merge(s0,dp_dllist(sp),nv);
         }
       }
       /* s0 : all the terms appeared in symbolic redunction */
       for ( s = s0, nred = 0; s; s = NEXT(s) ) {
         for ( r = gall;  r; r = NEXT(r) )
           if ( _dl_redble(BDY(ps[(int)BDY(r)])->dl,BDY(s),nv) )
             break;
         if ( r ) {
           dltod(BDY(s),nv,&tdp);
           dp_subd(tdp,ps[(int)BDY(r)],&sd);
           _dp_mod(sd,m,0,&sdm);
           mulmd_dup(m,sdm,ps[(int)BDY(r)],&f2);
           MKNODE(bt,f2,blist); blist = bt;
           s = symb_merge(s,dp_dllist(f2),nv);
           nred++;
         }
       }
   
       get_eg(&tmp1); add_eg(&eg_symb,&tmp0,&tmp1);
       init_eg(&eg_split_symb); add_eg(&eg_split_symb,&tmp0,&tmp1);
   
       /* the first nred polys in blist are reducers */
       /* row = the number of all the polys */
       for ( r = blist, row = 0; r; r = NEXT(r), row++ );
   
       /* head terms of reducers */
       ht = (DL *)MALLOC(nred*sizeof(DL));
       for ( r = blist, i = 0; i < nred; r = NEXT(r), i++ )
         ht[i] = BDY((DP)BDY(r))->dl;
   
       /* col = number of all terms */
       for ( s = s0, col = 0; s; s = NEXT(s), col++ );
   
       /* head terms of all terms */
       at = (DL *)MALLOC(col*sizeof(DL));
       for ( s = s0, i = 0; i < col; s = NEXT(s), i++ )
         at[i] = (DL)BDY(s);
   
       /* store coefficients separately in spmat and redmat */
       nsp = row-nred;
   
       /* reducer matrix */
       redmat = (int **)almat(nred,col);
       for ( i = 0, r = blist; i < nred; r = NEXT(r), i++ )
         _dpmod_to_vect(BDY(r),at,redmat[i]);
       /* XXX */
   /*    reduce_reducers_mod(redmat,nred,col,m); */
       /* register the position of the head term */
       indred = (int *)MALLOC(nred*sizeof(int));
       bzero(indred,nred*sizeof(int));
       isred = (int *)MALLOC(col*sizeof(int));
       bzero(isred,col*sizeof(int));
       for ( i = 0; i < nred; i++ ) {
         ri = redmat[i];
         for ( j = 0; j < col && !ri[j]; j++ );
         indred[i] = j;
         isred[j] = 1;
       }
   
       spcol = col-nred;
       /* head terms not in ht */
       st = (DL *)MALLOC(spcol*sizeof(DL));
       for ( j = 0, k = 0; j < col; j++ )
         if ( !isred[j] )
           st[k++] = at[j];
   
       /* spoly matrix; stored in reduced form; terms in ht[] are omitted */
       spmat = almat(nsp,spcol);
       w = (int *)MALLOC(col*sizeof(int));
       for ( ; i < row; r = NEXT(r), i++ ) {
         bzero(w,col*sizeof(int));
         _dpmod_to_vect(BDY(r),at,w);
         reduce_sp_by_red_mod(w,redmat,indred,nred,col,m);
         for ( j = 0, k = 0; j < col; j++ )
           if ( !isred[j] )
             spmat[i-nred][k++] = w[j];
       }
   
       get_eg(&tmp0); add_eg(&eg_elim1,&tmp1,&tmp0);
       init_eg(&eg_split_elim1); add_eg(&eg_split_elim1,&tmp1,&tmp0);
   
       colstat = (int *)MALLOC_ATOMIC(spcol*sizeof(int));
       for ( i = 0, nonzero=0; i < nsp; i++ )
         for ( j = 0; j < spcol; j++ )
           if ( spmat[i][j] )
             nonzero++;
       if ( DP_Print && nsp )
         fprintf(asir_out,"spmat : %d x %d (nonzero=%f%%)...",
           nsp,spcol,((double)nonzero*100)/(nsp*spcol));
       if ( nsp )
         rank = generic_gauss_elim_mod(spmat,nsp,spcol,m,colstat);
       else
         rank = 0;
       get_eg(&tmp1); add_eg(&eg_elim2,&tmp0,&tmp1);
       init_eg(&eg_split_elim2); add_eg(&eg_split_elim2,&tmp0,&tmp1);
   
       if ( DP_Print ) {
         fprintf(asir_out,"done rank = %d\n",rank,row,col);
         print_eg("Symb",&eg_split_symb);
         print_eg("Elim1",&eg_split_elim1);
         print_eg("Elim2",&eg_split_elim2);
         fprintf(asir_out,"\n");
       }
       for ( j = 0, i = 0; j < spcol; j++ )
         if ( colstat[j] ) {
           mp0 = 0;
           NEXTMP(mp0,mp); mp->dl = st[j]; mp->c = (Obj)STOI(1);
           for ( k = j+1; k < spcol; k++ )
             if ( !colstat[k] && spmat[i][k] ) {
               NEXTMP(mp0,mp); mp->dl = st[k];
               mp->c = (Obj)STOI(spmat[i][k]);
           }
           NEXT(mp) = 0;
           MKDP(nv,mp0,nf); nf->sugar = dm->sugar;
           nh = newps_mod(nf,m);
           d = updpairs(d,g,nh);
           g = updbase(g,nh);
           gall = append_one(gall,nh);
           i++;
         }
     }
     if ( DP_Print ) {
       print_eg("Symb",&eg_symb);
       print_eg("Elim1",&eg_elim1);
       print_eg("Elim2",&eg_elim2);
       fflush(asir_out);
     }
     return g;
 }  }
   
 void printdl(dl)  int DPPlength(DP_pairs n)
 DL dl;  
 {  {
         int i;    int i;
   
         fprintf(asir_out,"<<");    for ( i = 0; n; n = NEXT(n), i++ );
         for ( i = 0; i < CNVars-1; i++ )    return i;
                 fprintf(asir_out,"%d,",dl->d[i]);  
         fprintf(asir_out,"%d>>",dl->d[i]);  
 }  }
   
 void pltovl(l,vl)  void printdl(DL dl)
 LIST l;  
 VL *vl;  
 {  {
         NODE n;    int i;
         VL r,r0;  
   
         n = BDY(l);    fprintf(asir_out,"<<");
         for ( r0 = 0; n; n = NEXT(n) ) {    for ( i = 0; i < CNVars-1; i++ )
                 NEXTVL(r0,r); r->v = VR((P)BDY(n));      fprintf(asir_out,"%d,",dl->d[i]);
         }    fprintf(asir_out,"%d>>",dl->d[i]);
         if ( r0 ) NEXT(r) = 0;  
         *vl = r0;  
 }  }
   
 void makesubst(v,s)  void pltovl(LIST l,VL *vl)
 VL v;  
 NODE *s;  
 {  {
         NODE r,r0;    NODE n;
         Q q;    VL r,r0;
         unsigned int n;  
   
         for ( r0 = 0; v; v = NEXT(v) ) {    n = BDY(l);
                 NEXTNODE(r0,r); BDY(r) = (pointer)v->v;    for ( r0 = 0; n; n = NEXT(n) ) {
       NEXTVL(r0,r); r->v = VR((P)BDY(n));
     }
     if ( r0 ) NEXT(r) = 0;
     *vl = r0;
   }
   
   void vltopl(VL vl,LIST *l)
   {
     VL n;
     NODE r,r0;
     P p;
   
     n = vl;
     for ( r0 = 0; n; n = NEXT(n) ) {
       NEXTNODE(r0,r); MKV(n->v,p); BDY(r) = (pointer)p;
     }
     if ( r0 ) NEXT(r) = 0;
     MKLIST(*l,r0);
   }
   
   void makesubst(VL v,NODE *s)
   {
     NODE r,r0;
     Q q;
     unsigned int n;
   
     for ( r0 = 0; v; v = NEXT(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; UTOQ(n,q);
 #else  #else
                 n = random(); UTOQ(n,q);      n = random(); UTOQ(n,q);
 #endif  #endif
                 NEXTNODE(r0,r); BDY(r) = (pointer)q;      NEXTNODE(r0,r); BDY(r) = (pointer)q;
         }    }
         if ( r0 ) NEXT(r) = 0;    if ( r0 ) NEXT(r) = 0;
         *s = r0;    *s = r0;
 }  }
   
 void printsubst(s)  void printsubst(NODE s)
 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",QTOS((Q)BDY(s)));
                 if ( NEXT(s) ) {      if ( NEXT(s) ) {
                         fputc(',',asir_out); s = NEXT(s);        fputc(',',asir_out); s = NEXT(s);
                 } else      } else
                         break;        break;
         }    }
         fprintf(asir_out,"]\n"); return;    fprintf(asir_out,"]\n"); return;
 }  }
   
 void vlminus(v,w,d)  void vlminus(VL v,VL w,VL *d)
 VL v,w,*d;  
 {  {
         int i,j,n,m;    int i,j,n,m;
         V *va,*wa;    V *va,*wa;
         V a;    V a;
         VL r,r0;    VL r,r0;
         VL t;    VL t;
   
         for ( n = 0, t = v; t; t = NEXT(t), n++ );    for ( n = 0, t = v; t; t = NEXT(t), n++ );
         va = (V *)ALLOCA(n*sizeof(V));    va = (V *)ALLOCA(n*sizeof(V));
         for ( i = 0, t = v; t; t = NEXT(t), i++ )    for ( i = 0, t = v; t; t = NEXT(t), i++ )
                 va[i] = t->v;      va[i] = t->v;
         for ( m = 0, t = w; t; t = NEXT(t), m++ );    for ( m = 0, t = w; t; t = NEXT(t), m++ );
         wa = (V *)ALLOCA(m*sizeof(V));    wa = (V *)ALLOCA(m*sizeof(V));
         for ( i = 0, t = w; t; t = NEXT(t), i++ )    for ( i = 0, t = w; t; t = NEXT(t), i++ )
                 wa[i] = t->v;      wa[i] = t->v;
         for ( i = 0; i < n; i++ ) {    for ( i = 0; i < n; i++ ) {
                 a = va[i];      a = va[i];
                 for ( j = 0; j < m; j++ )      for ( j = 0; j < m; j++ )
                         if ( a == wa[j] )        if ( a == wa[j] )
                                 break;          break;
                 if ( j < m )      if ( j < m )
                         va[i] = 0;        va[i] = 0;
         }    }
         for ( r0 = 0, i = 0; i < n; i++ )    for ( r0 = 0, i = 0; i < n; i++ )
                 if ( va[i] ) { NEXTVL(r0,r); r->v = va[i];      }      if ( va[i] ) { NEXTVL(r0,r); r->v = va[i];  }
         if ( r0 ) NEXT(r) = 0;    if ( r0 ) NEXT(r) = 0;
         *d = r0;    *d = r0;
 }  }
   
 int validhc(a,m,s)  int validhc(P a,int m,NODE s)
 P a;  
 int m;  
 NODE s;  
 {  {
         P c,c1;    P c,c1;
         V v;    V v;
   
         if ( !a )    if ( !a )
                 return 0;      return 0;
         for ( c = a; s; s = NEXT(s) ) {    for ( c = a; s; s = NEXT(s) ) {
                 v = (V)BDY(s); s = NEXT(s);      v = (V)BDY(s); s = NEXT(s);
                 substp(CO,c,v,(P)BDY(s),&c1); c = c1;      substp(CO,c,v,(P)BDY(s),&c1); c = c1;
         }    }
         ptomp(m,c,&c1);    ptomp(m,c,&c1);
         return c1 ? 1 : 0;    return c1 ? 1 : 0;
 }  }
   
 void setup_arrays(f,m,r)  void setup_arrays(NODE f,int m,NODE *r)
 NODE f,*r;  
 int m;  
 {  {
         int i;    int i;
         NODE s,s0,f0;    NODE s,s0,f0;
   
         f0 = f = NODE_sortb(f,1);  #if 1
         psn = length(f); pslen = 2*psn;    f0 = f = NODE_sortb(f,1);
         ps = (DP *)MALLOC(pslen*sizeof(DP));  #else
         psh = (DL *)MALLOC(pslen*sizeof(DL));    f0 = f;
         pss = (int *)MALLOC(pslen*sizeof(int));  #endif
         psc = (P *)MALLOC(pslen*sizeof(P));    psn = length(f); pslen = 2*psn;
         for ( i = 0; i < psn; i++, f = NEXT(f) ) {    ps = (DP *)MALLOC(pslen*sizeof(DP));
                 prim_part((DP)BDY(f),m,&ps[i]);    psh = (DL *)MALLOC(pslen*sizeof(DL));
                 if ( Demand )    pss = (int *)MALLOC(pslen*sizeof(int));
                         dp_save(i,(Obj)ps[i],0);    psc = (P *)MALLOC(pslen*sizeof(P));
                 psh[i] = BDY(ps[i])->dl;    for ( i = 0; i < psn; i++, f = NEXT(f) ) {
                 pss[i] = ps[i]->sugar;      prim_part((DP)BDY(f),m,&ps[i]);
                 psc[i] = BDY(ps[i])->c;      if ( Demand )
         }        dp_save(i,(Obj)ps[i],0);
         if ( GenTrace && (OXCheck >= 0) ) {      psh[i] = BDY(ps[i])->dl;
                 Q q;      pss[i] = ps[i]->sugar;
                 STRING fname;      psc[i] = (P)BDY(ps[i])->c;
                 LIST input;    }
                 NODE arg;    if ( GenTrace ) {
                 Obj dmy;      Q q;
       STRING fname;
       LIST input;
       NODE arg,t,t1;
       Obj obj;
   
       t = 0;
       for ( i = psn-1; i >= 0; i-- ) {
         MKNODE(t1,ps[i],t);
         t = t1;
       }
       MKLIST(input,t);
   
                 STOQ(OXCheck,q);      if ( OXCheck >= 0 ) {
                 MKSTR(fname,"register_input");        STOQ(OXCheck,q);
                 MKLIST(input,f0);        MKSTR(fname,"register_input");
                 arg = mknode(3,q,fname,input);        arg = mknode(3,q,fname,input);
                 Pox_cmo_rpc(arg,&dmy);        Pox_cmo_rpc(arg,&obj);
         }      } else if ( OXCheck < 0 ) {
         for ( s0 = 0, i = 0; i < psn; i++ ) {        MKNODE(AllTraceList,input,0);
                 NEXTNODE(s0,s); BDY(s) = (pointer)i;      }
         }    }
         if ( s0 ) NEXT(s) = 0;    for ( s0 = 0, i = 0; i < psn; i++ ) {
         *r = s0;      NEXTNODE(s0,s); BDY(s) = (pointer)i;
     }
     if ( s0 ) NEXT(s) = 0;
     *r = s0;
 }  }
   
 void prim_part(f,m,r)  void prim_part(DP f,int m,DP *r)
 DP f,*r;  
 int m;  
 {  {
         P d,t;    P d,t;
   
         if ( m > 0 ) {    if ( m > 0 ) {
                 if ( PCoeffs )      if ( PCoeffs )
                         dp_prim_mod(f,m,r);        dp_prim_mod(f,m,r);
                 else      else
                         *r = f;        _dp_monic(f,m,r);
         } else {    } else {
                 if ( dp_fcoeffs )      if ( dp_fcoeffs || PCoeffs )
                         *r = f;        dp_prim(f,r);
                 else if ( PCoeffs )      else
                         dp_prim(f,r);        dp_ptozp(f,r);
                 else      if ( GenTrace && TraceList ) {
                         dp_ptozp(f,r);        /* adust the denominator according to the final
                 if ( GenTrace && TraceList ) {           content reduction */
                         divsp(CO,BDY(f)->c,BDY(*r)->c,&d);        divsp(CO,(P)BDY(f)->c,(P)BDY(*r)->c,&d);
                         mulp(CO,(P)ARG3(BDY((LIST)BDY(TraceList))),d,&t);        mulp(CO,(P)ARG3(BDY((LIST)BDY(TraceList))),d,&t);
                         ARG3(BDY((LIST)BDY(TraceList))) = t;        ARG3(BDY((LIST)BDY(TraceList))) = t;
                 }      }
         }    }
 }  }
   
 NODE /* of DP */ NODE_sortb_insert( newdp, nd, dec )  NODE /* of DP */ NODE_sortb_insert( DP newdp, NODE /* of DP */ nd, int dec )
 DP newdp;  
 NODE /* of DP */ nd;  
 int dec;  
 {  {
         register NODE last, p;    register NODE last, p;
         register DL newdl = BDY(newdp)->dl;    register DL newdl = BDY(newdp)->dl;
         register int (*cmpfun)() = cmpdl, nv = CNVars;    register int (*cmpfun)() = cmpdl, nv = CNVars;
         NODE newnd;    NODE newnd;
         int sgn = dec ? 1 : -1;    int sgn = dec ? 1 : -1;
         MKNODE( newnd, newdp, 0 );    MKNODE( newnd, newdp, 0 );
         if ( !(last = nd) || sgn*(*cmpfun)( nv, newdl, BDY((DP) BDY(last))->dl ) > 0 ) {    if ( !(last = nd) || sgn*(*cmpfun)( nv, newdl, BDY((DP) BDY(last))->dl ) > 0 ) {
                 NEXT(newnd) = last;      NEXT(newnd) = last;
                 return newnd;      return newnd;
         }    }
         for ( ; p = NEXT(last); last = p )    for ( ; p = NEXT(last); last = p )
                 if ( sgn*(*cmpfun)( nv, newdl, BDY((DP) BDY(p))->dl ) > 0 ) break;      if ( sgn*(*cmpfun)( nv, newdl, BDY((DP) BDY(p))->dl ) > 0 ) break;
         if ( p ) NEXT(NEXT(last) = newnd) = p;    if ( p ) NEXT(NEXT(last) = newnd) = p;
         else NEXT(last) = newnd;    else NEXT(last) = newnd;
         return nd;    return nd;
 }  }
   
 NODE NODE_sortb( node, dec )  NODE NODE_sortb( NODE node, int dec )
 NODE node;  
 int dec;  
 {  {
         register NODE nd, ans;    register NODE nd, ans;
   
         for ( ans = 0, nd = node; nd; nd = NEXT(nd) )    for ( ans = 0, nd = node; nd; nd = NEXT(nd) )
                 ans = NODE_sortb_insert( (DP) BDY(nd), ans, dec );      ans = NODE_sortb_insert( (DP) BDY(nd), ans, dec );
         return ans;    return ans;
 }  }
   
 NODE /* of index */ NODE_sortbi_insert( newdpi, nd, dec )  NODE /* of index */ NODE_sortbi_insert( int newdpi, NODE /* of index */ nd, int dec )
 int newdpi;  
 NODE /* of index */ nd;  
 int dec;  
 {  {
         register NODE last, p;    register NODE last, p;
         register DL newdl = psh[newdpi];    register DL newdl = psh[newdpi];
         register int (*cmpfun)() = cmpdl, nv = CNVars;    register int (*cmpfun)() = cmpdl, nv = CNVars;
         NODE newnd;    NODE newnd;
         int sgn = dec ? 1 : -1;    int sgn = dec ? 1 : -1;
         MKNODE( newnd, newdpi, 0 );    MKNODE( newnd, newdpi, 0 );
         if ( !(last = nd) || sgn*(*cmpfun)( nv, newdl, psh[(int)BDY(last)] ) > 0 ) {    if ( !(last = nd) || sgn*(*cmpfun)( nv, newdl, psh[(int)BDY(last)] ) > 0 ) {
                 NEXT(newnd) = last;      NEXT(newnd) = last;
                 return newnd;      return newnd;
         }    }
         for ( ; p = NEXT(last); last = p )    for ( ; p = NEXT(last); last = p )
                 if ( sgn*(*cmpfun)( nv, newdl, psh[(int)BDY(p)] ) > 0 ) break;      if ( sgn*(*cmpfun)( nv, newdl, psh[(int)BDY(p)] ) > 0 ) break;
         if ( p ) NEXT(NEXT(last) = newnd) = p;    if ( p ) NEXT(NEXT(last) = newnd) = p;
         else NEXT(last) = newnd;    else NEXT(last) = newnd;
         return nd;    return nd;
 }  }
   
 NODE NODE_sortbi( node, dec )  NODE NODE_sortbi( NODE node, int dec )
 NODE node;  
 int dec;  
 {  {
         register NODE nd, ans;    register NODE nd, ans;
   
         for ( ans = 0, nd = node; nd; nd = NEXT(nd) )    for ( ans = 0, nd = node; nd; nd = NEXT(nd) )
                 ans = NODE_sortbi_insert( (int) BDY(nd), ans, dec );      ans = NODE_sortbi_insert( (int) BDY(nd), ans, dec );
         return ans;    return ans;
 }  }
   
 void reduceall(in,h)  void reduceall(NODE in,NODE *h)
 NODE in;  
 NODE *h;  
 {  {
         NODE r,t,top;    NODE r,t,top;
         int n,i,j;    int n,i,j;
         int *w;    int *w;
         DP g,g1;    DP g,g1;
         struct oEGT tmp0,tmp1;    struct oEGT tmp0,tmp1;
   
         if ( NoRA ) {    if ( NoRA ) {
                 *h = in; return;      *h = in; return;
         }    }
         if ( DP_Print || DP_PrintShort ) {    if ( DP_Print || DP_PrintShort ) {
                 fprintf(asir_out,"reduceall\n"); fflush(asir_out);      fprintf(asir_out,"reduceall\n"); fflush(asir_out);
         }    }
         r = NODE_sortbi(in,0);    r = NODE_sortbi(in,0);
         n = length(r);    n = length(r);
         w = (int *)ALLOCA(n*sizeof(int));    w = (int *)ALLOCA(n*sizeof(int));
         for ( i = 0, t = r; i < n; i++, t = NEXT(t) )    for ( i = 0, t = r; i < n; i++, t = NEXT(t) )
                 w[i] = (int)BDY(t);      w[i] = (int)BDY(t);
         for ( i = 0; i < n; i++ ) {    /* w[i] < 0 : reduced to 0 */
                 for ( top = 0, j = n-1; j >= 0; j-- )    for ( i = 0; i < n; i++ ) {
                         if ( j != i ) {      for ( top = 0, j = n-1; j >= 0; j-- )
                                 MKNODE(t,(pointer)w[j],top); top = t;        if ( w[j] >= 0 && j != i ) {
                         }          MKNODE(t,(pointer)w[j],top); top = t;
                 get_eg(&tmp0);        }
                 dp_load(w[i],&ps[w[i]]);      get_eg(&tmp0);
       dp_load(w[i],&ps[w[i]]);
   
                 if ( GenTrace ) {      if ( GenTrace ) {
                         Q q;        Q q;
                         NODE node;        NODE node;
                         LIST hist;        LIST hist;
   
                         STOQ(w[i],q);        STOQ(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);
                 }      }
                 if ( !PtozpRA || !Multiple )      _dp_nf(top,ps[w[i]],ps,1,&g);
                         _dp_nf(top,ps[w[i]],ps,1,&g);      prim_part(g,0,&g1);
                 else      get_eg(&tmp1); add_eg(&eg_ra,&tmp0,&tmp1);
 #if MPI      if ( DP_Print || DP_PrintShort ) {
                         _dp_nf_ptozp_mpi(top,ps[w[i]],ps,1,Multiple,&g);        fprintf(asir_out,"."); fflush(asir_out);
 #else      }
                         _dp_nf_ptozp(top,ps[w[i]],ps,1,Multiple,&g);      if ( g1 ) {
 #endif        w[i] = newps(g1,0,(NODE)0);
                 prim_part(g,0,&g1);      } else {
                 get_eg(&tmp1); add_eg(&eg_ra,&tmp0,&tmp1);        w[i] = -1;
                 if ( DP_Print || DP_PrintShort ) {      }
                         fprintf(asir_out,"."); fflush(asir_out);    }
                 }    for ( top = 0, j = n-1; j >= 0; j-- ) {
                 w[i] = newps(g1,0,(NODE)0);      if ( w[j] >= 0 ) {
         }        MKNODE(t,(pointer)w[j],top); top = t;
         for ( top = 0, j = n-1; j >= 0; j-- ) {      }
                 MKNODE(t,(pointer)w[j],top); top = t;    }
         }    *h = top;
         *h = top;    if ( DP_Print || DP_PrintShort )
         if ( DP_Print || DP_PrintShort )      fprintf(asir_out,"\n");
                 fprintf(asir_out,"\n");  
 }  }
   
 void reduceall_mod(in,m,h)  void reduceall_mod(NODE in,int m,NODE *h)
 NODE in;  
 int m;  
 NODE *h;  
 {  {
         NODE r,t,top;    NODE r,t,top;
         int n,i,j;    int n,i,j;
         int *w;    int *w;
         DP g,p;    DP g,p;
         struct oEGT tmp0,tmp1;    struct oEGT tmp0,tmp1;
   
         if ( NoRA ) {    if ( NoRA ) {
                 *h = in; return;      *h = in; return;
         }    }
         if ( DP_Print || DP_PrintShort ) {    if ( DP_Print || DP_PrintShort ) {
                 fprintf(asir_out,"reduceall\n"); fflush(asir_out);      fprintf(asir_out,"reduceall\n"); fflush(asir_out);
         }    }
         r = NODE_sortbi(in,0);    r = NODE_sortbi(in,0);
         n = length(r);    n = length(r);
         w = (int *)ALLOCA(n*sizeof(int));    w = (int *)ALLOCA(n*sizeof(int));
         for ( i = 0, t = r; i < n; i++, t = NEXT(t) )    for ( i = 0, t = r; i < n; i++, t = NEXT(t) )
                 w[i] = (int)BDY(t);      w[i] = (int)BDY(t);
         for ( i = 0; i < n; i++ ) {    /* w[i] < 0 : reduced to 0 */
                 for ( top = 0, j = n-1; j >= 0; j-- )    for ( i = 0; i < n; i++ ) {
                         if ( j != i ) {      for ( top = 0, j = n-1; j >= 0; j-- )
                                 MKNODE(t,(pointer)w[j],top); top = t;        if ( w[j] >= 0 && j != i ) {
                         }          MKNODE(t,(pointer)w[j],top); top = t;
                 get_eg(&tmp0);        }
                 if ( PCoeffs )      get_eg(&tmp0);
                         dp_nf_mod(top,ps[w[i]],ps,m,1,&g);      if ( PCoeffs )
                 else {        dp_nf_mod(top,ps[w[i]],ps,m,1,&g);
                         dpto_dp(ps[w[i]],&p);      else {
                         _dp_nf_mod_destructive(top,p,ps,m,1,&g);        dpto_dp(ps[w[i]],&p);
                 }        _dp_nf_mod_destructive(top,p,ps,m,1,&g);
                 get_eg(&tmp1); add_eg(&eg_ra,&tmp0,&tmp1);      }
                 if ( DP_Print || DP_PrintShort ) {      get_eg(&tmp1); add_eg(&eg_ra,&tmp0,&tmp1);
                         fprintf(asir_out,"."); fflush(asir_out);      if ( DP_Print || DP_PrintShort ) {
                 }        fprintf(asir_out,"."); fflush(asir_out);
                 w[i] = newps_mod(g,m);      }
         }      if ( g ) {
         for ( top = 0, j = n-1; j >= 0; j-- ) {        w[i] = newps_mod(g,m);
                 MKNODE(t,(pointer)w[j],top); top = t;      } else {
         }        w[i] = -1;
         *h = top;      }
         if ( DP_Print || DP_PrintShort )    }
                 fprintf(asir_out,"\n");    for ( top = 0, j = n-1; j >= 0; j-- ) {
       if ( w[j] >= 0 ) {
         MKNODE(t,(pointer)w[j],top); top = t;
       }
     }
     *h = top;
     if ( DP_Print || DP_PrintShort )
       fprintf(asir_out,"\n");
 }  }
   
 int newps(a,m,subst)  int newps(DP a,int m,NODE subst)
 DP a;  
 int m;  
 NODE subst;  
 {  {
         if ( m && !validhc(!a?0:BDY(a)->c,m,subst) )    if ( m && !validhc(!a?0:(P)BDY(a)->c,m,subst) )
                 return -1;      return -1;
         if ( psn == pslen ) {    if ( psn == pslen ) {
                 pslen *= 2;      pslen *= 2;
                 ps = (DP *)REALLOC((char *)ps,pslen*sizeof(DP));      ps = (DP *)REALLOC((char *)ps,pslen*sizeof(DP));
                 psh = (DL *)REALLOC((char *)psh,pslen*sizeof(DL));      psh = (DL *)REALLOC((char *)psh,pslen*sizeof(DL));
                 pss = (int *)REALLOC((char *)pss,pslen*sizeof(int));      pss = (int *)REALLOC((char *)pss,pslen*sizeof(int));
                 psc = (P *)REALLOC((char *)psc,pslen*sizeof(P));      psc = (P *)REALLOC((char *)psc,pslen*sizeof(P));
                 if ( m )      if ( m )
                         psm = (DP *)REALLOC((char *)psm,pslen*sizeof(DP));        psm = (DP *)REALLOC((char *)psm,pslen*sizeof(DP));
         }    }
         if ( Demand ) {    if ( Demand ) {
                 if ( doing_f4 )      if ( doing_f4 )
                         ps[psn] = a;        ps[psn] = a;
                 else      else
                         ps[psn] = 0;        ps[psn] = 0;
                 dp_save(psn,(Obj)a,0);      dp_save(psn,(Obj)a,0);
         } else    } else
                 ps[psn] = a;      ps[psn] = a;
         psh[psn] = BDY(a)->dl;    psh[psn] = BDY(a)->dl;
         pss[psn] = a->sugar;    pss[psn] = a->sugar;
         psc[psn] = BDY(a)->c;    psc[psn] = (P)BDY(a)->c;
         if ( m )    if ( m )
                 _dp_mod(a,m,subst,&psm[psn]);      _dp_mod(a,m,subst,&psm[psn]);
         if ( GenTrace ) {    if ( GenTrace ) {
                 NODE tn,tr,tr1;      NODE tn,tr,tr1;
                 LIST trace;      LIST trace,trace1;
       NODE arg;
                 /* reverse the TraceList */      Q q1,q2;
                 tn = TraceList;      STRING fname;
                 for ( tr = 0; tn; tn = NEXT(tn) ) {      Obj obj;
                         MKNODE(tr1,BDY(tn),tr); tr = tr1;  
                 }      /* reverse the TraceList */
                 MKLIST(trace,tr);      tn = TraceList;
                 if ( OXCheck >= 0 ) {      for ( tr = 0; tn; tn = NEXT(tn) ) {
                         NODE arg;        MKNODE(tr1,BDY(tn),tr); tr = tr1;
                         Q q1,q2;      }
                         STRING fname;      MKLIST(trace,tr);
                         Obj dmy;      if ( OXCheck >= 0 ) {
         STOQ(OXCheck,q1);
                         STOQ(OXCheck,q1);        MKSTR(fname,"check_trace");
                         MKSTR(fname,"check_trace");        STOQ(psn,q2);
                         STOQ(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,&dmy);      } else if ( OXCheck < 0 ) {
                 } else        STOQ(psn,q1);
                         dp_save(psn,(Obj)trace,"t");        tn = mknode(2,q1,trace);
                 TraceList = 0;        MKLIST(trace1,tn);
         }        MKNODE(tr,trace1,AllTraceList);
         return psn++;        AllTraceList = tr;
       } else
         dp_save(psn,(Obj)trace,"t");
       TraceList = 0;
     }
     return psn++;
 }  }
   
 int newps_nosave(a,m,subst)  int newps_nosave(DP a,int m,NODE subst)
 DP a;  
 int m;  
 NODE subst;  
 {  {
         if ( m && !validhc(!a?0:BDY(a)->c,m,subst) )    if ( m && !validhc(!a?0:(P)BDY(a)->c,m,subst) )
                 return -1;      return -1;
         if ( psn == pslen ) {    if ( psn == pslen ) {
                 pslen *= 2;      pslen *= 2;
                 ps = (DP *)REALLOC((char *)ps,pslen*sizeof(DP));      ps = (DP *)REALLOC((char *)ps,pslen*sizeof(DP));
                 psh = (DL *)REALLOC((char *)psh,pslen*sizeof(DL));      psh = (DL *)REALLOC((char *)psh,pslen*sizeof(DL));
                 pss = (int *)REALLOC((char *)pss,pslen*sizeof(int));      pss = (int *)REALLOC((char *)pss,pslen*sizeof(int));
                 psc = (P *)REALLOC((char *)psc,pslen*sizeof(P));      psc = (P *)REALLOC((char *)psc,pslen*sizeof(P));
                 if ( m )      if ( m )
                         psm = (DP *)REALLOC((char *)psm,pslen*sizeof(DP));        psm = (DP *)REALLOC((char *)psm,pslen*sizeof(DP));
         }    }
         ps[psn] = 0;    ps[psn] = 0;
         psh[psn] = BDY(a)->dl;    psh[psn] = BDY(a)->dl;
         pss[psn] = a->sugar;    pss[psn] = a->sugar;
         psc[psn] = BDY(a)->c;    psc[psn] = (P)BDY(a)->c;
         if ( m )    if ( m )
                 _dp_mod(a,m,subst,&psm[psn]);      _dp_mod(a,m,subst,&psm[psn]);
         return psn++;    return psn++;
 }  }
   
 int newps_mod(a,m)  int newps_mod(DP a,int m)
 DP a;  
 int m;  
 {  {
         if ( psn == pslen ) {    if ( psn == pslen ) {
                 pslen *= 2;      pslen *= 2;
                 ps = (DP *)REALLOC((char *)ps,pslen*sizeof(DP));      ps = (DP *)REALLOC((char *)ps,pslen*sizeof(DP));
                 psh = (DL *)REALLOC((char *)psh,pslen*sizeof(DL));      psh = (DL *)REALLOC((char *)psh,pslen*sizeof(DL));
                 pss = (int *)REALLOC((char *)pss,pslen*sizeof(int));      pss = (int *)REALLOC((char *)pss,pslen*sizeof(int));
                 psc = (P *)REALLOC((char *)psc,pslen*sizeof(P)); /* XXX */      psc = (P *)REALLOC((char *)psc,pslen*sizeof(P)); /* XXX */
         }    }
         ps[psn] = a;    ps[psn] = a;
         psh[psn] = BDY(ps[psn])->dl;    psh[psn] = BDY(ps[psn])->dl;
         pss[psn] = ps[psn]->sugar;    pss[psn] = ps[psn]->sugar;
         return psn++;    return psn++;
 }  }
   
 void reducebase_dehomo(f,g)  void reducebase_dehomo(NODE f,NODE *g)
 NODE f,*g;  
 {  {
         int n,i,j,k;    int n,i,j,k;
         int *r;    int *r;
         DL *w,d;    DL *w,d;
         DP u;    DP u;
         NODE t,top;    NODE t,top;
   
         n = length(f);    n = length(f);
         w = (DL *)ALLOCA(n*sizeof(DL));    w = (DL *)ALLOCA(n*sizeof(DL));
         r = (int *)ALLOCA(n*sizeof(int));    r = (int *)ALLOCA(n*sizeof(int));
         for ( i = 0, t = f; i < n; i++, t = NEXT(t) ) {    for ( i = 0, t = f; i < n; i++, t = NEXT(t) ) {
                 r[i] = (int)BDY(t); w[i] = psh[r[i]];      r[i] = (int)BDY(t); w[i] = psh[r[i]];
         }    }
         for ( i = 0; i < n; i++ ) {    for ( i = 0; i < n; i++ ) {
                 for ( j = 0, d = w[i]; j < n; j++ ) {      for ( j = 0, d = w[i]; j < n; j++ ) {
                         if ( j != i ) {        if ( j != i ) {
                                 for ( k = 0; k < NVars; k++ )          for ( k = 0; k < NVars; k++ )
                                         if ( d->d[k] < w[j]->d[k] )            if ( d->d[k] < w[j]->d[k] )
                                                 break;              break;
                                 if ( k == NVars )          if ( k == NVars )
                                         break;            break;
                         }        }
                 }      }
                 if ( j != n )      if ( j != n )
                         r[i] = -1;        r[i] = -1;
         }    }
         for ( top = 0, i = n-1; i >= 0; i-- )    for ( top = 0, i = n-1; i >= 0; i-- )
                 if ( r[i] >= 0 ) {      if ( r[i] >= 0 ) {
                         dp_load(r[i],&ps[r[i]]); dp_dehomo(ps[r[i]],&u);        dp_load(r[i],&ps[r[i]]); dp_dehomo(ps[r[i]],&u);
                         if ( GenTrace ) {        if ( GenTrace ) {
                                 Q q;          Q q;
                                 LIST hist;          LIST hist;
                                 NODE node;          NODE node;
   
                                 STOQ(r[i],q);          STOQ(r[i],q);
                                 node = mknode(4,0,q,0,0);          node = mknode(4,NULLP,q,NULLP,NULLP);
                                 MKLIST(hist,node);          MKLIST(hist,node);
                                 MKNODE(TraceList,hist,0);          MKNODE(TraceList,hist,0);
                         }        }
                         j = newps(u,0,0);        j = newps(u,0,0);
                         MKNODE(t,j,top); top = t;        MKNODE(t,j,top); top = t;
                 }      }
         *g = top;    *g = top;
 }  }
   
 NODE append_one(f,n)  NODE append_one(NODE f,int n)
 NODE f;  
 int n;  
 {  {
         NODE t;    NODE t;
   
         if ( Reverse || !f ) {    if ( Reverse || !f ) {
                 MKNODE(t,(pointer)n,f); return t;      MKNODE(t,(pointer)n,f); return t;
         } else {    } else {
                 for ( t = f; NEXT(t); t = NEXT(t) );      for ( t = f; NEXT(t); t = NEXT(t) );
                 MKNODE(NEXT(t),(pointer)n,0);      MKNODE(NEXT(t),(pointer)n,0);
                 return f;      return f;
         }    }
 }  }
   
 DP_pairs minp( d, prest )  DP_pairs minp( DP_pairs d, DP_pairs *prest )
 DP_pairs d, *prest;  
 {  {
         register DP_pairs m, ml, p, l;    register DP_pairs m, ml, p, l;
         register DL lcm;    register DL lcm;
         register int s, nv = CNVars;    register int s, nv = CNVars;
         register int (*cmpfun)() = cmpdl;    register int (*cmpfun)() = cmpdl;
   
         if ( !(p = NEXT(m = d)) ) {    if ( !(p = NEXT(m = d)) ) {
                 *prest = p;      *prest = p;
                 NEXT(m) = 0;      NEXT(m) = 0;
                 return m;      return m;
         }    }
         for ( lcm = m->lcm, s = m->sugar, ml = 0, l = m; p; p = NEXT(l = p) )    for ( lcm = m->lcm, s = m->sugar, ml = 0, l = m; p; p = NEXT(l = p) )
                 if ( NoSugar ? (*cmpfun)( nv, lcm, p->lcm ) >= 0 :      if ( NoSugar ? (*cmpfun)( nv, lcm, p->lcm ) >= 0 :
                      (s > p->sugar || s == p->sugar && (*cmpfun)( nv, lcm, p->lcm ) >= 0) )           (s > p->sugar || s == p->sugar && (*cmpfun)( nv, lcm, p->lcm ) >= 0) )
                         ml = l,  lcm = (m = p)->lcm,  s = p->sugar;        ml = l,  lcm = (m = p)->lcm,  s = p->sugar;
         if ( !ml ) *prest = NEXT(m);    if ( !ml ) *prest = NEXT(m);
         else {    else {
                 NEXT(ml) = NEXT(m);      NEXT(ml) = NEXT(m);
                 *prest = d;      *prest = d;
         }    }
         NEXT(m) = 0;    NEXT(m) = 0;
         return m;    return m;
 }  }
   
 void minsugar(d,dm,dr)  void minsugar(DP_pairs d,DP_pairs *dm,DP_pairs *dr)
 DP_pairs d;  
 DP_pairs *dm,*dr;  
 {  {
         int msugar;    int msugar;
         DP_pairs t,dm0,dr0,dmt,drt;    DP_pairs t,dm0,dr0,dmt,drt;
   
         for ( msugar = d->sugar, t = NEXT(d); t; t = NEXT(t) )    for ( msugar = d->sugar, t = NEXT(d); t; t = NEXT(t) )
                 if ( t->sugar < msugar )      if ( t->sugar < msugar )
                         msugar = t->sugar;        msugar = t->sugar;
         dm0 = 0; dr0 = 0;    dm0 = 0; dr0 = 0;
         for ( t = d; t; t = NEXT(t) ) {    for ( t = d; t; t = NEXT(t) ) {
                 if ( t->sugar == msugar ) {      if ( t->sugar == msugar ) {
                         NEXTDPP(dm0,dmt);        NEXTDPP(dm0,dmt);
                         dmt->dp1 = t->dp1; dmt->dp2 = t->dp2;        dmt->dp1 = t->dp1; dmt->dp2 = t->dp2;
                         dmt->lcm = t->lcm; dmt->sugar = t->sugar;        dmt->lcm = t->lcm; dmt->sugar = t->sugar;
                 } else {      } else {
                         NEXTDPP(dr0,drt);        NEXTDPP(dr0,drt);
                         drt->dp1 = t->dp1; drt->dp2 = t->dp2;        drt->dp1 = t->dp1; drt->dp2 = t->dp2;
                         drt->lcm = t->lcm; drt->sugar = t->sugar;        drt->lcm = t->lcm; drt->sugar = t->sugar;
                 }      }
         }    }
         if ( dm0 ) NEXT(dmt) = 0;    if ( dm0 ) NEXT(dmt) = 0;
         if ( dr0 ) NEXT(drt) = 0;    if ( dr0 ) NEXT(drt) = 0;
         *dm = dm0; *dr = dr0;    *dm = dm0; *dr = dr0;
 }  }
   
 NODE gb(f,m,subst)  NODE gb(NODE f,int m,NODE subst)
 NODE f;  
 int m;  
 NODE subst;  
 {  {
         int i,nh,prev,mag;    int i,nh,prev,mag,mag0,magt;
         NODE r,g,gall;    NODE r,g,gall;
         DP_pairs d,d1;    DP_pairs d;
         DP_pairs l;    DP_pairs l;
         DP h,nf,nfm,dp1,dp2;    DP h,nf,nfm,dp1,dp2;
         MP mp;    MP mp;
         struct oEGT tnf0,tnf1,tnfm0,tnfm1,tpz0,tpz1,tsp0,tsp1,tspm0,tspm1,tnp0,tnp1,tmp0,tmp1;    struct oEGT tnf0,tnf1,tnfm0,tnfm1,tpz0,tpz1,tnp0,tnp1;
         int skip_nf_flag;    int skip_nf_flag;
         double t_0;    double t_0;
         Q q;    Q q;
         int new_sugar;    int new_sugar;
         static prev_sugar = -1;    static prev_sugar = -1;
   
         Max_mag = 0;    Max_mag = 0;
         prev = 1;    Max_coef = 0;
         doing_f4 = 0;    prev = 1;
         if ( m ) {    doing_f4 = 0;
                 psm = (DP *)MALLOC(pslen*sizeof(DP));    init_denomlist();
                 for ( i = 0; i < psn; i++ )    if ( m ) {
                         if ( psh[i] && !validhc(psc[i],m,subst) )      psm = (DP *)MALLOC(pslen*sizeof(DP));
                                 return 0;      for ( i = 0; i < psn; i++ )
                         else        if ( psh[i] && !validhc(psc[i],m,subst) )
                                 _dp_mod(ps[i],m,subst,&psm[i]);          return 0;
         }        else
         for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {          _dp_mod(ps[i],m,subst,&psm[i]);
                 i = (int)BDY(r);    }
                 d = updpairs(d,g,i);    for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {
                 g = updbase(g,i);      i = (int)BDY(r);
                 gall = append_one(gall,i);      d = updpairs(d,g,i);
         }      g = updbase(g,i);
         while ( d ) {      gall = append_one(gall,i);
                 get_eg(&tmp0);    }
                 l = minp(d,&d1); d = d1;    while ( d ) {
                 get_eg(&tmp1); add_eg(&eg_mp,&tmp0,&tmp1);      l = minp(d,&d);
                 if ( m ) {      if ( m ) {
                         get_eg(&tspm0);        _dp_sp_mod_dup(psm[l->dp1],psm[l->dp2],m,&h);
                         _dp_sp_mod_dup(psm[l->dp1],psm[l->dp2],m,&h);        if ( h )
                         if ( h )          new_sugar = h->sugar;
                                 new_sugar = h->sugar;        get_eg(&tnfm0);
                         get_eg(&tspm1); add_eg(&eg_spm,&tspm0,&tspm1);        _dp_nf_mod_destructive(gall,h,psm,m,0,&nfm);
                         get_eg(&tnfm0);        get_eg(&tnfm1); add_eg(&eg_nfm,&tnfm0,&tnfm1);
                         _dp_nf_mod_destructive(gall,h,psm,m,0,&nfm);      } else
                         get_eg(&tnfm1); add_eg(&eg_nfm,&tnfm0,&tnfm1);        nfm = (DP)1;
                 } else      if ( nfm ) {
                         nfm = (DP)1;        if ( Demand ) {
                 if ( nfm ) {          if ( dp_load_t(psn,&nf) ) {
                         get_eg(&tsp0);            skip_nf_flag = 1;
                         if ( Demand ) {            goto skip_nf;
                                 if ( dp_load_t(psn,&nf) ) {          } else {
                                         skip_nf_flag = 1;            skip_nf_flag = 0;
                                         get_eg(&tsp1); add_eg(&eg_nf,&tsp0,&tsp1);            dp_load(l->dp1,&dp1); dp_load(l->dp2,&dp2);
                                         tnf0=tsp0; tnf1=tsp1;            dp_sp(dp1,dp2,&h);
                                         goto skip_nf;          }
                                 } else {        } else
                                         skip_nf_flag = 0;          dp_sp(ps[l->dp1],ps[l->dp2],&h);
                                         dp_load(l->dp1,&dp1); dp_load(l->dp2,&dp2);        if ( GenTrace ) {
                                         dp_sp(dp1,dp2,&h);          STOQ(l->dp1,q); ARG1(BDY((LIST)BDY(NEXT(TraceList)))) = q;
                                 }          STOQ(l->dp2,q); ARG1(BDY((LIST)BDY(TraceList))) = q;
                         } else        }
                                 dp_sp(ps[l->dp1],ps[l->dp2],&h);        if ( h )
                         if ( GenTrace ) {          new_sugar = h->sugar;
                                 STOQ(l->dp1,q); ARG1(BDY((LIST)BDY(NEXT(TraceList)))) = q;        get_eg(&tnf0);
                                 STOQ(l->dp2,q); ARG1(BDY((LIST)BDY(TraceList))) = q;        t_0 = get_rtime();
                         }        if ( PCoeffs || dp_fcoeffs )
                         if ( h )          _dp_nf(gall,h,ps,!Top,&nf);
                                 new_sugar = h->sugar;        else
                         get_eg(&tsp1); add_eg(&eg_sp,&tsp0,&tsp1);          _dp_nf_z(gall,h,ps,!Top,DP_Multiple,&nf);
                         get_eg(&tnf0);        if ( DP_Print && nf )
                         t_0 = get_rtime();          fprintf(asir_out,"(%.3g)",get_rtime()-t_0);
                         if ( !Multiple )        get_eg(&tnf1); add_eg(&eg_nf,&tnf0,&tnf1);
                                 _dp_nf(gall,h,ps,!Top,&nf);      } else
                         else        nf = 0;
 #if MPI  
                                 _dp_nf_ptozp_mpi(gall,h,ps,!Top,Multiple,&nf);  
 #else  
                                 _dp_nf_ptozp(gall,h,ps,!Top,Multiple,&nf);  
 #endif  
                         if ( DP_Print )  
                                 fprintf(asir_out,"(%.3g)",get_rtime()-t_0);  
                         get_eg(&tnf1); add_eg(&eg_nf,&tnf0,&tnf1);  
                 } else  
                         nf = 0;  
 skip_nf:  skip_nf:
                 if ( nf ) {      if ( nf ) {
                         NZR++;        NZR++;
                         get_eg(&tpz0);        get_eg(&tpz0);
                         prim_part(nf,0,&h);        prim_part(nf,0,&h);
                         get_eg(&tpz1); add_eg(&eg_pz,&tpz0,&tpz1);        get_eg(&tpz1); add_eg(&eg_pz,&tpz0,&tpz1);
                         get_eg(&tnp0);        add_denomlist((P)BDY(h)->c);
                         if ( Demand && skip_nf_flag )        get_eg(&tnp0);
                                 nh = newps_nosave(h,m,subst);        if ( Demand && skip_nf_flag )
                         else          nh = newps_nosave(h,m,subst);
                                 nh = newps(h,m,subst);        else
                         get_eg(&tnp1); add_eg(&eg_np,&tnp0,&tnp1);          nh = newps(h,m,subst);
                         if ( nh < 0 )        get_eg(&tnp1); add_eg(&eg_np,&tnp0,&tnp1);
                                 return 0;        if ( nh < 0 )
                         d = updpairs(d,g,nh);          return 0;
                         g = updbase(g,nh);        d = updpairs(d,g,nh);
                         gall = append_one(gall,nh);        g = updbase(g,nh);
                         if ( !dp_fcoeffs && ShowMag ) {        gall = append_one(gall,nh);
                                 for ( mag = 0, mp = BDY(h); mp; mp = NEXT(mp) )        if ( !dp_fcoeffs && ShowMag ) {
                                         mag += p_mag((P)mp->c);          for ( mag = 0, mag0 = 0, mp = BDY(h); mp; mp = NEXT(mp) ) {
                                 Max_mag = MAX(Max_mag,mag);            magt = p_mag((P)mp->c);
                         }            mag0 = MAX(mag0,magt);
                         if ( DP_Print ) {            mag += magt;
                                 if ( !prev )          }
                                         fprintf(asir_out,"\n");          Max_coef = MAX(Max_coef,mag0);
                                 print_split_e(&tnf0,&tnf1); print_split_e(&tpz0,&tpz1);          Max_mag = MAX(Max_mag,mag);
                                 printdl(psh[nh]);        }
                                 fprintf(asir_out,"(%d,%d),nb=%d,nab=%d,rp=%d,sugar=%d",        if ( DP_Print ) {
                                         l->dp1,l->dp2,length(g),length(gall),DPPlength(d),          if ( !prev )
                                         pss[nh]);            fprintf(asir_out,"\n");
                                 if ( ShowMag )          print_split_e(&tnf0,&tnf1); print_split_e(&tpz0,&tpz1);
                                         fprintf(asir_out,",mag=%d",mag);          printdl(psh[nh]);
                                 fprintf(asir_out,"\n"); fflush(asir_out);          fprintf(asir_out,"(%d,%d),nb=%d,nab=%d,rp=%d,sugar=%d",
                         } else if ( DP_PrintShort ) {            l->dp1,l->dp2,length(g),length(gall),DPPlength(d),
                                 fprintf(asir_out,"+"); fflush(asir_out);            pss[nh]);
                         }          if ( ShowMag )
                         prev = 1;            fprintf(asir_out,",mag=(%d,%d)",mag,mag0);
                 } else {          fprintf(asir_out,"\n"); fflush(asir_out);
                         if ( m )        } else if ( DP_PrintShort ) {
                                 add_eg(&eg_znfm,&tnfm0,&tnfm1);          fprintf(asir_out,"+"); fflush(asir_out);
                         ZR++;        }
                         if ( DP_Print || DP_PrintShort ) {        prev = 1;
                                 if ( new_sugar != prev_sugar ) {      } else {
                                         fprintf(asir_out,"[%d]",new_sugar);        if ( m )
                                         prev_sugar = new_sugar;          add_eg(&eg_znfm,&tnfm0,&tnfm1);
                                 }        ZR++;
                                 fprintf(asir_out,"."); fflush(asir_out); prev = 0;        if ( DP_Print || DP_PrintShort ) {
                         }          if ( new_sugar != prev_sugar ) {
                 }            fprintf(asir_out,"[%d]",new_sugar);
         }            prev_sugar = new_sugar;
         if ( DP_Print || DP_PrintShort )          }
                 fprintf(asir_out,"gb done\n");          fprintf(asir_out,"."); fflush(asir_out); prev = 0;
         return g;        }
       }
     }
     if ( DP_Print || DP_PrintShort )
       fprintf(asir_out,"gb done\n");
     return g;
 }  }
   
 NODE gbd(f,m,subst,dlist)  NODE gb_mod(NODE f,int m)
 NODE f;  
 int m;  
 NODE subst;  
 NODE dlist;  
 {  {
         int i,nh,prev;    int i,nh,prev;
         NODE r,g,gall;    NODE r,g,gall;
         struct dp_pairs ol;    DP_pairs d;
         DP_pairs l;    DP_pairs l;
         NODE pair;    DP h,nf;
         DP h,nf,nfm,dp1,dp2;    struct oEGT tnfm0,tnfm1,tpz0,tpz1;
         struct oEGT tnf0,tnf1,tnfm0,tnfm1,tpz0,tpz1,tsp0,tsp1,tspm0,tspm1,tnp0,tnp1;  
   
         prev = 1;    prev = 1;
         l = &ol;    for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {
         if ( m ) {      i = (int)BDY(r);
                 psm = (DP *)MALLOC(pslen*sizeof(DP));      d = updpairs(d,g,i);
                 for ( i = 0; i < psn; i++ )      g = updbase(g,i);
                         if ( psh[i] && !validhc(psc[i],m,subst) )      gall = append_one(gall,i);
                                 return 0;    }
                         else    while ( d ) {
                                 _dp_mod(ps[i],m,subst,&psm[i]);      l = minp(d,&d);
         }      if ( PCoeffs ) {
         for ( gall = g = 0, r = f; r; r = NEXT(r) ) {        dp_sp_mod(ps[l->dp1],ps[l->dp2],m,&h);
                 i = (int)BDY(r);        get_eg(&tnfm0);
                 g = updbase(g,i);        dp_nf_mod(gall,h,ps,m,!Top,&nf);
                 gall = append_one(gall,i);      } else {
         }        _dp_sp_mod_dup(ps[l->dp1],ps[l->dp2],m,&h);
         while ( dlist ) {        get_eg(&tnfm0);
                 pair = BDY((LIST)BDY(dlist)); dlist = NEXT(dlist);        _dp_nf_mod_destructive(gall,h,ps,m,!Top,&nf);
                 l->dp1 = QTOS((Q)BDY(pair)); pair = NEXT(pair);      }
                 l->dp2 = QTOS((Q)BDY(pair));      get_eg(&tnfm1); add_eg(&eg_nfm,&tnfm0,&tnfm1);
                 if ( m ) {      if ( nf ) {
                         get_eg(&tspm0);        NZR++;
                         _dp_sp_mod_dup(ps[l->dp1],ps[l->dp2],m,&h);        get_eg(&tpz0);
                         get_eg(&tspm1); add_eg(&eg_spm,&tspm0,&tspm1);        prim_part(nf,m,&h);
                         get_eg(&tnfm0);        get_eg(&tpz1); add_eg(&eg_pz,&tpz0,&tpz1);
                         _dp_nf_mod_destructive(gall,h,ps,m,!Top,&nf);        nh = newps_mod(h,m);
                         get_eg(&tnfm1); add_eg(&eg_nfm,&tnfm0,&tnfm1);        if ( nh < 0 )
                 } else          return 0;
                         nfm = (DP)1;        d = updpairs(d,g,nh);
                 if ( nfm ) {        g = updbase(g,nh);
                         get_eg(&tsp0);        gall = append_one(gall,nh);
                         if ( Demand ) {        if ( DP_Print ) {
                                 dp_load(l->dp1,&dp1); dp_load(l->dp2,&dp2);          if ( !prev )
                                 dp_sp(dp1,dp2,&h);            fprintf(asir_out,"\n");
                         } else          print_split_eg(&tnfm0,&tnfm1); fflush(asir_out);
                                 dp_sp(ps[l->dp1],ps[l->dp2],&h);          fprintf(asir_out,"(%d,%d),nb=%d,nab=%d,rp=%d,sugar=%d",l->dp1,l->dp2,length(g),length(gall),DPPlength(d),pss[nh]);
                         get_eg(&tsp1); add_eg(&eg_sp,&tsp0,&tsp1);          printdl(psh[nh]); fprintf(asir_out,"\n"); fflush(asir_out);
                         get_eg(&tnfm0);        } else if ( DP_PrintShort ) {
                         get_eg(&tnf0);          fprintf(asir_out,"+"); fflush(asir_out);
                         if ( !Multiple )        }
                                 _dp_nf(gall,h,ps,!Top,&nf);        prev = 1;
                         else      } else {
                                 _dp_nf_ptozp(gall,h,ps,!Top,Multiple,&nf);        add_eg(&eg_znfm,&tnfm0,&tnfm1);
                         get_eg(&tnf1); add_eg(&eg_nf,&tnf0,&tnf1);        ZR++;
                 } else        if ( DP_Print || DP_PrintShort ) {
                         nf = 0;          fprintf(asir_out,"."); fflush(asir_out); prev = 0;
                 if ( nf ) {        }
                         NZR++;      }
                         get_eg(&tpz0);    }
                         prim_part(nf,0,&h);    if ( DP_Print || DP_PrintShort )
                         get_eg(&tpz1); add_eg(&eg_pz,&tpz0,&tpz1);      fprintf(asir_out,"gb_mod done\n");
                         get_eg(&tnp0);    return g;
                         nh = newps(h,m,subst);  
                         get_eg(&tnp1); add_eg(&eg_np,&tnp0,&tnp1);  
                         if ( nh < 0 )  
                                 return 0;  
                         g = updbase(g,nh);  
                         gall = append_one(gall,nh);  
                         if ( DP_Print ) {  
                                 if ( !prev )  
                                         fprintf(asir_out,"\n");  
                                 print_split_eg(&tnf0,&tnf1); fflush(asir_out);  
                                 fprintf(asir_out,"(%d,%d),nb=%d,nab=%d,rp=%d,sugar=%d",l->dp1,l->dp2,length(g),length(gall),length(dlist),pss[nh]);  
                                 printdl(psh[nh]); fprintf(asir_out,"\n"); fflush(asir_out);  
                         }  
                         prev = 1;  
                 } else {  
                         if ( m )  
                                 add_eg(&eg_znfm,&tnfm0,&tnfm1);  
                         ZR++;  
                         if ( DP_Print ) {  
                                 fprintf(asir_out,"."); fflush(asir_out); prev = 0;  
                         }  
                 }  
         }  
         if ( DP_Print )  
                 fprintf(asir_out,"gb done\n");  
         return g;  
 }  }
   
 NODE gb_mod(f,m)  DP_pairs updpairs( DP_pairs d, NODE /* of index */ g, int t)
 NODE f;  
 int m;  
 {  {
         int i,nh,prev;    register DP_pairs d1, dd, nd;
         NODE r,g,gall;    int dl,dl1;
         DP_pairs d,d1;  
         DP_pairs l;  
         DP h,nf;  
         struct oEGT tnfm0,tnfm1,tspm0,tspm1,tmp0,tmp1,tpz0,tpz1;  
   
         prev = 1;    if ( !g ) return d;
         for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {    if ( !NoCriB && d ) {
                 i = (int)BDY(r);      dl = DPPlength(d);
                 d = updpairs(d,g,i);      d = criterion_B( d, t );
                 g = updbase(g,i);      dl -= DPPlength(d); N_BP += dl;
                 gall = append_one(gall,i);    }
         }    d1 = newpairs( g, t );
         while ( d ) {    if ( NEXT(d1) ) {
                 get_eg(&tmp0);      dl = DPPlength(d1); TP += dl;
                 l = minp(d,&d1); d = d1;      d1 = criterion_M( d1 );
                 get_eg(&tmp1); add_eg(&eg_mp,&tmp0,&tmp1); get_eg(&tspm0);      dl1 = DPPlength(d1); NMP += (dl-dl1); dl = dl1;
                 if ( PCoeffs ) {      d1 = criterion_F( d1 );
                         dp_sp_mod(ps[l->dp1],ps[l->dp2],m,&h);      dl1 = DPPlength(d1); NFP += (dl-dl1); dl = dl1;
                         get_eg(&tspm1); add_eg(&eg_spm,&tspm0,&tspm1); get_eg(&tnfm0);    } else
                         dp_nf_mod(gall,h,ps,m,!Top,&nf);      dl = 1;
                 } else {    if ( !do_weyl )
                         _dp_sp_mod_dup(ps[l->dp1],ps[l->dp2],m,&h);      for ( dd = 0; d1; d1 = nd ) {
                         get_eg(&tspm1); add_eg(&eg_spm,&tspm0,&tspm1); get_eg(&tnfm0);        nd = NEXT(d1);
                         _dp_nf_mod_destructive(gall,h,ps,m,!Top,&nf);        if ( !criterion_2( d1->dp1, d1->dp2 ) ) {
                 }          NEXT(d1) = dd;
                 get_eg(&tnfm1); add_eg(&eg_nfm,&tnfm0,&tnfm1);          dd = d1;
                 if ( nf ) {        }
                         NZR++;      }
                         get_eg(&tpz0);    else
                         prim_part(nf,m,&h);      dd = d1;
                         get_eg(&tpz1); add_eg(&eg_pz,&tpz0,&tpz1);    dl1 = DPPlength(dd); NDP += (dl-dl1);
                         nh = newps_mod(h,m);    if ( !(nd = d) ) return dd;
                         if ( nh < 0 )    while ( nd = NEXT(d1 = nd) ) ;
                                 return 0;    NEXT(d1) = dd;
                         d = updpairs(d,g,nh);    return d;
                         g = updbase(g,nh);  
                         gall = append_one(gall,nh);  
                         if ( DP_Print ) {  
                                 if ( !prev )  
                                         fprintf(asir_out,"\n");  
                                 print_split_eg(&tnfm0,&tnfm1); fflush(asir_out);  
                                 fprintf(asir_out,"(%d,%d),nb=%d,nab=%d,rp=%d,sugar=%d",l->dp1,l->dp2,length(g),length(gall),DPPlength(d),pss[nh]);  
                                 printdl(psh[nh]); fprintf(asir_out,"\n"); fflush(asir_out);  
                         } else if ( DP_PrintShort ) {  
                                 fprintf(asir_out,"+"); fflush(asir_out);  
                         }  
                         prev = 1;  
                 } else {  
                         add_eg(&eg_znfm,&tnfm0,&tnfm1);  
                         ZR++;  
                         if ( DP_Print || DP_PrintShort ) {  
                                 fprintf(asir_out,"."); fflush(asir_out); prev = 0;  
                         }  
                 }  
         }  
         if ( DP_Print || DP_PrintShort )  
                 fprintf(asir_out,"gb_mod done\n");  
         return g;  
 }  }
   
 void _dp_nf(b,g,ps,full,rp)  DP_pairs newpairs( NODE /* of index */ g, int t )
 NODE b;  
 DP g;  
 DP *ps;  
 int full;  
 DP *rp;  
 {  {
         DP u,p,d,s,t,mult;    register NODE r;
         P coef;    register DL tdl = psh[t];
         NODE l;    register int ts;
         MP m,mr;    register DP_pairs p, last;
         int sugar,psugar;    int dp;
     register DL dl;
     register int s;
   
         if ( !g ) {    ts = pss[t] - tdl->td;
                 *rp = 0; return;    for ( last = 0, r = g; r; r = NEXT(r) ) {
         }      NEXT(p = NEWDP_pairs) = last;
         sugar = g->sugar;      last = p;
         for ( d = 0; g; ) {      dp = p->dp1 = (int)BDY(r);  p->dp2 = t;
                 for ( u = 0, l = b; l; l = NEXT(l) ) {      p->lcm = lcm_of_DL(CNVars, dl = psh[dp], tdl, (DL)0 );
                         if ( dl_redble(BDY(g)->dl,psh[(int)BDY(l)]) ) {  #if 0
                                 dp_load((int)BDY(l),&p);      if ( do_weyl )
                                 /* t+u = coef*(d+g) - mult*p (t = coef*d) */        p->sugar = dl_weyl_weight(p->lcm);
                                 dp_red(d,g,p,&t,&u,&coef,&mult);      else
                                 psugar = (BDY(g)->dl->td - BDY(p)->dl->td) + p->sugar;  #endif
                                 sugar = MAX(sugar,psugar);        p->sugar = (ts > (s = pss[dp] - dl->td) ? ts : s) + p->lcm->td;
                                 if ( GenTrace ) {    }
                                         LIST hist;    return last;
                                         Q cq;  
                                         NODE node,node0;  
   
                                         STOQ((int)BDY(l),cq);  
                                         node0 = mknode(4,coef,cq,mult,ONE);  
                                         MKLIST(hist,node0);  
                                         MKNODE(node,hist,TraceList); TraceList = node;  
                                 }  
                                 if ( !u ) {  
                                         if ( d )  
                                                 d->sugar = sugar;  
                                         *rp = d; return;  
                                 }  
                                 d = t;  
                                 break;  
                         }  
                 }  
                 if ( u )  
                         g = u;  
                 else if ( !full ) {  
                         if ( g ) {  
                                 MKDP(g->nv,BDY(g),t); t->sugar = sugar; g = t;  
                         }  
                         *rp = g; return;  
                 } else {  
                         m = BDY(g); NEWMP(mr); mr->dl = m->dl; mr->c = m->c;  
                         NEXT(mr) = 0; MKDP(g->nv,mr,t); t->sugar = mr->dl->td;  
                         addd(CO,d,t,&s); d = s;  
                         dp_rest(g,&t); g = t;  
                 }  
         }  
         if ( d )  
                 d->sugar = sugar;  
         *rp = d;  
 }  }
   
 #define SAFENM(q) ((q)?NM(q):0)  DP_pairs criterion_B( DP_pairs d, int s )
   {
     register DP_pairs dd, p;
     register DL tij, t = psh[s], dltmp;
   
 double pz_t_e, pz_t_d, pz_t_d1, pz_t_c, im_t_s, im_t_r;    if ( !d ) return 0;
     NEWDL( dltmp, CNVars );
     for ( dd = 0; d; d = p ) {
       p = NEXT(d),
       tij = d->lcm;
       if ( tij->td != lcm_of_DL(CNVars, tij, t, dltmp )->td
         || !dl_equal(CNVars, tij, dltmp )
         || (tij->td == lcm_of_DL(CNVars, psh[d->dp1], t, dltmp )->td
             && dl_equal(CNVars, dltmp, tij ))
         || (tij->td == lcm_of_DL(CNVars, psh[d->dp2], t, dltmp )->td
             && dl_equal(CNVars, dltmp, tij )) ) {
         NEXT(d) = dd;
         dd = d;
       }
     }
     return dd;
   }
   
 extern int GenTrace;  DP_pairs criterion_M( DP_pairs d1 )
 extern NODE TraceList;  
   
 void _dp_nf_ptozp(b,g,ps,full,multiple,r)  
 NODE b;  
 DP g;  
 DP *ps;  
 int full,multiple;  
 DP *r;  
 {  {
         DP u,dp,rp,t,t1,t2,red,shift;    register DP_pairs dd, e, d3, d2, p;
         Q dc,rc,dcq,rcq,cont,hr,hred,cr,cred,mcred,c,gcd,cq;    register DL itdl, jtdl;
         N gn,tn,cn;    register int itdltd, jtdltd;
         NODE l;  
         MP m,mr;  
         int hmag,denom;  
         int sugar,psugar;  
         NODE dist;  
         STRING imul;  
         int ndist;  
         int kara_bit;  
         extern int kara_mag;  
         double get_rtime();  
         double t_0,t_00,tt,ttt,t_p,t_m,t_m1,t_m2,t_s,t_g,t_a;  
         LIST hist;  
         NODE node;  
         Q rcred,mrcred;  
   
         if ( !g ) {    for ( dd = 0, e = d1; e; e = d3 ) {
                 *r = 0; return;      if ( !(d2 = NEXT(e)) ) {
         }        NEXT(e) = dd;
         pz_t_e = pz_t_d = pz_t_d1 = pz_t_c = 0;        return e;
         t_p = t_m = t_m1 = t_m2 = t_s = t_g = t_a = 0;      }
       itdltd = (itdl = e->lcm)->td;
       for ( d3 = 0; d2; d2 = p ) {
         p = NEXT(d2),
         jtdltd = (jtdl = d2->lcm)->td;
         if ( jtdltd == itdltd  )
           if ( dl_equal(CNVars, itdl, jtdl ) ) ;
           else if ( dl_redble( jtdl, itdl ) ) continue;
           else if ( dl_redble( itdl, jtdl ) ) goto delit;
           else ;
         else if ( jtdltd > itdltd )
           if ( dl_redble( jtdl, itdl ) ) continue;
           else ;
         else if ( dl_redble( itdl, jtdl ) ) goto delit;
         NEXT(d2) = d3;
         d3 = d2;
       }
       NEXT(e) = dd;
       dd = e;
       continue;
       /**/
     delit:  NEXT(d2) = d3;
       d3 = d2;
       for ( ; p; p = d2 ) {
         d2 = NEXT(p);
         NEXT(p) = d3;
         d3 = p;
       }
     }
     return dd;
   }
   
         denom = Denominator?Denominator:1;  static DP_pairs collect_pairs_of_hdlcm( DP_pairs d1, DP_pairs *prest )
         hmag = multiple*HMAG(g)/denom;  {
         kara_bit = kara_mag*27; /* XXX */    register DP_pairs w, p, r, s;
         if ( Dist ) {    register DL ti;
                 dist = BDY(Dist);    register int td;
                 ndist = length(dist);  
         }  
         sugar = g->sugar;  
   
         dc = 0; dp = 0; rc = ONE; rp = g;    td = (ti = (w = d1)->lcm)->td;
         MKSTR(imul,"dp_imul_index");    s = NEXT(w);
     NEXT(w) = 0;
     for ( r = 0; s; s = p ) {
       p = NEXT(s);
       if ( td == s->lcm->td && dl_equal(CNVars, ti, s->lcm ) )
       {
         NEXT(s) = w;
         w = s;
       } else {
         NEXT(s) = r;
         r = s;
       }
     }
     *prest = r;
     return w;
   }
   
         /* g = dc*dp+rc*rp */  int criterion_2( int dp1, int dp2 )
         for ( ; rp; ) {  {
                 for ( u = 0, l = b; l; l = NEXT(l) ) {    register int i, *d1, *d2;
                         if ( dl_redble(BDY(rp)->dl,psh[(int)BDY(l)]) ) {  
                                 t_0 = get_rtime();  
                                 dp_load((int)BDY(l),&red);  
                                 hr = (Q)BDY(rp)->c; hred = (Q)BDY(red)->c;  
                                 gcdn(NM(hr),NM(hred),&gn);  
                                 divsn(NM(hred),gn,&tn); NTOQ(tn,SGN(hred),cr);  
                                 divsn(NM(hr),gn,&tn); NTOQ(tn,SGN(hr),cred);  
                                 tt = get_rtime()-t_0; t_p += tt;  
   
                                 t_0 = get_rtime();    d1 = psh[dp1]->d,  d2 = psh[dp2]->d;
                                 dp_subd(rp,red,&shift);    for ( i = CNVars; --i >= 0; d1++, d2++ )
                                 if ( Dist && ndist && HMAG(red) > kara_bit ) {      if ( (*d1 <= *d2 ? *d1 : *d2) > 0  ) return 0;
                                         NODE n0,n1,n2,n3;    return 1;
                                         int i,s,id;  }
                                         Obj dmy;  
                                         Q ind;  
   
                                         if ( DP_Print )  DP_pairs criterion_F( DP_pairs d1 )
                                                 fprintf(asir_out,"d");  {
                                         i = (int)BDY(l); STOQ(i,ind);    DP_pairs rest, head;
                                         chsgnp((P)cred,(P *)&mcred);    register DP_pairs last, p, r, w;
     register int s;
   
                                         MKNODE(n3,ind,0); MKNODE(n2,mcred,n3);    for ( head = last = 0, p = d1; NEXT(p); ) {
                                         MKNODE(n1,imul,n2); MKNODE(n0,BDY(dist),n1);      s = (r = w = collect_pairs_of_hdlcm( p, &rest ))->sugar;
                                         Pox_rpc(n0,&dmy);      while ( w = NEXT(w) )
                                         muldc(CO,rp,(P)cr,&t);        if ( !do_weyl && criterion_2( w->dp1, w->dp2 ) ) {
                                         NEXT(n0)=0;          r = w;
                                         Pox_pop_local(n0,&t1);          break;
                                 } else {        } else if ( w->sugar < s ) s = (r = w)->sugar;
 /*      if ( last ) NEXT(last) = r;
                                         if ( DP_Print )      else head = r;
                                                 fprintf(asir_out,"l");      NEXT(last = r) = 0;
 */      if ( !(p = rest) ) return head;
                                         t_00 = get_rtime();    }
                                         muldc(CO,rp,(P)cr,&t);    if ( !last ) return p;
                                         ttt = get_rtime()-t_00; t_m1 += ttt/dp_nt(rp);    NEXT(last) = p;
                                         t_00 = get_rtime();    return head;
                                         chsgnp((P)cred,(P *)&mcred);  }
                                         muldc(CO,red,(P)mcred,&t1);  
                                         ttt = get_rtime()-t_00; t_m2 += ttt/dp_nt(red);  
                                 }  
                                 t_00 = get_rtime();  
                                 muld(CO,shift,t1,&t2);  
                                 addd(CO,t,t2,&u);  
                                 tt = get_rtime(); t_m += tt-t_0;  
                                 ttt = get_rtime(); t_s += ttt-t_00;  
   
                                 psugar = (BDY(rp)->dl->td - BDY(red)->dl->td) + red->sugar;  NODE updbase(NODE g,int t)
                                 sugar = MAX(sugar,psugar);  {
                                 if ( GenTrace ) {    g = remove_reducibles(g,t);
                                         /* u = cr*rp + (-cred)*shift*red */    g = append_one(g,t);
                                         STOQ((int)BDY(l),cq);    return g;
                                         node = mknode(4,cr,cq,0,0);  
                                         mulq(cred,rc,&rcred);  
                                         chsgnnum((Num)rcred,(Num *)&mrcred);  
                                         muldc(CO,shift,(P)mrcred,(DP *)&ARG2(node));  
                                         MKLIST(hist,node);  
                                 }  
                                 if ( !u ) {  
                                         if ( dp )  
                                                 dp->sugar = sugar;  
                                         *r = dp;  
                                         if ( GenTrace ) {  
                                                 ARG3(BDY(hist)) = ONE;  
                                                 MKNODE(node,hist,TraceList); TraceList = node;  
                                         }  
                                         goto final;  
                                 }  
                                 break;  
                         }  
                 }  
                 if ( u ) {  
                         if ( HMAG(u) > hmag ) {  
                                 t_0 = get_rtime();  
                                 if ( Dist && HMAG(u) > kara_bit ) {  
                                         if ( DP_Print )  
                                                 fprintf(asir_out,"D");  
                                         dp_ptozp_d(dist,ndist,u,&t);  
                                 } else {  
                                         if ( DP_Print )  
                                                 fprintf(asir_out,"L");  
                                         dp_ptozp_d(0,0,u,&t);  
                                 }  
                                 tt = get_rtime()-t_0; t_g += tt;  
                                 t_0 = get_rtime();  
                                 divsn(NM((Q)BDY(u)->c),NM((Q)BDY(t)->c),&cn); NTOQ(cn,1,cont);  
                                 if ( !dp_fcoeffs && DP_Print ) {  
                                         fprintf(asir_out,"(%d)",p_mag((P)cont)*100/p_mag((P)BDY(u)->c));  
                                         fflush(asir_out);  
                                 }  
                                 mulq(cr,dc,&dcq);  
                                 mulq(cont,rc,&rcq);  
                                 gcdn(SAFENM(dcq),SAFENM(rcq),&gn);  
                                 divsn(SAFENM(dcq),gn,&tn); NTOQ(tn,SGN(dcq),dc);  
                                 divsn(SAFENM(rcq),gn,&tn); NTOQ(tn,SGN(rcq),rc);  
                                 tt = get_rtime()-t_0; t_a += tt;  
                                 rp = t;  
                                 hmag = multiple*HMAG(rp)/denom;  
                                 if ( GenTrace ) {  
                                         NTOQ(gn,1,gcd);  
                                         ARG3(BDY(hist)) = (pointer)gcd;  
                                         MKNODE(node,hist,TraceList); TraceList = node;  
                                 }  
                         } else {  
                                 t_0 = get_rtime();  
                                 mulq(cr,dc,&dcq); dc = dcq;  
                                 tt = get_rtime()-t_0; t_a += tt;  
                                 rp = u;  
                                 if ( GenTrace ) {  
                                         ARG3(BDY(hist)) = (pointer)ONE;  
                                         MKNODE(node,hist,TraceList); TraceList = node;  
                                 }  
                         }  
                 } else if ( !full ) {  
                         if ( rp ) {  
                                 MKDP(rp->nv,BDY(rp),t); t->sugar = sugar; rp = t;  
                         }  
                         *r = rp;  
                         goto final;  
                 } else {  
                         t_0 = get_rtime();  
                         mulq((Q)BDY(rp)->c,rc,&c);  
                         gcdn(SAFENM(dc),SAFENM(c),&gn); NTOQ(gn,1,gcd);  
                         divsn(SAFENM(dc),gn,&tn); NTOQ(tn,SGN(dc),dcq);  
                         divsn(SAFENM(c),gn,&tn); NTOQ(tn,SGN(c),cq);  
                         muldc(CO,dp,(P)dcq,&t1);  
                         m = BDY(rp); NEWMP(mr); mr->dl = m->dl; mr->c = (P)cq;  
                         NEXT(mr) = 0; MKDP(rp->nv,mr,t); t->sugar = mr->dl->td;  
                         addd(CO,t,t1,&dp);  
                         dc = gcd;  
   
                         dp_rest(rp,&t); rp = t;  
                         tt = get_rtime()-t_0; t_a += tt;  
                 }  
         }  
         if ( GenTrace ) {  
                 mulq(ARG3(BDY((LIST)BDY(TraceList))),dc,&cq);  
                 ARG3(BDY((LIST)BDY(TraceList))) = (pointer)cq;  
         }  
         if ( dp )  
                 dp->sugar = sugar;  
         *r = dp;  
 final:  
         if ( DP_Print )  
                 fprintf(asir_out,"(%.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g)",  
                         t_p,t_m,t_m1,t_m2,t_s,  
                         t_g,t_a,  
                         pz_t_e, pz_t_d, pz_t_d1, pz_t_c);  
 }  }
   
 void dp_imul_d();  NODE /* of index */ remove_reducibles(NODE /* of index */ nd, int newdp )
 void imulv();  
   
 void _dp_nf_ptozp_mpi(b,g,ps,full,multiple,r)  
 NODE b;  
 DP g;  
 DP *ps;  
 int full,multiple;  
 DP *r;  
 {  {
         Obj dmy;    register DL dl, dln;
         DP u,dp,rp,t,t1,t2,red,shift;    register NODE last, p, head;
         Q dc,rc,dcq,rcq,cont,hr,hred,cr,cred,mcred,c,gcd,cq;    register int td;
         N gn,tn,cn;  
         NODE l,n0,n1,n2,n3;  
         MP m,mr;  
         int i,n;  
         int hmag,denom;  
         int sugar,psugar;  
         NODE dist;  
         STRING imul;  
         int ndist;  
         int kara_bit;  
         extern int kara_mag;  
         extern int mpi_mag;  
         Q ind;  
         int id;  
         double get_rtime();  
         double t_0,t_00,tt,ttt,t_p,t_m,t_m1,t_m2,t_s,t_g,t_a;  
   
         if ( !g ) {    dl = psh[newdp];
                 *r = 0; return;    td = dl->td;
         }    for ( head = last = 0, p = nd; p; ) {
         pz_t_e = pz_t_d = pz_t_d1 = pz_t_c = 0;      dln = psh[(int)BDY(p)];
         t_p = t_m = t_m1 = t_m2 = t_s = t_g = t_a = 0;      if ( dln->td >= td && dl_redble( dln, dl ) ) {
         p = NEXT(p);
         if ( last ) NEXT(last) = p;
       } else {
         if ( !last ) head = p;
         p = NEXT(last = p);
       }
     }
     return head;
   }
   
         denom = Denominator?Denominator:1;  int dl_redble(DL dl1,DL dl2)
         hmag = multiple*HMAG(g)/denom;  {
         kara_bit = kara_mag*27; /* XXX */    register int n, *d1, *d2;
         if ( Dist ) {  
                 dist = BDY(Dist);  
                 ndist = length(dist);  
         }  
         sugar = g->sugar;  
   
         dc = 0; dp = 0; rc = ONE; rp = g;    for ( d1 = dl1->d, d2 = dl2->d, n = CNVars; --n >= 0; d1++, d2++ )
         MKSTR(imul,"dp_imul_index");      if ( *d1 < *d2 ) return 0;
     return 1;
   }
   
         for ( ; rp; ) {  #if 0
                 for ( u = 0, l = b; l; l = NEXT(l) ) {  int dl_weyl_weight(DL dl)
                         if ( dl_redble(BDY(rp)->dl,psh[(int)BDY(l)]) ) {  {
                                 t_0 = get_rtime();    int n,w,i;
                                 dp_load((int)BDY(l),&red);  
                                 hr = (Q)BDY(rp)->c; hred = (Q)BDY(red)->c;  
                                 gcdn(NM(hr),NM(hred),&gn);  
                                 divsn(NM(hred),gn,&tn); NTOQ(tn,SGN(hred),cr);  
                                 divsn(NM(hr),gn,&tn); NTOQ(tn,SGN(hr),cred);  
                                 tt = get_rtime()-t_0; t_p += tt;  
   
                                 t_0 = get_rtime();    n = CNVars/2;
                                 dp_subd(rp,red,&shift);    for ( i = 0, w = 0; i < n; i++ )
                                 t_00 = get_rtime();      w += (-dl->d[i]+dl->d[n+i]);
                                 if ( Dist && ndist    return w;
                                         && HMAG(rp) > mpi_mag && p_mag((P)cr) > mpi_mag ) {  }
                                         if ( DP_Print ) fprintf(asir_out,"~");  #endif
                                         dp_imul_d(rp,cr,&t);  
                                 } else {  
                                         if ( DP_Print ) fprintf(asir_out,"_");  
                                         muldc(CO,rp,(P)cr,&t);  
                                 }  
                                 ttt = get_rtime()-t_00; t_m1 += ttt/dp_nt(rp);  
   
                                 t_00 = get_rtime();  int gbcheck(NODE f)
                                 chsgnp((P)cred,(P *)&mcred);  {
                                 if ( Dist && ndist    int i;
                                         && HMAG(red) > mpi_mag && p_mag((P)mcred) > mpi_mag ) {    NODE r,g,gall;
                                         if ( DP_Print ) fprintf(asir_out,"=");    DP_pairs d,l;
                                         dp_imul_d(red,mcred,&t1);    DP h,nf,dp1,dp2;
                                 } else {    struct oEGT tmp0,tmp1;
                                         if ( DP_Print ) fprintf(asir_out,"_");  
                                         muldc(CO,red,(P)mcred,&t1);  
                                 }  
                                 ttt = get_rtime()-t_00; t_m2 += ttt/dp_nt(red);  
   
                                 t_00 = get_rtime();    if ( NoGC )
                                 muld(CO,shift,t1,&t2);      return 1;
                                 addd(CO,t,t2,&u);    for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {
                                 tt = get_rtime(); t_m += tt-t_0;      i = (int)BDY(r);
                                 ttt = get_rtime(); t_s += ttt-t_00;      d = updpairs(d,g,i);
       g = updbase(g,i);
                                 psugar = (BDY(rp)->dl->td - BDY(red)->dl->td) + red->sugar;      gall = append_one(gall,i);
                                 sugar = MAX(sugar,psugar);    }
                                 if ( !u ) {    if ( DP_Print || DP_PrintShort ) {
                                         if ( dp )      fprintf(asir_out,"gbcheck total %d pairs\n",DPPlength(d)); fflush(asir_out);
                                                 dp->sugar = sugar;    }
                                         *r = dp;    while ( d ) {
                                         goto final;      l = d; d = NEXT(d);
                                 }      get_eg(&tmp0);
                                 break;      dp_load(l->dp1,&dp1); dp_load(l->dp2,&dp2);
                         }      dp_sp(dp1,dp2,&h);
                 }  /* fprintf(stderr,"{%d,%d}",l->dp1,l->dp2); */
                 if ( u ) {      _dp_nf(gall,h,ps,1,&nf);
                         if ( HMAG(u) > hmag ) {      get_eg(&tmp1); add_eg(&eg_gc,&tmp0,&tmp1);
                                 t_0 = get_rtime();      if ( DP_Print || DP_PrintShort ) {
                                 if ( Dist && HMAG(u) > mpi_mag ) {        fprintf(asir_out,"."); fflush(asir_out);
                                         if ( DP_Print )      }
                                                 fprintf(asir_out,"D");      if ( nf )
                                         dp_ptozp_d(dist,ndist,u,&t);        return 0;
                                 } else {    }
                                         if ( DP_Print )    if ( DP_Print || DP_PrintShort )
                                                 fprintf(asir_out,"L");      fprintf(asir_out,"\n");
                                         dp_ptozp_d(0,0,u,&t);    return 1;
                                 }  
                                 tt = get_rtime()-t_0; t_g += tt;  
                                 t_0 = get_rtime();  
                                 divsn(NM((Q)BDY(u)->c),NM((Q)BDY(t)->c),&cn); NTOQ(cn,1,cont);  
                                 if ( DP_Print ) {  
                                         fprintf(asir_out,"(%d)",p_mag((P)cont)*100/p_mag((P)BDY(u)->c));  
                                         fflush(asir_out);  
                                 }  
                                 mulq(cr,dc,&dcq);  
                                 mulq(cont,rc,&rcq);  
                                 gcdn(SAFENM(dcq),SAFENM(rcq),&gn);  
                                 divsn(SAFENM(dcq),gn,&tn); NTOQ(tn,SGN(dcq),dc);  
                                 divsn(SAFENM(rcq),gn,&tn); NTOQ(tn,SGN(rcq),rc);  
                                 tt = get_rtime()-t_0; t_a += tt;  
                                 rp = t;  
                                 hmag = multiple*HMAG(rp)/denom;  
                         } else {  
                                 t_0 = get_rtime();  
                                 mulq(cr,dc,&dcq); dc = dcq;  
                                 tt = get_rtime()-t_0; t_a += tt;  
                                 rp = u;  
                         }  
                 } else if ( !full ) {  
                         if ( rp ) {  
                                 MKDP(rp->nv,BDY(rp),t); t->sugar = sugar; rp = t;  
                         }  
                         *r = rp;  
                         goto final;  
                 } else {  
                         t_0 = get_rtime();  
                         mulq((Q)BDY(rp)->c,(Q)rc,(Q *)&c);  
                         gcdn(SAFENM(dc),SAFENM(c),&gn); NTOQ(gn,1,gcd);  
                         divsn(SAFENM(dc),gn,&tn); NTOQ(tn,SGN(dc),dcq);  
                         divsn(SAFENM(c),gn,&tn); NTOQ(tn,SGN(c),cq);  
                         muldc(CO,dp,(P)dcq,&t1);  
                         m = BDY(rp); NEWMP(mr); mr->dl = m->dl; mr->c = (P)cq;  
                         NEXT(mr) = 0; MKDP(rp->nv,mr,t); t->sugar = mr->dl->td;  
                         addd(CO,t,t1,&dp);  
                         dc = gcd;  
   
                         dp_rest(rp,&t); rp = t;  
                         tt = get_rtime()-t_0; t_a += tt;  
                 }  
         }  
         if ( dp )  
                 dp->sugar = sugar;  
         *r = dp;  
 final:  
         if ( DP_Print )  
                 fprintf(asir_out,"(%.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g)",  
                         t_p,t_m,t_m1,t_m2,t_s,  
                         t_g,t_a,  
                         pz_t_e, pz_t_d, pz_t_d1, pz_t_c);  
 }  }
   
 void dp_imul_d(p,q,rp)  void gbcheck_list(NODE f,int n,VECT *gp,LIST *pp)
 DP p;  
 Q q;  
 DP *rp;  
 {  {
         int nsep,ndist,i,j,k,l,n;    int i;
         double t0,t1,t2;    NODE r,g,gall,u,u0,t;
         Q *s;    VECT vect;
         pointer *b;    LIST pair;
         VECT c,cs,ri;    DP_pairs d,l;
         VECT *r;    Q q1,q2;
         MP m;  
         NODE tn,dist,n0,n1,n2;  
         Obj dmy;  
         STRING imul;  
   
         extern LIST Dist;    /* we need the following settings */
     NVars = CNVars = n;
     setup_arrays(f,0,&r);
     for ( gall = g = 0, d = 0; r; r = NEXT(r) ) {
       i = (int)BDY(r);
       d = updpairs(d,g,i);
       g = updbase(g,i);
       gall = append_one(gall,i);
     }
     NEWVECT(vect); vect->len = psn; vect->body = (pointer)ps;
     *gp = vect;
   
         if ( !p || !q ) {    for ( u0 = 0, l = d; l; l = NEXT(l) ) {
                 *rp = 0; return;      NEXTNODE(u0,u);
         }      STOQ(l->dp1,q1);
         dist = BDY(Dist);      STOQ(l->dp2,q2);
         for ( tn = dist, ndist = 0; tn; tn = NEXT(tn), ndist++ );      t = mknode(2,q1,q2);
         nsep = ndist + 1;      MKLIST(pair,t);
         for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ );      BDY(u) = (pointer)pair;
         if ( n <= nsep ) {    }
                 muldc(CO,p,(P)q,rp); return;    if ( u0 )
         }      NEXT(u) = 0;
         MKSTR(imul,"imulv");    MKLIST(*pp,u0);
         t0 = get_rtime();  
         dp_dtov(p,&c);  
         sepvect(c,nsep,&cs);  
         r = (VECT *)CALLOC(nsep,sizeof(VECT *));  
         for ( i = 0, tn = dist, b = BDY(cs); i < ndist; i++, tn = NEXT(tn) ) {  
                 n0 = mknode(4,BDY(tn),imul,b[i],q);  
                 Pox_rpc(n0,&dmy);  
         }  
         t1 = get_rtime();  
         im_t_s += t1 - t0;  
         imulv(b[i],q,&r[i]);  
         t1 = get_rtime();  
         for ( i = 0, tn = dist; i < ndist; i++, tn = NEXT(tn) ) {  
                 MKNODE(n0,BDY(tn),0);  
                 Pox_pop_local(n0,&r[i]);  
                 if ( OID(r[i]) == O_ERR ) {  
                         printexpr(CO,(Obj)r[i]);  
                         error("dp_imul_d : aborted");  
                 }  
         }  
         t2 = get_rtime();  
         im_t_r += t2 - t1;  
         s = (Q *)CALLOC(n,sizeof(Q));  
         for ( i = j = 0; i < nsep; i++ ) {  
                 for ( k = 0, ri = r[i], l = ri->len; k < l; k++, j++ ) {  
                         s[j] = (Q)BDY(ri)[k];  
                 }  
         }  
         dp_vtod(s,p,rp);  
 }  }
   
 void imulv(w,c,rp)  int membercheck(NODE f,NODE x)
 VECT w;  
 Q c;  
 VECT *rp;  
 {  {
         int n,i;    DP g;
         VECT r;    struct oEGT tmp0,tmp1;
   
         n = w->len;    if ( NoMC )
         MKVECT(r,n); *rp = r;      return 1;
         for ( i = 0; i < n; i++ )    if ( DP_Print || DP_PrintShort ) {
                 mulq((Q)BDY(w)[i],(Q)c,(Q *)&BDY(r)[i]);      fprintf(asir_out,"membercheck\n"); fflush(asir_out);
     }
     for ( ; f; f = NEXT(f) ) {
       get_eg(&tmp0);
       _dp_nf(x,(DP)BDY(f),ps,1,&g);
       get_eg(&tmp1); add_eg(&eg_mc,&tmp0,&tmp1);
       if ( DP_Print ) {
         print_split_eg(&tmp0,&tmp1); fflush(asir_out);
       } else if ( DP_PrintShort ) {
         fprintf(asir_out,"."); fflush(asir_out);
       }
       if ( g )
         return 0;
     }
     if ( DP_Print || DP_PrintShort )
       fprintf(asir_out,"\n");
     return 1;
 }  }
   
 DP_pairs updpairs( d, g, t)  void dp_set_flag(Obj name,Obj value)
 DP_pairs d;  
 NODE /* of index */ g;  
 int t;  
 {  {
         register DP_pairs d1, dd, nd;    char *n;
         struct oEGT tup0,tup1;    int v;
         int dl,dl1;    Q ratio;
   
         if ( !g ) return d;    if ( OID(name) != O_STR )
         get_eg(&tup0);      return;
         if ( !NoCriB && d ) {    n = BDY((STRING)name);
                 dl = DPPlength(d);    if ( !strcmp(n,"Demand") ) {
                 d = criterion_B( d, t );      Demand = value ? BDY((STRING)value) : 0; return;
                 dl -= DPPlength(d); NBP += dl;    }
         }    if ( !strcmp(n,"Dist") ) {
         d1 = newpairs( g, t );      Dist = (LIST)value; return;
         if ( NEXT(d1) ) {    }
                 dl = DPPlength(d1); TP += dl;    if ( !strcmp(n,"Content") ) {
                 d1 = criterion_M( d1 );      ratio = (Q)value;
                 dl1 = DPPlength(d1); NMP += (dl-dl1); dl = dl1;      if ( ratio ) {
                 d1 = criterion_F( d1 );        DP_Multiple = BD(NM(ratio))[0];
                 dl1 = DPPlength(d1); NFP += (dl-dl1); dl = dl1;        Denominator = INT(ratio) ? 1 : BD(DN(ratio))[0];
         } else      } else {
                 dl = 1;        DP_Multiple = 0;
         if ( !do_weyl )        Denominator = 1;
                 for ( dd = 0; d1; d1 = nd ) {      }
                         nd = NEXT(d1);    }
                         if ( !criterion_2( d1->dp1, d1->dp2 ) ) {    if ( value && OID(value) != O_N )
                                 NEXT(d1) = dd;      return;
                                 dd = d1;    v = QTOS((Q)value);
                         }    if ( !strcmp(n,"NoSugar") )
                 }      NoSugar = v;
         else    else if ( !strcmp(n,"NoCriB") )
                 dd = d1;      NoCriB = v;
         dl1 = DPPlength(dd); NDP += (dl-dl1);    else if ( !strcmp(n,"NoGC") )
         get_eg(&tup1);      NoGC = v;
         add_eg(&eg_up,&tup0,&tup1);    else if ( !strcmp(n,"NoMC") )
         if ( !(nd = d) ) return dd;      NoMC = v;
         while ( nd = NEXT(d1 = nd) ) ;    else if ( !strcmp(n,"NoRA") )
         NEXT(d1) = dd;      NoRA = v;
         return d;    else if ( !strcmp(n,"NoGCD") )
       NoGCD = v;
     else if ( !strcmp(n,"Top") )
       Top = v;
     else if ( !strcmp(n,"ShowMag") )
       ShowMag = v;
     else if ( !strcmp(n,"PrintShort") )
       DP_PrintShort = v;
     else if ( !strcmp(n,"Print") )
       DP_Print = v;
     else if ( !strcmp(n,"NFStat") )
       DP_NFStat = v;
     else if ( !strcmp(n,"Stat") )
       Stat = v;
     else if ( !strcmp(n,"Reverse") )
       Reverse = v;
     else if ( !strcmp(n,"Multiple") )
       DP_Multiple = v;
     else if ( !strcmp(n,"Denominator") )
       Denominator = v;
     else if ( !strcmp(n,"PtozpRA") )
       PtozpRA = v;
     else if ( !strcmp(n,"GenTrace") )
       GenTrace = v;
     else if ( !strcmp(n,"OXCheck") )
       OXCheck = v;
     else if ( !strcmp(n,"GenSyz") )
       GenSyz = v;
     else if ( !strcmp(n,"OneZeroHomo") )
       OneZeroHomo = v;
     else if ( !strcmp(n,"MaxDeg") )
       MaxDeg = v;
 }  }
   
 DP_pairs newpairs( g, t )  void dp_make_flaglist(LIST *list)
 NODE /* of index */ g;  
 register int t;  
 {  {
         register NODE r;    Q v,nm,dn;
         register DL tdl = psh[t];    STRING name,path;
         register int ts;    NODE n,n1;
         register DP_pairs p, last;  
         int dp;  
         register DL dl;  
         register int s;  
   
         ts = pss[t] - tdl->td;  
         for ( last = 0, r = g; r; r = NEXT(r) ) {  
                 NEXT(p = NEWDP_pairs) = last;  
                 last = p;  
                 dp = p->dp1 = (int)BDY(r);  p->dp2 = t;  
                 p->lcm = lcm_of_DL(CNVars, dl = psh[dp], tdl, (DL)0 );  
 #if 0  #if 0
                 if ( do_weyl )    STOQ(DP_Multiple,v); MKNODE(n,v,0); MKSTR(name,"DP_Multiple"); MKNODE(n1,name,n); n = n1;
                         p->sugar = dl_weight(p->lcm);    STOQ(Denominator,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Denominator"); MKNODE(n1,name,n); n = n1;
                 else  #else
     if ( DP_Multiple ) {
       STOQ(DP_Multiple,nm); STOQ(Denominator,dn); divq(nm,dn,&v);
     } else
       v = 0;
     MKNODE(n,v,0); MKSTR(name,"Content"); MKNODE(n1,name,n); n = n1;
 #endif  #endif
                         p->sugar = (ts > (s = pss[dp] - dl->td) ? ts : s) + p->lcm->td;    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;
         return last;    STOQ(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;
     STOQ(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;
     STOQ(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;
     STOQ(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;
     STOQ(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;
     STOQ(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;
     STOQ(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;
     STOQ(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;
     STOQ(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;
     if ( Demand )
       MKSTR(path,Demand);
     else
       path = 0;
     MKNODE(n1,path,n); n = n1; MKSTR(name,"Demand"); MKNODE(n1,name,n); n = n1;
     MKLIST(*list,n);
 }  }
   
 DP_pairs criterion_B( d, s )  #define DELIM '/'
 register DP_pairs d;  
 int s;  void dp_save(int index,Obj p,char *prefix)
 {  {
         register DP_pairs dd, p;    FILE *fp;
         register DL tij, t = psh[s], dltmp;    char path[BUFSIZ];
   
         if ( !d ) return 0;    if ( prefix )
         NEWDL( dltmp, CNVars );      sprintf(path,"%s%c%s%d",Demand,DELIM,prefix,index);
         for ( dd = 0; d; d = p ) {    else
                 p = NEXT(d),      sprintf(path,"%s%c%d",Demand,DELIM,index);
                 tij = d->lcm;    if ( !(fp = fopen(path,"wb") ) )
                 if ( tij->td != lcm_of_DL(CNVars, tij, t, dltmp )->td      error("dp_save : cannot open a file");
                         || !dl_equal(CNVars, tij, dltmp )    savevl(fp,VC); saveobj(fp,p); fclose(fp);
                         || (tij->td == lcm_of_DL(CNVars, psh[d->dp1], t, dltmp )->td  
                             && dl_equal(CNVars, dltmp, tij ))  
                         || (tij->td == lcm_of_DL(CNVars, psh[d->dp2], t, dltmp )->td  
                             && dl_equal(CNVars, dltmp, tij )) ) {  
                         NEXT(d) = dd;  
                         dd = d;  
                 }  
         }  
         return dd;  
 }  }
   
 DP_pairs criterion_M( d1 )  void dp_load(int index,DP *p)
 DP_pairs d1;  
 {  {
         register DP_pairs dd, e, d3, d2, p;    FILE *fp;
         register DL itdl, jtdl;    char path[BUFSIZ];
         register int itdltd, jtdltd;  
   
         for ( dd = 0, e = d1; e; e = d3 ) {    if ( !Demand || ps[index] )
                 if ( !(d2 = NEXT(e)) ) {      *p = ps[index];
                         NEXT(e) = dd;    else {
                         return e;      sprintf(path,"%s%c%d",Demand,DELIM,index);
                 }      if ( !(fp = fopen(path,"rb") ) )
                 itdltd = (itdl = e->lcm)->td;        error("dp_load : cannot open a file");
                 for ( d3 = 0; d2; d2 = p ) {      if ( PCoeffs )
                         p = NEXT(d2),        loadvl(fp);
                         jtdltd = (jtdl = d2->lcm)->td;      else
                         if ( jtdltd == itdltd  )        skipvl(fp);
                                 if ( dl_equal(CNVars, itdl, jtdl ) ) ;      loadobj(fp,(Obj *)p); fclose(fp);
                                 else if ( dl_redble( jtdl, itdl ) ) continue;    }
                                 else if ( dl_redble( itdl, jtdl ) ) goto delit;  
                                 else ;  
                         else if ( jtdltd > itdltd )  
                                 if ( dl_redble( jtdl, itdl ) ) continue;  
                                 else ;  
                         else if ( dl_redble( itdl, jtdl ) ) goto delit;  
                         NEXT(d2) = d3;  
                         d3 = d2;  
                 }  
                 NEXT(e) = dd;  
                 dd = e;  
                 continue;  
                 /**/  
         delit:  NEXT(d2) = d3;  
                 d3 = d2;  
                 for ( ; p; p = d2 ) {  
                         d2 = NEXT(p);  
                         NEXT(p) = d3;  
                         d3 = p;  
                 }  
         }  
         return dd;  
 }  }
   
 static DP_pairs collect_pairs_of_hdlcm( d1, prest )  int dp_load_t(int index,DP *p)
 DP_pairs d1, *prest;  
 {  {
         register DP_pairs w, p, r, s;    FILE *fp;
         register DL ti;    char path[BUFSIZ];
         register int td;  
   
         td = (ti = (w = d1)->lcm)->td;    sprintf(path,"%s%c%d",Demand,DELIM,index);
         s = NEXT(w);    if ( !(fp = fopen(path,"rb") ) )
         NEXT(w) = 0;      return 0;
         for ( r = 0; s; s = p ) {    else {
                 p = NEXT(s);      if ( PCoeffs )
                 if ( td == s->lcm->td && dl_equal(CNVars, ti, s->lcm ) )        loadvl(fp);
                 {      else
                         NEXT(s) = w;        skipvl(fp);
                         w = s;      loadobj(fp,(Obj *)p); fclose(fp); return 1;
                 } else {    }
                         NEXT(s) = r;  
                         r = s;  
                 }  
         }  
         *prest = r;  
         return w;  
 }  }
   
 int criterion_2( dp1, dp2 )  void init_stat() {
 int dp1, dp2;    init_eg(&eg_nf); init_eg(&eg_nfm); init_eg(&eg_znfm);
 {    init_eg(&eg_pz); init_eg(&eg_np);
         register int i, *d1, *d2;    init_eg(&eg_ra); init_eg(&eg_mc); init_eg(&eg_gc);
     ZR = NZR = TP = NMP = N_BP = NFP = NDP = 0;
   }
   
         d1 = psh[dp1]->d,  d2 = psh[dp2]->d;  void print_stat() {
         for ( i = CNVars; --i >= 0; d1++, d2++ )    if ( !DP_Print && !Stat )
                 if ( (*d1 <= *d2 ? *d1 : *d2) > 0  ) return 0;      return;
         return 1;    print_eg("NF",&eg_nf); print_eg("NFM",&eg_nfm); print_eg("ZNFM",&eg_znfm);
     print_eg("PZ",&eg_pz); print_eg("NP",&eg_np);
     print_eg("RA",&eg_ra); print_eg("MC",&eg_mc); print_eg("GC",&eg_gc);
     fprintf(asir_out,"T=%d,B=%d M=%d F=%d D=%d ZR=%d NZR=%d\n",TP,N_BP,NMP,NFP,NDP,ZR,NZR);
 }  }
   
 DP_pairs criterion_F( d1 )  /*
 DP_pairs d1;   * dp_nf used in gb()
 {   *
         DP_pairs rest, head;   */
         register DP_pairs last, p, r, w;  
         register int s;  
   
         for ( head = last = 0, p = d1; NEXT(p); ) {  double pz_t_e, pz_t_d, pz_t_d1, pz_t_c, im_t_s, im_t_r;
                 s = (r = w = collect_pairs_of_hdlcm( p, &rest ))->sugar;  
                 while ( w = NEXT(w) )  
                         if ( criterion_2( w->dp1, w->dp2 ) ) {  
                                 r = w;  
                                 break;  
                         } else if ( w->sugar < s ) s = (r = w)->sugar;  
                 if ( last ) NEXT(last) = r;  
                 else head = r;  
                 NEXT(last = r) = 0;  
                 if ( !(p = rest) ) return head;  
         }  
         if ( !last ) return p;  
         NEXT(last) = p;  
         return head;  
 }  
   
 NODE updbase(g,t)  extern int GenTrace;
 NODE g;  extern NODE TraceList;
 int t;  extern int mpi_mag;
 {  
         g = remove_reducibles(g,t);  
         g = append_one(g,t);  
         return g;  
 }  
   
 NODE /* of index */ remove_reducibles( nd, newdp )  void dp_mulc_d(DP p,P c,DP *r)
 NODE /* of index */ nd;  
 int newdp;  
 {  {
         register DL dl, dln;    if ( Dist && BDY(Dist)
         register NODE last, p, head;      && HMAG(p) > mpi_mag
         register int td;      && p_mag((P)c) > mpi_mag ) {
       if ( DP_NFStat ) fprintf(asir_out,"~");
         dl = psh[newdp];      dp_imul_d(p,(Q)c,r);
         td = dl->td;    } else {
         for ( head = last = 0, p = nd; p; ) {      if ( DP_NFStat ) fprintf(asir_out,"_");
                 dln = psh[(int)BDY(p)];      muldc(CO,p,(Obj)c,r);
                 if ( dln->td >= td && dl_redble( dln, dl ) ) {    }
                         p = NEXT(p);  
                         if ( last ) NEXT(last) = p;  
                 } else {  
                         if ( !last ) head = p;  
                         p = NEXT(last = p);  
                 }  
         }  
         return head;  
 }  }
   
 int dl_redble( dl1, dl2 )  void _dp_nf(NODE b,DP g,DP *ps,int full,DP *rp)
 DL dl1, dl2;  
 {  {
         register int n, *d1, *d2;    DP u,p,d,s,t,mult;
     P coef;
     NODE l;
     MP m,mr;
     int sugar,psugar;
   
         for ( d1 = dl1->d, d2 = dl2->d, n = CNVars; --n >= 0; d1++, d2++ )    if ( !g ) {
                 if ( *d1 < *d2 ) return 0;      *rp = 0; return;
         return 1;    }
     sugar = g->sugar;
     for ( d = 0; g; ) {
       for ( u = 0, l = b; l; l = NEXT(l) ) {
         if ( dl_redble(BDY(g)->dl,psh[(int)BDY(l)]) ) {
           dp_load((int)BDY(l),&p);
           /* t+u = coef*(d+g) - mult*p (t = coef*d) */
           dp_red(d,g,p,&t,&u,&coef,&mult);
           psugar = (BDY(g)->dl->td - BDY(p)->dl->td) + p->sugar;
           sugar = MAX(sugar,psugar);
           if ( GenTrace ) {
             LIST hist;
             Q cq;
             NODE node,node0;
   
             STOQ((int)BDY(l),cq);
             node0 = mknode(4,coef,cq,mult,ONE);
             MKLIST(hist,node0);
             MKNODE(node,hist,TraceList); TraceList = node;
           }
           if ( !u ) {
             if ( d )
               d->sugar = sugar;
             *rp = d; return;
           }
           d = t;
           break;
         }
       }
       if ( u )
         g = u;
       else if ( !full ) {
         if ( g ) {
           MKDP(g->nv,BDY(g),t); t->sugar = sugar; g = t;
         }
         *rp = g; return;
       } else {
         m = BDY(g); NEWMP(mr); mr->dl = m->dl; mr->c = m->c;
         NEXT(mr) = 0; MKDP(g->nv,mr,t); t->sugar = mr->dl->td;
         addd(CO,d,t,&s); d = s;
         dp_rest(g,&t); g = t;
       }
     }
     if ( d )
       d->sugar = sugar;
     *rp = d;
 }  }
   
 int dl_weight(dl)  void _dp_nf_z(NODE b,DP g,DP *ps,int full,int multiple,DP *r)
 DL dl;  
 {  {
         int n,w,i;    DP u,dp,rp,t,t1,red,shift;
     Q dc,rc,dcq,rcq,cont,hr,hred,cr,cred,mcred,c,gcd,cq;
     NODE l;
     int hmag,denom;
     int sugar,psugar;
     STRING imul;
     double t_0,tt,t_p,t_m,t_g,t_a;
     LIST hist;
     NODE node;
     Q rcred,mrcred;
   
         n = CNVars/2;    if ( !g ) {
         for ( i = 0, w = 0; i < n; i++ )      *r = 0; return;
                 w += (-dl->d[i]+dl->d[n+i]);    }
         return w;    pz_t_e = pz_t_d = pz_t_d1 = pz_t_c = 0;
 }    t_p = t_m = t_g = t_a = 0;
   
 int gbcheck(f)    denom = Denominator?Denominator:1;
 NODE f;    hmag = multiple*HMAG(g)/denom;
 {    sugar = g->sugar;
         int i;  
         NODE r,g,gall;  
         DP_pairs d,l;  
         DP h,nf,dp1,dp2;  
         struct oEGT tmp0,tmp1;  
   
         if ( NoGC )    dc = 0; dp = 0; rc = ONE; rp = g;
                 return 1;    MKSTR(imul,"dp_imul_index");
         for ( gall = g = 0, d = 0, r = f; r; r = NEXT(r) ) {  
                 i = (int)BDY(r);  
                 d = updpairs(d,g,i);  
                 g = updbase(g,i);  
                 gall = append_one(gall,i);  
         }  
         if ( DP_Print || DP_PrintShort ) {  
                 fprintf(asir_out,"gbcheck total %d pairs\n",DPPlength(d)); fflush(asir_out);  
         }  
         while ( d ) {  
                 l = d; d = NEXT(d);  
                 get_eg(&tmp0);  
                 dp_load(l->dp1,&dp1); dp_load(l->dp2,&dp2); dp_sp(dp1,dp2,&h);  
                 _dp_nf(gall,h,ps,1,&nf);  
                 get_eg(&tmp1); add_eg(&eg_gc,&tmp0,&tmp1);  
                 if ( DP_Print || DP_PrintShort ) {  
                         fprintf(asir_out,"."); fflush(asir_out);  
                 }  
                 if ( nf )  
                         return 0;  
         }  
         if ( DP_Print || DP_PrintShort )  
                 fprintf(asir_out,"\n");  
         return 1;  
 }  
   
 int membercheck(f,x)    /* g = dc*dp+rc*rp */
 NODE f,x;    for ( ; rp; ) {
 {      for ( u = 0, l = b; l; l = NEXT(l) ) {
         DP g;        if ( dl_redble(BDY(rp)->dl,psh[(int)BDY(l)]) ) {
         struct oEGT tmp0,tmp1;          t_0 = get_rtime();
           dp_load((int)BDY(l),&red);
           hr = (Q)BDY(rp)->c; hred = (Q)BDY(red)->c;
           igcd_cofactor((Q)BDY(rp)->c,(Q)BDY(red)->c,&gcd,&cred,&cr);
           tt = get_rtime(); t_p += tt-t_0;
   
         if ( NoMC )          dp_subd(rp,red,&shift);
                 return 1;          dp_mulc_d(rp,(P)cr,&t);
         if ( DP_Print || DP_PrintShort ) {          chsgnp((P)cred,(P *)&mcred);
                 fprintf(asir_out,"membercheck\n"); fflush(asir_out);          dp_mulc_d(red,(P)mcred,&t1);
         }          muld(CO,shift,t1,&t1);
         for ( ; f; f = NEXT(f) ) {          addd(CO,t,t1,&u);
                 get_eg(&tmp0);          t_m += get_rtime()-tt;
                 _dp_nf(x,(DP)BDY(f),ps,1,&g);  
                 get_eg(&tmp1); add_eg(&eg_mc,&tmp0,&tmp1);  
                 if ( DP_Print ) {  
                         print_split_eg(&tmp0,&tmp1); fflush(asir_out);  
                 } else if ( DP_PrintShort ) {  
                         fprintf(asir_out,"."); fflush(asir_out);  
                 }  
                 if ( g )  
                         return 0;  
         }  
         if ( DP_Print || DP_PrintShort )  
                 fprintf(asir_out,"\n");  
         return 1;  
 }  
   
 void dp_set_flag(name,value)          psugar = (BDY(rp)->dl->td - BDY(red)->dl->td) + red->sugar;
 Obj name,value;          sugar = MAX(sugar,psugar);
 {  
         char *n;  
         int v;  
   
         if ( OID(name) != O_STR )          if ( GenTrace ) {
                 return;            /* u = cr*rp + (-cred)*shift*red */
         n = BDY((STRING)name);            STOQ((int)BDY(l),cq);
         if ( !strcmp(n,"Demand") ) {            node = mknode(4,cr,cq,NULLP,NULLP);
                 Demand = value ? BDY((STRING)value) : 0; return;            mulq(cred,rc,&rcred);
         }            chsgnnum((Num)rcred,(Num *)&mrcred);
         if ( !strcmp(n,"Dist") ) {            muldc(CO,shift,(Obj)mrcred,(DP *)&ARG2(node));
                 Dist = (LIST)value; return;            MKLIST(hist,node);
         }          }
         if ( value && OID(value) != O_N )  
                 return;  
         v = QTOS((Q)value);  
         if ( !strcmp(n,"NoSugar") )  
                 NoSugar = v;  
         else if ( !strcmp(n,"NoCriB") )  
                 NoCriB = v;  
         else if ( !strcmp(n,"NoGC") )  
                 NoGC = v;  
         else if ( !strcmp(n,"NoMC") )  
                 NoMC = v;  
         else if ( !strcmp(n,"NoRA") )  
                 NoRA = v;  
         else if ( !strcmp(n,"NoGCD") )  
                 NoGCD = v;  
         else if ( !strcmp(n,"Top") )  
                 Top = v;  
         else if ( !strcmp(n,"ShowMag") )  
                 ShowMag = v;  
         else if ( !strcmp(n,"DP_PrintShort") )  
                 DP_PrintShort = v;  
         else if ( !strcmp(n,"DP_Print") )  
                 DP_Print = v;  
         else if ( !strcmp(n,"Stat") )  
                 Stat = v;  
         else if ( !strcmp(n,"Reverse") )  
                 Reverse = v;  
         else if ( !strcmp(n,"Multiple") )  
                 Multiple = v;  
         else if ( !strcmp(n,"Denominator") )  
                 Denominator = v;  
         else if ( !strcmp(n,"PtozpRA") )  
                 PtozpRA = v;  
         else if ( !strcmp(n,"GenTrace") )  
                 GenTrace = v;  
         else if ( !strcmp(n,"OXCheck") )  
                 OXCheck = v;  
 }  
   
 void dp_make_flaglist(list)          if ( !u ) {
 LIST *list;            if ( dp )
 {              dp->sugar = sugar;
         Q v;            *r = dp;
         STRING name,path;            if ( GenTrace ) {
         NODE n,n1;              ARG3(BDY(hist)) = ONE;
               MKNODE(node,hist,TraceList); TraceList = node;
             }
             goto final;
           }
           break;
         }
       }
       if ( u ) {
         if ( multiple && HMAG(u) > hmag ) {
           t_0 = get_rtime();
           dp_ptozp_d(u,&rp);
           tt = get_rtime(); t_g += tt-t_0;
   
         STOQ(Multiple,v); MKNODE(n,v,0); MKSTR(name,"Multiple"); MKNODE(n1,name,n); n = n1;          divsq((Q)BDY(u)->c,(Q)BDY(rp)->c,&cont);
         STOQ(Denominator,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Denominator"); MKNODE(n1,name,n); n = n1;          if ( !dp_fcoeffs && DP_NFStat ) {
         MKNODE(n1,Dist,n); n = n1; MKSTR(name,"Dist"); MKNODE(n1,name,n); n = n1;            fprintf(asir_out,
         STOQ(Reverse,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Reverse"); MKNODE(n1,name,n); n = n1;              "(%d)",p_mag((P)cont)*100/p_mag((P)BDY(u)->c));
         STOQ(Stat,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"Stat"); MKNODE(n1,name,n); n = n1;            fflush(asir_out);
         STOQ(DP_Print,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"DP_Print"); MKNODE(n1,name,n); n = n1;          }
         STOQ(OXCheck,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"OXCheck"); MKNODE(n1,name,n); n = n1;          mulq(cr,dc,&dcq); mulq(cont,rc,&rcq);
         STOQ(GenTrace,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"GenTrace"); MKNODE(n1,name,n); n = n1;          igcd_cofactor(dcq,rcq,&gcd,&dc,&rc);
         STOQ(PtozpRA,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"PtozpRA"); MKNODE(n1,name,n); n = n1;          t_a = get_rtime()-tt;
         STOQ(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;          hmag = multiple*HMAG(rp)/denom;
         STOQ(NoGCD,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoGCD"); MKNODE(n1,name,n); n = n1;          if ( GenTrace ) {
         STOQ(NoRA,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoRA"); MKNODE(n1,name,n); n = n1;            ARG3(BDY(hist)) = (pointer)gcd;
         STOQ(NoMC,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoMC"); MKNODE(n1,name,n); n = n1;            MKNODE(node,hist,TraceList); TraceList = node;
         STOQ(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;        } else {
         STOQ(NoSugar,v); MKNODE(n1,v,n); n = n1; MKSTR(name,"NoSugar"); MKNODE(n1,name,n); n = n1;          rp = u;
         if ( Demand )          t_0 = get_rtime();
                 MKSTR(path,Demand);          mulq(cr,dc,&dc);
         else          t_a += get_rtime()-t_0;
                 path = 0;          if ( GenTrace ) {
         MKNODE(n1,path,n); n = n1; MKSTR(name,"Demand"); MKNODE(n1,name,n); n = n1;            ARG3(BDY(hist)) = (pointer)ONE;
         MKLIST(*list,n);            MKNODE(node,hist,TraceList); TraceList = node;
           }
         }
       } else if ( !full ) {
         if ( rp ) {
           MKDP(rp->nv,BDY(rp),t); t->sugar = sugar; rp = t;
         }
         *r = rp;
         goto final;
       } else {
         t_0 = get_rtime();
         mulq((Q)BDY(rp)->c,rc,&c);
         igcd_cofactor(dc,c,&dc,&dcq,&cq);
         muldc(CO,dp,(Obj)dcq,&t);
         dp_hm(rp,&t1); BDY(t1)->c = (Obj)cq;  addd(CO,t,t1,&dp);
         dp_rest(rp,&rp);
         t_a += get_rtime()-t_0;
       }
     }
     if ( GenTrace ) {
       mulq(ARG3(BDY((LIST)BDY(TraceList))),dc,&cq);
       ARG3(BDY((LIST)BDY(TraceList))) = (pointer)cq;
     }
     if ( dp )
       dp->sugar = sugar;
     *r = dp;
   final:
     if ( DP_NFStat )
       fprintf(asir_out,
         "(%.3g %.3g %.3g %.3g %.3g %.3g %.3g %.3g)",
         t_p,t_m,t_g,t_a,
         pz_t_e, pz_t_d, pz_t_d1, pz_t_c);
 }  }
   
 #define DELIM '/'  void imulv();
   
 void dp_save(index,p,prefix)  void dp_imul_d(DP p,Q q,DP *rp)
 int index;  
 Obj p;  
 char *prefix;  
 {  {
         FILE *fp;    int nsep,ndist,i,j,k,l,n;
         char path[BUFSIZ];    double t0,t1,t2;
     Q *s;
     pointer *b;
     VECT c,cs,ri;
     VECT *r;
     MP m;
     NODE tn,dist,n0;
     Obj obj;
     STRING imul;
     extern LIST Dist;
   
         if ( prefix )    if ( !p || !q ) {
                 sprintf(path,"%s%c%s%d",Demand,DELIM,prefix,index);      *rp = 0; return;
         else    }
                 sprintf(path,"%s%c%d",Demand,DELIM,index);    dist = BDY(Dist);
         if ( !(fp = fopen(path,"wb") ) )    for ( tn = dist, ndist = 0; tn; tn = NEXT(tn), ndist++ );
                 error("dp_save : cannot open a file");    nsep = ndist + 1;
         savevl(fp,VC); saveobj(fp,p); fclose(fp);    for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ );
     if ( n <= nsep ) {
       muldc(CO,p,(Obj)q,rp); return;
     }
     MKSTR(imul,"imulv");
     t0 = get_rtime();
     dp_dtov(p,&c);
     sepvect(c,nsep,&cs);
     r = (VECT *)CALLOC(nsep,sizeof(VECT *));
     for ( i = 0, tn = dist, b = BDY(cs); i < ndist; i++, tn = NEXT(tn) ) {
       n0 = mknode(4,BDY(tn),imul,b[i],q);
       Pox_rpc(n0,&obj);
     }
     t1 = get_rtime();
     im_t_s += t1 - t0;
     imulv(b[i],q,&r[i]);
     t1 = get_rtime();
     for ( i = 0, tn = dist; i < ndist; i++, tn = NEXT(tn) ) {
       MKNODE(n0,BDY(tn),0);
       Pox_pop_local(n0,&obj); r[i] = (VECT)obj;
       if ( OID(r[i]) == O_ERR ) {
         printexpr(CO,(Obj)r[i]);
         error("dp_imul_d : aborted");
       }
     }
     t2 = get_rtime();
     im_t_r += t2 - t1;
     s = (Q *)CALLOC(n,sizeof(Q));
     for ( i = j = 0; i < nsep; i++ ) {
       for ( k = 0, ri = r[i], l = ri->len; k < l; k++, j++ ) {
         s[j] = (Q)BDY(ri)[k];
       }
     }
     dp_vtod(s,p,rp);
 }  }
   
 void dp_load(index,p)  void imulv(VECT w,Q c,VECT *rp)
 int index;  
 DP *p;  
 {  {
         FILE *fp;    int n,i;
         char path[BUFSIZ];    VECT r;
   
         if ( !Demand || ps[index] )    n = w->len;
                 *p = ps[index];    MKVECT(r,n); *rp = r;
         else {    for ( i = 0; i < n; i++ )
                 sprintf(path,"%s%c%d",Demand,DELIM,index);      mulq((Q)BDY(w)[i],(Q)c,(Q *)&BDY(r)[i]);
                 if ( !(fp = fopen(path,"rb") ) )  
                         error("dp_load : cannot open a file");  
                 skipvl(fp); loadobj(fp,(Obj *)p); fclose(fp);  
         }  
 }  }
   
 int dp_load_t(index,p)  void dptoca(DP p,unsigned int **rp)
 int index;  
 DP *p;  
 {  {
         FILE *fp;    int i;
         char path[BUFSIZ];    MP m;
     unsigned int *r;
   
         sprintf(path,"%s%c%d",Demand,DELIM,index);    if ( !p )
         if ( !(fp = fopen(path,"rb") ) )      *rp = 0;
                 return 0;    else {
         else {      for ( m = BDY(p), i = 0; m; m = NEXT(m), i++ );
                 skipvl(fp); loadobj(fp,(Obj *)p); fclose(fp); return 1;      *rp = r = (unsigned int *)MALLOC_ATOMIC(i*sizeof(unsigned int));
         }      for ( m = BDY(p), i = 0; m; m = NEXT(m), i++ )
         r[i] = ITOS(C(m));
     }
 }  }
   
 void init_stat() {  
         init_eg(&eg_up); init_eg(&eg_sp); init_eg(&eg_spm);  
         init_eg(&eg_nf); init_eg(&eg_nfm); init_eg(&eg_znfm);  
         init_eg(&eg_pz); init_eg(&eg_np); init_eg(&eg_mp);  
         init_eg(&eg_ra); init_eg(&eg_mc); init_eg(&eg_gc);  
         ZR = NZR = TP = NBP = NMP = NFP = NDP = 0;  
 }  
   
 void print_stat() {  
         if ( !DP_Print && !Stat )  
                 return;  
         print_eg("UP",&eg_up); print_eg("SP",&eg_sp); print_eg("SPM",&eg_spm);  
         print_eg("NF",&eg_nf); print_eg("NFM",&eg_nfm); print_eg("ZNFM",&eg_znfm);  
         print_eg("PZ",&eg_pz); print_eg("NP",&eg_np); print_eg("MP",&eg_mp);  
         print_eg("RA",&eg_ra); print_eg("MC",&eg_mc); print_eg("GC",&eg_gc);  
         fprintf(asir_out,"T=%d,B=%d M=%d F=%d D=%d ZR=%d NZR=%d\n",TP,NBP,NMP,NFP,NDP,ZR,NZR);  
 }  
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.74
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>