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

Diff for /OpenXM_contrib2/asir2018/engine/nd.c between version 1.21 and 1.66

version 1.21, 2019/09/19 06:29:48 version 1.66, 2022/04/03 06:45:47
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.20 2019/09/15 08:46:12 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2018/engine/nd.c,v 1.65 2022/04/03 00:39:12 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
 int Nnd_add,Nf4_red;  void print_siglist(NODE l);
   
   NODE nd_hpdata;
   int Nnd_add,Nf4_red,NcriB,NcriMF,Ncri2,Npairs,Nnewpair;
 struct oEGT eg_search,f4_symb,f4_conv,f4_elim1,f4_elim2;  struct oEGT eg_search,f4_symb,f4_conv,f4_elim1,f4_elim2;
   
 int diag_period = 6;  int diag_period = 6;
Line 18  NM _nm_free_list;
Line 21  NM _nm_free_list;
 ND _nd_free_list;  ND _nd_free_list;
 ND_pairs _ndp_free_list;  ND_pairs _ndp_free_list;
 NODE nd_hcf;  NODE nd_hcf;
   int Nsyz,Nsamesig;
   
 Obj nd_top_weight;  Obj nd_top_weight;
   
Line 49  static NDV *nd_ps_trace;
Line 53  static NDV *nd_ps_trace;
 static NDV *nd_ps_sym;  static NDV *nd_ps_sym;
 static NDV *nd_ps_trace_sym;  static NDV *nd_ps_trace_sym;
 static RHist *nd_psh;  static RHist *nd_psh;
 static int nd_psn,nd_pslen;  static int nd_psn,nd_pslen,nd_nbase;
 static RHist *nd_red;  static RHist *nd_red;
 static int *nd_work_vector;  static int *nd_work_vector;
 static int **nd_matrix;  static int **nd_matrix;
Line 64  static int nd_module_rank,nd_poly_weight_len;
Line 68  static int nd_module_rank,nd_poly_weight_len;
 static int *nd_poly_weight,*nd_module_weight;  static int *nd_poly_weight,*nd_module_weight;
 static NODE nd_tracelist;  static NODE nd_tracelist;
 static NODE nd_alltracelist;  static NODE nd_alltracelist;
 static int nd_gentrace,nd_gensyz,nd_nora,nd_newelim,nd_intersect,nd_lf;  static int nd_gentrace,nd_gensyz,nd_nora,nd_newelim,nd_intersect,nd_lf,nd_norb;
   static int nd_f4_td,nd_sba_f4step,nd_sba_pot,nd_sba_largelcm,nd_sba_dontsort,nd_sba_redundant_check;
   static int nd_top,nd_sba_syz,nd_sba_inputisgb,nd_sba_heu;
 static int *nd_gbblock;  static int *nd_gbblock;
 static NODE nd_nzlist,nd_check_splist;  static NODE nd_nzlist,nd_check_splist;
 static int nd_splist;  static int nd_splist;
 static int *nd_sugarweight;  static int *nd_sugarweight;
 static int nd_f4red,nd_rank0,nd_last_nonzero;  static int nd_f4red,nd_rank0,nd_last_nonzero;
   static DL *nd_sba_hm;
   static NODE *nd_sba_pos;
   
   struct comp_sig_spec {
     int n;
     // current_i <-> oldv[i]
     int *oldv;
     int *weight;
     struct order_pair *order_pair;
     int block_length;
     int **matrix;
     int row;
     int (*cmpdl)(int n,DL d1,DL d2);
   };
   
   struct comp_sig_spec *nd_sba_modord;
   
   DL ndltodl(int n,UINT *ndl);
 NumberField get_numberfield();  NumberField get_numberfield();
 UINT *nd_det_compute_bound(NDV **dm,int n,int j);  UINT *nd_det_compute_bound(NDV **dm,int n,int j);
 void nd_det_reconstruct(NDV **dm,int n,int j,NDV d);  void nd_det_reconstruct(NDV **dm,int n,int j,NDV d);
Line 88  P ndc_div(int mod,union oNDC a,union oNDC b);
Line 111  P ndc_div(int mod,union oNDC a,union oNDC b);
 P ndctop(int mod,union oNDC c);  P ndctop(int mod,union oNDC c);
 void finalize_tracelist(int i,P cont);  void finalize_tracelist(int i,P cont);
 void conv_ilist(int demand,int trace,NODE g,int **indp);  void conv_ilist(int demand,int trace,NODE g,int **indp);
 void parse_nd_option(NODE opt);  void parse_nd_option(VL vl,NODE opt);
 void dltondl(int n,DL dl,UINT *r);  void dltondl(int n,DL dl,UINT *r);
 DP ndvtodp(int mod,NDV p);  DP ndvtodp(int mod,NDV p);
 DP ndtodp(int mod,ND p);  DP ndtodp(int mod,ND p);
 DPM ndvtodpm(int mod,NDV p);  DPM ndvtodpm(int mod,NDV p);
   NDV dptondv(int mod,DP p);
 NDV dpmtondv(int mod,DPM p);  NDV dpmtondv(int mod,DPM p);
   int dp_getdeg(DP p);
 int dpm_getdeg(DPM p,int *rank);  int dpm_getdeg(DPM p,int *rank);
 void dpm_ptozp(DPM p,Z *cont,DPM *r);  void dpm_ptozp(DPM p,Z *cont,DPM *r);
 int compdmm(int nv,DMM a,DMM b);  int compdmm(int nv,DMM a,DMM b);
   DPM sigtodpm(SIG s);
   SIG dup_sig(SIG sig);
   
 void Pdp_set_weight(NODE,VECT *);  void Pdp_set_weight(NODE,VECT *);
 void Pox_cmo_rpc(NODE,Obj *);  void Pox_cmo_rpc(NODE,Obj *);
Line 111  NODE nd_f4_red_mod64_main(int m,ND_pairs sp0,int nsp,U
Line 138  NODE nd_f4_red_mod64_main(int m,ND_pairs sp0,int nsp,U
 NODE nd_f4_red_lf_main(int m,ND_pairs sp0,int nsp,int trace,UINT *s0vect,int col,  NODE nd_f4_red_lf_main(int m,ND_pairs sp0,int nsp,int trace,UINT *s0vect,int col,
         NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred);          NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred);
 int nd_gauss_elim_lf(mpz_t **mat0,int *sugar,int row,int col,int *colstat);  int nd_gauss_elim_lf(mpz_t **mat0,int *sugar,int row,int col,int *colstat);
   int nd_gauss_elim_mod_s(UINT **mat,int *sugar,ND_pairs *spactive,int row,int col,int md,int *colstat,SIG *sig);
 NODE nd_f4_lf_trace_main(int m,int **indp);  NODE nd_f4_lf_trace_main(int m,int **indp);
 void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp);  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp);
   
Line 150  void _NM_alloc()
Line 178  void _NM_alloc()
 matrix alloc_matrix(int row,int col)  matrix alloc_matrix(int row,int col)
 {  {
   unsigned long **a;    unsigned long **a;
   int i,len,blen;  
   matrix mat;    matrix mat;
   
   mat = (matrix)MALLOC(sizeof(struct matrix));    mat = (matrix)MALLOC(sizeof(struct matrix));
Line 355  void ndl_dehomogenize(UINT *d)
Line 382  void ndl_dehomogenize(UINT *d)
 void ndl_lcm(UINT *d1,unsigned *d2,UINT *d)  void ndl_lcm(UINT *d1,unsigned *d2,UINT *d)
 {  {
     UINT t1,t2,u,u1,u2;      UINT t1,t2,u,u1,u2;
     int i,j,l;      int i,j;
   
     if ( nd_module && (MPOS(d1) != MPOS(d2)) )      if ( nd_module && (MPOS(d1) != MPOS(d2)) )
         error("ndl_lcm : inconsistent monomials");          error("ndl_lcm : inconsistent monomials");
Line 453  void ndl_lcm(UINT *d1,unsigned *d2,UINT *d)
Line 480  void ndl_lcm(UINT *d1,unsigned *d2,UINT *d)
 void ndl_max(UINT *d1,unsigned *d2,UINT *d)  void ndl_max(UINT *d1,unsigned *d2,UINT *d)
 {  {
     UINT t1,t2,u,u1,u2;      UINT t1,t2,u,u1,u2;
     int i,j,l;      int i,j;
   
     for ( i = nd_exporigin; i < nd_wpd; i++ ) {      for ( i = nd_exporigin; i < nd_wpd; i++ ) {
         u1 = d1[i]; u2 = d2[i];          u1 = d1[i]; u2 = d2[i];
Line 493  int ndl_weight(UINT *d)
Line 520  int ndl_weight(UINT *d)
 int ndl_weight2(UINT *d)  int ndl_weight2(UINT *d)
 {  {
     int t,u;      int t,u;
     int i,j;      int i;
   
     for ( i = 0, t = 0; i < nd_nvar; i++ ) {      for ( i = 0, t = 0; i < nd_nvar; i++ ) {
         u = GET_EXP(d,i);          u = GET_EXP(d,i);
Line 551  int ndl_lex_compare(UINT *d1,UINT *d2)
Line 578  int ndl_lex_compare(UINT *d1,UINT *d2)
   
 int ndl_block_compare(UINT *d1,UINT *d2)  int ndl_block_compare(UINT *d1,UINT *d2)
 {  {
     int i,l,j,ord_o,ord_l;      int i,l,j,ord_o;
     struct order_pair *op;      struct order_pair *op;
     UINT t1,t2,m;      UINT t1,t2,m;
     UINT *mask;      UINT *mask;
Line 720  int ndl_module_glex_compare(UINT *d1,UINT *d2)
Line 747  int ndl_module_glex_compare(UINT *d1,UINT *d2)
       break;        break;
   
     case 1:      case 1:
       if ( nd_pot_nelim && MPOS(d1)>=nd_pot_nelim+1 && MPOS(d2) >= nd_pot_nelim+1 ) {        if ( nd_pot_nelim && MPOS(d1)>=(UINT)(nd_pot_nelim+1) && MPOS(d2) >= (UINT)(nd_pot_nelim+1) ) {
          if ( TD(d1) > TD(d2) ) return 1;           if ( TD(d1) > TD(d2) ) return 1;
          else if ( TD(d1) < TD(d2) ) return -1;           else if ( TD(d1) < TD(d2) ) return -1;
          if ( (c = ndl_lex_compare(d1,d2)) != 0 ) return c;           if ( (c = ndl_lex_compare(d1,d2)) != 0 ) return c;
Line 744  int ndl_module_glex_compare(UINT *d1,UINT *d2)
Line 771  int ndl_module_glex_compare(UINT *d1,UINT *d2)
   
     default:      default:
       error("ndl_module_glex_compare : invalid module_ordtype");        error("ndl_module_glex_compare : invalid module_ordtype");
         return 0;
   }    }
 }  }
   
Line 754  int ndl_module_compare(UINT *d1,UINT *d2)
Line 782  int ndl_module_compare(UINT *d1,UINT *d2)
   
   switch ( nd_module_ordtype ) {    switch ( nd_module_ordtype ) {
     case 0:      case 0:
       if ( (c = ndl_lex_compare(d1,d2)) != 0 ) return c;        if ( (c = (*ndl_base_compare_function)(d1,d2)) != 0 ) return c;
       else if ( MPOS(d1) > MPOS(d2) ) return -1;        else if ( MPOS(d1) > MPOS(d2) ) return -1;
       else if ( MPOS(d1) < MPOS(d2) ) return 1;        else if ( MPOS(d1) < MPOS(d2) ) return 1;
       else return 0;        else return 0;
Line 776  int ndl_module_compare(UINT *d1,UINT *d2)
Line 804  int ndl_module_compare(UINT *d1,UINT *d2)
   
     default:      default:
       error("ndl_module_compare : invalid module_ordtype");        error("ndl_module_compare : invalid module_ordtype");
         return 0;
   }    }
 }  }
   
 extern DMMstack dmm_stack;  extern DMMstack dmm_stack;
 void _addtodl(int n,DL d1,DL d2);  void _addtodl(int n,DL d1,DL d2);
   void _adddl(int n,DL d1,DL d2,DL d3);
 int _eqdl(int n,DL d1,DL d2);  int _eqdl(int n,DL d1,DL d2);
   
 int ndl_module_schreyer_compare(UINT *m1,UINT *m2)  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
 {  {
   int pos1,pos2,t,j;    int pos1,pos2,t,j,retpot;
   DMM *in;    DMM *in;
   DMMstack s;    DMMstack s;
   static DL d1=0;    static DL d1=0;
   static DL d2=0;    static DL d2=0;
   static int dlen=0;    static int dlen=0;
     extern int ReversePOT;
   
     if ( ReversePOT ) retpot = -1;
     else retpot = 1;
   pos1 = MPOS(m1); pos2 = MPOS(m2);    pos1 = MPOS(m1); pos2 = MPOS(m2);
   if ( pos1 == pos2 ) return (*ndl_base_compare_function)(m1,m2);    if ( pos1 == pos2 ) return (*ndl_base_compare_function)(m1,m2);
   if ( nd_nvar > dlen ) {    if ( nd_nvar > dlen ) {
Line 808  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
Line 841  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
     _addtodl(nd_nvar,in[pos1]->dl,d1);      _addtodl(nd_nvar,in[pos1]->dl,d1);
     _addtodl(nd_nvar,in[pos2]->dl,d2);      _addtodl(nd_nvar,in[pos2]->dl,d2);
     if ( in[pos1]->pos == in[pos2]->pos && _eqdl(nd_nvar,d1,d2)) {      if ( in[pos1]->pos == in[pos2]->pos && _eqdl(nd_nvar,d1,d2)) {
       if ( pos1 < pos2 ) return 1;        if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return 0;        else return 0;
     }      }
     pos1 = in[pos1]->pos;      pos1 = in[pos1]->pos;
Line 817  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
Line 850  int ndl_module_schreyer_compare(UINT *m1,UINT *m2)
     if ( pos1 == pos2 ) return (*dl_base_compare_function)(nd_nvar,d1,d2);      if ( pos1 == pos2 ) return (*dl_base_compare_function)(nd_nvar,d1,d2);
   }    }
   // comparison by the bottom order    // comparison by the bottom order
 LAST:  
   switch ( nd_base_ordtype ) {    switch ( nd_base_ordtype ) {
     case 0:      case 0:
       t = (*dl_base_compare_function)(nd_nvar,d1,d2);        t = (*dl_base_compare_function)(nd_nvar,d1,d2);
       if ( t ) return t;        if ( t ) return t;
       else if ( pos1 < pos2 ) return 1;        else if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return 0;        else return 0;
       break;        break;
     case 1:      case 1:
       if ( pos1 < pos2 ) return 1;        if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return (*dl_base_compare_function)(nd_nvar,d1,d2);        else return (*dl_base_compare_function)(nd_nvar,d1,d2);
       break;        break;
     case 2:      case 2:
       if ( d1->td > d2->td  ) return 1;        if ( d1->td > d2->td  ) return 1;
       else if ( d1->td < d2->td ) return -1;        else if ( d1->td < d2->td ) return -1;
       else if ( pos1 < pos2 ) return 1;        else if ( pos1 < pos2 ) return retpot;
       else if ( pos1 > pos2 ) return -1;        else if ( pos1 > pos2 ) return -retpot;
       else return (*dl_base_compare_function)(nd_nvar,d1,d2);        else return (*dl_base_compare_function)(nd_nvar,d1,d2);
       break;        break;
     default:      default:
       error("ndl_schreyer_compare : invalid base ordtype");        error("ndl_schreyer_compare : invalid base ordtype");
         return 0;
   }    }
 }  }
   
Line 960  INLINE void ndl_sub(UINT *d1,UINT *d2,UINT *d)
Line 993  INLINE void ndl_sub(UINT *d1,UINT *d2,UINT *d)
   
 int ndl_disjoint(UINT *d1,UINT *d2)  int ndl_disjoint(UINT *d1,UINT *d2)
 {  {
     UINT t1,t2,u,u1,u2;      UINT t1,t2,u1,u2;
     int i,j;      int i,j;
   
     if ( nd_module && (MPOS(d1) == MPOS(d2)) ) return 0;      if ( nd_module && (MPOS(d1) == MPOS(d2)) ) return 0;
Line 1197  INLINE int ndl_find_reducer(UINT *dg)
Line 1230  INLINE int ndl_find_reducer(UINT *dg)
     return -1;      return -1;
 }  }
   
   INLINE int ndl_find_reducer_nonsig(UINT *dg)
   {
       RHist r;
       int i;
   
       for ( i = 0; i < nd_psn; i++ ) {
         r = nd_psh[i];
         if ( ndl_reducible(dg,DL(r)) ) return i;
       }
       return -1;
   }
   
   // ret=0,...,nd_psn-1 => reducer found
   // ret=nd_psn => reducer not found
   // ret=-1 => singular top reducible
   
   int comp_sig(SIG s1,SIG s2);
   void _ndltodl(UINT *ndl,DL dl);
   
   void print_sig(SIG s)
   {
     int i;
   
     fprintf(asir_out,"<<");
     for ( i = 0; i < nd_nvar; i++ ) {
       fprintf(asir_out,"%d",s->dl->d[i]);
       if ( i != nd_nvar-1 ) fprintf(asir_out,",");
     }
     fprintf(asir_out,">>*e%d",s->pos);
   }
   
   void print_siglist(NODE l)
   {
     for ( ; l; l = NEXT(l) )
      print_sig((SIG)l->body);
   }
   
   
   // assuming increasing order wrt signature
   
   INLINE int ndl_find_reducer_s(UINT *dg,SIG sig)
   {
     RHist r;
     int i,singular,ret,d,k;
     static int wpd,nvar;
     static SIG quo;
     static UINT *tmp;
   
     if ( !quo || nvar != nd_nvar ) NEWSIG(quo);
     if ( wpd != nd_wpd ) {
       wpd = nd_wpd;
       tmp = (UINT *)MALLOC(wpd*sizeof(UINT));
     }
     d = ndl_hash_value(dg);
   #if 1
     for ( r = nd_red[d], k = 0; r; r = NEXT(r), k++ ) {
       if ( ndl_equal(dg,DL(r)) ) {
         return r->index;
       }
     }
   #endif
     singular = 0;
     for ( i = 0; i < nd_psn; i++ ) {
       r = nd_psh[i];
       if ( ndl_reducible(dg,DL(r)) ) {
         ndl_sub(dg,DL(r),tmp);
         _ndltodl(tmp,DL(quo));
         _addtodl(nd_nvar,DL(nd_psh[i]->sig),DL(quo));
         quo->pos = nd_psh[i]->sig->pos;
         _adddl(nd_nvar,DL(quo),nd_sba_hm[quo->pos],DL2(quo));
         ret = comp_sig(sig,quo);
   //      if ( ret >= 0 ) { singular = 0; break; }
         if ( ret > 0 ) { singular = 0; break; }
         if ( ret == 0 ) { /* fprintf(asir_out,"s"); fflush(asir_out); */ singular = 1; }
       }
     }
     if ( singular ) return -1;
     else if ( i < nd_psn )
       nd_append_red(dg,i);
     return i;
   }
   
 ND nd_merge(ND p1,ND p2)  ND nd_merge(ND p1,ND p2)
 {  {
     int n,c;      int c;
     int t,can,td1,td2;      int can;
     ND r;  
     NM m1,m2,mr0,mr,s;      NM m1,m2,mr0,mr,s;
   
     if ( !p1 ) return p2;      if ( !p1 ) return p2;
     else if ( !p2 ) return p1;      else if ( !p2 ) return p1;
     else {      else {
         can = 0;          can = 0;
         for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {          for (  m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {
             c = DL_COMPARE(DL(m1),DL(m2));              c = DL_COMPARE(DL(m1),DL(m2));
             switch ( c ) {              switch ( c ) {
                 case 0:                  case 0:
Line 1241  ND nd_merge(ND p1,ND p2)
Line 1355  ND nd_merge(ND p1,ND p2)
   
 ND nd_add(int mod,ND p1,ND p2)  ND nd_add(int mod,ND p1,ND p2)
 {  {
     int n,c;      int c;
     int t,can,td1,td2;      int t,can;
     ND r;  
     NM m1,m2,mr0,mr,s;      NM m1,m2,mr0,mr,s;
   
     Nnd_add++;      Nnd_add++;
Line 1254  ND nd_add(int mod,ND p1,ND p2)
Line 1367  ND nd_add(int mod,ND p1,ND p2)
     else if ( !mod ) return nd_add_q(p1,p2);      else if ( !mod ) return nd_add_q(p1,p2);
     else {      else {
         can = 0;          can = 0;
         for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {          for ( m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {
             c = DL_COMPARE(DL(m1),DL(m2));              c = DL_COMPARE(DL(m1),DL(m2));
             switch ( c ) {              switch ( c ) {
                 case 0:                  case 0:
Line 1315  int nm_destructive_add_q(NM *m1,NM *m2,NM *mr0,NM *mr)
Line 1428  int nm_destructive_add_q(NM *m1,NM *m2,NM *mr0,NM *mr)
   
 ND nd_add_q(ND p1,ND p2)  ND nd_add_q(ND p1,ND p2)
 {  {
     int n,c,can;      int c,can;
     ND r;  
     NM m1,m2,mr0,mr,s;      NM m1,m2,mr0,mr,s;
     P t;      P t;
   
Line 1324  ND nd_add_q(ND p1,ND p2)
Line 1436  ND nd_add_q(ND p1,ND p2)
     else if ( !p2 ) return p1;      else if ( !p2 ) return p1;
     else {      else {
         can = 0;          can = 0;
         for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {          for ( m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {
             c = DL_COMPARE(DL(m1),DL(m2));              c = DL_COMPARE(DL(m1),DL(m2));
             switch ( c ) {              switch ( c ) {
                 case 0:                  case 0:
Line 1366  ND nd_add_q(ND p1,ND p2)
Line 1478  ND nd_add_q(ND p1,ND p2)
   
 ND nd_add_sf(ND p1,ND p2)  ND nd_add_sf(ND p1,ND p2)
 {  {
     int n,c,can;      int c,can;
     ND r;  
     NM m1,m2,mr0,mr,s;      NM m1,m2,mr0,mr,s;
     int t;      int t;
   
Line 1375  ND nd_add_sf(ND p1,ND p2)
Line 1486  ND nd_add_sf(ND p1,ND p2)
     else if ( !p2 ) return p1;      else if ( !p2 ) return p1;
     else {      else {
         can = 0;          can = 0;
         for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {          for ( m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {
             c = DL_COMPARE(DL(m1),DL(m2));              c = DL_COMPARE(DL(m1),DL(m2));
             switch ( c ) {              switch ( c ) {
                 case 0:                  case 0:
Line 1457  ND nd_reduce2(int mod,ND d,ND g,NDV p,NM mul,NDC dn,Ob
Line 1568  ND nd_reduce2(int mod,ND d,ND g,NDV p,NM mul,NDC dn,Ob
 /* ret=1 : success, ret=0 : overflow */  /* ret=1 : success, ret=0 : overflow */
 int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)  int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
 {  {
     NM m,mrd,tail;      NM m,tail;
     NM mul;      NM mul;
     int n,sugar,psugar,sugar0,stat,index;      int n,sugar,index;
     int c,c1,c2,dummy;  
     RHist h;      RHist h;
     NDV p,red;      NDV p;
     Q cg,cred,gcd,tq,qq;  
     Z iq;      Z iq;
     DP dmul;      DP dmul;
     NODE node;      NODE node;
     LIST hist;      LIST hist;
     double hmag;      double hmag;
     P tp,tp1;      Obj div;
     Obj tr,tr1,div;  
     union oNDC hg;      union oNDC hg;
     P cont;      P cont;
   
Line 1480  int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
Line 1588  int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
     }      }
     if ( !mod ) hmag = ((double)p_mag(HCP(g)))*nd_scale;      if ( !mod ) hmag = ((double)p_mag(HCP(g)))*nd_scale;
   
     sugar0 = sugar = SG(g);      sugar = SG(g);
     n = NV(g);      n = NV(g);
     mul = (NM)MALLOC(sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT));      mul = (NM)MALLOC(sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT));
     if ( d )      if ( d )
Line 1503  int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
Line 1611  int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
                 nmtodp(mod,mul,&dmul);                  nmtodp(mod,mul,&dmul);
                 node = mknode(4,div,iq,dmul,ONE);                  node = mknode(4,div,iq,dmul,ONE);
             }              }
             sugar = MAX(sugar,SG(p)+TD(DL(mul)));              sugar = MAX(sugar,(int)(SG(p)+TD(DL(mul))));
             if ( !mod && g && !nd_vc && ((double)(p_mag(HCP(g))) > hmag) ) {              if ( !mod && g && !nd_vc && ((double)(p_mag(HCP(g))) > hmag) ) {
                 hg = HCU(g);                  hg = HCU(g);
                 nd_removecont2(d,g);                  nd_removecont2(d,g);
Line 1539  int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
Line 1647  int nd_nf(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
     return 1;      return 1;
 }  }
   
   // ret=1 => success
   // ret=0 => overflow
   // ret=-1 => singular top reducible
   
   int nd_nf_s(int mod,ND d,ND g,NDV *ps,int full,ND *rp)
   {
       NM m,tail;
       NM mul;
       int n,sugar,index;
       RHist h;
       NDV p;
       Z iq;
       DP dmul;
       NODE node;
       LIST hist;
       double hmag;
       Obj div;
   //    union oNDC hg;
       SIG sig;
   
       if ( !g ) {
           *rp = d;
           return 1;
       }
       if ( !mod ) hmag = ((double)p_mag(HCP(g)))*nd_scale;
   
       sugar = SG(g);
       n = NV(g);
       mul = (NM)MALLOC(sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT));
       if ( d )
           for ( tail = BDY(d); NEXT(tail); tail = NEXT(tail) );
       sig = g->sig;
       for ( ; g; ) {
           index = ndl_find_reducer_s(HDL(g),sig);
           if ( index >= 0 && index < nd_psn ) {
               // reducer found
               h = nd_psh[index];
               ndl_sub(HDL(g),DL(h),DL(mul));
               if ( ndl_check_bound2(index,DL(mul)) ) {
                   nd_free(g); nd_free(d);
                   return 0;
               }
               p = ps[index];
               /* d+g -> div*(d+g)+mul*p */
               g = nd_reduce2(mod,d,g,p,mul,0,&div);
               sugar = MAX(sugar,(int)(SG(p)+TD(DL(mul))));
               if ( !mod && g && ((double)(p_mag(HCP(g))) > hmag) ) {
   //                hg = HCU(g);
                   nd_removecont2(d,g);
                   hmag = ((double)p_mag(HCP(g)))*nd_scale;
               }
               if ( nd_gentrace ) {
                   /* Trace=[div,index,mul,ONE] */
                   STOZ(index,iq);
                   nmtodp(mod,mul,&dmul);
                   node = mknode(4,div,iq,dmul,ONE);
               }
               MKLIST(hist,node);
               MKNODE(node,hist,nd_tracelist); nd_tracelist = node;
           } else if ( index == -1 ) {
             // singular top reducible
             return -1;
           } else if ( !full ) {
               *rp = g;
               g->sig = sig;
               return 1;
           } else {
               m = BDY(g);
               if ( NEXT(m) ) {
                   BDY(g) = NEXT(m); NEXT(m) = 0; LEN(g)--;
               } else {
                   FREEND(g); g = 0;
               }
               if ( d ) {
                   NEXT(tail)=m; tail=m; LEN(d)++;
               } else {
                   MKND(n,m,1,d); tail = BDY(d);
               }
           }
       }
       if ( d ) {
         SG(d) = sugar;
         d->sig = sig;
       }
       *rp = d;
       return 1;
   }
   
 int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *rp)  int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *rp)
 {  {
     int hindex,index;      int hindex,index;
     NDV p;      NDV p;
     ND u,d,red;      ND d,red;
     NODE l;      NM mul,m,tail;
     NM mul,m,mrd,tail;      int sugar,psugar,n;
     int sugar,psugar,n,h_reducible;  
     PGeoBucket bucket;      PGeoBucket bucket;
     int c,c1,c2;      int c,c1,c2;
     Z cg,cred,gcd,zzz;      Z cg,cred,gcd;
     RHist h;      RHist h;
     double hmag,gmag;      double hmag,gmag;
     int count = 0;      int count = 0;
Line 1651  int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *rp
Line 1846  int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *rp
     }      }
 }  }
   
   int nd_nf_pbucket_s(int mod,ND g,NDV *ps,int full,ND *rp)
   {
     int hindex,index;
     NDV p;
     ND d,red;
     NM mul,m,tail;
     int sugar,psugar,n;
     PGeoBucket bucket;
     int c,c1,c2;
     Z cg,cred,gcd;
     RHist h;
     double hmag,gmag;
     int count = 0;
     int hcount = 0;
     SIG sig;
   
     if ( !g ) {
       *rp = 0;
       return 1;
     }
     sugar = SG(g);
     n = NV(g);
     if ( !mod ) hmag = ((double)p_mag((P)HCZ(g)))*nd_scale;
     bucket = create_pbucket();
     add_pbucket(mod,bucket,g);
     d = 0;
     mul = (NM)MALLOC(sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT));
     sig = g->sig;
     while ( 1 ) {
       if ( mod > 0 || mod == -1 )
         hindex = head_pbucket(mod,bucket);
       else if ( mod == -2 )
         hindex = head_pbucket_lf(bucket);
       else
         hindex = head_pbucket_q(bucket);
       if ( hindex < 0 ) {
         if ( DP_Print > 3 ) printf("(%d %d)",count,hcount);
         if ( d ) {
           SG(d) = sugar;
           d->sig = sig;
         }
         *rp = d;
         return 1;
       }
       g = bucket->body[hindex];
       index = ndl_find_reducer_s(HDL(g),sig);
       if ( index >= 0 && index < nd_psn ) {
         count++;
         if ( !d ) hcount++;
         h = nd_psh[index];
         ndl_sub(HDL(g),DL(h),DL(mul));
         if ( ndl_check_bound2(index,DL(mul)) ) {
           nd_free(d);
           free_pbucket(bucket);
           *rp = 0;
           return 0;
         }
         p = ps[index];
         if ( mod == -1 )
           CM(mul) = _mulsf(_invsf(HCM(p)),_chsgnsf(HCM(g)));
         else if ( mod == -2 ) {
           Z inv,t;
           divlf(ONE,HCZ(p),&inv);
           chsgnlf(HCZ(g),&t);
           mullf(inv,t,&CZ(mul));
         } else if ( mod ) {
           c1 = invm(HCM(p),mod); c2 = mod-HCM(g);
           DMAR(c1,c2,0,mod,c); CM(mul) = c;
         } else {
           igcd_cofactor(HCZ(g),HCZ(p),&gcd,&cg,&cred);
           chsgnz(cg,&CZ(mul));
           nd_mul_c_q(d,(P)cred);
           mulq_pbucket(bucket,cred);
           g = bucket->body[hindex];
           gmag = (double)p_mag((P)HCZ(g));
         }
         red = ndv_mul_nm(mod,mul,p);
         bucket->body[hindex] = nd_remove_head(g);
         red = nd_remove_head(red);
         add_pbucket(mod,bucket,red);
         psugar = SG(p)+TD(DL(mul));
         sugar = MAX(sugar,psugar);
         if ( !mod && hmag && (gmag > hmag) ) {
            g = normalize_pbucket(mod,bucket);
            if ( !g ) {
              if ( d ) {
                SG(d) = sugar;
                d->sig = sig;
              }
              *rp = d;
              return 1;
            }
            nd_removecont2(d,g);
            hmag = ((double)p_mag((P)HCZ(g)))*nd_scale;
            add_pbucket(mod,bucket,g);
         }
       } else if ( index == -1 ) {
         // singular top reducible
         return -1;
       } else if ( !full ) {
         g = normalize_pbucket(mod,bucket);
         if ( g ) {
           SG(g) = sugar;
           g->sig = sig;
         }
         *rp = g;
         return 1;
       } else {
         m = BDY(g);
         if ( NEXT(m) ) {
           BDY(g) = NEXT(m); NEXT(m) = 0; LEN(g)--;
         } else {
           FREEND(g); g = 0;
         }
         bucket->body[hindex] = g;
         NEXT(m) = 0;
         if ( d ) {
           NEXT(tail)=m; tail=m; LEN(d)++;
         } else {
           MKND(n,m,1,d); tail = BDY(d);
         }
       }
     }
   }
   
 /* input : list of NDV, cand : list of NDV */  /* input : list of NDV, cand : list of NDV */
   
 int ndv_check_membership(int m,NODE input,int obpe,int oadv,EPOS oepos,NODE cand)  int ndv_check_membership(int m,NODE input,int obpe,int oadv,EPOS oepos,NODE cand)
 {  {
     int n,i,stat;      int i,stat;
     ND nf,d;      ND nf,d;
     NDV r;      NDV r;
     NODE t,s;      NODE t,s;
     union oNDC dn;  
     Z q;      Z q;
     LIST list;      LIST list;
   
     ndv_setup(m,0,cand,nd_gentrace?1:0,1);      ndv_setup(m,0,cand,nd_gentrace?1:0,1,0);
     n = length(cand);  //    n = length(cand);
   
   if ( nd_gentrace ) { nd_alltracelist = 0; nd_tracelist = 0; }    if ( nd_gentrace ) { nd_alltracelist = 0; nd_tracelist = 0; }
     /* membercheck : list is a subset of Id(cand) ? */      /* membercheck : list is a subset of Id(cand) ? */
Line 1750  void free_pbucket(PGeoBucket b) {
Line 2069  void free_pbucket(PGeoBucket b) {
     GCFREE(b);      GCFREE(b);
 }  }
   
   #if 0
 void add_pbucket_symbolic(PGeoBucket g,ND d)  void add_pbucket_symbolic(PGeoBucket g,ND d)
 {  {
     int l,i,k,m;      int l,i,k,m;
Line 1767  void add_pbucket_symbolic(PGeoBucket g,ND d)
Line 2087  void add_pbucket_symbolic(PGeoBucket g,ND d)
     g->body[k] = d;      g->body[k] = d;
     g->m = MAX(g->m,k);      g->m = MAX(g->m,k);
 }  }
   #else
   void add_pbucket_symbolic(PGeoBucket g,ND d)
   {
     int l,k,m,m0;
   
     if ( !d )
       return;
     m0 = g->m;
     while ( 1 ) {
       l = LEN(d);
       for ( k = 0, m = 1; l > m; k++, m <<= 1 );
       /* 2^(k-1) < l <= 2^k (=m) */
       if ( g->body[k] == 0 ) {
         g->body[k] = d;
         m0 = MAX(k,m0);
         break;
       } else {
         d = nd_merge(g->body[k],d);
         g->body[k] = 0;
       }
     }
     g->m = m0;
   }
   #endif
   
   #if 0
 void add_pbucket(int mod,PGeoBucket g,ND d)  void add_pbucket(int mod,PGeoBucket g,ND d)
 {  {
     int l,i,k,m;      int l,i,k,m;
Line 1785  void add_pbucket(int mod,PGeoBucket g,ND d)
Line 2130  void add_pbucket(int mod,PGeoBucket g,ND d)
     g->body[k] = d;      g->body[k] = d;
     g->m = MAX(g->m,k);      g->m = MAX(g->m,k);
 }  }
   #else
   void add_pbucket(int mod,PGeoBucket g,ND d)
   {
     int l,k,m,m0;
   
     m0 = g->m;
     while ( d != 0 ) {
       l = LEN(d);
       for ( k = 0, m = 1; l > m; k++, m <<= 1 );
       /* 2^(k-1) < l <= 2^k (=m) */
       if ( g->body[k] == 0 ) {
         g->body[k] = d;
         m0 = MAX(k,m0);
         break;
       } else {
         d = nd_add(mod,g->body[k],d);
         g->body[k] = 0;
       }
     }
     g->m = m0;
   }
   #endif
   
 void mulq_pbucket(PGeoBucket g,Z c)  void mulq_pbucket(PGeoBucket g,Z c)
 {  {
     int k;      int k;
Line 1828  NM remove_head_pbucket_symbolic(PGeoBucket g)
Line 2195  NM remove_head_pbucket_symbolic(PGeoBucket g)
   
 int head_pbucket(int mod,PGeoBucket g)  int head_pbucket(int mod,PGeoBucket g)
 {  {
     int j,i,c,k,nv,sum;      int j,i,c,k,sum;
     UINT *di,*dj;      UINT *dj;
     ND gi,gj;      ND gi,gj;
   
     k = g->m;      k = g->m;
Line 1874  int head_pbucket(int mod,PGeoBucket g)
Line 2241  int head_pbucket(int mod,PGeoBucket g)
   
 int head_pbucket_q(PGeoBucket g)  int head_pbucket_q(PGeoBucket g)
 {  {
     int j,i,c,k,nv;      int j,i,c,k;
     Z sum,t;      Z sum,t;
     ND gi,gj;      ND gi,gj;
   
Line 1888  int head_pbucket_q(PGeoBucket g)
Line 2255  int head_pbucket_q(PGeoBucket g)
                 gj = g->body[j];                  gj = g->body[j];
                 sum = HCZ(gj);                  sum = HCZ(gj);
             } else {              } else {
                 nv = NV(gi);  //               nv = NV(gi);
                 c = DL_COMPARE(HDL(gi),HDL(gj));                  c = DL_COMPARE(HDL(gi),HDL(gj));
                 if ( c > 0 ) {                  if ( c > 0 ) {
                     if ( sum ) HCZ(gj) = sum;                      if ( sum ) HCZ(gj) = sum;
Line 1914  int head_pbucket_q(PGeoBucket g)
Line 2281  int head_pbucket_q(PGeoBucket g)
   
 int head_pbucket_lf(PGeoBucket g)  int head_pbucket_lf(PGeoBucket g)
 {  {
     int j,i,c,k,nv;      int j,i,c,k;
     Z sum,t;      Z sum,t;
     ND gi,gj;      ND gi,gj;
   
Line 1928  int head_pbucket_lf(PGeoBucket g)
Line 2295  int head_pbucket_lf(PGeoBucket g)
                 gj = g->body[j];                  gj = g->body[j];
                 sum = HCZ(gj);                  sum = HCZ(gj);
             } else {              } else {
                 nv = NV(gi);  //              nv = NV(gi);
                 c = DL_COMPARE(HDL(gi),HDL(gj));                  c = DL_COMPARE(HDL(gi),HDL(gj));
                 if ( c > 0 ) {                  if ( c > 0 ) {
                     if ( sum ) HCZ(gj) = sum;                      if ( sum ) HCZ(gj) = sum;
Line 1955  int head_pbucket_lf(PGeoBucket g)
Line 2322  int head_pbucket_lf(PGeoBucket g)
 ND normalize_pbucket(int mod,PGeoBucket g)  ND normalize_pbucket(int mod,PGeoBucket g)
 {  {
     int i;      int i;
     ND r,t;      ND r;
   
     r = 0;      r = 0;
     for ( i = 0; i <= g->m; i++ ) {      for ( i = 0; i <= g->m; i++ ) {
Line 1996  void register_hcf(NDV p)
Line 2363  void register_hcf(NDV p)
 #else  #else
 void register_hcf(NDV p)  void register_hcf(NDV p)
 {  {
     DCP dc,t;  
     P hc,h,q;      P hc,h,q;
     Q dmy;      Q dmy;
     int c;      int c;
Line 2030  void register_hcf(NDV p)
Line 2396  void register_hcf(NDV p)
   
 int do_diagonalize(int sugar,int m)  int do_diagonalize(int sugar,int m)
 {  {
   int i,nh,stat;    int i,stat;
   NODE r,g,t;    ND nf,s,head;
   ND h,nf,s,head;  
   NDV nfv;    NDV nfv;
   Q q;    P cont;
   P nm,nmp,dn,mnp,dnp,cont,cont1;  
   union oNDC hc;    union oNDC hc;
   NODE node;    NODE node;
   LIST l;    LIST l;
   Z iq;    Z iq;
   
     if ( diag_period == 0 ) return 1;
   for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) {    for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) {
     if ( nd_gentrace ) {      if ( nd_gentrace ) {
       /* Trace = [1,index,1,1] */        /* Trace = [1,index,1,1] */
Line 2096  LIST compute_splist()
Line 2461  LIST compute_splist()
   return l0;    return l0;
 }  }
   
   typedef struct oHPDATA {
     int n;
     P hn; // HP(t)=hn(t)/((1-t^w0)*...*(1-t^w(n-1)))
     VECT x; // BDY(x)[i] = <<0,...,1,...,0>>
     int *w;
   } *HPDATA;
   
   void make_reduced(VECT b,int nv);
   void mhp_rec(VECT b,VECT x,P t,P *r);
   P mhp_ctop(P *r,P *plist,int n);
   DL monomial_colon(DL a,DL b,int n);
   LIST dp_monomial_hilbert_poincare(VECT b,VECT x);
   LIST dp_monomial_hilbert_poincare_weight(VECT b,VECT x,int *w);
   
   void setup_hpdata(HPDATA final,HPDATA current)
   {
     int n,i;
     DL *p;
     VECT b,x;
     DL dl;
     LIST ret;
     int *w;
   
     final->n = n = nd_nvar;
     final->hn = (P)BDY(nd_hpdata);
   #if 0
     if ( NEXT(nd_hpdata) != 0 && (weight=(LIST)BDY(NEXT(nd_hpdata))) != 0 ) {
       wlen = length(BDY(weight));
       if ( n != wlen )
         error("setup_hpdata : inconsistent weight length");
       w = (int *)MALLOC(n*sizeof(int));
       for ( i = 0, nd = BDY((LIST)weight); i < n; i++, nd = NEXT(nd) )
         w[i] = ZTOS((Z)BDY(nd));
     } else
       w = 0;
   #else
     w = current_dl_weight_vector;
   #endif
     MKVECT(x,n);
     for ( i = 0; i < n; i++ ) {
       NEWDL(dl,n); dl->d[i] = 1; dl->td = 1; BDY(x)[i] = dl;
     }
     final->x = x;
     final->w = w;
   
     MKVECT(b,nd_psn); p = (DL *)BDY(b);
     for ( i = 0; i < nd_psn; i++ ) {
       p[i] = ndltodl(n,nd_psh[i]->dl);
     }
     if ( w ) {
       ret = dp_monomial_hilbert_poincare_weight(b,x,w);
     } else
       ret = dp_monomial_hilbert_poincare(b,x);
     current->n = n;
     current->hn = (P)BDY(BDY(ret));
     current->x = x;
     current->w = w;
   }
   
   int comp_hn(P a, P b)
   {
     P s;
     DCP dc;
   
     subp(CO,a,b,&s);
     if ( !s ) return -1;
     else if ( OID(s) == 1 ) return 0;
     else {
       for ( dc = DC(s); NEXT(dc); dc = NEXT(dc) );
       return (int)ZTOS((Z)dc->d);
     }
   }
   
   void update_hpdata(HPDATA current,int nh)
   {
     NODE data1;
     DL new;
     int i,n;
     Z dz;
     DL *p;
     VECT b;
     P tv,td,s,hn;
     LIST list1;
   
     n = nd_nvar;
     new = ndltodl(n,nd_psh[nh]->dl);
     MKVECT(b,nh); p = (DL *)BDY(b);
     for ( i = 0; i < nh; i++ ) {
       p[i] = monomial_colon(ndltodl(n,nd_psh[i]->dl),new,n);
     }
     // compute HP(I:new)
     if ( current->w )
       list1 = dp_monomial_hilbert_poincare_weight(b,current->x,current->w);
     else
       list1 = dp_monomial_hilbert_poincare(b,current->x);
     data1 = BDY((LIST)list1);
     // HP(I+<new>) = H(I)-t^d*H(I:new), d=tdeg(new)
     makevar("t",&tv); UTOZ(new->td,dz);
     pwrp(CO,tv,dz,&td);
     mulp(CO,(P)ARG0(data1),td,&s);
     subp(CO,current->hn,s,&hn);
     current->hn = hn;
   }
   
   ND_pairs nd_remove_same_sugar( ND_pairs d, int sugar)
   {
     struct oND_pairs root;
     ND_pairs prev,cur;
     int i;
   
     root.next = d;
     prev = &root; cur = d;
     i = 0;
     while ( cur ) {
       if ( SG(cur) == sugar ) {
         prev->next = cur->next;
         i++;
       } else
         prev = cur;
       cur = cur->next;
     }
     if ( DP_Print && i ) fprintf(asir_out,"[%d]",i);
     return root.next;
   }
   
 /* return value = 0 => input is not a GB */  /* return value = 0 => input is not a GB */
   
 NODE nd_gb(int m,int ishomo,int checkonly,int gensyz,int **indp)  NODE nd_gb(int m,int ishomo,int checkonly,int gensyz,int **indp)
 {  {
   int i,nh,sugar,stat;    int i,nh,sugar,stat;
   NODE r,g,t;    NODE g,t;
   ND_pairs d;    ND_pairs d;
   ND_pairs l;    ND_pairs l;
   ND h,nf,s,head,nf1;    ND h,nf;
   NDV nfv;    NDV nfv;
   Z q;    Z q;
   union oNDC dn,hc;    union oNDC hc;
   int diag_count = 0;    int diag_count = 0;
     int Nnfnz = 0,Nnfz = 0;
   P cont;    P cont;
   LIST list;    LIST list;
     struct oHPDATA current_hpdata,final_hpdata;
   
   struct oEGT eg1,eg2,eg_update;
   
   init_eg(&eg_update);
   Nnd_add = 0;    Nnd_add = 0;
   g = 0; d = 0;    g = 0; d = 0;
   for ( i = 0; i < nd_psn; i++ ) {    for ( i = 0; i < nd_psn; i++ ) {
Line 2119  NODE nd_gb(int m,int ishomo,int checkonly,int gensyz,i
Line 2614  NODE nd_gb(int m,int ishomo,int checkonly,int gensyz,i
     g = update_base(g,i);      g = update_base(g,i);
   }    }
   sugar = 0;    sugar = 0;
     if ( nd_hpdata ) {
       if ( DP_Print ) fprintf(asir_out,"Hilbert driven algorithm.\n");
       setup_hpdata(&final_hpdata,&current_hpdata);
     }
   while ( d ) {    while ( d ) {
 again:  again:
     l = nd_minp(d,&d);      l = nd_minp(d,&d);
Line 2143  again:
Line 2642  again:
       goto again;        goto again;
     }      }
 #if USE_GEOBUCKET  #if USE_GEOBUCKET
     stat = (m&&!nd_gentrace)?nd_nf_pbucket(m,h,nd_ps,!Top,&nf)      stat = (m&&!nd_gentrace)?nd_nf_pbucket(m,h,nd_ps,!nd_top&&!Top,&nf)
       :nd_nf(m,0,h,nd_ps,!Top,&nf);        :nd_nf(m,0,h,nd_ps,!nd_top&&!Top,&nf);
 #else  #else
     stat = nd_nf(m,0,h,nd_ps,!Top,&nf);      stat = nd_nf(m,0,h,nd_ps,!nd_top&&!Top,&nf);
 #endif  #endif
     if ( !stat ) {      if ( !stat ) {
       NEXT(l) = d; d = l;        NEXT(l) = d; d = l;
       d = nd_reconstruct(0,d);        d = nd_reconstruct(0,d);
       goto again;        goto again;
     } else if ( nf ) {      } else if ( nf ) {
         Nnfnz++;
       if ( checkonly || gensyz ) return 0;        if ( checkonly || gensyz ) return 0;
       if ( nd_newelim ) {        if ( nd_newelim ) {
         if ( nd_module ) {          if ( nd_module ) {
Line 2176  again:
Line 2676  again:
         }          }
       }        }
       nfv = ndtondv(m,nf); nd_free(nf);        nfv = ndtondv(m,nf); nd_free(nf);
       nh = ndv_newps(m,nfv,0,0);        nh = ndv_newps(m,nfv,0);
       if ( !m && (ishomo && ++diag_count == diag_period) ) {        if ( !m && (ishomo && ++diag_count == diag_period) ) {
         diag_count = 0;          diag_count = 0;
         stat = do_diagonalize(sugar,m);          stat = do_diagonalize(sugar,m);
Line 2186  again:
Line 2686  again:
           goto again;            goto again;
         }          }
       }        }
   get_eg(&eg1);
       d = update_pairs(d,g,nh,0);        d = update_pairs(d,g,nh,0);
   get_eg(&eg2); add_eg(&eg_update,&eg1,&eg2);
       g = update_base(g,nh);        g = update_base(g,nh);
       FREENDP(l);        FREENDP(l);
    } else {        if ( nd_hpdata ) {
      if ( nd_gentrace && gensyz ) {          int dg,sugar0;
        nd_tracelist = reverse_node(nd_tracelist);  
        MKLIST(list,nd_tracelist);          update_hpdata(&current_hpdata,nh);
        STOZ(-1,q); t = mknode(2,q,list); MKLIST(list,t);          dg = comp_hn(final_hpdata.hn,current_hpdata.hn);
        MKNODE(t,list,nd_alltracelist);          if ( dg < 0 ) {
        nd_alltracelist = t; nd_tracelist = 0;             int d_len;
      }             for ( d_len = 0; d; d = d->next, d_len++);
      if ( DP_Print ) { printf("."); fflush(stdout); }             fprintf(asir_out,"[%d] We found a gb\n",d_len);
        FREENDP(l);          }
    }          sugar0 = sugar;
  }          while ( d && dg > sugar0 ) {
  conv_ilist(nd_demand,0,g,indp);            d = nd_remove_same_sugar(d,sugar0);
     if ( !checkonly && DP_Print ) { printf("nd_gb done. Number of nd_add=%d\n",Nnd_add); fflush(stdout); }            sugar0++;
     return g;          }
         }
       } else {
         Nnfz++;
         if ( nd_gentrace && gensyz ) {
           nd_tracelist = reverse_node(nd_tracelist);
           MKLIST(list,nd_tracelist);
           STOZ(-1,q); t = mknode(2,q,list); MKLIST(list,t);
           MKNODE(t,list,nd_alltracelist);
           nd_alltracelist = t; nd_tracelist = 0;
         }
         if ( DP_Print ) { fprintf(asir_out,"."); fflush(asir_out); }
           FREENDP(l);
       }
     }
     conv_ilist(nd_demand,0,g,indp);
     if ( !checkonly && DP_Print ) {
       fprintf(asir_out,"\nnd_gb done. Nbase=%d,Nnd_add=%d,Npairs=%d, Nnfnz=%d,Nnfz=%d,",nd_psn,Nnd_add,Npairs,Nnfnz,Nnfz);
       fprintf(asir_out,"Nremoved=%d\n",NcriB+NcriMF+Ncri2);
       fflush(asir_out);
     }
     if ( DP_Print ) {
       print_eg("update",&eg_update); fprintf(asir_out,"\n");
     }
     return g;
 }  }
   
   ND_pairs update_pairs_s(ND_pairs d,int t,NODE *syz);
   int update_pairs_array_s(ND_pairs *d,int t,NODE *syz);
   ND_pairs nd_newpairs_s(int t ,NODE *syz);
   ND_pairs *nd_newpairs_array_s(int t ,NODE *syz);
   
   int nd_nf_pbucket_s(int mod,ND g,NDV *ps,int full,ND *nf);
   int nd_nf_s(int mod,ND d,ND g,NDV *ps,int full,ND *nf);
   
   void _copydl(int n,DL d1,DL d2);
   void _subfromdl(int n,DL d1,DL d2);
   extern int (*cmpdl)(int n,DL d1,DL d2);
   int _dl_redble_ext(DL,DL,DL,int);
   
   int primitive_irred(ND p,SIG sig)
   {
     static int wpd=0,dlen=0;
     static DL dquo,squo;
     static UINT *quo;
     int i;
   
     if ( dlen < nd_nvar ) {
       NEWDL(dquo,nd_nvar);
       NEWDL(squo,nd_nvar);
       dlen = nd_nvar;
     }
     if ( wpd != nd_wpd ) {
       wpd = nd_wpd;
       quo = (UINT *)MALLOC(wpd*sizeof(UINT));
     }
     for ( i = 0; i < nd_psn; i++ ) {
       if ( sig->pos == nd_psh[i]->sig->pos &&
         _dl_redble_ext(DL(nd_psh[i]->sig),DL(sig),squo,nd_nvar) )
         if ( ndl_reducible(HDL(p),DL(nd_psh[i])) ) {
           if ( DP_Print ) fprintf(asir_out,"D");
           ndl_sub(HDL(p),DL(nd_psh[i]),quo);
           _ndltodl(quo,dquo);
           if ( _eqdl(nd_nvar,squo,dquo) )
             return 0;
         }
     }
     return 1;
   }
   
   NODE insert_sig(NODE l,SIG s)
   {
     int pos;
     DL sig;
     struct oNODE root;
     NODE p,prev,r;
     SIG t;
   
     pos = s->pos; sig = DL(s);
     root.next = l; prev = &root;
     for ( p = l; p; p = p->next ) {
       t = (SIG)p->body;
       if ( t->pos == pos ) {
         if ( _dl_redble(DL(t),sig,nd_nvar) )
           return root.next;
         else if ( _dl_redble(sig,DL(t),nd_nvar) )
           // remove p
           prev->next = p->next;
         else
           prev = p;
       } else
         prev = p;
     }
     NEWNODE(r); r->body = (pointer)s; r->next = 0;
     for ( p = &root; p->next; p = p->next );
     p->next = r;
   //  r->next = root.next;
   //  return r;
     return root.next;
   }
   
   ND_pairs remove_spair_s(ND_pairs d,SIG sig)
   {
     struct oND_pairs root;
     ND_pairs prev,p;
     SIG spsig;
   
     root.next = d;
     prev = &root; p = d;
     while ( p ) {
       spsig = p->sig;
       if ( sig->pos == spsig->pos && _dl_redble(DL(sig),DL(spsig),nd_nvar) ) {
         // remove p
         prev->next = p->next;
         Nsyz++;
       } else
         prev = p;
       p = p->next;
     }
     return (ND_pairs)root.next;
   }
   
   int small_lcm(ND_pairs l)
   {
     SIG sig;
     int i;
     NODE t;
     static DL lcm,mul,quo;
     static int nvar = 0;
   
     if ( nd_sba_largelcm ) return 0;
     if ( nvar < nd_nvar ) {
       nvar = nd_nvar; NEWDL(lcm,nvar); NEWDL(quo,nvar); NEWDL(mul,nvar);
     }
     sig = l->sig;
     _ndltodl(l->lcm,lcm);
   #if 0
     for ( i = 0; i < nd_psn; i++ ) {
       if ( sig->pos == nd_psh[i]->sig->pos &&
         _dl_redble_ext(DL(nd_psh[i]->sig),DL(sig),quo,nd_nvar) ) {
         _ndltodl(DL(nd_psh[i]),mul);
         _addtodl(nd_nvar,quo,mul);
         if ( (*cmpdl)(nd_nvar,lcm,mul) > 0 )
           break;
       }
     }
     if ( i < nd_psn ) return 1;
     else return 0;
   #else
     for ( t = nd_sba_pos[sig->pos]; t; t = t->next ) {
       i = (long)BDY(t);
       if ( _dl_redble_ext(DL(nd_psh[i]->sig),DL(sig),quo,nd_nvar) ) {
         _ndltodl(DL(nd_psh[i]),mul);
         _addtodl(nd_nvar,quo,mul);
         if ( (*cmpdl)(nd_nvar,lcm,mul) > 0 )
           break;
       }
     }
     if ( t ) return 1;
     else return 0;
   #endif
   }
   
   ND_pairs find_smallest_lcm(ND_pairs l)
   {
     SIG sig;
     int i,minindex;
     NODE t;
     ND_pairs r;
     struct oSIG sig1;
     static DL mul,quo,quo2,minlm;
     static int nvar = 0;
   
     if ( nvar < nd_nvar ) {
       nvar = nd_nvar;
       NEWDL(quo,nvar); NEWDL(quo2,nvar); NEWDL(mul,nvar);
       NEWDL(minlm,nvar);
     }
     sig = l->sig;
     // find mg s.t. m*s(g)=sig and m*lm(g) is minimal
     _ndltodl(l->lcm,minlm); minindex = -1;
     for ( t = nd_sba_pos[sig->pos]; t; t = t->next ) {
       i = (long)BDY(t);
       if ( _dl_redble_ext(DL(nd_psh[i]->sig),DL(sig),quo,nd_nvar) ) {
         _ndltodl(DL(nd_psh[i]),mul);
         _addtodl(nd_nvar,quo,mul);
         if ( (*cmpdl)(nd_nvar,minlm,mul) > 0 ) {
           minindex = i;
           break;
           _copydl(nd_nvar,mul,minlm);
         }
       }
     }
     // l->lcm is minimal; return l itself
     if ( minindex < 0 ) return l;
     else return 0;
     for ( i = 0; i < nd_psn; i++ ) {
       if ( i == minindex ) continue;
       _ndltodl(DL(nd_psh[i]),mul);
       if ( _dl_redble_ext(mul,minlm,quo,nd_nvar) ) {
         _addtodl(nd_nvar,nd_ps[i]->sig->dl,quo);
         sig1.pos = nd_ps[i]->sig->pos;
         sig1.dl = quo;
         sig1.dl2 = quo2;
         _adddl(nd_nvar,sig1.dl,nd_sba_hm[sig1.pos],sig1.dl2);
         if ( comp_sig(sig,&sig1) > 0 ) {
   //        printf("X");
           NEWND_pairs(r);
           r->sig = sig;
           r->i1 = minindex;
           r->i2 = i;
           dltondl(nd_nvar,minlm,r->lcm);
           r->next = 0;
           return r;
         }
       }
     }
     // there is no suitable spair
     return 0;
   }
   
   ND_pairs remove_large_lcm(ND_pairs d)
   {
     struct oND_pairs root;
     ND_pairs prev,p;
   
     root.next = d;
     prev = &root; p = d;
     while ( p ) {
   #if 0
       if ( small_lcm(p) ) {
         // remove p
         prev->next = p->next;
       } else
   #else
       if ( find_smallest_lcm(p) == 0 ) {
         // remove p
         prev->next = p->next;
       } else
   #endif
         prev = p;
       p = p->next;
     }
     return (ND_pairs)root.next;
   }
   
   struct oEGT eg_create,eg_newpairs,eg_merge;
   
   NODE conv_ilist_s(int demand,int trace,int **indp);
   
   // S(fj*ei-fi*ej)
   
   void _subdl(int,DL,DL,DL);
   
   SIG trivial_sig(int i,int j)
   {
     static DL lcm;
     static struct oSIG sigi,sigj;
     static int nvar = 0;
     SIG sig;
   
     if ( nvar != nd_nvar ) {
       nvar = nd_nvar; NEWDL(lcm,nvar); NEWDL(sigi.dl,nvar); NEWDL(sigj.dl,nvar);
       NEWDL(sigi.dl2,nvar); NEWDL(sigj.dl2,nvar);
     }
     if ( nd_sba_inputisgb != 0 ) {
       lcm_of_DL(nd_nvar,nd_sba_hm[i],nd_sba_hm[j],lcm);
       sigi.pos = i; _subdl(nd_nvar,lcm,nd_sba_hm[i],sigi.dl);
       _copydl(nd_nvar,lcm,sigi.dl2);
       sigj.pos = j; _subdl(nd_nvar,lcm,nd_sba_hm[j],sigj.dl);
       _copydl(nd_nvar,lcm,sigj.dl2);
       if ( comp_sig(&sigi,&sigj) > 0 ) sig = dup_sig(&sigi);
       else sig = dup_sig(&sigj);
     } else {
       sigi.pos = i; _copydl(nd_nvar,nd_sba_hm[j],sigi.dl);
       _adddl(nd_nvar,sigi.dl,nd_sba_hm[i],sigi.dl2);
       sigj.pos = j; _copydl(nd_nvar,nd_sba_hm[i],sigj.dl);
       _adddl(nd_nvar,sigj.dl,nd_sba_hm[j],sigj.dl2);
       if ( comp_sig(&sigi,&sigj) > 0 ) sig = dup_sig(&sigi);
       else sig = dup_sig(&sigj);
     }
     return sig;
   }
   
   int nd_minsig(ND_pairs *d)
   {
     int min,i,ret;
   
     min = -1;
     for ( i = 0; i < nd_nbase; i++ ) {
       if ( d[i] != 0 ) {
         if ( min < 0 ) min = i;
         else {
           ret = comp_sig(d[i]->sig,d[min]->sig);
           if ( ret < 0 ) min = i;
         }
       }
     }
     return min;
   }
   
   int dlength(ND_pairs d)
   {
     int i;
     for ( i = 0; d; d = d->next, i++ );
     return i;
   }
   
   NODE nd_sba_buch(int m,int ishomo,int **indp,NODE *syzp)
   {
     int i,j,nh,sugar,stat,pos;
     NODE g;
     ND_pairs *d;
     ND_pairs l,l1;
     ND h,nf;
     NDV nfv;
   //  union oNDC hc;
     SIG sig;
     NODE *syzlist;
     int ngen,ind;
     int Nnominimal,Nredundant;
     DL lcm,quo,mul;
     struct oHPDATA final_hpdata,current_hpdata;
     struct oEGT eg1,eg2,eg3,eg4,eg_remove,eg_large,eg_nf,eg_nfzero;
     struct oEGT eg_minsig,eg_smallest,eg_removecont,eg_hpdata,eg_updatepairs,eg_sbabuch,eg_sp;
     int Nnfs=0,Nnfz=0,Nnfnz=0,dlen,nsyz;
   
   init_eg(&eg_remove);
     syzlist = (NODE *)MALLOC(nd_psn*sizeof(NODE));
     d = (ND_pairs *)MALLOC(nd_psn*sizeof(ND_pairs));
     nd_nbase = nd_psn;
     Nsyz = 0;
     Nnd_add = 0;
     Nnominimal = 0;
     Nredundant = 0;
     ngen = nd_psn;
     if ( !do_weyl || nd_sba_inputisgb ) {
       for ( i = 0; i < nd_psn; i++ )
         for ( j = i+1; j < nd_psn; j++ ) {
           sig = trivial_sig(i,j);
           syzlist[sig->pos] = insert_sig(syzlist[sig->pos],sig);
         }
       }
     dlen = 0;
     for ( i = 0; i < nd_psn; i++ ) {
       dlen += update_pairs_array_s(d,i,syzlist);
     }
     sugar = 0;
     pos = 0;
     if ( nd_hpdata ) {
       setup_hpdata(&final_hpdata,&current_hpdata);
     }
     NEWDL(lcm,nd_nvar); NEWDL(quo,nd_nvar); NEWDL(mul,nd_nvar);
   init_eg(&eg_sp);
   init_eg(&eg_create);
   init_eg(&eg_merge);
   init_eg(&eg_minsig);
   init_eg(&eg_smallest);
   init_eg(&eg_large);
   init_eg(&eg_nf);
   init_eg(&eg_nfzero);
   init_eg(&eg_removecont);
   init_eg(&eg_updatepairs);
   init_eg(&eg_hpdata);
   init_eg(&eg_sbabuch);
   get_eg(&eg3);
     while ( 1 ) {
       if ( DP_Print && !nd_hpdata && dlen%1000 == 0 ) fprintf(asir_out,"(%d)",dlen);
   again :
   get_eg(&eg1);
       ind = nd_minsig(d);
   get_eg(&eg2); add_eg(&eg_minsig,&eg1,&eg2);
       if ( ind < 0 ) break;
       l = d[ind];
   //    printf("(%d,%d)",l->i1,l->i2); print_sig(l->sig); printf("\n");
   get_eg(&eg1);
       l1 = find_smallest_lcm(l);
   get_eg(&eg2); add_eg(&eg_smallest,&eg1,&eg2);
       if ( l1 == 0 ) {
         d[ind] = d[ind]->next; dlen--;
   //      if ( DP_Print && !nd_hpdata ) fprintf(asir_out,"M");
         Nnominimal++;
         continue;
       }
       if ( SG(l1) != sugar ) {
         sugar = SG(l1);
         if ( DP_Print ) fprintf(asir_out,"%d",sugar);
       }
       sig = l1->sig;
       if ( DP_Print && nd_sba_pot ) {
         if ( sig->pos != pos ) {
           fprintf(asir_out,"[%d]",sig->pos);
           pos = sig->pos;
         }
       }
   get_eg(&eg1);
       stat = nd_sp(m,0,l1,&h);
   get_eg(&eg2); add_eg(&eg_sp,&eg1,&eg2);
       if ( !stat ) {
         nd_reconstruct_s(0,d);
         goto again;
       }
   get_eg(&eg1);
   #if USE_GEOBUCKET
       stat = (m&&!nd_gentrace)?nd_nf_pbucket_s(m,h,nd_ps,!nd_top&&!Top,&nf):nd_nf_s(m,0,h,nd_ps,!nd_top&&!Top,&nf);
   #else
       stat = nd_nf_s(m,0,h,nd_ps,!nd_top&&!Top,&nf);
   #endif
   get_eg(&eg2);
       if ( !stat ) {
         nd_reconstruct_s(0,d);
         goto again;
       } else if ( stat == -1 ) {
         d[ind] = d[ind]->next; dlen--;
         Nnfs++;
         if ( DP_Print ) { printf("S"); fflush(stdout); }
       } else if ( nf ) {
         d[ind] = d[ind]->next; dlen--;
         Nnfnz++;
         if ( DP_Print ) {
           if ( nd_sba_redundant_check ) {
             if ( ndl_find_reducer_nonsig(HDL(nf)) >= 0 ) {
               Nredundant++;
               printf("R");
             } else
               printf("+");
           } else
             printf("+");
           fflush(stdout);
         }
         add_eg(&eg_nf,&eg1,&eg2);
   //      hc = HCU(nf);
         get_eg(&eg1);
         nd_removecont(m,nf);
         get_eg(&eg2); add_eg(&eg_removecont,&eg1,&eg2);
         nfv = ndtondv(m,nf); nd_free(nf);
         nh = ndv_newps(m,nfv,0);
   
         get_eg(&eg1);
         dlen += update_pairs_array_s(d,nh,syzlist);
         get_eg(&eg2); add_eg(&eg_updatepairs,&eg1,&eg2);
         nd_sba_pos[sig->pos] = append_one(nd_sba_pos[sig->pos],nh);
         if ( nd_hpdata ) {
           int dg,sugar0;
   
           get_eg(&eg1);
           update_hpdata(&current_hpdata,nh);
           get_eg(&eg2); add_eg(&eg_hpdata,&eg1,&eg2);
           dg = comp_hn(final_hpdata.hn,current_hpdata.hn);
   //        if ( !compp(CO,final_hpdata.hn,current_hpdata.hn) )
           if ( dg < 0 ) {
             if ( DP_Print ) { printf("\nWe found a gb.\n"); }
               break;
           } else if ( nd_sba_heu == 1 ) {
             for ( i = 0; i < ngen; i++ ) {
               sugar0 = sugar;
               while ( d[i] && dg > sugar0 ) {
                 d[i] = nd_remove_same_sugar(d[i],sugar0);
                 sugar0++;
               }
             }
           }
         }
       } else {
         d[ind] = d[ind]->next; dlen--;
         Nnfz++;
         add_eg(&eg_nfzero,&eg1,&eg2);
       // syzygy
   get_eg(&eg1);
         nsyz = Nsyz;
         d[sig->pos] = remove_spair_s(d[sig->pos],sig);
         dlen -= Nsyz-nsyz;
   get_eg(&eg2); add_eg(&eg_remove,&eg1,&eg2);
         syzlist[sig->pos] = insert_sig(syzlist[sig->pos],sig);
         if ( DP_Print ) { printf("."); fflush(stdout); }
       }
     }
     get_eg(&eg4); add_eg(&eg_sbabuch,&eg3,&eg4);
     g = conv_ilist_s(nd_demand,0,indp);
     if ( DP_Print ) {
       printf("\ndlen=%d,nd_sba done. nd_add=%d,Nsyz=%d,Nsamesig=%d,Nnominimal=%d\n",dlen,Nnd_add,Nsyz,Nsamesig,Nnominimal);
       printf("Nbase=%d,Nnfnz=%d,Nnfz=%d,Nnfsingular=%d\n",nd_psn,Nnfnz,Nnfz,Nnfs);
       fflush(stdout);
       if ( nd_sba_redundant_check )
         printf("Nredundant=%d\n",Nredundant);
       fflush(stdout);
        print_eg("sp",&eg_sp);
        print_eg("create",&eg_create);
        print_eg("merge",&eg_merge);
        print_eg("minsig",&eg_minsig);
        print_eg("smallest",&eg_smallest);
        print_eg("remove",&eg_remove);
        printf("\n");
        print_eg("nf",&eg_nf);
        print_eg("nfzero",&eg_nfzero);
        print_eg("removecont",&eg_removecont);
        print_eg("updatepairs",&eg_updatepairs);
        print_eg("hpdata",&eg_hpdata);
        print_eg("total",&eg_sbabuch);
        printf("\n");
     }
     if ( nd_sba_syz ) {
       print_eg("remove",&eg_remove);
       print_eg("nf",&eg_nf);
       print_eg("nfzero",&eg_nfzero);
       printf("\n");
     }
     if ( nd_sba_syz ) {
       NODE hsyz,tsyz,prev;
   
       hsyz = 0;
       for ( i = 0; i < ngen; i++ ) {
         tsyz = syzlist[i];
         for ( prev = 0; tsyz != 0; prev = tsyz, tsyz = NEXT(tsyz))
           BDY(tsyz) = (pointer)sigtodpm((SIG)BDY(tsyz));
         if ( prev != 0 ) {
           prev->next = hsyz; hsyz = syzlist[i];
         }
       }
       *syzp = hsyz;
     } else *syzp = 0;
     return g;
   }
   
 /* splist = [[i1,i2],...] */  /* splist = [[i1,i2],...] */
   
 int check_splist(int m,NODE splist)  int check_splist(int m,NODE splist)
Line 2218  int check_splist(int m,NODE splist)
Line 3241  int check_splist(int m,NODE splist)
   for ( d = 0, t = splist; t; t = NEXT(t) ) {    for ( d = 0, t = splist; t; t = NEXT(t) ) {
     p = BDY((LIST)BDY(t));      p = BDY((LIST)BDY(t));
     NEXTND_pairs(d,r);      NEXTND_pairs(d,r);
     r->i1 = ZTOS((Q)ARG0(p)); r->i2 = ZTOS((Q)ARG1(p));      r->i1 = (int)ZTOS((Q)ARG0(p)); r->i2 = (int)ZTOS((Q)ARG1(p));
     ndl_lcm(DL(nd_psh[r->i1]),DL(nd_psh[r->i2]),r->lcm);      ndl_lcm(DL(nd_psh[r->i1]),DL(nd_psh[r->i2]),r->lcm);
     SG(r) = TD(LCM(r)); /* XXX */      SG(r) = TD(LCM(r)); /* XXX */
   }    }
Line 2233  again:
Line 3256  again:
       d = nd_reconstruct(0,d);        d = nd_reconstruct(0,d);
       goto again;        goto again;
     }      }
     stat = nd_nf(m,0,h,nd_ps,!Top,&nf);      stat = nd_nf(m,0,h,nd_ps,!nd_top&&!Top,&nf);
     if ( !stat ) {      if ( !stat ) {
       NEXT(l) = d; d = l;        NEXT(l) = d; d = l;
       d = nd_reconstruct(0,d);        d = nd_reconstruct(0,d);
Line 2256  int check_splist_f4(int m,NODE splist)
Line 3279  int check_splist_f4(int m,NODE splist)
   for ( d = 0, t = splist; t; t = NEXT(t) ) {    for ( d = 0, t = splist; t; t = NEXT(t) ) {
     p = BDY((LIST)BDY(t));      p = BDY((LIST)BDY(t));
     NEXTND_pairs(d,r);      NEXTND_pairs(d,r);
     r->i1 = ZTOS((Q)ARG0(p)); r->i2 = ZTOS((Q)ARG1(p));      r->i1 = (int)ZTOS((Q)ARG0(p)); r->i2 = (int)ZTOS((Q)ARG1(p));
     ndl_lcm(DL(nd_psh[r->i1]),DL(nd_psh[r->i2]),r->lcm);      ndl_lcm(DL(nd_psh[r->i1]),DL(nd_psh[r->i2]),r->lcm);
     SG(r) = TD(LCM(r)); /* XXX */      SG(r) = TD(LCM(r)); /* XXX */
   }    }
Line 2287  int check_splist_f4(int m,NODE splist)
Line 3310  int check_splist_f4(int m,NODE splist)
   
 int do_diagonalize_trace(int sugar,int m)  int do_diagonalize_trace(int sugar,int m)
 {  {
   int i,nh,stat;    int i,stat;
   NODE r,g,t;    ND nf,s,head;
   ND h,nf,nfq,s,head;    NDV nfv;
   NDV nfv,nfqv;  
   Q q,den,num;  
   union oNDC hc;    union oNDC hc;
   NODE node;    NODE node;
   LIST l;    LIST l;
   Z iq;    Z iq;
   P cont,cont1;    P cont;
   
     if ( diag_period == 0 ) return 1;
   for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) {    for ( i = nd_psn-1; i >= 0 && SG(nd_psh[i]) == sugar; i-- ) {
     if ( nd_gentrace ) {      if ( nd_gentrace ) {
         /* Trace = [1,index,1,1] */          /* Trace = [1,index,1,1] */
Line 2358  void nd_subst_vector(VL vl,P p,NODE subst,P *r)
Line 3380  void nd_subst_vector(VL vl,P p,NODE subst,P *r)
 NODE nd_gb_trace(int m,int ishomo,int **indp)  NODE nd_gb_trace(int m,int ishomo,int **indp)
 {  {
   int i,nh,sugar,stat;    int i,nh,sugar,stat;
   NODE r,g,t;    NODE g,t;
   ND_pairs d;    ND_pairs d;
   ND_pairs l;    ND_pairs l;
   ND h,nf,nfq,s,head;    ND h,nf,nfq;
   NDV nfv,nfqv;    NDV nfv,nfqv;
   Z q,den,num;    Z q;
   P hc;    P hc;
   union oNDC dn,hnfq;    union oNDC hnfq;
   struct oEGT eg_monic,egm0,egm1;    struct oEGT eg_monic,egm0,egm1;
   int diag_count = 0;    int diag_count = 0;
   P cont;    P cont;
   LIST list;    LIST list;
     struct oHPDATA current_hpdata,final_hpdata;
   
   init_eg(&eg_monic);    init_eg(&eg_monic);
   init_eg(&eg_invdalg);    init_eg(&eg_invdalg);
Line 2380  NODE nd_gb_trace(int m,int ishomo,int **indp)
Line 3403  NODE nd_gb_trace(int m,int ishomo,int **indp)
     g = update_base(g,i);      g = update_base(g,i);
   }    }
   sugar = 0;    sugar = 0;
     if ( nd_hpdata ) {
       if ( DP_Print ) fprintf(asir_out,"Hilbert driven algorithm.\n");
       setup_hpdata(&final_hpdata,&current_hpdata);
     }
   
   while ( d ) {    while ( d ) {
 again:  again:
     l = nd_minp(d,&d);      l = nd_minp(d,&d);
Line 2408  again:
Line 3436  again:
       goto again;        goto again;
     }      }
 #if USE_GEOBUCKET  #if USE_GEOBUCKET
     stat = nd_nf_pbucket(m,h,nd_ps,!Top,&nf);      stat = nd_nf_pbucket(m,h,nd_ps,!nd_top&&!Top,&nf);
 #else  #else
     stat = nd_nf(m,0,h,nd_ps,!Top,&nf);      stat = nd_nf(m,0,h,nd_ps,!nd_top&&!Top,&nf);
 #endif  #endif
     if ( !stat ) {      if ( !stat ) {
       NEXT(l) = d; d = l;        NEXT(l) = d; d = l;
Line 2423  again:
Line 3451  again:
       } else        } else
         nfq = 0;          nfq = 0;
       if ( !nfq ) {        if ( !nfq ) {
         if ( !nd_sp(0,1,l,&h) || !nd_nf(0,0,h,nd_ps_trace,!Top,&nfq) ) {          if ( !nd_sp(0,1,l,&h) || !nd_nf(0,0,h,nd_ps_trace,!nd_top&&!Top,&nfq) ) {
           NEXT(l) = d; d = l;            NEXT(l) = d; d = l;
           d = nd_reconstruct(1,d);            d = nd_reconstruct(1,d);
           goto again;            goto again;
Line 2459  again:
Line 3487  again:
             nd_tracelist = t;              nd_tracelist = t;
           }            }
         }          }
         nh = ndv_newps(0,nfv,nfqv,0);          nh = ndv_newps(0,nfv,nfqv);
         if ( ishomo && ++diag_count == diag_period ) {          if ( ishomo && ++diag_count == diag_period ) {
           diag_count = 0;            diag_count = 0;
           if ( DP_Print > 2 ) fprintf(asir_out,"|");            if ( DP_Print > 2 ) fprintf(asir_out,"|");
Line 2473  again:
Line 3501  again:
         }          }
         d = update_pairs(d,g,nh,0);          d = update_pairs(d,g,nh,0);
         g = update_base(g,nh);          g = update_base(g,nh);
           if ( nd_hpdata ) {
             int dg,sugar0;
   
             update_hpdata(&current_hpdata,nh);
             dg = comp_hn(final_hpdata.hn,current_hpdata.hn);
             if ( dg < 0 ) {
                int d_len;
                for ( d_len = 0; d; d = d->next, d_len++);
                fprintf(asir_out,"[%d] We found a gb\n",d_len);
             }
             sugar0 = sugar;
             while ( d && dg > sugar0 ) {
               d = nd_remove_same_sugar(d,sugar0);
               sugar0++;
             }
           }
       } else {        } else {
         if ( DP_Print ) { printf("*"); fflush(stdout); }          if ( DP_Print ) { fprintf(asir_out,"*"); fflush(asir_out); }
       }        }
     } else {      } else {
       if ( DP_Print ) { printf("."); fflush(stdout); }        if ( DP_Print ) { fprintf(asir_out,"."); fflush(asir_out); }
     }      }
     FREENDP(l);      FREENDP(l);
   }    }
Line 2489  again:
Line 3533  again:
     }      }
   }    }
   conv_ilist(nd_demand,1,g,indp);    conv_ilist(nd_demand,1,g,indp);
   if ( DP_Print ) { printf("nd_gb_trace done.\n"); fflush(stdout); }    if ( DP_Print ) { fprintf(asir_out,"\nnd_gb_trace done. Nbase=%d\n",nd_psn); fflush(asir_out); }
   return g;    return g;
 }  }
   
Line 2518  NODE ndv_reduceall(int m,NODE f)
Line 3562  NODE ndv_reduceall(int m,NODE f)
   int i,j,n,stat;    int i,j,n,stat;
   ND nf,g,head;    ND nf,g,head;
   NODE t,a0,a;    NODE t,a0,a;
   union oNDC dn;  
   Q q,num,den;  
   NODE node;    NODE node;
   LIST l;    LIST l;
   Z iq,jq;    Z iq,jq;
   int *perm;    int *perm;
   union oNDC hc;    union oNDC hc;
   P cont,cont1;    P cont;
   
   if ( nd_nora ) return f;    if ( nd_nora ) return f;
   n = length(f);    n = length(f);
   ndv_setup(m,0,f,0,1);    ndv_setup(m,0,f,0,1,0);
   perm = (int *)MALLOC(n*sizeof(int));    perm = (int *)MALLOC(n*sizeof(int));
   if ( nd_gentrace ) {    if ( nd_gentrace ) {
     for ( t = nd_tracelist, i = 0; i < n; i++, t = NEXT(t) )      for ( t = nd_tracelist, i = 0; i < n; i++, t = NEXT(t) )
       perm[i] = ZTOS((Q)ARG1(BDY((LIST)BDY(t))));        perm[i] = (int)ZTOS((Q)ARG1(BDY((LIST)BDY(t))));
   }    }
   for ( i = 0; i < n; ) {    for ( i = 0; i < n; ) {
     if ( nd_gentrace ) {      if ( nd_gentrace ) {
Line 2552  NODE ndv_reduceall(int m,NODE f)
Line 3594  NODE ndv_reduceall(int m,NODE f)
       hc = HCU(nf); nd_removecont(m,nf);        hc = HCU(nf); nd_removecont(m,nf);
       if ( nd_gentrace ) {        if ( nd_gentrace ) {
         for ( t = nd_tracelist; t; t = NEXT(t) ) {          for ( t = nd_tracelist; t; t = NEXT(t) ) {
           jq = ARG1(BDY((LIST)BDY(t))); j = ZTOS(jq);            jq = ARG1(BDY((LIST)BDY(t))); j = (int)ZTOS(jq);
           STOZ(perm[j],jq); ARG1(BDY((LIST)BDY(t))) = jq;            STOZ(perm[j],jq); ARG1(BDY((LIST)BDY(t))) = jq;
         }          }
         /* exact division */          /* exact division */
Line 2577  NODE ndv_reduceall(int m,NODE f)
Line 3619  NODE ndv_reduceall(int m,NODE f)
   return a0;    return a0;
 }  }
   
   int ndplength(ND_pairs d)
   {
     int i;
     for ( i = 0; d; i++ ) d = NEXT(d);
     return i;
   }
   
 ND_pairs update_pairs( ND_pairs d, NODE /* of index */ g, int t, int gensyz)  ND_pairs update_pairs( ND_pairs d, NODE /* of index */ g, int t, int gensyz)
 {  {
   ND_pairs d1,nd,cur,head,prev,remove;    ND_pairs d1,nd,cur,head,prev,remove;
     int len0;
   
   if ( !g ) return d;    if ( !g ) return d;
   /* for testing */    /* for testing */
Line 2596  ND_pairs update_pairs( ND_pairs d, NODE /* of index */
Line 3646  ND_pairs update_pairs( ND_pairs d, NODE /* of index */
   }    }
   d = crit_B(d,t);    d = crit_B(d,t);
   d1 = nd_newpairs(g,t);    d1 = nd_newpairs(g,t);
     len0 = ndplength(d1);
   d1 = crit_M(d1);    d1 = crit_M(d1);
   d1 = crit_F(d1);    d1 = crit_F(d1);
     NcriMF += len0-ndplength(d1);
   if ( gensyz || do_weyl )    if ( gensyz || do_weyl )
     head = d1;      head = d1;
   else {    else {
Line 2607  ND_pairs update_pairs( ND_pairs d, NODE /* of index */
Line 3659  ND_pairs update_pairs( ND_pairs d, NODE /* of index */
         remove = cur;          remove = cur;
         if ( !prev ) head = cur = NEXT(cur);          if ( !prev ) head = cur = NEXT(cur);
         else cur = NEXT(prev) = NEXT(cur);          else cur = NEXT(prev) = NEXT(cur);
         FREENDP(remove);          FREENDP(remove); Ncri2++;
       } else {        } else {
         prev = cur; cur = NEXT(cur);          prev = cur; cur = NEXT(cur);
       }        }
Line 2623  ND_pairs update_pairs( ND_pairs d, NODE /* of index */
Line 3675  ND_pairs update_pairs( ND_pairs d, NODE /* of index */
   }    }
 }  }
   
   ND_pairs merge_pairs_s(ND_pairs d,ND_pairs d1);
   
   ND_pairs update_pairs_s( ND_pairs d, int t,NODE *syz)
   {
     ND_pairs d1;
     struct oEGT eg1,eg2,eg3;
   
     if ( !t ) return d;
   get_eg(&eg1);
     d1 = nd_newpairs_s(t,syz);
   get_eg(&eg2); add_eg(&eg_create,&eg1,&eg2);
     d = merge_pairs_s(d,d1);
   get_eg(&eg3); add_eg(&eg_merge,&eg2,&eg3);
     return d;
   }
   
   int update_pairs_array_s( ND_pairs *d, int t,NODE *syz)
   {
     ND_pairs *d1;
     struct oEGT eg1,eg2,eg3;
     int i;
   
     if ( !t ) return 0;
   get_eg(&eg1);
     Nnewpair = 0;
     d1 = nd_newpairs_array_s(t,syz);
   get_eg(&eg2); add_eg(&eg_create,&eg1,&eg2);
     for ( i = 0; i < nd_nbase; i++ )
       d[i] = merge_pairs_s(d[i],d1[i]);
   get_eg(&eg3); add_eg(&eg_merge,&eg2,&eg3);
     return Nnewpair;
   }
   
 ND_pairs nd_newpairs( NODE g, int t )  ND_pairs nd_newpairs( NODE g, int t )
 {  {
   NODE h;    NODE h;
Line 2633  ND_pairs nd_newpairs( NODE g, int t )
Line 3717  ND_pairs nd_newpairs( NODE g, int t )
   
   dl = DL(nd_psh[t]);    dl = DL(nd_psh[t]);
   ts = SG(nd_psh[t]) - TD(dl);    ts = SG(nd_psh[t]) - TD(dl);
   if ( nd_module && nd_intersect && (MPOS(dl) > nd_intersect) ) return 0;    if ( nd_module && nd_intersect && (MPOS(dl) > (UINT)nd_intersect) ) return 0;
   for ( r0 = 0, h = g; h; h = NEXT(h) ) {    for ( r0 = 0, h = g; h; h = NEXT(h) ) {
     if ( nd_module && (MPOS(DL(nd_psh[(long)BDY(h)])) != MPOS(dl)) )      if ( nd_module && (MPOS(DL(nd_psh[(long)BDY(h)])) != MPOS(dl)) )
       continue;        continue;
Line 2647  ND_pairs nd_newpairs( NODE g, int t )
Line 3731  ND_pairs nd_newpairs( NODE g, int t )
       if ( nd_gbblock[i] >= 0 )        if ( nd_gbblock[i] >= 0 )
         continue;          continue;
     }      }
     NEXTND_pairs(r0,r);      NEXTND_pairs(r0,r); Npairs++;
     r->i1 = (long)BDY(h);      r->i1 = (long)BDY(h);
     r->i2 = t;      r->i2 = t;
     ndl_lcm(DL(nd_psh[r->i1]),dl,r->lcm);      ndl_lcm(DL(nd_psh[r->i1]),dl,r->lcm);
Line 2661  ND_pairs nd_newpairs( NODE g, int t )
Line 3745  ND_pairs nd_newpairs( NODE g, int t )
   return r0;    return r0;
 }  }
   
   int sig_cmpdl_op(int n,DL d1,DL d2)
   {
     int e1,e2,i,j,l;
     int *t1,*t2;
     int len;
     struct order_pair *pair;
   
     len = nd_sba_modord->block_length;
     pair = nd_sba_modord->order_pair;
   
     for ( i = 0, t1 = d1->d, t2 = d2->d; i < len; i++ ) {
       l = pair[i].length;
       switch ( pair[i].order ) {
         case 0:
           for ( j = 0, e1 = e2 = 0; j < l; j++ ) {
             e1 += t1[j];
             e2 += t2[j];
           }
           if ( e1 > e2 )
             return 1;
           else if ( e1 < e2 )
             return -1;
           else {
             for ( j = l - 1; j >= 0 && t1[j] == t2[j]; j-- );
             if ( j >= 0 )
               return t1[j] < t2[j] ? 1 : -1;
           }
           break;
         case 1:
           for ( j = 0, e1 = e2 = 0; j < l; j++ ) {
             e1 += t1[j];
             e2 += t2[j];
           }
           if ( e1 > e2 )
             return 1;
           else if ( e1 < e2 )
             return -1;
           else {
             for ( j = 0; j < l && t1[j] == t2[j]; j++ );
             if ( j < l )
               return t1[j] > t2[j] ? 1 : -1;
           }
           break;
         case 2:
           for ( j = 0; j < l && t1[j] == t2[j]; j++ );
           if ( j < l )
             return t1[j] > t2[j] ? 1 : -1;
           break;
         default:
           error("sig_cmpdl_op : invalid order"); break;
       }
       t1 += l; t2 += l; // head += l;
     }
     return 0;
   }
   
   int sig_cmpdl_mat(int n,DL d1,DL d2)
   {
     int *v,*t1,*t2;
     int s,i,j,len;
     int **matrix;
     static int *w;
     static int nvar = 0;
   
     if ( nvar != n ) {
       nvar = n; w = (int *)MALLOC(n*sizeof(int));
     }
     for ( i = 0, t1 = d1->d, t2 = d2->d; i < n; i++ )
       w[i] = t1[i]-t2[i];
     len = nd_sba_modord->row;
     matrix = nd_sba_modord->matrix;
     for ( j = 0; j < len; j++ ) {
       v = matrix[j];
       for ( i = 0, s = 0; i < n; i++ )
         s += v[i]*w[i];
       if ( s > 0 )
         return 1;
       else if ( s < 0 )
         return -1;
     }
     return 0;
   }
   
   struct comp_sig_spec *create_comp_sig_spec(VL current_vl,VL old_vl,Obj ord,Obj weight)
   {
     struct comp_sig_spec *spec;
     VL ovl,vl;
     V ov;
     int i,j,n,nvar,s;
     NODE node,t,tn;
     struct order_pair *l;
     MAT m;
     Obj **b;
     int **w;
     int *a;
   
     spec = (struct comp_sig_spec *)MALLOC(sizeof(struct comp_sig_spec));
     for ( i = 0, vl = current_vl; vl; vl = NEXT(vl), i++ );
     spec->n = nvar = i;
     if ( old_vl != 0 ) {
       spec->oldv = (int *)MALLOC(nvar*sizeof(int));
       for ( i = 0, ovl = old_vl; i < nvar; ovl = NEXT(ovl), i++ ) {
         ov = ovl->v;
         for ( j = 0, vl = current_vl; vl; vl = NEXT(vl), j++ )
           if ( ov == vl->v ) break;
         spec->oldv[i] = j;
       }
     } else
       spec->oldv = 0;
     if ( !ord || NUM(ord) ) {
       switch ( ZTOS((Z)ord) ) {
         case 0:
           spec->cmpdl = cmpdl_revgradlex; break;
         case 1:
           spec->cmpdl = cmpdl_gradlex; break;
         case 2:
           spec->cmpdl = cmpdl_lex; break;
         default:
           error("create_comp_sig_spec : invalid spec"); break;
       }
     } else if ( OID(ord) == O_LIST ) {
       node = BDY((LIST)ord);
       for ( n = 0, t = node; t; t = NEXT(t), n++ );
       l = (struct order_pair *)MALLOC_ATOMIC(n*sizeof(struct order_pair));
       for ( i = 0, t = node, s = 0; i < n; t = NEXT(t), i++ ) {
         tn = BDY((LIST)BDY(t)); l[i].order = (int)ZTOS((Q)BDY(tn));
         tn = NEXT(tn); l[i].length = (int)ZTOS((Q)BDY(tn));
         s += l[i].length;
       }
       if ( s != nvar )
         error("create_comp_sig_spec : invalid spec");
       spec->order_pair = l;
       spec->block_length = n;
       spec->cmpdl = sig_cmpdl_op;
     } else if ( OID(ord) == O_MAT ) {
       m = (MAT)ord; b = (Obj **)BDY(m);
       if ( m->col != nvar )
         error("create_comp_sig_spec : invalid spec");
       w = almat(m->row,m->col);
       for ( i = 0; i < m->row; i++ )
         for ( j = 0; j < m->col; j++ )
           w[i][j] = (int)ZTOS((Q)b[i][j]);
       spec->row = m->row;
       spec->matrix = w;
       spec->cmpdl = sig_cmpdl_mat;
     } else
       error("create_comp_sig_spec : invalid spec");
     if ( weight != 0 ) {
       node = BDY((LIST)weight);
       a = (int *)MALLOC(nvar*sizeof(int));
       for ( i = 0; i < nvar; i++, node = NEXT(node) )
         a[i] = (int)ZTOS((Z)BDY(node));
       spec->weight = a;
     }
     return spec;
   }
   
   #define SIG_MUL_WEIGHT(a,i) (weight?(a)*weight[i]:(a))
   
   int comp_sig_monomial(int n,DL d1,DL d2)
   {
     static DL m1,m2;
     static int nvar = 0;
     int *oldv,*weight;
     int i,w1,w2;
   
     if ( nvar != n ) {
       nvar = n; NEWDL(m1,nvar); NEWDL(m2,nvar);
     }
     if ( !nd_sba_modord )
       return (*cmpdl)(n,d1,d2);
     else if ( !nd_sba_modord->weight && !nd_sba_modord->oldv ) {
       return (*nd_sba_modord->cmpdl)(n,d1,d2);
     } else {
       weight = nd_sba_modord->weight;
       oldv = nd_sba_modord->oldv;
       if ( oldv ) {
         for ( i = 0; i < n; i++ ) {
           m1->d[i] = d1->d[oldv[i]]; m2->d[i] = d2->d[oldv[i]];
         }
       } else {
         for ( i = 0; i < n; i++ ) {
           m1->d[i] = d1->d[i]; m2->d[i] = d2->d[i];
         }
       }
       for ( i = 0, w1 = w2 = 0; i < n; i++ ) {
         w1 += SIG_MUL_WEIGHT(m1->d[i],i);
         w2 += SIG_MUL_WEIGHT(m2->d[i],i);
       }
       m1->td = w1; m2->td = w2;
       return (*nd_sba_modord->cmpdl)(n,m1,m2);
     }
   }
   
   int comp_sig(SIG s1,SIG s2)
   {
     if ( nd_sba_pot ) {
       if ( s1->pos > s2->pos ) return 1;
       else if ( s1->pos < s2->pos ) return -1;
       else return comp_sig_monomial(nd_nvar,DL(s1),DL(s2));
     } else {
       static DL m1,m2;
       static int nvar = 0;
       int ret;
   
       if ( nvar != nd_nvar ) {
         nvar = nd_nvar; NEWDL(m1,nvar); NEWDL(m2,nvar);
       }
       if ( !nd_sba_modord )
         ret = (*cmpdl)(nd_nvar,DL2(s1),DL2(s2));
       else
         ret = comp_sig_monomial(nd_nvar,DL2(s1),DL2(s2));
       if ( ret != 0 ) return ret;
       else if ( s1->pos > s2->pos ) return 1;
       else if ( s1->pos < s2->pos ) return -1;
       else return 0;
     }
   }
   
   int _create_spair_s(int i1,int i2,ND_pairs sp,SIG sig1,SIG sig2)
   {
     int ret,s1,s2;
     RHist p1,p2;
     static int wpd;
     static UINT *lcm;
   
     sp->i1 = i1;
     sp->i2 = i2;
     p1 = nd_psh[i1];
     p2 = nd_psh[i2];
     ndl_lcm(DL(p1),DL(p2),sp->lcm);
   #if 0
     s1 = SG(p1)-TD(DL(p1));
     s2 = SG(p2)-TD(DL(p2));
     SG(sp) = MAX(s1,s2) + TD(sp->lcm);
   #endif
   
     if ( wpd != nd_wpd ) {
       wpd = nd_wpd;
       lcm = (UINT *)MALLOC(wpd*sizeof(UINT));
     }
     // DL(sig1) <- sp->lcm
     // DL(sig1) -= DL(p1)
     // DL(sig1) += DL(p1->sig)
     ndl_sub(sp->lcm,DL(p1),lcm);
     _ndltodl(lcm,DL(sig1));
     _addtodl(nd_nvar,DL(p1->sig),DL(sig1));
     sig1->pos = p1->sig->pos;
     _adddl(nd_nvar,DL(sig1),nd_sba_hm[sig1->pos],DL2(sig1));
   
     // DL(sig2) <- sp->lcm
     // DL(sig2) -= DL(p2)
     // DL(sig2) += DL(p2->sig)
     ndl_sub(sp->lcm,DL(p2),lcm);
     _ndltodl(lcm,DL(sig2));
     _addtodl(nd_nvar,DL(p2->sig),DL(sig2));
     sig2->pos = p2->sig->pos;
     _adddl(nd_nvar,DL(sig2),nd_sba_hm[sig2->pos],DL2(sig2));
   
     ret = comp_sig(sig1,sig2);
     if ( ret == 0 ) return 0;
     else if ( ret > 0 ) sp->sig = sig1;
     else sp->sig = sig2;
   
     s1 = DL(sig1)->td+nd_sba_hm[p1->sig->pos]->td;
     s2 = DL(sig2)->td+nd_sba_hm[p2->sig->pos]->td;
     SG(sp) = MAX(s1,s2);
   
     return 1;
   }
   
   SIG dup_sig(SIG sig)
   {
     SIG r;
   
     if ( !sig ) return 0;
     else {
       NEWSIG(r);
       _copydl(nd_nvar,DL(sig),DL(r));
       _copydl(nd_nvar,DL2(sig),DL2(r));
       r->pos = sig->pos;
       return r;
     }
   }
   
   void dup_ND_pairs(ND_pairs to,ND_pairs from)
   {
     to->i1 = from->i1;
     to->i2 = from->i2;
     to->sugar = from->sugar;
     to->sugar2 = from->sugar2;
     ndl_copy(from->lcm,to->lcm);
     to->sig = dup_sig(from->sig);
   }
   
   ND_pairs merge_pairs_s(ND_pairs p1,ND_pairs p2)
   {
     struct oND_pairs root;
     ND_pairs q1,q2,r;
     int ret;
   
     r = &root;
     for ( q1 = p1, q2 = p2; q1 != 0 && q2 != 0; ) {
       ret = comp_sig(q1->sig,q2->sig);
       if ( ret < 0 ) {
         r->next = q1; r = q1; q1 = q1->next;
       } else if ( ret > 0 ) {
         r->next = q2; r = q2; q2 = q2->next;
       } else {
         Nnewpair--;
         ret = DL_COMPARE(q1->lcm,q2->lcm);
         Nsamesig++;
         if ( ret < 0 ) {
           r->next = q1; r = q1; q1 = q1->next;
           q2 = q2->next;
         } else {
           r->next = q2; r = q2; q2 = q2->next;
           q1 = q1->next;
         }
       }
     }
     if ( q1 ) {
       r->next = q1;
     } else {
       r->next = q2;
     }
     return root.next;
   }
   
   ND_pairs insert_pair_s(ND_pairs l,ND_pairs s)
   {
     ND_pairs p,prev;
     int ret=1;
   
     for ( p = l, prev = 0; p != 0; prev = p, p = p->next ) {
       if ( (ret = comp_sig(s->sig,p->sig)) <= 0 )
         break;
     }
     if ( ret == 0 ) {
       ret = DL_COMPARE(s->lcm,p->lcm);
       if ( ret < 0 ) {
         // replace p with s
         s->next = p->next;
         if ( prev == 0 ) {
           return s;
         } else {
           prev->next = s;
           return l;
         }
       } else
         return l;
     } else {
       // insert s between prev and p
       Nnewpair++;
       s->next = p;
       if ( prev == 0 ) {
         return s;
       } else {
         prev->next = s;
         return l;
       }
     }
   }
   
   INLINE int __dl_redble(DL d1,DL d2,int nvar)
   {
     int i;
   
     if ( d1->td > d2->td )
       return 0;
     for ( i = nvar-1; i >= 0; i-- )
       if ( d1->d[i] > d2->d[i] )
         break;
     if ( i >= 0 )
       return 0;
     else
       return 1;
   }
   
   ND_pairs nd_newpairs_s(int t, NODE *syz)
   {
     NODE s;
     int ret,i;
     ND_pairs r0,_sp,sp;
     SIG spsig,tsig;
     static int nvar = 0;
     static SIG _sig1,_sig2;
   
     NEWND_pairs(_sp);
     if ( !_sig1 || nvar != nd_nvar ) {
       nvar = nd_nvar; NEWSIG(_sig1); NEWSIG(_sig2);
     }
     r0 = 0;
     for ( i = 0; i < t; i++ ) {
       ret = _create_spair_s(i,t,_sp,_sig1,_sig2);
       if ( ret ) {
         spsig = _sp->sig;
         for ( s = syz[spsig->pos]; s; s = s->next ) {
           tsig = (SIG)s->body;
           if ( _dl_redble(DL(tsig),DL(spsig),nd_nvar) )
             break;
         }
         if ( s == 0 ) {
           NEWND_pairs(sp);
           dup_ND_pairs(sp,_sp);
           r0 = insert_pair_s(r0,sp);
         } else
           Nsyz++;
       }
     }
     return r0;
   }
   
   ND_pairs *nd_newpairs_array_s(int t, NODE *syz)
   {
     NODE s;
     int ret,i;
     ND_pairs _sp,sp;
     ND_pairs *d;
     SIG spsig,tsig;
     static int nvar = 0;
     static SIG _sig1,_sig2;
   
     NEWND_pairs(_sp);
     if ( !_sig1 || nvar != nd_nvar ) {
       nvar = nd_nvar; NEWSIG(_sig1); NEWSIG(_sig2);
     }
     d = (ND_pairs *)MALLOC(nd_nbase*sizeof(ND_pairs));
     Nnewpair = 0;
     for ( i = 0; i < t; i++ ) {
       ret = _create_spair_s(i,t,_sp,_sig1,_sig2);
       if ( ret ) {
         spsig = _sp->sig;
         for ( s = syz[spsig->pos]; s; s = s->next ) {
           tsig = (SIG)s->body;
           if ( _dl_redble(DL(tsig),DL(spsig),nd_nvar) )
             break;
         }
         if ( s == 0 ) {
           NEWND_pairs(sp);
           dup_ND_pairs(sp,_sp);
           d[spsig->pos] = insert_pair_s(d[spsig->pos],sp);
         } else
           Nsyz++;
       }
     }
     return d;
   }
   
 /* ipair = [i1,i2],[i1,i2],... */  /* ipair = [i1,i2],[i1,i2],... */
 ND_pairs nd_ipairtospair(NODE ipair)  ND_pairs nd_ipairtospair(NODE ipair)
 {  {
Line 2671  ND_pairs nd_ipairtospair(NODE ipair)
Line 4204  ND_pairs nd_ipairtospair(NODE ipair)
   for ( r0 = 0, t = ipair; t; t = NEXT(t) ) {    for ( r0 = 0, t = ipair; t; t = NEXT(t) ) {
     NEXTND_pairs(r0,r);      NEXTND_pairs(r0,r);
     tn = BDY((LIST)BDY(t));      tn = BDY((LIST)BDY(t));
     r->i1 = ZTOS((Q)ARG0(tn));      r->i1 = (int)ZTOS((Q)ARG0(tn));
     r->i2 = ZTOS((Q)ARG1(tn));      r->i2 = (int)ZTOS((Q)ARG1(tn));
     ndl_lcm(DL(nd_psh[r->i1]),DL(nd_psh[r->i2]),r->lcm);      ndl_lcm(DL(nd_psh[r->i1]),DL(nd_psh[r->i2]),r->lcm);
     s1 = SG(nd_psh[r->i1])-TD(DL(nd_psh[r->i1]));      s1 = SG(nd_psh[r->i1])-TD(DL(nd_psh[r->i1]));
     s2 = SG(nd_psh[r->i2])-TD(DL(nd_psh[r->i2]));      s2 = SG(nd_psh[r->i2])-TD(DL(nd_psh[r->i2]));
Line 2691  ND_pairs crit_B( ND_pairs d, int s )
Line 4224  ND_pairs crit_B( ND_pairs d, int s )
 {  {
   ND_pairs cur,head,prev,remove;    ND_pairs cur,head,prev,remove;
   UINT *t,*tl,*lcm;    UINT *t,*tl,*lcm;
   int td,tdl;  
   
   if ( !d ) return 0;    if ( !d ) return 0;
   t = DL(nd_psh[s]);    t = DL(nd_psh[s]);
Line 2711  ND_pairs crit_B( ND_pairs d, int s )
Line 4243  ND_pairs crit_B( ND_pairs d, int s )
           } else {            } else {
             cur = NEXT(prev) = NEXT(cur);              cur = NEXT(prev) = NEXT(cur);
           }            }
           FREENDP(remove);            FREENDP(remove); NcriB++;
         } else {          } else {
           prev = cur; cur = NEXT(cur);            prev = cur; cur = NEXT(cur);
         }          }
Line 2857  NODE update_base(NODE nd,int ndp)
Line 4389  NODE update_base(NODE nd,int ndp)
 ND_pairs nd_minp( ND_pairs d, ND_pairs *prest )  ND_pairs nd_minp( ND_pairs d, ND_pairs *prest )
 {  {
     ND_pairs m,ml,p,l;      ND_pairs m,ml,p,l;
     UINT *lcm;      int s;
     int s,td,len,tlen,c,c1;  
   
     if ( !(p = NEXT(m = d)) ) {      if ( !(p = NEXT(m = d)) ) {
         *prest = p;          *prest = p;
Line 2938  ND_pairs nd_minsugarp( ND_pairs d, ND_pairs *prest )
Line 4469  ND_pairs nd_minsugarp( ND_pairs d, ND_pairs *prest )
     return dm0;      return dm0;
 }  }
   
   ND_pairs nd_minsugarp_s( ND_pairs d, ND_pairs *prest )
   {
     int msugar;
     ND_pairs t,last;
   
     for ( msugar = SG(d), t = d; t; t = NEXT(t) )
       if ( SG(t) == msugar ) last = t;
     *prest = last->next;
     last->next = 0;
     return d;
   }
   
 int nd_tdeg(NDV c)  int nd_tdeg(NDV c)
 {  {
   int wmax = 0;    int wmax = 0;
Line 2946  int nd_tdeg(NDV c)
Line 4489  int nd_tdeg(NDV c)
   
   len = LEN(c);    len = LEN(c);
   for ( a = BDY(c), i = 0; i < len; i++, NMV_ADV(a) )    for ( a = BDY(c), i = 0; i < len; i++, NMV_ADV(a) )
     wmax = MAX(TD(DL(a)),wmax);      wmax = MAX((int)TD(DL(a)),wmax);
   return wmax;    return wmax;
 }  }
   
 int ndv_newps(int m,NDV a,NDV aq,int f4)  int ndv_newps(int m,NDV a,NDV aq)
 {  {
     int len;  
     RHist r;      RHist r;
     NDV b;  
     NODE tn;      NODE tn;
     LIST l;      LIST l;
     Z iq;      Z iq;
Line 2986  int ndv_newps(int m,NDV a,NDV aq,int f4)
Line 4527  int ndv_newps(int m,NDV a,NDV aq,int f4)
         SG(r) = nd_tdeg(aq);          SG(r) = nd_tdeg(aq);
 #endif  #endif
         ndl_copy(HDL(aq),DL(r));          ndl_copy(HDL(aq),DL(r));
           r->sig = dup_sig(aq->sig);
     } else {      } else {
         if ( !m ) register_hcf(a);          if ( !m ) register_hcf(a);
         nd_bound[nd_psn] = ndv_compute_bound(a);          nd_bound[nd_psn] = ndv_compute_bound(a);
Line 2995  int ndv_newps(int m,NDV a,NDV aq,int f4)
Line 4537  int ndv_newps(int m,NDV a,NDV aq,int f4)
         SG(r) = nd_tdeg(a);          SG(r) = nd_tdeg(a);
 #endif  #endif
         ndl_copy(HDL(a),DL(r));          ndl_copy(HDL(a),DL(r));
           r->sig = dup_sig(a->sig);
     }      }
     if ( nd_demand ) {      if ( nd_demand ) {
         if ( aq ) {          if ( aq ) {
Line 3017  int ndv_newps(int m,NDV a,NDV aq,int f4)
Line 4560  int ndv_newps(int m,NDV a,NDV aq,int f4)
     return nd_psn++;      return nd_psn++;
 }  }
   
   // find LM wrt the specified modord
   void ndv_lm_modord(NDV p,DL d)
   {
     NMV m;
     DL tmp;
     int len,i,ret;
   
     NEWDL(tmp,nd_nvar);
     m = BDY(p); len = LEN(p);
     _ndltodl(DL(m),d); // printdl(d); printf("->");
     for ( i = 1, NMV_ADV(m); i < len; i++, NMV_ADV(m) ) {
       _ndltodl(DL(m),tmp);
       ret = comp_sig_monomial(nd_nvar,tmp,d);
       if ( ret > 0 ) _copydl(nd_nvar,tmp,d);
     }
   //   printdl(d); printf("\n");
   }
   
 /* nd_tracelist = [[0,index,div],...,[nd_psn-1,index,div]] */  /* nd_tracelist = [[0,index,div],...,[nd_psn-1,index,div]] */
 /* return 1 if success, 0 if failure (HC(a mod p)) */  /* return 1 if success, 0 if failure (HC(a mod p)) */
   
 int ndv_setup(int mod,int trace,NODE f,int dont_sort,int dont_removecont)  int ndv_setup(int mod,int trace,NODE f,int dont_sort,int dont_removecont,int sba)
 {  {
   int i,j,td,len,max;    int i,j;
   NODE s,s0,f0,tn;    NODE s,tn;
   UINT *d;  
   RHist r;    RHist r;
   NDVI w;    NDVI w;
   NDV a,am;    NDV a,am;
   union oNDC hc;    union oNDC hc;
   NODE node;    NODE node;
   P hcp;  
   Z iq,jq;    Z iq,jq;
   LIST l;    LIST l;
   
Line 3044  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
Line 4603  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
     if ( BDY(s) ) { w[i].p = BDY(s); w[i].i = j; i++; }      if ( BDY(s) ) { w[i].p = BDY(s); w[i].i = j; i++; }
   if ( !dont_sort ) {    if ( !dont_sort ) {
     /* XXX heuristic */      /* XXX heuristic */
     if ( !nd_ord->id && (nd_ord->ord.simple<2) )      if ( !sba && !nd_ord->id && (nd_ord->ord.simple<2) )
       qsort(w,nd_psn,sizeof(struct oNDVI),        qsort(w,nd_psn,sizeof(struct oNDVI),
         (int (*)(const void *,const void *))ndvi_compare_rev);          (int (*)(const void *,const void *))ndvi_compare_rev);
     else      else
Line 3079  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
Line 4638  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
         if ( !dont_removecont) ndv_removecont(-2,a);          if ( !dont_removecont) ndv_removecont(-2,a);
         am = nd_ps[i] = ndv_dup(trace,w[i].p);          am = nd_ps[i] = ndv_dup(trace,w[i].p);
         ndv_mod(trace,am);          ndv_mod(trace,am);
       if ( DL_COMPARE(HDL(am),HDL(a)) )          if ( DL_COMPARE(HDL(am),HDL(a)) ) return 0;
         return 0;  
         ndv_removecont(trace,am);          ndv_removecont(trace,am);
       } else {        } else {
         a = nd_ps_trace[i] = ndv_dup(0,w[i].p);          a = nd_ps_trace[i] = ndv_dup(0,w[i].p);
Line 3088  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
Line 4646  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
         register_hcf(a);          register_hcf(a);
         am = nd_ps[i] = ndv_dup(mod,a);          am = nd_ps[i] = ndv_dup(mod,a);
         ndv_mod(mod,am);          ndv_mod(mod,am);
       if ( DL_COMPARE(HDL(am),HDL(a)) )          if ( DL_COMPARE(HDL(am),HDL(a)) ) return 0;
         return 0;  
         ndv_removecont(mod,am);          ndv_removecont(mod,am);
       }        }
     } else {      } else {
Line 3120  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
Line 4677  int ndv_setup(int mod,int trace,NODE f,int dont_sort,i
       }        }
     }      }
   }    }
     if ( sba ) {
       nd_sba_hm = (DL *)MALLOC(nd_psn*sizeof(DL));
      // setup signatures
       for ( i = 0; i < nd_psn; i++ ) {
         SIG sig;
   
         NEWSIG(sig); sig->pos = i;
         nd_ps[i]->sig = sig;
         if ( nd_demand ) nd_ps_sym[i]->sig = sig;
           nd_psh[i]->sig = sig;
         if ( trace ) {
           nd_ps_trace[i]->sig = sig;
           if ( nd_demand ) nd_ps_trace_sym[i]->sig = sig;
         }
         NEWDL(nd_sba_hm[i],nd_nvar);
         if ( nd_sba_modord )
           ndv_lm_modord(nd_ps[i],nd_sba_hm[i]);
         else
           _ndltodl(DL(nd_psh[i]),nd_sba_hm[i]);
         _adddl(nd_nvar,DL(sig),nd_sba_hm[i],DL2(sig));
       }
       nd_sba_pos = (NODE *)MALLOC(nd_psn*sizeof(NODE));
       for ( i = 0; i < nd_psn; i++ ) {
         j = nd_psh[i]->sig->pos;
         nd_sba_pos[j] = append_one(nd_sba_pos[j],i);
       }
     }
   if ( nd_gentrace && nd_tracelist ) NEXT(tn) = 0;    if ( nd_gentrace && nd_tracelist ) NEXT(tn) = 0;
   return 1;    return 1;
 }  }
Line 3133  static VECT prev_weight_vector_obj;
Line 4717  static VECT prev_weight_vector_obj;
 void preprocess_algcoef(VL vv,VL av,struct order_spec *ord,LIST f,  void preprocess_algcoef(VL vv,VL av,struct order_spec *ord,LIST f,
     struct order_spec **ord1p,LIST *f1p,NODE *alistp)      struct order_spec **ord1p,LIST *f1p,NODE *alistp)
 {  {
     NODE alist,t,s,r0,r,arg;      NODE alist,t,s,arg;
     VL tv;      VL tv;
     P poly;      P poly;
     DP d;      DP d;
Line 3142  void preprocess_algcoef(VL vv,VL av,struct order_spec 
Line 4726  void preprocess_algcoef(VL vv,VL av,struct order_spec 
     MP m;      MP m;
     int i,nvar,nalg,n;      int i,nvar,nalg,n;
     NumberField nf;      NumberField nf;
     LIST f1,f2;      LIST f1;
     struct order_spec *current_spec;  
     VECT obj,obj0;      VECT obj,obj0;
     VECT tmp;      VECT tmp;
   
Line 3232  NODE postprocess_algcoef(VL av,NODE alist,NODE r)
Line 4815  NODE postprocess_algcoef(VL av,NODE alist,NODE r)
   
 void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int f4,struct order_spec *ord,LIST *rp)  void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int f4,struct order_spec *ord,LIST *rp)
 {  {
     VL tv,fv,vv,vc,av;    VL tv,fv,vv,av;
     NODE fd,fd0,r,r0,t,x,s,xx,alist;    NODE fd,fd0,r0,t,x,s,xx,alist,r;
     int e,max,nvar,i;    int e,max,nvar,i;
     NDV b;    NDV b;
     int ishomo,nalg,mrank,trank,wmax,len;    int ishomo,nalg,mrank,trank,wmax,len;
   NMV a;    NMV a;
     Alg alpha,dp;    P zp;
     P p,zp;    Q dmy;
     Q dmy;    LIST f1,zpl;
     LIST f1,f2,zpl;    struct order_spec *ord1;
     Obj obj;    NODE tr,tl1,tl2,tl3,tl4,nzlist;
     NumberField nf;    LIST l1,l2,l3,l4,l5;
     struct order_spec *ord1;  
     NODE tr,tl1,tl2,tl3,tl4,nzlist;  
     LIST l1,l2,l3,l4,l5;  
   int j;    int j;
   Z jq,bpe,last_nonzero;    Z jq,bpe,last_nonzero;
     int *perm;    int *perm;
     EPOS oepos;    EPOS oepos;
     int obpe,oadv,ompos,cbpe;    int obpe,oadv,ompos,cbpe;
     VECT hvect;    VECT hvect;
   
     nd_module = 0;    NcriB = NcriMF = Ncri2 = 0;
     if ( !m && Demand ) nd_demand = 1;    nd_module = 0;
     else nd_demand = 0;    if ( !m && Demand ) nd_demand = 1;
     parse_nd_option(current_option);    else nd_demand = 0;
   
     if ( DP_Multiple )    if ( DP_Multiple )
         nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);      nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);
 #if 0  #if 0
     ndv_alloc = 0;    ndv_alloc = 0;
 #endif  #endif
     get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
     if ( m && nd_vc )    parse_nd_option(vv,current_option);
        error("nd_{gr,f4} : computation over Fp(X) is unsupported. Use dp_gr_mod_main().");    if ( m && nd_vc )
     for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );      error("nd_{gr,f4} : computation over Fp(X) is unsupported. Use dp_gr_mod_main().");
     switch ( ord->id ) {    for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
         case 1:    switch ( ord->id ) {
             if ( ord->nv != nvar )      case 1:
                 error("nd_{gr,f4} : invalid order specification");        if ( ord->nv != nvar )
             break;          error("nd_{gr,f4} : invalid order specification");
         default:        break;
             break;      default:
         break;
     }
     nd_nalg = 0;
     av = 0;
     if ( !m ) {
       get_algtree((Obj)f,&av);
       for ( nalg = 0, tv = av; tv; tv = NEXT(tv), nalg++ );
       nd_ntrans = nvar;
       nd_nalg = nalg;
       /* #i -> t#i */
       if ( nalg ) {
         preprocess_algcoef(vv,av,ord,f,&ord1,&f1,&alist);
         ord = ord1;
         f = f1;
     }      }
     nd_nalg = 0;      nvar += nalg;
     av = 0;    }
     if ( !m ) {    nd_init_ord(ord);
         get_algtree((Obj)f,&av);    mrank = 0;
         for ( nalg = 0, tv = av; tv; tv = NEXT(tv), nalg++ );    for ( t = BDY(f), max = 1; t; t = NEXT(t) )
         nd_ntrans = nvar;      for ( tv = vv; tv; tv = NEXT(tv) ) {
         nd_nalg = nalg;  
         /* #i -> t#i */  
         if ( nalg ) {  
             preprocess_algcoef(vv,av,ord,f,&ord1,&f1,&alist);  
             ord = ord1;  
             f = f1;  
         }  
         nvar += nalg;  
     }  
     nd_init_ord(ord);  
     mrank = 0;  
     for ( t = BDY(f), max = 1; t; t = NEXT(t) )  
         for ( tv = vv; tv; tv = NEXT(tv) ) {  
             if ( nd_module ) {  
                 if ( OID(BDY(t)) == O_DPM ) {  
                   e = dpm_getdeg((DPM)BDY(t),&trank);  
                   max = MAX(e,max);  
                   mrank = MAX(mrank,trank);  
                 } else {  
                   s = BDY((LIST)BDY(t));  
                   trank = length(s);  
                   mrank = MAX(mrank,trank);  
                   for ( ; s; s = NEXT(s) ) {  
                       e = getdeg(tv->v,(P)BDY(s));  
                       max = MAX(e,max);  
                   }  
                 }  
             } else {  
                 e = getdeg(tv->v,(P)BDY(t));  
                 max = MAX(e,max);  
             }  
         }  
     nd_setup_parameters(nvar,nd_nzlist?0:max);  
     obpe = nd_bpe; oadv = nmv_adv; oepos = nd_epos; ompos = nd_mpos;  
     ishomo = 1;  
     for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {  
       if ( nd_module ) {        if ( nd_module ) {
         if ( OID(BDY(t)) == O_DPM ) {          if ( OID(BDY(t)) == O_DPM ) {
           Z cont;            e = dpm_getdeg((DPM)BDY(t),&trank);
           DPM zdpm;            max = MAX(e,max);
             mrank = MAX(mrank,trank);
           if ( !m && !nd_gentrace ) dpm_ptozp((DPM)BDY(t),&cont,&zdpm);  
           else zdpm = (DPM)BDY(t);  
           b = (pointer)dpmtondv(m,zdpm);  
         } else {          } else {
           if ( !m && !nd_gentrace ) pltozpl((LIST)BDY(t),&dmy,&zpl);            s = BDY((LIST)BDY(t));
           else zpl = (LIST)BDY(t);            trank = length(s);
           b = (pointer)pltondv(CO,vv,zpl);            mrank = MAX(mrank,trank);
             for ( ; s; s = NEXT(s) ) {
               e = getdeg(tv->v,(P)BDY(s));
               max = MAX(e,max);
             }
         }          }
       } else {        } else {
           if ( OID(BDY(t)) == O_DP ) {
             e = dp_getdeg((DP)BDY(t));
             max = MAX(e,max);
           } else {
             e = getdeg(tv->v,(P)BDY(t));
             max = MAX(e,max);
           }
         }
       }
     nd_setup_parameters(nvar,nd_nzlist?0:max);
     obpe = nd_bpe; oadv = nmv_adv; oepos = nd_epos; ompos = nd_mpos;
     ishomo = 1;
     for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {
       if ( nd_module ) {
         if ( OID(BDY(t)) == O_DPM ) {
           Z cont;
           DPM zdpm;
   
           if ( !m && !nd_gentrace ) dpm_ptozp((DPM)BDY(t),&cont,&zdpm);
           else zdpm = (DPM)BDY(t);
           b = (pointer)dpmtondv(m,zdpm);
         } else {
           if ( !m && !nd_gentrace ) pltozpl((LIST)BDY(t),&dmy,&zpl);
           else zpl = (LIST)BDY(t);
           b = (pointer)pltondv(CO,vv,zpl);
         }
       } else {
         if ( OID(BDY(t)) == O_DP ) {
           DP zdp;
   
           if ( !m && !nd_gentrace ) dp_ptozp((DP)BDY(t),&zdp);
           else zdp = (DP)BDY(t);
           b = (pointer)dptondv(m,zdp);
         } else {
         if ( !m && !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);          if ( !m && !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);
         else zp = (P)BDY(t);          else zp = (P)BDY(t);
         b = (pointer)ptondv(CO,vv,zp);          b = (pointer)ptondv(CO,vv,zp);
       }        }
       if ( ishomo )  
         ishomo = ishomo && ndv_ishomo(b);  
       if ( m ) ndv_mod(m,b);  
       if ( b ) { NEXTNODE(fd0,fd); BDY(fd) = (pointer)b; }  
     }      }
     if ( fd0 ) NEXT(fd) = 0;      if ( ishomo )
         ishomo = ishomo && ndv_ishomo(b);
       if ( m ) ndv_mod(m,b);
       if ( b ) { NEXTNODE(fd0,fd); BDY(fd) = (pointer)b; }
     }
     if ( fd0 ) NEXT(fd) = 0;
   
   if ( !ishomo && homo ) {    if ( !ishomo && homo ) {
         for ( t = fd0, wmax = max; t; t = NEXT(t) ) {          for ( t = fd0, wmax = max; t; t = NEXT(t) ) {
             b = (NDV)BDY(t); len = LEN(b);              b = (NDV)BDY(t); len = LEN(b);
             for ( a = BDY(b), i = 0; i < len; i++, NMV_ADV(a) )              for ( a = BDY(b), i = 0; i < len; i++, NMV_ADV(a) )
                 wmax = MAX(TD(DL(a)),wmax);                  wmax = MAX((int)TD(DL(a)),wmax);
         }          }
         homogenize_order(ord,nvar,&ord1);          homogenize_order(ord,nvar,&ord1);
         nd_init_ord(ord1);          nd_init_ord(ord1);
         nd_setup_parameters(nvar+1,nd_nzlist?0:wmax);          nd_setup_parameters(nvar+1,nd_nzlist?0:wmax);
         for ( t = fd0; t; t = NEXT(t) )          for ( t = fd0; t; t = NEXT(t) )
             ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos,ompos);              ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos,ompos);
     }    }
   
     ndv_setup(m,0,fd0,(nd_gbblock||nd_splist||nd_check_splist)?1:0,0);    ndv_setup(m,0,fd0,(nd_gbblock||nd_splist||nd_check_splist)?1:0,0,0);
     if ( nd_gentrace ) {    if ( nd_gentrace ) {
         MKLIST(l1,nd_tracelist); MKNODE(nd_alltracelist,l1,0);        MKLIST(l1,nd_tracelist); MKNODE(nd_alltracelist,l1,0);
     }    }
   if ( nd_splist ) {    if ( nd_splist ) {
     *rp = compute_splist();      *rp = compute_splist();
     return;      return;
   }    }
   if ( nd_check_splist ) {    if ( nd_check_splist ) {
         if ( f4 ) {      if ( f4 ) {
             if ( check_splist_f4(m,nd_check_splist) ) *rp = (LIST)ONE;        if ( check_splist_f4(m,nd_check_splist) ) *rp = (LIST)ONE;
             else *rp = 0;        else *rp = 0;
         } else {      } else {
             if ( check_splist(m,nd_check_splist) ) *rp = (LIST)ONE;        if ( check_splist(m,nd_check_splist) ) *rp = (LIST)ONE;
             else *rp = 0;        else *rp = 0;
         }      }
     return;      return;
   }    }
     x = f4?nd_f4(m,0,&perm):nd_gb(m,ishomo || homo,0,0,&perm);    x = f4?nd_f4(m,0,&perm):nd_gb(m,ishomo || homo,0,0,&perm);
   if ( !x ) {    if ( !x ) {
     *rp = 0; return;      *rp = 0; return;
   }    }
   if ( nd_gentrace ) {    if ( nd_gentrace ) {
     MKVECT(hvect,nd_psn);      MKVECT(hvect,nd_psn);
     for ( i = 0; i < nd_psn; i++ )      for ( i = 0; i < nd_psn; i++ )
        ndltodp(nd_psh[i]->dl,(DP *)&BDY(hvect)[i]);        ndltodp(nd_psh[i]->dl,(DP *)&BDY(hvect)[i]);
   }    }
   if ( !ishomo && homo ) {    if ( !ishomo && homo ) {
        /* dehomogenization */         /* dehomogenization */
Line 3389  void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int 
Line 4983  void nd_gr(LIST f,LIST v,int m,int homo,int retdp,int 
     nd_init_ord(ord);      nd_init_ord(ord);
     nd_setup_parameters(nvar,0);      nd_setup_parameters(nvar,0);
   }    }
     nd_demand = 0;    nd_demand = 0;
   if ( nd_module && nd_intersect ) {    if ( nd_module && nd_intersect ) {
     for ( j = nd_psn-1, x = 0; j >= 0; j-- )      for ( j = nd_psn-1, x = 0; j >= 0; j-- )
       if ( MPOS(DL(nd_psh[j])) > nd_intersect ) {        if ( MPOS(DL(nd_psh[j])) > (UINT)nd_intersect ) {
         MKNODE(xx,(pointer)((unsigned long)j),x); x = xx;          MKNODE(xx,(pointer)((unsigned long)j),x); x = xx;
       }        }
     conv_ilist(nd_demand,0,x,0);      conv_ilist(nd_demand,0,x,0);
     goto FINAL;      goto FINAL;
   }    }
     if ( nd_gentrace  && f4 ) { nzlist = nd_alltracelist; }    if ( nd_gentrace  && f4 ) { nzlist = nd_alltracelist; }
     x = ndv_reducebase(x,perm);    x = ndv_reducebase(x,perm);
     if ( nd_gentrace  && !f4 ) { tl1 = nd_alltracelist; nd_alltracelist = 0; }    if ( nd_gentrace  && !f4 ) { tl1 = nd_alltracelist; nd_alltracelist = 0; }
     x = ndv_reduceall(m,x);    x = ndv_reduceall(m,x);
     cbpe = nd_bpe;    cbpe = nd_bpe;
     if ( nd_gentrace && !f4 ) {    if ( nd_gentrace && !f4 ) {
         tl2 = nd_alltracelist; nd_alltracelist = 0;      tl2 = nd_alltracelist; nd_alltracelist = 0;
         ndv_check_membership(m,fd0,obpe,oadv,oepos,x);      ndv_check_membership(m,fd0,obpe,oadv,oepos,x);
         tl3 = nd_alltracelist; nd_alltracelist = 0;      tl3 = nd_alltracelist; nd_alltracelist = 0;
         if ( nd_gensyz ) {      if ( nd_gensyz ) {
           nd_gb(m,0,1,1,0);        nd_gb(m,0,1,1,0);
             tl4 = nd_alltracelist; nd_alltracelist = 0;        tl4 = nd_alltracelist; nd_alltracelist = 0;
         } else tl4 = 0;      } else tl4 = 0;
     }    }
     nd_bpe = cbpe;    nd_bpe = cbpe;
     nd_setup_parameters(nd_nvar,0);    nd_setup_parameters(nd_nvar,0);
 FINAL:  FINAL:
     for ( r0 = 0, t = x; t; t = NEXT(t) ) {    for ( r0 = 0, t = x; t; t = NEXT(t) ) {
       NEXTNODE(r0,r);      NEXTNODE(r0,r);
       if ( nd_module ) {      if ( nd_module ) {
         if ( retdp ) BDY(r) = ndvtodpm(m,BDY(t));        if ( retdp ) BDY(r) = ndvtodpm(m,BDY(t));
         else BDY(r) = ndvtopl(m,CO,vv,BDY(t),mrank);        else BDY(r) = ndvtopl(m,CO,vv,BDY(t),mrank);
       } else if ( retdp ) BDY(r) = ndvtodp(m,BDY(t));      } else if ( retdp ) BDY(r) = ndvtodp(m,BDY(t));
       else BDY(r) = ndvtop(m,CO,vv,BDY(t));      else BDY(r) = ndvtop(m,CO,vv,BDY(t));
     }    }
     if ( r0 ) NEXT(r) = 0;    if ( r0 ) NEXT(r) = 0;
     if ( !m && nd_nalg )    if ( !m && nd_nalg )
         r0 = postprocess_algcoef(av,alist,r0);          r0 = postprocess_algcoef(av,alist,r0);
     MKLIST(*rp,r0);    MKLIST(*rp,r0);
     if ( nd_gentrace ) {    if ( nd_gentrace ) {
   if ( f4 ) {    if ( f4 ) {
             STOZ(16,bpe);              STOZ(16,bpe);
             STOZ(nd_last_nonzero,last_nonzero);              STOZ(nd_last_nonzero,last_nonzero);
Line 3459  FINAL:
Line 5053  FINAL:
 #endif  #endif
 }  }
   
   NODE nd_sba_f4(int m,int **indp);
   
   void nd_sba(LIST f,LIST v,int m,int homo,int retdp,int f4,struct order_spec *ord,LIST *rp)
   {
     VL tv,fv,vv;
     NODE fd,fd0,r,r0,t,x,nd,nd1,syz;
     int e,max,nvar,i;
     NDV b;
     int ishomo,wmax,len;
     NMV a;
     P zp;
     Q dmy;
     struct order_spec *ord1;
     int *perm;
     EPOS oepos;
     int obpe,oadv,ompos;
     struct oEGT eg0,eg1,egconv,egintred;
     LIST l1,redind;
     Z z;
   
     nd_module = 0;
     nd_demand = 0;
     Nsamesig = 0;
     if ( DP_Multiple )
       nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);
     get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
     parse_nd_option(vv,current_option);
     if ( m && nd_vc )
       error("nd_sba : computation over Fp(X) is unsupported. Use dp_gr_mod_main().");
     for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
     switch ( ord->id ) {
       case 1:
         if ( ord->nv != nvar )
           error("nd_sba : invalid order specification");
         break;
       default:
         break;
     }
     nd_nalg = 0;
     nd_init_ord(ord);
     // for SIG comparison
     initd(ord);
     for ( t = BDY(f), max = 1; t; t = NEXT(t) ) {
       for ( tv = vv; tv; tv = NEXT(tv) ) {
         if ( OID(BDY(t)) == O_DP ) {
           e = dp_getdeg((DP)BDY(t));
           max = MAX(e,max);
         } else {
           e = getdeg(tv->v,(P)BDY(t));
           max = MAX(e,max);
         }
       }
     }
     nd_setup_parameters(nvar,max);
     obpe = nd_bpe; oadv = nmv_adv; oepos = nd_epos; ompos = nd_mpos;
     ishomo = 1;
     for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) {
       if ( OID(BDY(t)) == O_DP ) {
         DP zdp;
   
         if ( !m ) dp_ptozp((DP)BDY(t),&zdp);
         else zdp = (DP)BDY(t);
         b = (pointer)dptondv(m,zdp);
       } else {
         if ( !m ) ptozp((P)BDY(t),1,&dmy,&zp);
         else zp = (P)BDY(t);
         b = (pointer)ptondv(CO,vv,zp);
       }
       if ( ishomo )
         ishomo = ishomo && ndv_ishomo(b);
       if ( m ) ndv_mod(m,b);
       if ( b ) { NEXTNODE(fd0,fd); BDY(fd) = (pointer)b; }
     }
     if ( fd0 ) NEXT(fd) = 0;
   
     if ( !ishomo && homo ) {
       for ( t = fd0, wmax = max; t; t = NEXT(t) ) {
         b = (NDV)BDY(t); len = LEN(b);
         for ( a = BDY(b), i = 0; i < len; i++, NMV_ADV(a) )
           wmax = MAX((int)TD(DL(a)),wmax);
         }
         homogenize_order(ord,nvar,&ord1);
         nd_init_ord(ord1);
         // for SIG comparison
         initd(ord1);
         nd_setup_parameters(nvar+1,nd_nzlist?0:wmax);
         for ( t = fd0; t; t = NEXT(t) )
           ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos,ompos);
     }
   
     ndv_setup(m,0,fd0,nd_sba_dontsort,0,1);
     if ( nd_gentrace ) {
       MKLIST(l1,nd_tracelist); MKNODE(nd_alltracelist,l1,0);
     }
     x = f4 ? nd_sba_f4(m,&perm) : nd_sba_buch(m,ishomo || homo,&perm,&syz);
     if ( !x ) {
       *rp = 0; return;
     }
     if ( !ishomo && homo ) {
          /* dehomogenization */
       for ( t = x; t; t = NEXT(t) ) ndv_dehomogenize((NDV)BDY(t),ord);
       nd_init_ord(ord);
       // for SIG comparison
       initd(ord);
       nd_setup_parameters(nvar,0);
     }
     nd_demand = 0;
     get_eg(&eg0);
     x = ndv_reducebase(x,perm);
     for ( nd = 0, i = length(x)-1; i >= 0; i-- ) {
       STOZ(perm[i],z); MKNODE(nd1,z,nd); nd = nd1;
     }
     MKLIST(redind,nd);
     x = ndv_reduceall(m,x);
     get_eg(&eg1); init_eg(&egintred); add_eg(&egintred,&eg0,&eg1);
     nd_setup_parameters(nd_nvar,0);
     get_eg(&eg0);
     for ( r0 = 0, t = x; t; t = NEXT(t) ) {
       NEXTNODE(r0,r);
       if ( retdp ) BDY(r) = ndvtodp(m,BDY(t));
       else BDY(r) = ndvtop(m,CO,vv,BDY(t));
     }
     if ( r0 ) NEXT(r) = 0;
     if ( nd_sba_syz ) {
       LIST gb,hsyz;
   
       MKLIST(gb,r0);
       MKLIST(hsyz,syz);
       nd = mknode(2,gb,hsyz);
       MKLIST(*rp,nd);
     } else if ( nd_gentrace ) {
       LIST gb,trace;
   
       MKLIST(trace,nd_alltracelist);
       MKLIST(gb,r0);
       nd = mknode(3,gb,redind,trace);
       MKLIST(*rp,nd);
     } else
       MKLIST(*rp,r0);
     get_eg(&eg1); init_eg(&egconv); add_eg(&egconv,&eg0,&eg1);
     print_eg("intred",&egintred); fprintf(asir_out,"\n");
     print_eg("conv",&egconv); fprintf(asir_out,"\n");
   }
   
 void nd_gr_postproc(LIST f,LIST v,int m,struct order_spec *ord,int do_check,LIST *rp)  void nd_gr_postproc(LIST f,LIST v,int m,struct order_spec *ord,int do_check,LIST *rp)
 {  {
     VL tv,fv,vv,vc,av;      VL tv,fv,vv,av;
     NODE fd,fd0,r,r0,t,x,s,xx,alist;      NODE fd,fd0,r,r0,t,x,alist;
     int e,max,nvar,i;      int e,max,nvar,i,nalg;
     NDV b;      NDV b;
     int ishomo,nalg;      int ishomo;
     Alg alpha,dp;      P zp;
     P p,zp;  
     Q dmy;      Q dmy;
     LIST f1,f2;      LIST f1;
     Obj obj;  
     NumberField nf;  
     struct order_spec *ord1;      struct order_spec *ord1;
     int *perm;      int *perm;
   
     parse_nd_option(current_option);  
     get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);      get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
       parse_nd_option(vv,current_option);
     for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );      for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
     switch ( ord->id ) {      switch ( ord->id ) {
         case 1:          case 1:
Line 3518  void nd_gr_postproc(LIST f,LIST v,int m,struct order_s
Line 5253  void nd_gr_postproc(LIST f,LIST v,int m,struct order_s
         if ( b ) { NEXTNODE(fd0,fd); BDY(fd) = (pointer)b; }          if ( b ) { NEXTNODE(fd0,fd); BDY(fd) = (pointer)b; }
     }      }
     if ( fd0 ) NEXT(fd) = 0;      if ( fd0 ) NEXT(fd) = 0;
     ndv_setup(m,0,fd0,0,1);      ndv_setup(m,0,fd0,0,1,0);
     for ( x = 0, i = 0; i < nd_psn; i++ )      for ( x = 0, i = 0; i < nd_psn; i++ )
         x = update_base(x,i);          x = update_base(x,i);
     if ( do_check ) {      if ( do_check ) {
Line 3607  NDV recompute_trace(NODE ti,NDV *p,int mod)
Line 5342  NDV recompute_trace(NODE ti,NDV *p,int mod)
   
 void nd_gr_recompute_trace(LIST f,LIST v,int m,struct order_spec *ord,LIST tlist,LIST *rp)  void nd_gr_recompute_trace(LIST f,LIST v,int m,struct order_spec *ord,LIST tlist,LIST *rp)
 {  {
     VL tv,fv,vv,vc,av;    VL tv,fv,vv;
     NODE fd,fd0,r,r0,t,x,s,xx,alist;    NODE r,r0,t;
     int e,max,nvar,i;    int nvar,i;
     NDV b;    NDV b;
     int ishomo,nalg;    P zp;
     Alg alpha,dp;    Q dmy;
     P p,zp;  
     Q dmy;  
     LIST f1,f2;  
     Obj obj;  
     NumberField nf;  
     struct order_spec *ord1;  
   NODE permtrace,intred,ind,perm,trace,ti;    NODE permtrace,intred,ind,perm,trace,ti;
   int len,n,j;    int len,n,j;
   NDV *db,*pb;    NDV *db,*pb;
   
     parse_nd_option(current_option);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
     get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);    parse_nd_option(vv,current_option);
     for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );    for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
     switch ( ord->id ) {    switch ( ord->id ) {
         case 1:      case 1:
             if ( ord->nv != nvar )        if ( ord->nv != nvar )
                 error("nd_check : invalid order specification");          error("nd_check : invalid order specification");
             break;        break;
         default:      default:
             break;        break;
     }    }
     nd_init_ord(ord);    nd_init_ord(ord);
   nd_bpe = ZTOS((Q)ARG7(BDY(tlist)));    nd_bpe = (int)ZTOS((Q)ARG7(BDY(tlist)));
     nd_setup_parameters(nvar,0);    nd_setup_parameters(nvar,0);
   
   len = length(BDY(f));    len = length(BDY(f));
   db = (NDV *)MALLOC(len*sizeof(NDV *));    db = (NDV *)MALLOC(len*sizeof(NDV *));
   for ( i = 0, t = BDY(f); t; i++, t = NEXT(t) ) {    for ( i = 0, t = BDY(f); t; i++, t = NEXT(t) ) {
       ptozp((P)BDY(t),1,&dmy,&zp);      ptozp((P)BDY(t),1,&dmy,&zp);
       b = ptondv(CO,vv,zp);      b = ptondv(CO,vv,zp);
         ndv_mod(m,b);      ndv_mod(m,b);
     ndv_mul_c(m,b,invm(HCM(b),m));      ndv_mul_c(m,b,invm(HCM(b),m));
     db[i] = b;      db[i] = b;
     }    }
   
   permtrace = BDY((LIST)ARG2(BDY(tlist)));    permtrace = BDY((LIST)ARG2(BDY(tlist)));
   intred = BDY((LIST)ARG3(BDY(tlist)));    intred = BDY((LIST)ARG3(BDY(tlist)));
Line 3655  void nd_gr_recompute_trace(LIST f,LIST v,int m,struct 
Line 5384  void nd_gr_recompute_trace(LIST f,LIST v,int m,struct 
   trace = NEXT(permtrace);    trace = NEXT(permtrace);
   
   for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {    for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {
     j = ZTOS((Q)ARG0(BDY((LIST)BDY(t))));      j = (int)ZTOS((Q)ARG0(BDY((LIST)BDY(t))));
     if ( j > i ) i = j;      if ( j > i ) i = j;
   }    }
   n = i+1;    n = i+1;
Line 3677  void nd_gr_recompute_trace(LIST f,LIST v,int m,struct 
Line 5406  void nd_gr_recompute_trace(LIST f,LIST v,int m,struct 
     pb[ZTOS((Q)ARG0(ti))] = recompute_trace(BDY((LIST)ARG1(ti)),pb,m);      pb[ZTOS((Q)ARG0(ti))] = recompute_trace(BDY((LIST)ARG1(ti)),pb,m);
     if ( !pb[ZTOS((Q)ARG0(ti))] ) { *rp = 0; return; }      if ( !pb[ZTOS((Q)ARG0(ti))] ) { *rp = 0; return; }
       if ( DP_Print ) {        if ( DP_Print ) {
            fprintf(asir_out,"*"); fflush(asir_out);          fprintf(asir_out,"*"); fflush(asir_out);
       }      }
   }    }
     for ( r0 = 0, t = ind; t; t = NEXT(t) ) {    for ( r0 = 0, t = ind; t; t = NEXT(t) ) {
         NEXTNODE(r0,r);      NEXTNODE(r0,r);
     b = pb[ZTOS((Q)BDY(t))];      b = pb[ZTOS((Q)BDY(t))];
         ndv_mul_c(m,b,invm(HCM(b),m));      ndv_mul_c(m,b,invm(HCM(b),m));
 #if 0  #if 0
         BDY(r) = ndvtop(m,CO,vv,pb[ZTOS((Q)BDY(t))]);      BDY(r) = ndvtop(m,CO,vv,pb[ZTOS((Q)BDY(t))]);
 #else  #else
         BDY(r) = ndvtodp(m,pb[ZTOS((Q)BDY(t))]);      BDY(r) = ndvtodp(m,pb[ZTOS((Q)BDY(t))]);
 #endif  #endif
     }    }
     if ( r0 ) NEXT(r) = 0;    if ( r0 ) NEXT(r) = 0;
     MKLIST(*rp,r0);    MKLIST(*rp,r0);
     if ( DP_Print ) fprintf(asir_out,"\n");    if ( DP_Print ) fprintf(asir_out,"\n");
 }  }
   
 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int retdp,int f4,struct order_spec *ord,LIST *rp)  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int retdp,int f4,struct order_spec *ord,LIST *rp)
 {  {
     VL tv,fv,vv,vc,av;      VL tv,fv,vv,av;
     NODE fd,fd0,in0,in,r,r0,t,s,cand,alist;      NODE fd,fd0,in0,in,r,t,s,cand,alist;
     int m,nocheck,nvar,mindex,e,max;      int m,nocheck,nvar,mindex,e,max;
     NDV c;      NDV c;
     NMV a;      NMV a;
     P p,zp;      P zp;
     Q dmy;      Q dmy;
     EPOS oepos;      EPOS oepos;
     int obpe,oadv,wmax,i,len,cbpe,ishomo,nalg,mrank,trank,ompos;      int obpe,oadv,wmax,i,len,cbpe,ishomo,nalg,mrank,trank,ompos;
     Alg alpha,dp;      LIST f1,zpl;
     P poly;  
     LIST f1,f2,zpl;  
     Obj obj;  
     NumberField nf;  
     struct order_spec *ord1;      struct order_spec *ord1;
     struct oEGT eg_check,eg0,eg1;      struct oEGT eg_check,eg_gb,eg_intred,eg0,eg1;
     NODE tr,tl1,tl2,tl3,tl4;      NODE tr,tl1,tl2,tl3,tl4;
     LIST l1,l2,l3,l4,l5;      LIST l1,l2,l3,l4,l5;
     int *perm;      int *perm;
Line 3720  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5445  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
     Z jq,bpe;      Z jq,bpe;
     VECT hvect;      VECT hvect;
   
       NcriB = NcriMF = Ncri2 = 0;
     nd_module = 0;      nd_module = 0;
     nd_lf = 0;      nd_lf = 0;
     parse_nd_option(current_option);      get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
       parse_nd_option(vv,current_option);
     if ( nd_lf ) {      if ( nd_lf ) {
       if ( f4 )        if ( f4 )
         nd_f4_lf_trace(f,v,trace,homo,ord,rp);          nd_f4_lf_trace(f,v,trace,homo,ord,rp);
Line 3733  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5460  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
     if ( DP_Multiple )      if ( DP_Multiple )
         nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);          nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);
   
     get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);  
     for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );      for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
     switch ( ord->id ) {      switch ( ord->id ) {
         case 1:          case 1:
Line 3787  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5513  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
                 }                  }
               }                }
             } else {              } else {
                 if ( OID(BDY(t)) == O_DP ) {
                   e = dp_getdeg((DP)BDY(t));
                   max = MAX(e,max);
                 } else {
                 e = getdeg(tv->v,(P)BDY(t));                  e = getdeg(tv->v,(P)BDY(t));
                 max = MAX(e,max);                  max = MAX(e,max);
                 }
             }              }
         }          }
     nd_setup_parameters(nvar,max);      nd_setup_parameters(nvar,max);
Line 3809  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5540  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
             c = (pointer)pltondv(CO,vv,zpl);              c = (pointer)pltondv(CO,vv,zpl);
           }            }
         } else {          } else {
           if ( !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);            if ( OID(BDY(t)) == O_DP ) {
           else zp = (P)BDY(t);              DP zdp;
           c = (pointer)ptondv(CO,vv,zp);  
               if ( !nd_gentrace ) dp_ptozp((DP)BDY(t),&zdp);
               else zdp = (DP)BDY(t);
               c = (pointer)dptondv(m,zdp);
             } else {
               if ( !nd_gentrace ) ptozp((P)BDY(t),1,&dmy,&zp);
               else zp = (P)BDY(t);
               c = (pointer)ptondv(CO,vv,zp);
             }
         }          }
         if ( ishomo )          if ( ishomo )
             ishomo = ishomo && ndv_ishomo(c);              ishomo = ishomo && ndv_ishomo(c);
Line 3826  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5565  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
         for ( t = in0, wmax = max; t; t = NEXT(t) ) {          for ( t = in0, wmax = max; t; t = NEXT(t) ) {
             c = (NDV)BDY(t); len = LEN(c);              c = (NDV)BDY(t); len = LEN(c);
             for ( a = BDY(c), i = 0; i < len; i++, NMV_ADV(a) )              for ( a = BDY(c), i = 0; i < len; i++, NMV_ADV(a) )
                 wmax = MAX(TD(DL(a)),wmax);                  wmax = MAX((int)TD(DL(a)),wmax);
         }          }
         homogenize_order(ord,nvar,&ord1);          homogenize_order(ord,nvar,&ord1);
         nd_init_ord(ord1);          nd_init_ord(ord1);
Line 3836  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5575  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
     }      }
     if ( MaxDeg > 0 ) nocheck = 1;      if ( MaxDeg > 0 ) nocheck = 1;
     while ( 1 ) {      while ( 1 ) {
     tl1 = tl2 = tl3 = tl4 = 0;          get_eg(&eg0);
           tl1 = tl2 = tl3 = tl4 = 0;
         if ( Demand )          if ( Demand )
             nd_demand = 1;              nd_demand = 1;
         ret = ndv_setup(m,1,fd0,nd_gbblock?1:0,0);          ret = ndv_setup(m,1,fd0,nd_gbblock?1:0,0,0);
         if ( nd_gentrace ) {          if ( nd_gentrace ) {
             MKLIST(l1,nd_tracelist); MKNODE(nd_alltracelist,l1,0);              MKLIST(l1,nd_tracelist); MKNODE(nd_alltracelist,l1,0);
         }          }
Line 3862  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5602  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
             nd_init_ord(ord);              nd_init_ord(ord);
             nd_setup_parameters(nvar,0);              nd_setup_parameters(nvar,0);
         }          }
           get_eg(&eg1); init_eg(&eg_gb); add_eg(&eg_gb,&eg0,&eg1);
         nd_demand = 0;          nd_demand = 0;
         cand = ndv_reducebase(cand,perm);          cand = ndv_reducebase(cand,perm);
         if ( nd_gentrace ) { tl1 = nd_alltracelist; nd_alltracelist = 0; }          if ( nd_gentrace ) { tl1 = nd_alltracelist; nd_alltracelist = 0; }
           get_eg(&eg0);
         cand = ndv_reduceall(0,cand);          cand = ndv_reduceall(0,cand);
           get_eg(&eg1); init_eg(&eg_intred); add_eg(&eg_intred,&eg0,&eg1);
         cbpe = nd_bpe;          cbpe = nd_bpe;
         if ( nd_gentrace ) { tl2 = nd_alltracelist; nd_alltracelist = 0; }          if ( nd_gentrace ) { tl2 = nd_alltracelist; nd_alltracelist = 0; }
         get_eg(&eg0);          get_eg(&eg0);
Line 3876  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5619  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
           tl3 = nd_alltracelist; nd_alltracelist = 0;            tl3 = nd_alltracelist; nd_alltracelist = 0;
         } else tl3 = 0;          } else tl3 = 0;
             /* gbcheck : cand is a GB of Id(cand) ? */              /* gbcheck : cand is a GB of Id(cand) ? */
             if ( nd_vc || nd_gentrace || nd_gensyz )              if ( nd_vc || nd_gentrace || nd_gensyz || do_weyl )
               ret = nd_gb(0,0,1,nd_gensyz?1:0,0)!=0;                ret = nd_gb(0,0,1,nd_gensyz?1:0,0)!=0;
             else              else
               ret = nd_f4(0,1,0)!=0;                ret = nd_f4(0,1,0)!=0;
Line 3903  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 5646  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
     }      }
     get_eg(&eg1); init_eg(&eg_check); add_eg(&eg_check,&eg0,&eg1);      get_eg(&eg1); init_eg(&eg_check); add_eg(&eg_check,&eg0,&eg1);
     if ( DP_Print )      if ( DP_Print )
         fprintf(asir_out,"check=%.3fsec\n",eg_check.exectime);          fprintf(asir_out,"gb=%.3fsec,check=%.3fsec,intred=%.3fsec\n",
             eg_gb.exectime,eg_check.exectime,eg_intred.exectime);
     /* dp->p */      /* dp->p */
     nd_bpe = cbpe;      nd_bpe = cbpe;
     nd_setup_parameters(nd_nvar,0);      nd_setup_parameters(nd_nvar,0);
Line 3972  DL ndltodl(int n,UINT *ndl)
Line 5716  DL ndltodl(int n,UINT *ndl)
     int i,j,l,s,ord_l;      int i,j,l,s,ord_l;
     struct order_pair *op;      struct order_pair *op;
   
     NEWDL(dl,n);      NEWDL_NOINIT(dl,n);
     dl->td = TD(ndl);      dl->td = TD(ndl);
     d = dl->d;      d = dl->d;
     if ( nd_blockmask ) {      if ( nd_blockmask ) {
Line 3988  DL ndltodl(int n,UINT *ndl)
Line 5732  DL ndltodl(int n,UINT *ndl)
     return dl;      return dl;
 }  }
   
   void _ndltodl(UINT *ndl,DL dl)
   {
       int *d;
       int i,j,l,s,ord_l,n;
       struct order_pair *op;
   
       n = nd_nvar;
       dl->td = TD(ndl);
       d = dl->d;
       if ( nd_blockmask ) {
           l = nd_blockmask->n;
           op = nd_blockmask->order_pair;
           for ( j = 0, s = 0; j < l; j++ ) {
               ord_l = op[j].length;
               for ( i = 0; i < ord_l; i++, s++ ) d[s] = GET_EXP(ndl,s);
           }
       } else {
           for ( i = 0; i < n; i++ ) d[i] = GET_EXP(ndl,i);
       }
   }
   
 void nmtodp(int mod,NM m,DP *r)  void nmtodp(int mod,NM m,DP *r)
 {  {
     DP dp;      DP dp;
Line 4045  void nd_print(ND p)
Line 5810  void nd_print(ND p)
     else {      else {
         for ( m = BDY(p); m; m = NEXT(m) ) {          for ( m = BDY(p); m; m = NEXT(m) ) {
             if ( CM(m) & 0x80000000 ) printf("+@_%d*",IFTOF(CM(m)));              if ( CM(m) & 0x80000000 ) printf("+@_%d*",IFTOF(CM(m)));
             else printf("+%d*",CM(m));              else printf("+%ld*",CM(m));
             ndl_print(DL(m));              ndl_print(DL(m));
         }          }
         printf("\n");          printf("\n");
Line 4132  void ndv_removecont(int mod,NDV p)
Line 5897  void ndv_removecont(int mod,NDV p)
     int i,len,all_p;      int i,len,all_p;
     Z *c;      Z *c;
     P *w;      P *w;
     Z dvr,t;      Z dvr;
     P g,cont,tp;      P g,cont,tp;
     NMV m;      NMV m;
   
Line 4177  void ndv_homogenize(NDV p,int obpe,int oadv,EPOS oepos
Line 5942  void ndv_homogenize(NDV p,int obpe,int oadv,EPOS oepos
   
     len = p->len;      len = p->len;
     for ( m = BDY(p), i = 0, max = 0; i < len; NMV_OADV(m), i++ )      for ( m = BDY(p), i = 0, max = 0; i < len; NMV_OADV(m), i++ )
         max = MAX(max,TD(DL(m)));          max = MAX(max,(int)TD(DL(m)));
     mr0 = nmv_adv>oadv?(NMV)REALLOC(BDY(p),len*nmv_adv):BDY(p);      mr0 = nmv_adv>oadv?(NMV)REALLOC(BDY(p),len*nmv_adv):BDY(p);
     m = (NMV)((char *)mr0+(len-1)*oadv);      m = (NMV)((char *)mr0+(len-1)*oadv);
     mr = (NMV)((char *)mr0+(len-1)*nmv_adv);      mr = (NMV)((char *)mr0+(len-1)*nmv_adv);
Line 4195  void ndv_dehomogenize(NDV p,struct order_spec *ord)
Line 5960  void ndv_dehomogenize(NDV p,struct order_spec *ord)
 {  {
     int i,j,adj,len,newnvar,newwpd,newadv,newexporigin,newmpos;      int i,j,adj,len,newnvar,newwpd,newadv,newexporigin,newmpos;
     int pos;      int pos;
     Q *w;  
     Q dvr,t;  
     NMV m,r;      NMV m,r;
   
     len = p->len;      len = p->len;
Line 4256  void nd_heu_nezgcdnpz(VL vl,P *pl,int m,int full,P *pr
Line 6019  void nd_heu_nezgcdnpz(VL vl,P *pl,int m,int full,P *pr
   
 void removecont_array(P *p,int n,int full)  void removecont_array(P *p,int n,int full)
 {  {
     int all_p,all_q,i;      int all_q,i;
     Z *c;      Z *c;
     P *w;      P *w;
     P t,s;      P t,s;
Line 4285  void removecont_array_q(Z *c,int n)
Line 6048  void removecont_array_q(Z *c,int n)
 {  {
   struct oVECT v;    struct oVECT v;
   Z d0,d1,a,u,u1,gcd;    Z d0,d1,a,u,u1,gcd;
   int i,j;    int i;
   Z *q,*r;    Z *q,*r;
   
   q = (Z *)MALLOC(n*sizeof(Z));    q = (Z *)MALLOC(n*sizeof(Z));
Line 4429  void ndv_free(NDV p)
Line 6192  void ndv_free(NDV p)
   
 void nd_append_red(UINT *d,int i)  void nd_append_red(UINT *d,int i)
 {  {
     RHist m,m0;      RHist m;
     int h;      int h;
   
     NEWRHist(m);      NEWRHist(m);
Line 4473  UINT *nd_compute_bound(ND p)
Line 6236  UINT *nd_compute_bound(ND p)
 {  {
     UINT *d1,*d2,*t;      UINT *d1,*d2,*t;
     UINT u;      UINT u;
     int i,j,k,l,len,ind;      int i,j,k,l,ind;
     NM m;      NM m;
   
     if ( !p )      if ( !p )
         return 0;          return 0;
     d1 = (UINT *)MALLOC(nd_wpd*sizeof(UINT));      d1 = (UINT *)MALLOC(nd_wpd*sizeof(UINT));
     d2 = (UINT *)MALLOC(nd_wpd*sizeof(UINT));      d2 = (UINT *)MALLOC(nd_wpd*sizeof(UINT));
     len = LEN(p);  
     m = BDY(p); ndl_copy(DL(m),d1); m = NEXT(m);      m = BDY(p); ndl_copy(DL(m),d1); m = NEXT(m);
     for ( m = NEXT(m); m; m = NEXT(m) ) {      for ( m = NEXT(m); m; m = NEXT(m) ) {
         ndl_lcm(DL(m),d1,d2);          ndl_lcm(DL(m),d1,d2);
Line 4524  int nd_get_exporigin(struct order_spec *ord)
Line 6286  int nd_get_exporigin(struct order_spec *ord)
 }  }
   
 void nd_setup_parameters(int nvar,int max) {  void nd_setup_parameters(int nvar,int max) {
     int i,j,n,elen,ord_o,ord_l,l,s,wpd;      int i,j,elen,wpd;
     struct order_pair *op;      extern int CNVars;
   
     nd_nvar = nvar;      nd_nvar = nvar;
       CNVars = nvar;
     if ( max ) {      if ( max ) {
         /* XXX */          /* XXX */
         if ( do_weyl ) nd_bpe = 32;          if ( do_weyl ) nd_bpe = 32;
Line 4544  void nd_setup_parameters(int nvar,int max) {
Line 6307  void nd_setup_parameters(int nvar,int max) {
     }      }
     if ( !do_weyl && weight_check && (current_dl_weight_vector || nd_matrix) ) {      if ( !do_weyl && weight_check && (current_dl_weight_vector || nd_matrix) ) {
         UINT t;          UINT t;
     int st;          int st;
         int *v;          int *v;
   /* t = max(weights) */    /* t = max(weights) */
         t = 0;          t = 0;
         if ( current_dl_weight_vector )          if ( current_dl_weight_vector )
             for ( i = 0, t = 0; i < nd_nvar; i++ ) {              for ( i = 0, t = 0; i < nd_nvar; i++ ) {
                 if ( (st=current_dl_weight_vector[i]) < 0 ) st = -st;                  if ( (st=current_dl_weight_vector[i]) < 0 ) st = -st;
                 if ( t < st ) t = st;                  if ( t < (UINT)st ) t = st;
             }              }
         if ( nd_matrix )          if ( nd_matrix )
             for ( i = 0; i < nd_matrix_len; i++ )              for ( i = 0; i < nd_matrix_len; i++ )
                 for ( j = 0, v = nd_matrix[i]; j < nd_nvar; j++ ) {                  for ( j = 0, v = nd_matrix[i]; j < nd_nvar; j++ ) {
                     if ( (st=v[j]) < 0 ) st = -st;                      if ( (st=v[j]) < 0 ) st = -st;
                     if ( t < st ) t = st;                      if ( t < (UINT)st ) t = st;
         }          }
         /* i = bitsize of t */          /* i = bitsize of t */
         for ( i = 0; t; t >>=1, i++ );          for ( i = 0; t; t >>=1, i++ );
Line 4601  ND_pairs nd_reconstruct(int trace,ND_pairs d)
Line 6364  ND_pairs nd_reconstruct(int trace,ND_pairs d)
     int i,obpe,oadv,h;      int i,obpe,oadv,h;
     static NM prev_nm_free_list;      static NM prev_nm_free_list;
     static ND_pairs prev_ndp_free_list;      static ND_pairs prev_ndp_free_list;
     RHist mr0,mr;      RHist mr;
     RHist r;      RHist r;
     RHist *old_red;      RHist *old_red;
     ND_pairs s0,s,t;      ND_pairs s0,s,t;
Line 4640  ND_pairs nd_reconstruct(int trace,ND_pairs d)
Line 6403  ND_pairs nd_reconstruct(int trace,ND_pairs d)
         NEXTND_pairs(s0,s);          NEXTND_pairs(s0,s);
         s->i1 = t->i1;          s->i1 = t->i1;
         s->i2 = t->i2;          s->i2 = t->i2;
           s->sig = t->sig;
         SG(s) = SG(t);          SG(s) = SG(t);
         ndl_reconstruct(LCM(t),LCM(s),obpe,oepos);          ndl_reconstruct(LCM(t),LCM(s),obpe,oepos);
     }      }
Line 4658  ND_pairs nd_reconstruct(int trace,ND_pairs d)
Line 6422  ND_pairs nd_reconstruct(int trace,ND_pairs d)
             h = ndl_hash_value(DL(mr));              h = ndl_hash_value(DL(mr));
             NEXT(mr) = nd_red[h];              NEXT(mr) = nd_red[h];
             nd_red[h] = mr;              nd_red[h] = mr;
               mr->sig = r->sig;
         }          }
     for ( i = 0; i < REDTAB_LEN; i++ ) old_red[i] = 0;      for ( i = 0; i < REDTAB_LEN; i++ ) old_red[i] = 0;
     old_red = 0;      old_red = 0;
     for ( i = 0; i < nd_psn; i++ ) {      for ( i = 0; i < nd_psn; i++ ) {
         NEWRHist(r); SG(r) = SG(nd_psh[i]);          NEWRHist(r); SG(r) = SG(nd_psh[i]);
         ndl_reconstruct(DL(nd_psh[i]),DL(r),obpe,oepos);          ndl_reconstruct(DL(nd_psh[i]),DL(r),obpe,oepos);
           r->sig = nd_psh[i]->sig;
         nd_psh[i] = r;          nd_psh[i] = r;
     }      }
     if ( s0 ) NEXT(s) = 0;      if ( s0 ) NEXT(s) = 0;
Line 4675  ND_pairs nd_reconstruct(int trace,ND_pairs d)
Line 6441  ND_pairs nd_reconstruct(int trace,ND_pairs d)
     return s0;      return s0;
 }  }
   
   void nd_reconstruct_s(int trace,ND_pairs *d)
   {
       int i,obpe,oadv,h;
       static NM prev_nm_free_list;
       static ND_pairs prev_ndp_free_list;
       RHist mr;
       RHist r;
       RHist *old_red;
       ND_pairs s0,s,t;
       EPOS oepos;
   
       obpe = nd_bpe;
       oadv = nmv_adv;
       oepos = nd_epos;
       if ( obpe < 2 ) nd_bpe = 2;
       else if ( obpe < 3 ) nd_bpe = 3;
       else if ( obpe < 4 ) nd_bpe = 4;
       else if ( obpe < 5 ) nd_bpe = 5;
       else if ( obpe < 6 ) nd_bpe = 6;
       else if ( obpe < 8 ) nd_bpe = 8;
       else if ( obpe < 10 ) nd_bpe = 10;
       else if ( obpe < 16 ) nd_bpe = 16;
       else if ( obpe < 32 ) nd_bpe = 32;
       else error("nd_reconstruct_s : exponent too large");
   
       nd_setup_parameters(nd_nvar,0);
       prev_nm_free_list = _nm_free_list;
       prev_ndp_free_list = _ndp_free_list;
       _nm_free_list = 0;
       _ndp_free_list = 0;
       for ( i = nd_psn-1; i >= 0; i-- ) {
           ndv_realloc(nd_ps[i],obpe,oadv,oepos);
           ndv_realloc(nd_ps_sym[i],obpe,oadv,oepos);
       }
       if ( trace )
           for ( i = nd_psn-1; i >= 0; i-- ) {
               ndv_realloc(nd_ps_trace[i],obpe,oadv,oepos);
               ndv_realloc(nd_ps_trace_sym[i],obpe,oadv,oepos);
           }
   
       for ( i = 0; i < nd_nbase; i++ ) {
         s0 = 0;
         for ( t = d[i]; t; t = NEXT(t) ) {
             NEXTND_pairs(s0,s);
             s->i1 = t->i1;
             s->i2 = t->i2;
             s->sig = t->sig;
             SG(s) = SG(t);
             ndl_reconstruct(LCM(t),LCM(s),obpe,oepos);
         }
         if ( s0 ) NEXT(s) = 0;
         d[i] = s0;
       }
   
       old_red = (RHist *)MALLOC(REDTAB_LEN*sizeof(RHist));
       for ( i = 0; i < REDTAB_LEN; i++ ) {
           old_red[i] = nd_red[i];
           nd_red[i] = 0;
       }
       for ( i = 0; i < REDTAB_LEN; i++ )
           for ( r = old_red[i]; r; r = NEXT(r) ) {
               NEWRHist(mr);
               mr->index = r->index;
               SG(mr) = SG(r);
               ndl_reconstruct(DL(r),DL(mr),obpe,oepos);
               h = ndl_hash_value(DL(mr));
               NEXT(mr) = nd_red[h];
               nd_red[h] = mr;
               mr->sig = r->sig;
           }
       for ( i = 0; i < REDTAB_LEN; i++ ) old_red[i] = 0;
       old_red = 0;
       for ( i = 0; i < nd_psn; i++ ) {
           NEWRHist(r); SG(r) = SG(nd_psh[i]);
           ndl_reconstruct(DL(nd_psh[i]),DL(r),obpe,oepos);
           r->sig = nd_psh[i]->sig;
           nd_psh[i] = r;
       }
       if ( s0 ) NEXT(s) = 0;
       prev_nm_free_list = 0;
       prev_ndp_free_list = 0;
   #if 0
       GC_gcollect();
   #endif
   }
   
 void ndl_reconstruct(UINT *d,UINT *r,int obpe,EPOS oepos)  void ndl_reconstruct(UINT *d,UINT *r,int obpe,EPOS oepos)
 {  {
     int n,i,ei,oepw,omask0,j,s,ord_l,l;      int n,i,ei,omask0,j,s,ord_l,l;
     struct order_pair *op;      struct order_pair *op;
   
     n = nd_nvar;      n = nd_nvar;
     oepw = (sizeof(UINT)*8)/obpe;  //    oepw = (sizeof(UINT)*8)/obpe;
     omask0 = (1<<obpe)-1;      omask0 = (1<<obpe)-1;
     TD(r) = TD(d);      TD(r) = TD(d);
     for ( i = nd_exporigin; i < nd_wpd; i++ ) r[i] = 0;      for ( i = nd_exporigin; i < nd_wpd; i++ ) r[i] = 0;
Line 4709  void ndl_reconstruct(UINT *d,UINT *r,int obpe,EPOS oep
Line 6561  void ndl_reconstruct(UINT *d,UINT *r,int obpe,EPOS oep
 ND nd_copy(ND p)  ND nd_copy(ND p)
 {  {
     NM m,mr,mr0;      NM m,mr,mr0;
     int c,n;  
     ND r;      ND r;
   
     if ( !p )      if ( !p )
Line 4736  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
Line 6587  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
     P gp,tp;      P gp,tp;
     Z g,t;      Z g,t;
     Z iq;      Z iq;
     int td;  
     LIST hist;      LIST hist;
     NODE node;      NODE node;
     DP d;      DP d;
Line 4783  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
Line 6633  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
         MKLIST(hist,node); MKNODE(node,hist,nd_tracelist);          MKLIST(hist,node); MKNODE(node,hist,nd_tracelist);
         nd_tracelist = node;          nd_tracelist = node;
     }      }
       if ( *rp )
         (*rp)->sig = p->sig;
     FREENM(m1); FREENM(m2);      FREENM(m1); FREENM(m2);
     return 1;      return 1;
 }  }
Line 4987  ND ndv_mul_nm_symbolic(NM m0,NDV p)
Line 6839  ND ndv_mul_nm_symbolic(NM m0,NDV p)
 {  {
     NM mr,mr0;      NM mr,mr0;
     NMV m;      NMV m;
     UINT *d,*dt,*dm;      UINT *d;
     int c,n,td,i,c1,c2,len;      int i,len;
     Q q;  
     ND r;      ND r;
   
     if ( !p ) return 0;      if ( !p ) return 0;
     else {      else {
         n = NV(p); m = BDY(p);          m = BDY(p);
         d = DL(m0);          d = DL(m0);
         len = LEN(p);          len = LEN(p);
         mr0 = 0;          mr0 = 0;
         td = TD(d);  
         c = CM(m0);  
         for ( i = 0; i < len; i++, NMV_ADV(m) ) {          for ( i = 0; i < len; i++, NMV_ADV(m) ) {
             NEXTNM(mr0,mr);              NEXTNM(mr0,mr);
             CM(mr) = 1;              CM(mr) = 1;
Line 5016  ND ndv_mul_nm(int mod,NM m0,NDV p)
Line 6865  ND ndv_mul_nm(int mod,NM m0,NDV p)
 {  {
     NM mr,mr0;      NM mr,mr0;
     NMV m;      NMV m;
     UINT *d,*dt,*dm;      UINT *d;
     int c,n,td,i,c1,c2,len;      int c,i,c1,c2,len;
     P q;      P q;
     ND r;      ND r;
   
Line 5029  ND ndv_mul_nm(int mod,NM m0,NDV p)
Line 6878  ND ndv_mul_nm(int mod,NM m0,NDV p)
         } else          } else
             return weyl_ndv_mul_nm(mod,m0,p);              return weyl_ndv_mul_nm(mod,m0,p);
     } else {      } else {
         n = NV(p); m = BDY(p);          m = BDY(p);
         d = DL(m0);          d = DL(m0);
         len = LEN(p);          len = LEN(p);
         mr0 = 0;          mr0 = 0;
         td = TD(d);  
         if ( mod == -1 ) {          if ( mod == -1 ) {
             c = CM(m0);              c = CM(m0);
             for ( i = 0; i < len; i++, NMV_ADV(m) ) {              for ( i = 0; i < len; i++, NMV_ADV(m) ) {
Line 5077  ND nd_quo(int mod,PGeoBucket bucket,NDV d)
Line 6925  ND nd_quo(int mod,PGeoBucket bucket,NDV d)
 {  {
     NM mq0,mq;      NM mq0,mq;
     NMV tm;      NMV tm;
     Q q;      int i,nv,c,c1,c2,hindex;
     int i,nv,sg,c,c1,c2,hindex;  
     ND p,t,r;      ND p,t,r;
   
     if ( bucket->m < 0 ) return 0;      if ( bucket->m < 0 ) return 0;
Line 5161  NDV ndv_dup_realloc(NDV p,int obpe,int oadv,EPOS oepos
Line 7008  NDV ndv_dup_realloc(NDV p,int obpe,int oadv,EPOS oepos
     }      }
     MKNDV(NV(p),mr0,len,r);      MKNDV(NV(p),mr0,len,r);
     SG(r) = SG(p);      SG(r) = SG(p);
       r->sig = p->sig;
     return r;      return r;
 }  }
   
Line 5224  ND nd_dup(ND p)
Line 7072  ND nd_dup(ND p)
 void ndv_mod(int mod,NDV p)  void ndv_mod(int mod,NDV p)
 {  {
     NMV t,d;      NMV t,d;
     int r,s,u;      int r;
     int i,len,dlen;      int i,len,dlen;
     P cp;      P cp;
     Z c;      Z c;
Line 5336  NDV pltondv(VL vl,VL dvl,LIST p)
Line 7184  NDV pltondv(VL vl,VL dvl,LIST p)
   
 ND ptond(VL vl,VL dvl,P p)  ND ptond(VL vl,VL dvl,P p)
 {  {
     int n,i,j,k,e;      int j,k,e,i;
     VL tvl;      VL tvl;
     V v;      V v;
     DCP dc;      DCP dc;
     DCP *w;      DCP *w;
     ND r,s,t,u;      ND r,t,s;
     P x;      P x;
     int c;  
     UINT *d;      UINT *d;
     NM m,m0;      NM m,m0;
   
Line 5352  ND ptond(VL vl,VL dvl,P p)
Line 7199  ND ptond(VL vl,VL dvl,P p)
     else if ( NUM(p) ) {      else if ( NUM(p) ) {
         NEWNM(m);          NEWNM(m);
         ndl_zero(DL(m));          ndl_zero(DL(m));
         if ( !INT((Q)p) )          if ( RATN(p) && !INT((Q)p) )
           error("ptond : input must be integer-coefficient");            error("ptond : input must be integer-coefficient");
         CZ(m) = (Z)p;          CZ(m) = (Z)p;
         NEXT(m) = 0;          NEXT(m) = 0;
Line 5375  ND ptond(VL vl,VL dvl,P p)
Line 7222  ND ptond(VL vl,VL dvl,P p)
         } else {          } else {
             NEWNM(m0); d = DL(m0);              NEWNM(m0); d = DL(m0);
             for ( j = k-1, s = 0; j >= 0; j-- ) {              for ( j = k-1, s = 0; j >= 0; j-- ) {
                 ndl_zero(d); e = ZTOS(DEG(w[j])); PUT_EXP(d,i,e);                  ndl_zero(d); e = (int)ZTOS(DEG(w[j])); PUT_EXP(d,i,e);
                 TD(d) = MUL_WEIGHT(e,i);                  TD(d) = MUL_WEIGHT(e,i);
                 if ( nd_blockmask) ndl_weight_mask(d);                  if ( nd_blockmask) ndl_weight_mask(d);
                 if ( nd_module ) MPOS(d) = 0;                  if ( nd_module ) MPOS(d) = 0;
Line 5435  LIST ndvtopl(int mod,VL vl,VL dvl,NDV p,int rank)
Line 7282  LIST ndvtopl(int mod,VL vl,VL dvl,NDV p,int rank)
     Z q;      Z q;
     P c;      P c;
     UINT *d;      UINT *d;
     P s,r,u,t,w;      P r,u,t,w;
     GFS gfs;      GFS gfs;
     P *a;      P *a;
     LIST l;      LIST l;
Line 5493  NDV ndtondv(int mod,ND p)
Line 7340  NDV ndtondv(int mod,ND p)
     }      }
     MKNDV(NV(p),m0,len,d);      MKNDV(NV(p),m0,len,d);
     SG(d) = SG(p);      SG(d) = SG(p);
       d->sig = p->sig;
     return d;      return d;
 }  }
   
Line 5538  void dpm_sort(DPM p,DPM *rp)
Line 7386  void dpm_sort(DPM p,DPM *rp)
   
 int dpm_comp(DPM *a,DPM *b)  int dpm_comp(DPM *a,DPM *b)
 {  {
   return compdpm(CO,*a,*b);    return -compdpm(CO,*a,*b);
 }  }
   
 NODE dpm_sort_list(NODE l)  NODE dpm_sort_list(NODE l)
Line 5565  int nmv_comp(NMV a,NMV b)
Line 7413  int nmv_comp(NMV a,NMV b)
   return -t;    return -t;
 }  }
   
   NDV dptondv(int mod,DP p)
   {
     NDV d;
     NMV m,m0;
     MP t;
     int i,len,n;
   
     if ( !p ) return 0;
     for ( t = BDY(p), len = 0; t; t = NEXT(t), len++ );
     n = p->nv;
     if ( mod > 0 || mod == -1 )
       m0 = m = (NMV)MALLOC_ATOMIC_IGNORE_OFF_PAGE(len*nmv_adv);
     else
       m0 = m = MALLOC(len*nmv_adv);
     for ( i = 0, t = BDY(p); i < len; i++, NMV_ADV(m), t = NEXT(t) ) {
       dltondl(n,t->dl,DL(m));
       TD(DL(m)) = ndl_weight(DL(m));
       CZ(m) = (Z)t->c;
     }
     qsort(m0,len,nmv_adv,(int (*)(const void *,const void *))nmv_comp);
     MKNDV(NV(p),m0,len,d);
     SG(d) = SG(p);
     return d;
   }
   
 NDV dpmtondv(int mod,DPM p)  NDV dpmtondv(int mod,DPM p)
 {  {
   NDV d;    NDV d;
Line 5616  ND ndvtond(int mod,NDV p)
Line 7489  ND ndvtond(int mod,NDV p)
     NEXT(m) = 0;      NEXT(m) = 0;
     MKND(NV(p),m0,len,d);      MKND(NV(p),m0,len,d);
     SG(d) = SG(p);      SG(d) = SG(p);
       d->sig = p->sig;
     return d;      return d;
 }  }
   
Line 5640  DP ndvtodp(int mod,NDV p)
Line 7514  DP ndvtodp(int mod,NDV p)
     return d;      return d;
 }  }
   
   DPM sigtodpm(SIG s)
   {
     DMM m;
     DPM d;
   
     NEWDMM(m);
     m->c = (Obj)ONE;
     m->dl = s->dl;
     m->pos = s->pos+1;
     m->next = 0;
     MKDPM(nd_nvar,m,d);
     SG(d) = s->dl->td;
     return d;
   }
   
 DPM ndvtodpm(int mod,NDV p)  DPM ndvtodpm(int mod,NDV p)
 {  {
   DMM m,m0;    DMM m,m0;
Line 5665  DPM ndvtodpm(int mod,NDV p)
Line 7554  DPM ndvtodpm(int mod,NDV p)
   
 DP ndtodp(int mod,ND p)  DP ndtodp(int mod,ND p)
 {  {
     MP m,m0;    MP m,m0;
   DP d;    DP d;
     NM t;    NM t;
     int i,len;  
   
     if ( !p ) return 0;    if ( !p ) return 0;
     m0 = 0;    m0 = 0;
     len = p->len;    for ( t = BDY(p); t; t = NEXT(t) ) {
     for ( t = BDY(p); t; t = NEXT(t) ) {      NEXTMP(m0,m);
         NEXTMP(m0,m);      m->dl = ndltodl(nd_nvar,DL(t));
       m->dl = ndltodl(nd_nvar,DL(t));      m->c = (Obj)ndctop(mod,t->c);
       m->c = (Obj)ndctop(mod,t->c);    }
     }    NEXT(m) = 0;
     NEXT(m) = 0;  
   MKDP(nd_nvar,m0,d);    MKDP(nd_nvar,m0,d);
     SG(d) = SG(p);    SG(d) = SG(p);
     return d;    return d;
 }  }
   
 void ndv_print(NDV p)  void ndv_print(NDV p)
Line 5694  void ndv_print(NDV p)
Line 7581  void ndv_print(NDV p)
         len = LEN(p);          len = LEN(p);
         for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) {          for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) {
             if ( CM(m) & 0x80000000 ) printf("+@_%d*",IFTOF(CM(m)));              if ( CM(m) & 0x80000000 ) printf("+@_%d*",IFTOF(CM(m)));
             else printf("+%d*",CM(m));              else printf("+%ld*",CM(m));
             ndl_print(DL(m));              ndl_print(DL(m));
         }          }
         printf("\n");          printf("\n");
Line 5725  NODE ndv_reducebase(NODE x,int *perm)
Line 7612  NODE ndv_reducebase(NODE x,int *perm)
     NDVI w;      NDVI w;
     NODE t,t0;      NODE t,t0;
   
       if ( nd_norb ) return x;
     len = length(x);      len = length(x);
     w = (NDVI)MALLOC(len*sizeof(struct oNDVI));      w = (NDVI)MALLOC(len*sizeof(struct oNDVI));
     for ( i = 0, t = x; i < len; i++, t = NEXT(t) ) {      for ( i = 0, t = x; i < len; i++, t = NEXT(t) ) {
Line 5980  void nd_nf_p(Obj f,LIST g,LIST v,int m,struct order_sp
Line 7868  void nd_nf_p(Obj f,LIST g,LIST v,int m,struct order_sp
     NDV ndvf;      NDV ndvf;
     VL vv,tv;      VL vv,tv;
     int stat,nvar,max,mrank;      int stat,nvar,max,mrank;
     union oNDC dn;  
     Q cont;      Q cont;
     P pp;      P pp;
     LIST ppl;      LIST ppl;
Line 6027  void nd_nf_p(Obj f,LIST g,LIST v,int m,struct order_sp
Line 7914  void nd_nf_p(Obj f,LIST g,LIST v,int m,struct order_sp
     ndf = (pointer)ndvtond(m,ndvf);      ndf = (pointer)ndvtond(m,ndvf);
   
     /* dont sort, dont removecont */      /* dont sort, dont removecont */
     ndv_setup(m,0,in0,1,1);      ndv_setup(m,0,in0,1,1,0);
     nd_scale=2;      nd_scale=2;
     stat = nd_nf(m,0,ndf,nd_ps,1,&nf);      stat = nd_nf(m,0,ndf,nd_ps,1,&nf);
     if ( !stat )      if ( !stat )
Line 6149  IndArray nm_ind_pair_to_vect_compress(int trace,UINT *
Line 8036  IndArray nm_ind_pair_to_vect_compress(int trace,UINT *
 {  {
     NM m;      NM m;
     NMV mr;      NMV mr;
     UINT *d,*t,*s,*u;      UINT *d,*t,*u;
     NDV p;      NDV p;
     unsigned char *ivc;      unsigned char *ivc;
     unsigned short *ivs;      unsigned short *ivs;
     UINT *v,*ivi,*s0v;      UINT *v,*ivi;
     int i,j,len,prev,diff,cdiff,h,st,ed,md,c;      int i,j,len,prev,diff,cdiff,st,ed,md,c;
     IndArray r;      IndArray r;
   
     m = pair->mul;      m = pair->mul;
Line 6315  int ndv_reduce_vect_q(Z *svect0,int trace,int col,IndA
Line 8202  int ndv_reduce_vect_q(Z *svect0,int trace,int col,IndA
     unsigned int *ivi;      unsigned int *ivi;
     NDV redv;      NDV redv;
     NMV mr;      NMV mr;
     NODE rp;  
     int maxrs;      int maxrs;
     double hmag;      double hmag;
     int l;  
     static mpz_t *svect;      static mpz_t *svect;
     static int svect_len=0;      static int svect_len=0;
   
Line 6400  int ndv_reduce_vect_q(Z *svect0,int trace,int col,IndA
Line 8285  int ndv_reduce_vect_q(Z *svect0,int trace,int col,IndA
 }  }
 #endif  #endif
   
 int ndv_reduce_vect(int m,UINT *svect,int col,IndArray *imat,NM_ind_pair *rp0,int nred)  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray *imat,NM_ind_pair *rp0,int nred,SIG sig)
 {  {
     int i,j,k,len,pos,prev;      int i,j,k,len,pos,prev;
     UINT c,c1,c2,c3,up,lo,dmy;      UINT c,c1,c2,c3,up,lo,dmy;
Line 6410  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
Line 8295  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
     unsigned int *ivi;      unsigned int *ivi;
     NDV redv;      NDV redv;
     NMV mr;      NMV mr;
     NODE rp;  
     int maxrs;      int maxrs;
   
     maxrs = 0;      maxrs = 0;
     for ( i = 0; i < nred; i++ ) {      for ( i = 0; i < nred; i++ ) {
         ivect = imat[i];          ivect = imat[i];
         k = ivect->head; svect[k] %= m;          k = ivect->head; svect[k] %= m;
         if ( (c = svect[k]) != 0 ) {          if ( (c = svect[k]) != 0 && (sig == 0 || comp_sig(sig,rp0[i]->sig) > 0 ) ) {
             maxrs = MAX(maxrs,rp0[i]->sugar);              maxrs = MAX(maxrs,rp0[i]->sugar);
             c = m-c; redv = nd_ps[rp0[i]->index];              c = m-c; redv = nd_ps[rp0[i]->index];
             len = LEN(redv); mr = BDY(redv);              len = LEN(redv); mr = BDY(redv);
Line 6427  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
Line 8311  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
                     ivc = ivect->index.c;                      ivc = ivect->index.c;
                     for ( j = 1, NMV_ADV(mr); j < len; j++, NMV_ADV(mr) ) {                      for ( j = 1, NMV_ADV(mr); j < len; j++, NMV_ADV(mr) ) {
                         pos = prev+ivc[j]; c1 = CM(mr); prev = pos;                          pos = prev+ivc[j]; c1 = CM(mr); prev = pos;
             if ( c1 ) {                          if ( c1 ) {
               c2 = svect[pos];                            c2 = svect[pos];
                           DMA(c1,c,c2,up,lo);                            DMA(c1,c,c2,up,lo);
                           if ( up ) { DSAB(m,up,lo,dmy,c3); svect[pos] = c3;                            if ( up ) { DSAB(m,up,lo,dmy,c3); svect[pos] = c3;
                           } else svect[pos] = lo;                            } else svect[pos] = lo;
             }                          }
                     }                      }
                     break;                      break;
                 case 2:                  case 2:
Line 6440  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
Line 8324  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
                     for ( j = 1, NMV_ADV(mr); j < len; j++, NMV_ADV(mr) ) {                      for ( j = 1, NMV_ADV(mr); j < len; j++, NMV_ADV(mr) ) {
                         pos = prev+ivs[j]; c1 = CM(mr);                          pos = prev+ivs[j]; c1 = CM(mr);
                         prev = pos;                          prev = pos;
             if ( c1 ) {                          if ( c1 ) {
               c2 = svect[pos];                            c2 = svect[pos];
                           DMA(c1,c,c2,up,lo);                            DMA(c1,c,c2,up,lo);
                           if ( up ) { DSAB(m,up,lo,dmy,c3); svect[pos] = c3;                            if ( up ) { DSAB(m,up,lo,dmy,c3); svect[pos] = c3;
                           } else svect[pos] = lo;                            } else svect[pos] = lo;
             }                          }
                     }                      }
                     break;                      break;
                 case 4:                  case 4:
Line 6453  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
Line 8337  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
                     for ( j = 1, NMV_ADV(mr); j < len; j++, NMV_ADV(mr) ) {                      for ( j = 1, NMV_ADV(mr); j < len; j++, NMV_ADV(mr) ) {
                         pos = prev+ivi[j]; c1 = CM(mr);                          pos = prev+ivi[j]; c1 = CM(mr);
                         prev = pos;                          prev = pos;
             if ( c1 ) {                          if ( c1 ) {
               c2 = svect[pos];                            c2 = svect[pos];
                           DMA(c1,c,c2,up,lo);                            DMA(c1,c,c2,up,lo);
                           if ( up ) { DSAB(m,up,lo,dmy,c3); svect[pos] = c3;                            if ( up ) { DSAB(m,up,lo,dmy,c3); svect[pos] = c3;
                           } else svect[pos] = lo;                            } else svect[pos] = lo;
             }                          }
                     }                      }
                     break;                      break;
             }              }
Line 6472  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
Line 8356  int ndv_reduce_vect(int m,UINT *svect,int col,IndArray
 int ndv_reduce_vect_sf(int m,UINT *svect,int col,IndArray *imat,NM_ind_pair *rp0,int nred)  int ndv_reduce_vect_sf(int m,UINT *svect,int col,IndArray *imat,NM_ind_pair *rp0,int nred)
 {  {
     int i,j,k,len,pos,prev;      int i,j,k,len,pos,prev;
     UINT c,c1,c2,c3,up,lo,dmy;      UINT c;
     IndArray ivect;      IndArray ivect;
     unsigned char *ivc;      unsigned char *ivc;
     unsigned short *ivs;      unsigned short *ivs;
     unsigned int *ivi;      unsigned int *ivi;
     NDV redv;      NDV redv;
     NMV mr;      NMV mr;
     NODE rp;  
     int maxrs;      int maxrs;
   
     maxrs = 0;      maxrs = 0;
Line 6521  int ndv_reduce_vect_sf(int m,UINT *svect,int col,IndAr
Line 8404  int ndv_reduce_vect_sf(int m,UINT *svect,int col,IndAr
   
 ND nd_add_lf(ND p1,ND p2)  ND nd_add_lf(ND p1,ND p2)
 {  {
     int n,c,can;      int c,can;
     ND r;  
     NM m1,m2,mr0,mr,s;      NM m1,m2,mr0,mr,s;
     Z t;      Z t;
   
Line 6530  ND nd_add_lf(ND p1,ND p2)
Line 8412  ND nd_add_lf(ND p1,ND p2)
     else if ( !p2 ) return p1;      else if ( !p2 ) return p1;
     else {      else {
         can = 0;          can = 0;
         for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {          for (  m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {
             c = DL_COMPARE(DL(m1),DL(m2));              c = DL_COMPARE(DL(m1),DL(m2));
             switch ( c ) {              switch ( c ) {
                 case 0:                  case 0:
Line 6569  ND nd_add_lf(ND p1,ND p2)
Line 8451  ND nd_add_lf(ND p1,ND p2)
 int ndv_reduce_vect_lf(mpz_t *svect,int trace,int col,IndArray *imat,NM_ind_pair *rp0,int nred)  int ndv_reduce_vect_lf(mpz_t *svect,int trace,int col,IndArray *imat,NM_ind_pair *rp0,int nred)
 {  {
     int i,j,k,len,pos,prev;      int i,j,k,len,pos,prev;
     mpz_t c,mc,c1;  
     IndArray ivect;      IndArray ivect;
     unsigned char *ivc;      unsigned char *ivc;
     unsigned short *ivs;      unsigned short *ivs;
     unsigned int *ivi;      unsigned int *ivi;
     NDV redv;      NDV redv;
     NMV mr;      NMV mr;
     NODE rp;  
     int maxrs;      int maxrs;
   
     maxrs = 0;      maxrs = 0;
Line 6627  int ndv_reduce_vect_lf(mpz_t *svect,int trace,int col,
Line 8507  int ndv_reduce_vect_lf(mpz_t *svect,int trace,int col,
 int nd_gauss_elim_lf(mpz_t **mat0,int *sugar,int row,int col,int *colstat)  int nd_gauss_elim_lf(mpz_t **mat0,int *sugar,int row,int col,int *colstat)
 {  {
     int i,j,k,l,rank,s;      int i,j,k,l,rank,s;
     mpz_t a,a1,inv;      mpz_t a,inv;
     mpz_t *t,*pivot,*pk;      mpz_t *t,*pivot,*pk;
     mpz_t **mat;      mpz_t **mat;
     struct oEGT eg0,eg1,eg_forward,eg_mod,eg_back;  
     int size,size1;  
   
     mpz_init(inv);      mpz_init(inv);
     mpz_init(a);      mpz_init(a);
     mat = (mpz_t **)mat0;      mat = (mpz_t **)mat0;
         size = 0;  
     for ( rank = 0, j = 0; j < col; j++ ) {      for ( rank = 0, j = 0; j < col; j++ ) {
         for ( i = rank; i < row; i++ ) {          for ( i = rank; i < row; i++ ) {
             mpz_mod(mat[i][j],mat[i][j],BDY(current_mod_lf));              mpz_mod(mat[i][j],mat[i][j],BDY(current_mod_lf));
Line 6726  NDV vect_to_ndv(UINT *vect,int spcol,int col,int *rhea
Line 8603  NDV vect_to_ndv(UINT *vect,int spcol,int col,int *rhea
     }      }
 }  }
   
 NDV vect_to_ndv_2(unsigned long *vect,int col,UINT *s0vect)  NDV vect_to_ndv_s(UINT *vect,int col,UINT *s0vect)
 {  {
     int j,k,len;      int j,k,len;
     UINT *p;      UINT *p;
       UINT c;
     NDV r;      NDV r;
     NMV mr0,mr;      NMV mr0,mr;
   
       for ( j = 0, len = 0; j < col; j++ ) if ( vect[j] ) len++;
       if ( !len ) return 0;
       else {
           mr0 = (NMV)MALLOC_ATOMIC_IGNORE_OFF_PAGE(nmv_adv*len);
           mr = mr0;
           p = s0vect;
           for ( j = k = 0; j < col; j++, p += nd_wpd )
             if ( (c = vect[k++]) != 0 ) {
               ndl_copy(p,DL(mr)); CM(mr) = c; NMV_ADV(mr);
             }
           MKNDV(nd_nvar,mr0,len,r);
           return r;
       }
   }
   
   NDV vect_to_ndv_2(unsigned long *vect,int col,UINT *s0vect)
   {
       int j,len;
       UINT *p;
       NDV r;
       NMV mr0,mr;
   
     for ( j = 0, len = 0; j < col; j++ ) if ( vect[j/BLEN] & (1L<<(j%BLEN)) ) len++;      for ( j = 0, len = 0; j < col; j++ ) if ( vect[j/BLEN] & (1L<<(j%BLEN)) ) len++;
     if ( !len ) return 0;      if ( !len ) return 0;
     else {      else {
Line 6772  NDV vect_to_ndv_q(Z *vect,int spcol,int col,int *rhead
Line 8672  NDV vect_to_ndv_q(Z *vect,int spcol,int col,int *rhead
         if ( (c = vect[k++]) != 0 ) {          if ( (c = vect[k++]) != 0 ) {
           if ( !INT(c) )            if ( !INT(c) )
             error("vect_to_ndv_q : components must be integers");              error("vect_to_ndv_q : components must be integers");
             ndl_copy(p,DL(mr)); CZ(mr) = c; NMV_ADV(mr);            ndl_copy(p,DL(mr)); CZ(mr) = c; NMV_ADV(mr);
         }          }
       }        }
     }      }
Line 6843  NDV plain_vect_to_ndv_q(Z *vect,int col,UINT *s0vect)
Line 8743  NDV plain_vect_to_ndv_q(Z *vect,int col,UINT *s0vect)
 int nd_sp_f4(int m,int trace,ND_pairs l,PGeoBucket bucket)  int nd_sp_f4(int m,int trace,ND_pairs l,PGeoBucket bucket)
 {  {
     ND_pairs t;      ND_pairs t;
     NODE sp0,sp;  
     int stat;      int stat;
     ND spol;      ND spol;
   
Line 6887  int nd_symbolic_preproc(PGeoBucket bucket,int trace,UI
Line 8786  int nd_symbolic_preproc(PGeoBucket bucket,int trace,UI
             if ( ndl_check_bound2(index,DL(mul)) )              if ( ndl_check_bound2(index,DL(mul)) )
                 return 0;                  return 0;
             sugar = TD(DL(mul))+SG(ps[index]);              sugar = TD(DL(mul))+SG(ps[index]);
             MKNM_ind_pair(pair,mul,index,sugar);              MKNM_ind_pair(pair,mul,index,sugar,0);
             red = ndv_mul_nm_symbolic(mul,ps[index]);              red = ndv_mul_nm_symbolic(mul,ps[index]);
             add_pbucket_symbolic(bucket,nd_remove_head(red));              add_pbucket_symbolic(bucket,nd_remove_head(red));
             NEXTNODE(rp0,rp); BDY(rp) = (pointer)pair;              NEXTNODE(rp0,rp); BDY(rp) = (pointer)pair;
Line 6915  void print_ndp(ND_pairs l)
Line 8814  void print_ndp(ND_pairs l)
   
 NODE nd_f4(int m,int checkonly,int **indp)  NODE nd_f4(int m,int checkonly,int **indp)
 {  {
     int i,nh,stat,index,f4red;    int i,nh,stat,f4red;
     NODE r,g,tn0,tn,node;    NODE r,g,tn0,tn,node;
     ND_pairs d,l,t,ll0,ll,lh;    ND_pairs d,l,t,ll,lh;
   LIST l0,l1;    LIST l0,l1;
     ND spol,red;    NDV nf;
     NDV nf,redv;    NODE rp0,nflist,nzlist,nzlist_t;
     NM s0,s;    int col;
     NODE rp0,srp0,nflist,nzlist,nzlist_t;    UINT *s0vect;
     int nsp,nred,col,rank,len,k,j,a,i1s,i2s;    int sugar;
     UINT c;    PGeoBucket bucket;
     UINT **spmat;    struct oEGT eg0,eg1,eg_f4;
     UINT *s0vect,*svect,*p,*v;    Z i1,i2,sugarq;
     int *colstat;  
     IndArray *imat;  
     int *rhead;  
     int spcol,sprow;  
     int sugar,sugarh;  
     PGeoBucket bucket;  
     struct oEGT eg0,eg1,eg_f4;  
     Z i1,i2,sugarq;  
   
     init_eg(&f4_symb); init_eg(&f4_conv); init_eg(&f4_conv); init_eg(&f4_elim1); init_eg(&f4_elim2);    init_eg(&f4_symb); init_eg(&f4_conv); init_eg(&f4_conv); init_eg(&f4_elim1); init_eg(&f4_elim2);
 #if 0  #if 0
     ndv_alloc = 0;    ndv_alloc = 0;
 #endif  #endif
     Nf4_red=0;    Nf4_red=0;
     g = 0; d = 0;    g = 0; d = 0;
     for ( i = 0; i < nd_psn; i++ ) {    for ( i = 0; i < nd_psn; i++ ) {
         d = update_pairs(d,g,i,0);      d = update_pairs(d,g,i,0);
         g = update_base(g,i);      g = update_base(g,i);
     }    }
   nzlist = 0;    nzlist = 0;
     nzlist_t = nd_nzlist;    nzlist_t = nd_nzlist;
     f4red = 1;    f4red = 1;
     nd_last_nonzero = 0;    nd_last_nonzero = 0;
     while ( d ) {    while ( d ) {
         get_eg(&eg0);      get_eg(&eg0);
         l = nd_minsugarp(d,&d);      l = nd_minsugarp(d,&d);
         sugar = nd_sugarweight?l->sugar2:SG(l);      sugar = nd_sugarweight?l->sugar2:SG(l);
         if ( MaxDeg > 0 && sugar > MaxDeg ) break;      if ( MaxDeg > 0 && sugar > MaxDeg ) break;
         if ( nzlist_t ) {      if ( nzlist_t ) {
             node = BDY((LIST)BDY(nzlist_t));        node = BDY((LIST)BDY(nzlist_t));
             sugarh = ZTOS((Q)ARG0(node));  //      sugarh = (int)ZTOS((Q)ARG0(node));
             tn = BDY((LIST)ARG1(node));        tn = BDY((LIST)ARG1(node));
             if ( !tn ) {        if ( !tn ) {
               nzlist_t = NEXT(nzlist_t);          nzlist_t = NEXT(nzlist_t);
               continue;          continue;
             }        }
             /* tn = [[i1,i2],...] */        /* tn = [[i1,i2],...] */
             lh = nd_ipairtospair(tn);            lh = nd_ipairtospair(tn);
         }      }
         bucket = create_pbucket();      bucket = create_pbucket();
         stat = nd_sp_f4(m,0,l,bucket);      stat = nd_sp_f4(m,0,l,bucket);
         if ( !stat ) {      if ( !stat ) {
             for ( t = l; NEXT(t); t = NEXT(t) );        for ( t = l; NEXT(t); t = NEXT(t) );
             NEXT(t) = d; d = l;        NEXT(t) = d; d = l;
             d = nd_reconstruct(0,d);        d = nd_reconstruct(0,d);
             continue;        continue;
         }      }
         if ( bucket->m < 0 ) continue;      if ( bucket->m < 0 ) continue;
         col = nd_symbolic_preproc(bucket,0,&s0vect,&rp0);      col = nd_symbolic_preproc(bucket,0,&s0vect,&rp0);
         if ( !col ) {      if ( !col ) {
             for ( t = l; NEXT(t); t = NEXT(t) );        for ( t = l; NEXT(t); t = NEXT(t) );
             NEXT(t) = d; d = l;        NEXT(t) = d; d = l;
             d = nd_reconstruct(0,d);        d = nd_reconstruct(0,d);
             continue;        continue;
         }      }
         get_eg(&eg1); init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg1); add_eg(&f4_symb,&eg0,&eg1);      get_eg(&eg1); init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg1); add_eg(&f4_symb,&eg0,&eg1);
         if ( DP_Print )      if ( DP_Print )
             fprintf(asir_out,"sugar=%d,symb=%.3fsec,",        fprintf(asir_out,"sugar=%d,symb=%.3fsec,",sugar,eg_f4.exectime);
                 sugar,eg_f4.exectime);      nflist = nd_f4_red(m,nd_nzlist?lh:l,0,s0vect,col,rp0,nd_gentrace?&ll:0);
         nflist = nd_f4_red(m,nd_nzlist?lh:l,0,s0vect,col,rp0,nd_gentrace?&ll:0);      if ( checkonly && nflist ) return 0;
         if ( checkonly && nflist ) return 0;      /* adding new bases */
         /* adding new bases */      if ( nflist ) nd_last_nonzero = f4red;
         if ( nflist ) nd_last_nonzero = f4red;      for ( r = nflist; r; r = NEXT(r) ) {
         for ( r = nflist; r; r = NEXT(r) ) {        nf = (NDV)BDY(r);
             nf = (NDV)BDY(r);        if ( nd_f4_td ) SG(nf) = nd_tdeg(nf);
             ndv_removecont(m,nf);        ndv_removecont(m,nf);
             if ( !m && nd_nalg ) {        if ( !m && nd_nalg ) {
                 ND nf1;          ND nf1;
   
                 nf1 = ndvtond(m,nf);          nf1 = ndvtond(m,nf);
                 nd_monic(0,&nf1);          nd_monic(0,&nf1);
                 nd_removecont(m,nf1);          nd_removecont(m,nf1);
                 nf = ndtondv(m,nf1);          nf = ndtondv(m,nf1);
             }        }
             nh = ndv_newps(m,nf,0,1);        nh = ndv_newps(m,nf,0);
             d = update_pairs(d,g,nh,0);        d = update_pairs(d,g,nh,0);
             g = update_base(g,nh);        g = update_base(g,nh);
         }      }
         if ( DP_Print ) {      if ( DP_Print ) {
           fprintf(asir_out,"f4red=%d,gblen=%d\n",f4red,length(g)); fflush(asir_out);        fprintf(asir_out,"f4red=%d,gblen=%d\n",f4red,length(g)); fflush(asir_out);
         }      }
         if ( nd_gentrace ) {      if ( nd_gentrace ) {
       for ( t = ll, tn0 = 0; t; t = NEXT(t) ) {        for ( t = ll, tn0 = 0; t; t = NEXT(t) ) {
         NEXTNODE(tn0,tn);          NEXTNODE(tn0,tn);
                 STOZ(t->i1,i1); STOZ(t->i2,i2);          STOZ(t->i1,i1); STOZ(t->i2,i2);
                 node = mknode(2,i1,i2); MKLIST(l0,node);          node = mknode(2,i1,i2); MKLIST(l0,node);
         BDY(tn) = l0;          BDY(tn) = l0;
       }        }
       if ( tn0 ) NEXT(tn) = 0; MKLIST(l0,tn0);        if ( tn0 ) NEXT(tn) = 0; MKLIST(l0,tn0);
             STOZ(sugar,sugarq); node = mknode(2,sugarq,l0); MKLIST(l1,node);          STOZ(sugar,sugarq); node = mknode(2,sugarq,l0); MKLIST(l1,node);
             MKNODE(node,l1,nzlist); nzlist = node;          MKNODE(node,l1,nzlist); nzlist = node;
         }  
         if ( nd_nzlist ) nzlist_t = NEXT(nzlist_t);  
         f4red++;  
         if ( nd_f4red && f4red > nd_f4red ) break;  
         if ( nd_rank0 && !nflist ) break;  
     }      }
     if ( nd_gentrace ) {      if ( nd_nzlist ) nzlist_t = NEXT(nzlist_t);
     MKLIST(l0,reverse_node(nzlist));      f4red++;
         MKNODE(nd_alltracelist,l0,0);      if ( nd_f4red && f4red > nd_f4red ) break;
     }      if ( nd_rank0 && !nflist ) break;
     }
     if ( nd_gentrace ) {
     MKLIST(l0,reverse_node(nzlist));
         MKNODE(nd_alltracelist,l0,0);
     }
 #if 0  #if 0
     fprintf(asir_out,"ndv_alloc=%d\n",ndv_alloc);    fprintf(asir_out,"ndv_alloc=%d\n",ndv_alloc);
 #endif  #endif
   if ( DP_Print ) {    if ( DP_Print ) {
     fprintf(asir_out,"number of red=%d,",Nf4_red);    fprintf(asir_out,"number of red=%d,",Nf4_red);
     fprintf(asir_out,"symb=%.3fsec,conv=%.3fsec,elim1=%.3fsec,elim2=%.3fsec\n",    fprintf(asir_out,"symb=%.3fsec,conv=%.3fsec,elim1=%.3fsec,elim2=%.3fsec\n",
       f4_symb.exectime,f4_conv.exectime,f4_elim1.exectime,f4_elim2.exectime);      f4_symb.exectime,f4_conv.exectime,f4_elim1.exectime,f4_elim2.exectime);
     fprintf(asir_out,"number of removed pairs=%d\n,",NcriB+NcriMF+Ncri2);
   }    }
   conv_ilist(nd_demand,0,g,indp);    conv_ilist(nd_demand,0,g,indp);
     return g;    return g;
 }  }
   
 NODE nd_f4_trace(int m,int **indp)  NODE nd_f4_trace(int m,int **indp)
 {  {
     int i,nh,stat,index;      int i,nh,stat;
     NODE r,g;      NODE r,g;
     ND_pairs d,l,l0,t;      ND_pairs d,l,l0,t;
     ND spol,red;      NDV nfqv,nfv;
     NDV nf,redv,nfqv,nfv;      NODE rp0,nflist;
     NM s0,s;      int col;
     NODE rp0,srp0,nflist;      UINT *s0vect;
     int nsp,nred,col,rank,len,k,j,a;  
     UINT c;  
     UINT **spmat;  
     UINT *s0vect,*svect,*p,*v;  
     int *colstat;  
     IndArray *imat;  
     int *rhead;  
     int spcol,sprow;  
     int sugar;      int sugar;
     PGeoBucket bucket;      PGeoBucket bucket;
     struct oEGT eg0,eg1,eg_f4;      struct oEGT eg0,eg1,eg_f4;
Line 7129  NODE nd_f4_trace(int m,int **indp)
Line 9013  NODE nd_f4_trace(int m,int **indp)
             nfv = ndv_dup(0,nfqv);              nfv = ndv_dup(0,nfqv);
             ndv_mod(m,nfv);              ndv_mod(m,nfv);
             ndv_removecont(m,nfv);              ndv_removecont(m,nfv);
             nh = ndv_newps(0,nfv,nfqv,1);              nh = ndv_newps(0,nfv,nfqv);
             d = update_pairs(d,g,nh,0);              d = update_pairs(d,g,nh,0);
             g = update_base(g,nh);              g = update_base(g,nh);
         }          }
Line 7243  void red_by_vect_2(matrix mat,int *sugar,unsigned long
Line 9127  void red_by_vect_2(matrix mat,int *sugar,unsigned long
   
 NODE nd_f4_red_2(ND_pairs sp0,UINT *s0vect,int col,NODE rp0,ND_pairs *nz)  NODE nd_f4_red_2(ND_pairs sp0,UINT *s0vect,int col,NODE rp0,ND_pairs *nz)
 {  {
     int nsp,nred,i,i0,k,rank,row;    int nsp,i,rank;
     NODE r0,rp;    NODE r0,rp;
     ND_pairs sp;    ND_pairs sp;
   ND spol;    ND spol;
   NM_ind_pair rt;    NM_ind_pair rt;
     int *s0hash;    int *s0hash;
   UINT *s;    UINT *s;
   int *pivot,*sugar,*head;    int *sugar;
   matrix mat;    matrix mat;
     NM m;    NODE r;
     NODE r;  
   struct oEGT eg0,eg1,eg2,eg_elim1,eg_elim2;    struct oEGT eg0,eg1,eg2,eg_elim1,eg_elim2;
   int rhead,rsugar,size;    int rhead,rsugar,size;
     unsigned long *v;    unsigned long *v;
   
     get_eg(&eg0);    get_eg(&eg0);
     for ( sp = sp0, nsp = 0; sp; sp = NEXT(sp), nsp++ );    for ( sp = sp0, nsp = 0; sp; sp = NEXT(sp), nsp++ );
     nred = length(rp0);  //  nred = length(rp0);
     mat = alloc_matrix(nsp,col);    mat = alloc_matrix(nsp,col);
     s0hash = (int *)MALLOC(col*sizeof(int));    s0hash = (int *)MALLOC(col*sizeof(int));
     for ( i = 0, s = s0vect; i < col; i++, s += nd_wpd )    for ( i = 0, s = s0vect; i < col; i++, s += nd_wpd )
         s0hash[i] = ndl_hash_value(s);      s0hash[i] = ndl_hash_value(s);
   
   sugar = (int *)MALLOC(nsp*sizeof(int));    sugar = (int *)MALLOC(nsp*sizeof(int));
   for ( i = 0, sp = sp0; sp; sp = NEXT(sp) ) {    for ( i = 0, sp = sp0; sp; sp = NEXT(sp) ) {
Line 7318  NODE nd_f4_red(int m,ND_pairs sp0,int trace,UINT *s0ve
Line 9201  NODE nd_f4_red(int m,ND_pairs sp0,int trace,UINT *s0ve
     NODE r0,rp;      NODE r0,rp;
     ND_pairs sp;      ND_pairs sp;
     NM_ind_pair *rvect;      NM_ind_pair *rvect;
     UINT *s;  
     int *s0hash;  
     struct oEGT eg0,eg1,eg_conv;      struct oEGT eg0,eg1,eg_conv;
   
     if ( m == 2 && nd_rref2 )      if ( m == 2 && nd_rref2 )
Line 7370  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
Line 9251  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
         NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred,ND_pairs *nz)          NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred,ND_pairs *nz)
 {  {
     int spcol,sprow,a;      int spcol,sprow,a;
     int i,j,k,l,rank;      int i,j,k,rank;
     NODE r0,r;      NODE r0,r;
     ND_pairs sp;      ND_pairs sp;
     ND spol;      ND spol;
Line 7396  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
Line 9277  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
         if ( m == -1 )          if ( m == -1 )
             maxrs = ndv_reduce_vect_sf(m,svect,col,imat,rvect,nred);              maxrs = ndv_reduce_vect_sf(m,svect,col,imat,rvect,nred);
         else          else
             maxrs = ndv_reduce_vect(m,svect,col,imat,rvect,nred);              maxrs = ndv_reduce_vect(m,svect,col,imat,rvect,nred,0);
         for ( i = 0; i < col; i++ ) if ( svect[i] ) break;          for ( i = 0; i < col; i++ ) if ( svect[i] ) break;
         if ( i < col ) {          if ( i < col ) {
             spmat[sprow] = v = (UINT *)MALLOC_ATOMIC(spcol*sizeof(UINT));              spmat[sprow] = v = (UINT *)MALLOC_ATOMIC(spcol*sizeof(UINT));
Line 7452  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
Line 9333  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
     return r0;      return r0;
 }  }
   
   NODE nd_f4_red_main_s(int m,ND_pairs sp0,int nsp,UINT *s0vect,int col,
           NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred,NODE *syzlistp)
   {
       int sprow,a;
       int i,rank;
       NODE r0,r;
       ND_pairs sp;
       ND spol;
       UINT **spmat;
       UINT *svect;
       int *colstat;
       struct oEGT eg0,eg1,eg2,eg_f4,eg_f4_1,eg_f4_2;
       int maxrs;
       int *spsugar;
       SIG *spsig;
   
       get_eg(&eg0);
       /* elimination (1st step) */
       spmat = (UINT **)MALLOC(nsp*sizeof(UINT *));
       spsugar = (int *)MALLOC(nsp*sizeof(int));
       spsig = (SIG *)MALLOC(nsp*sizeof(SIG));
       for ( a = sprow = 0, sp = sp0; a < nsp; a++, sp = NEXT(sp) ) {
           nd_sp(m,0,sp,&spol);
           if ( !spol ) {
             syzlistp[sp->sig->pos] = insert_sig(syzlistp[sp->sig->pos],sp->sig);
             continue;
           }
           svect = (UINT *)MALLOC(col*sizeof(UINT));
           nd_to_vect(m,s0vect,col,spol,svect);
           maxrs = ndv_reduce_vect(m,svect,col,imat,rvect,nred,spol->sig);
           for ( i = 0; i < col; i++ ) if ( svect[i] ) break;
           if ( i < col ) {
               spmat[sprow] = svect;
               spsugar[sprow] = MAX(maxrs,SG(spol));
               spsig[sprow] = sp->sig;
               sprow++;
           } else {
             syzlistp[sp->sig->pos] = insert_sig(syzlistp[sp->sig->pos],sp->sig);
           }
           nd_free(spol);
       }
       get_eg(&eg1); init_eg(&eg_f4_1); add_eg(&eg_f4_1,&eg0,&eg1); add_eg(&f4_elim1,&eg0,&eg1);
       if ( DP_Print ) {
           fprintf(asir_out,"elim1=%.3fsec,",eg_f4_1.exectime);
           fflush(asir_out);
       }
       /* free index arrays */
       for ( i = 0; i < nred; i++ ) GCFREE(imat[i]->index.c);
   
       /* elimination (2nd step) */
       colstat = (int *)MALLOC(col*sizeof(int));
       rank = nd_gauss_elim_mod_s(spmat,spsugar,0,sprow,col,m,colstat,spsig);
       r0 = 0;
       for ( i = 0; i < sprow; i++ ) {
           if ( spsugar[i] >= 0 ) {
             NEXTNODE(r0,r);
             BDY(r) = vect_to_ndv_s(spmat[i],col,s0vect);
             SG((NDV)BDY(r)) = spsugar[i];
             ((NDV)BDY(r))->sig = spsig[i];
           } else
             syzlistp[spsig[i]->pos] = insert_sig(syzlistp[spsig[i]->pos],spsig[i]);
           GCFREE(spmat[i]);
       }
       if ( r0 ) NEXT(r) = 0;
       get_eg(&eg2); init_eg(&eg_f4_2); add_eg(&eg_f4_2,&eg1,&eg2); add_eg(&f4_elim2,&eg1,&eg2);
       init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg2);
       if ( DP_Print ) {
           fprintf(asir_out,"elim2=%.3fsec,",eg_f4_2.exectime);
           fprintf(asir_out,"nsp=%d,nred=%d,spmat=(%d,%d),rank=%d ",
               nsp,nred,sprow,col,rank);
           fprintf(asir_out,"%.3fsec,",eg_f4.exectime);
       }
       return r0;
   }
   
   
 /* for small finite fields */  /* for small finite fields */
   
 NODE nd_f4_red_sf_main(int m,ND_pairs sp0,int nsp,UINT *s0vect,int col,  NODE nd_f4_red_sf_main(int m,ND_pairs sp0,int nsp,UINT *s0vect,int col,
         NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred,ND_pairs *nz)          NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred,ND_pairs *nz)
 {  {
     int spcol,sprow,a;      int spcol,sprow,a;
     int i,j,k,l,rank;      int i,j,k,rank;
     NODE r0,r;      NODE r0,r;
     ND_pairs sp;      ND_pairs sp;
     ND spol;      ND spol;
Line 7539  NODE nd_f4_red_lf_main(int m,ND_pairs sp0,int nsp,int 
Line 9495  NODE nd_f4_red_lf_main(int m,ND_pairs sp0,int nsp,int 
         NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred)          NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred)
 {  {
     int spcol,sprow,a;      int spcol,sprow,a;
     int i,j,k,l,rank;      int i,j,k,rank;
     NODE r0,r;      NODE r0,r;
     ND_pairs sp;      ND_pairs sp;
     ND spol;      ND spol;
Line 7621  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
Line 9577  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
         NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred)          NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred)
 {  {
     int spcol,sprow,a;      int spcol,sprow,a;
     int i,j,k,l,rank;      int i,j,k,rank;
     NODE r0,r;      NODE r0,r;
     ND_pairs sp;      ND_pairs sp;
     ND spol;      ND spol;
Line 7775  NDV nd_recv_ndv()
Line 9731  NDV nd_recv_ndv()
   
 int nd_gauss_elim_q(Z **mat0,int *sugar,int row,int col,int *colstat)  int nd_gauss_elim_q(Z **mat0,int *sugar,int row,int col,int *colstat)
 {  {
     int i,j,t,c,rank,inv;      int i,j,c,rank;
     int *ci,*ri;      int *ci,*ri;
     Z dn;      Z dn;
     MAT m,nm;      MAT m,nm;
Line 7864  int nd_gauss_elim_mod(UINT **mat0,int *sugar,ND_pairs 
Line 9820  int nd_gauss_elim_mod(UINT **mat0,int *sugar,ND_pairs 
     return rank;      return rank;
 }  }
   
   int nd_gauss_elim_mod_s(UINT **mat,int *sugar,ND_pairs *spactive,int row,int col,int md,int *colstat,SIG *sig)
   {
     int i,j,k,rank,s;
     UINT inv;
     UINT a;
     UINT *pk;
     int *used;
   
     used = (int *)MALLOC(row*sizeof(int));
     for ( j = 0; j < col; j++ ) {
       for ( i = 0; i < row; i++ )
         a = mat[i][j] %= md;
       for ( i = 0; i < row; i++ )
         if ( !used[i] && mat[i][j] ) break;
       if ( i == row ) {
         colstat[j] = 0;
         continue;
       } else {
         colstat[j] = 1;
         used[i] = 1;
       }
       /* column j is normalized */
       s = sugar[i];
       inv = invm(mat[i][j],md);
       /* normalize pivot row */
       for ( k = j, pk = mat[i]+j; k < col; k++, pk++ ) {
         DMAR(*pk,inv,0,md,*pk);
       }
       for ( k = i+1; k < row; k++ ) {
         if ( (a = mat[k][j]) != 0 ) {
           sugar[k] = MAX(sugar[k],s);
           red_by_vect(md,mat[k]+j,mat[i]+j,(int)(md-a),col-j);
           Nf4_red++;
         }
       }
     }
     rank = 0;
     for ( i = 0; i < row; i++ ) {
       for ( j = 0; j < col; j++ )
         if ( mat[i][j] ) break;
       if ( j == col ) sugar[i] = -1;
       else rank++;
     }
     return rank;
   }
   
   
 int nd_gauss_elim_sf(UINT **mat0,int *sugar,int row,int col,int md,int *colstat)  int nd_gauss_elim_sf(UINT **mat0,int *sugar,int row,int col,int md,int *colstat)
 {  {
     int i,j,k,l,inv,a,rank,s;      int i,j,k,l,inv,a,rank,s;
Line 7938  void ndv_save(NDV p,int index)
Line 9940  void ndv_save(NDV p,int index)
     FILE *s;      FILE *s;
     char name[BUFSIZ];      char name[BUFSIZ];
     short id;      short id;
     int nv,sugar,len,n,i,td,e,j;      int nv,sugar,len,i,td,e,j;
     NMV m;      NMV m;
     unsigned int *dl;      unsigned int *dl;
     int mpos;      int mpos;
Line 8004  NDV ndv_load(int index)
Line 10006  NDV ndv_load(int index)
     FILE *s;      FILE *s;
     char name[BUFSIZ];      char name[BUFSIZ];
     short id;      short id;
     int nv,sugar,len,n,i,td,e,j;      int nv,sugar,len,i,td,e,j;
     NDV d;      NDV d;
     NMV m0,m;      NMV m0,m;
     unsigned int *dl;      unsigned int *dl;
Line 8237  ND ndv_mul_nmv_trunc(int mod,NMV m0,NDV p,UINT *d)
Line 10239  ND ndv_mul_nmv_trunc(int mod,NMV m0,NDV p,UINT *d)
     NM mr,mr0;      NM mr,mr0;
     NM tnm;      NM tnm;
     NMV m;      NMV m;
     UINT *d0,*dt,*dm;      UINT *d0;
     int c,n,td,i,c1,c2,len;      int c,i,c1,c2,len;
     Z q;      Z q;
     ND r;      ND r;
   
     if ( !p ) return 0;      if ( !p ) return 0;
     else {      else {
         n = NV(p); m = BDY(p); len = LEN(p);          m = BDY(p); len = LEN(p);
         d0 = DL(m0);          d0 = DL(m0);
         td = TD(d);  
         mr0 = 0;          mr0 = 0;
         NEWNM(tnm);          NEWNM(tnm);
         if ( mod ) {          if ( mod ) {
Line 8284  ND ndv_mul_nmv_trunc(int mod,NMV m0,NDV p,UINT *d)
Line 10285  ND ndv_mul_nmv_trunc(int mod,NMV m0,NDV p,UINT *d)
   
 void nd_det_reconstruct(NDV **dm,int n,int j,NDV d)  void nd_det_reconstruct(NDV **dm,int n,int j,NDV d)
 {  {
     int i,obpe,oadv,h,k,l;      int obpe,oadv,k,l;
     static NM prev_nm_free_list;      static NM prev_nm_free_list;
     EPOS oepos;      EPOS oepos;
   
Line 8319  void nd_det_reconstruct(NDV **dm,int n,int j,NDV d)
Line 10320  void nd_det_reconstruct(NDV **dm,int n,int j,NDV d)
   
 UINT *nd_det_compute_bound(NDV **dm,int n,int j)  UINT *nd_det_compute_bound(NDV **dm,int n,int j)
 {  {
     UINT *d0,*d1,*d,*t,*r;      UINT *d0,*d;
     int k,l,i;      int k,l,i;
   
     d0 = (UINT *)MALLOC(nd_nvar*sizeof(UINT));      d0 = (UINT *)MALLOC(nd_nvar*sizeof(UINT));
Line 8336  UINT *nd_det_compute_bound(NDV **dm,int n,int j)
Line 10337  UINT *nd_det_compute_bound(NDV **dm,int n,int j)
   
 DL nd_separate_d(UINT *d,UINT *trans)  DL nd_separate_d(UINT *d,UINT *trans)
 {  {
     int n,td,i,e,j;      int td,i,e,j;
     DL a;      DL a;
   
     ndl_zero(trans);      ndl_zero(trans);
Line 8369  DL nd_separate_d(UINT *d,UINT *trans)
Line 10370  DL nd_separate_d(UINT *d,UINT *trans)
   
 int nd_monic(int mod,ND *p)  int nd_monic(int mod,ND *p)
 {  {
     UINT *trans,*t;  
     DL alg;      DL alg;
     MP mp0,mp;      MP mp0,mp;
     NM m,m0,m1,ma0,ma,mb,mr0,mr;      NM m,ma0,ma,mb,mr0,mr;
     ND r;      ND r;
     DL dl;      DL dl;
     DP nm;      DP nm;
     NDV ndv;  
     DAlg inv,cd;      DAlg inv,cd;
     ND s,c;  
     Z l,mul;      Z l,mul;
     Z ln;      Z ln;
     int n,ntrans,i,e,td,is_lc,len;      int i,e,td,is_lc,len;
     NumberField nf;      NumberField nf;
     struct oEGT eg0,eg1;      struct oEGT eg0,eg1;
   
Line 8540  void conv_ilist(int demand,int trace,NODE g,int **indp
Line 10538  void conv_ilist(int demand,int trace,NODE g,int **indp
   if ( indp ) *indp = ind;    if ( indp ) *indp = ind;
 }  }
   
 void parse_nd_option(NODE opt)  NODE conv_ilist_s(int demand,int trace,int **indp)
 {  {
     NODE t,p,u;    int n,i;
     int *ind;
     NODE g0,g;
   
     n = nd_psn;
     ind = (int *)MALLOC(n*sizeof(int));
     g0 = 0;
     for ( i = 0; i < n; i++ ) {
       ind[i] = i;
       NEXTNODE(g0,g);
       BDY(g) = (pointer)(demand?ndv_load(i):(trace?nd_ps_trace[i]:nd_ps[i]));
     }
     if ( g0 ) NEXT(g) = 0;
     if ( indp ) *indp = ind;
     return g0;
   }
   
   void parse_nd_option(VL vl,NODE opt)
   {
     NODE t,p,u;
   int i,s,n;    int i,s,n;
     char *key;    char *key;
     Obj value;    Obj value;
     VL oldvl;
   
     nd_gentrace = 0; nd_gensyz = 0; nd_nora = 0; nd_gbblock = 0;    nd_gentrace = 0; nd_gensyz = 0; nd_nora = 0; nd_norb = 0; nd_gbblock = 0;
   nd_newelim = 0; nd_intersect = 0; nd_nzlist = 0;    nd_newelim = 0; nd_intersect = 0; nd_nzlist = 0;
   nd_splist = 0; nd_check_splist = 0;    nd_splist = 0; nd_check_splist = 0;
     nd_sugarweight = 0;    nd_sugarweight = 0; nd_f4red =0; nd_rank0 = 0;
     nd_f4red =0;    nd_f4_td = 0; nd_sba_f4step = 2; nd_sba_pot = 0; nd_sba_largelcm = 0;
     nd_rank0 = 0;    nd_sba_dontsort = 0; nd_top = 0; nd_sba_redundant_check = 0;
     for ( t = opt; t; t = NEXT(t) ) {    nd_sba_syz = 0; nd_sba_modord = 0; nd_sba_inputisgb = 0;
         p = BDY((LIST)BDY(t));    nd_hpdata = 0; nd_sba_heu = 0;
         key = BDY((STRING)BDY(p));  
         value = (Obj)BDY(NEXT(p));    for ( t = opt; t; t = NEXT(t) ) {
         if ( !strcmp(key,"gentrace") )      p = BDY((LIST)BDY(t));
             nd_gentrace = value?1:0;      key = BDY((STRING)BDY(p));
         else if ( !strcmp(key,"gensyz") )      value = (Obj)BDY(NEXT(p));
             nd_gensyz = value?1:0;      if ( !strcmp(key,"gentrace") )
         else if ( !strcmp(key,"nora") )        nd_gentrace = value?1:0;
             nd_nora = value?1:0;      else if ( !strcmp(key,"gensyz") )
         else if ( !strcmp(key,"gbblock") ) {        nd_gensyz = value?1:0;
             if ( value && OID(value) == O_LIST ) {      else if ( !strcmp(key,"nora") )
         nd_nora = value?1:0;
       else if ( !strcmp(key,"norb") )
         nd_norb = value?1:0;
       else if ( !strcmp(key,"gbblock") ) {
         if ( value && OID(value) == O_LIST ) {
         u = BDY((LIST)value);          u = BDY((LIST)value);
               nd_gbblock = MALLOC((2*length(u)+1)*sizeof(int));          nd_gbblock = MALLOC((2*length(u)+1)*sizeof(int));
         for ( i = 0; u; u = NEXT(u) ) {          for ( i = 0; u; u = NEXT(u) ) {
           p = BDY((LIST)BDY(u));            p = BDY((LIST)BDY(u));
           s = nd_gbblock[i++] = ZTOS((Q)BDY(p));            s = nd_gbblock[i++] = (int)ZTOS((Q)BDY(p));
           nd_gbblock[i++] = s+ZTOS((Q)BDY(NEXT(p)))-1;            nd_gbblock[i++] = (int)(s+ZTOS((Q)BDY(NEXT(p)))-1);
         }          }
         nd_gbblock[i] = -1;          nd_gbblock[i] = -1;
             } else        } else
               nd_gbblock = 0;          nd_gbblock = 0;
     } else if ( !strcmp(key,"newelim") )      } else if ( !strcmp(key,"newelim") )
             nd_newelim = value?1:0;              nd_newelim = value?1:0;
     else if ( !strcmp(key,"intersect") )      else if ( !strcmp(key,"intersect") )
             nd_intersect = value?1:0;              nd_intersect = value?1:0;
     else if ( !strcmp(key,"syzgen") )      else if ( !strcmp(key,"syzgen") )
             nd_intersect = ZTOS((Q)value);              nd_intersect = (int)ZTOS((Q)value);
     else if ( !strcmp(key,"lf") )      else if ( !strcmp(key,"lf") )
             nd_lf = value?1:0;              nd_lf = value?1:0;
     else if ( !strcmp(key,"trace") ) {      else if ( !strcmp(key,"trace") ) {
            if ( value ) {        if ( value ) {
                u = BDY((LIST)value);          u = BDY((LIST)value);
            nd_nzlist = BDY((LIST)ARG2(u));          nd_nzlist = BDY((LIST)ARG2(u));
            nd_bpe = ZTOS((Q)ARG3(u));          nd_bpe = (int)ZTOS((Q)ARG3(u));
            }        }
     } else if ( !strcmp(key,"f4red") ) {      } else if ( !strcmp(key,"f4red") ) {
        nd_f4red = ZTOS((Q)value);        nd_f4red = (int)ZTOS((Q)value);
     } else if ( !strcmp(key,"rank0") ) {      } else if ( !strcmp(key,"rank0") ) {
             nd_rank0 = value?1:0;        nd_rank0 = value?1:0;
     } else if ( !strcmp(key,"splist") ) {      } else if ( !strcmp(key,"splist") ) {
             nd_splist = value?1:0;        nd_splist = value?1:0;
     } else if ( !strcmp(key,"check_splist") ) {      } else if ( !strcmp(key,"check_splist") ) {
       nd_check_splist = BDY((LIST)value);        nd_check_splist = BDY((LIST)value);
       } else if ( !strcmp(key,"hpdata") ) {
         if ( value )
           nd_hpdata = BDY((LIST)value);
     } else if ( !strcmp(key,"sugarweight") ) {      } else if ( !strcmp(key,"sugarweight") ) {
       u = BDY((LIST)value);        u = BDY((LIST)value);
             n = length(u);        n = length(u);
             nd_sugarweight = MALLOC(n*sizeof(int));        nd_sugarweight = MALLOC(n*sizeof(int));
       for ( i = 0; i < n; i++, u = NEXT(u) )        for ( i = 0; i < n; i++, u = NEXT(u) )
                 nd_sugarweight[i] = ZTOS((Q)BDY(u));          nd_sugarweight[i] = (int)ZTOS((Q)BDY(u));
       } else if ( !strcmp(key,"f4_td") ) {
         nd_f4_td = value?1:0;
       } else if ( !strcmp(key,"sba_f4step") ) {
         nd_sba_f4step = value?(int)ZTOS((Q)value):0;
       } else if ( !strcmp(key,"sba_pot") ) {
         nd_sba_pot = (int)ZTOS((Q)value);
       } else if ( !strcmp(key,"sba_largelcm") ) {
         nd_sba_largelcm = value?1:0;
       } else if ( !strcmp(key,"sba_dontsort") ) {
         nd_sba_dontsort = value?1:0;
       } else if ( !strcmp(key,"sba_syz") ) {
         nd_sba_syz = value?1:0;
       } else if ( !strcmp(key,"sba_heu") ) {
         nd_sba_heu = value?1:0;
       } else if ( !strcmp(key,"sba_modord") ) {
         // value=[vlist,ordspec,weight]
         u = BDY((LIST)value);
         pltovl((LIST)ARG0(u),&oldvl);
         nd_sba_modord = create_comp_sig_spec(vl,oldvl,(Obj)ARG1(u),argc(u)==3?ARG2(u):0);
       } else if ( !strcmp(key,"sba_gbinput") ) {
         nd_sba_inputisgb = value?1:0;
         if ( nd_sba_inputisgb != 0 ) {
           // value=[vlist,ordspec,weight]
           u = BDY((LIST)value);
           pltovl((LIST)ARG0(u),&oldvl);
           nd_sba_modord = create_comp_sig_spec(vl,oldvl,(Obj)ARG1(u),argc(u)==3?ARG2(u):0);
         }
       } else if ( !strcmp(key,"sba_redundant_check") ) {
         nd_sba_redundant_check = value?1:0;
       } else if ( !strcmp(key,"top") ) {
         nd_top = value?1:0;
     }      }
     }    }
     if ( nd_sba_syz ) nd_sba_dontsort = 1;
 }  }
   
 ND mdptond(DP d);  ND mdptond(DP d);
Line 8657  ND nd_mul_nm(int mod,NM m0,ND p)
Line 10715  ND nd_mul_nm(int mod,NM m0,ND p)
 ND nd_mul_nm_lf(NM m0,ND p)  ND nd_mul_nm_lf(NM m0,ND p)
 {  {
   UINT *d0;    UINT *d0;
   Z c0,c1,c;    Z c0,c1;
   NM tm,mr,mr0;    NM tm,mr,mr0;
   ND r;    ND r;
   
Line 8667  ND nd_mul_nm_lf(NM m0,ND p)
Line 10725  ND nd_mul_nm_lf(NM m0,ND p)
   mr0 = 0;    mr0 = 0;
   for ( tm = BDY(p); tm; tm = NEXT(tm) ) {    for ( tm = BDY(p); tm; tm = NEXT(tm) ) {
     NEXTNM(mr0,mr);      NEXTNM(mr0,mr);
     c = CZ(tm); mullf(c0,CZ(tm),&c1); CZ(mr) = c1;      mullf(c0,CZ(tm),&c1); CZ(mr) = c1;
     ndl_add(d0,DL(tm),DL(mr));      ndl_add(d0,DL(tm),DL(mr));
   }    }
   NEXT(mr) = 0;    NEXT(mr) = 0;
Line 8756  ND *btog_lf(NODE ti,ND **p,int nb)
Line 10814  ND *btog_lf(NODE ti,ND **p,int nb)
 ND btog_one(NODE ti,ND *p,int nb,int mod)  ND btog_one(NODE ti,ND *p,int nb,int mod)
 {  {
   PGeoBucket r;    PGeoBucket r;
   int i,ci,j;    int ci,j;
   NODE t,s;    NODE t,s;
   ND m,tp;    ND m,tp;
   ND pi,rd;    ND pi,rd;
Line 8770  ND btog_one(NODE ti,ND *p,int nb,int mod)
Line 10828  ND btog_one(NODE ti,ND *p,int nb,int mod)
     ptomp(mod,(P)HCZ(m),&c);      ptomp(mod,(P)HCZ(m),&c);
     if ( (ci = ((MQ)c)->cont) != 0 ) {      if ( (ci = ((MQ)c)->cont) != 0 ) {
       HCM(m) = ci;        HCM(m) = ci;
       pi = p[j=ZTOS((Q)ARG1(s))];        pi = p[j=(int)ZTOS((Q)ARG1(s))];
     if ( !pi ) {      if ( !pi ) {
       pi = nd_load_mod(j);        pi = nd_load_mod(j);
       tp = nd_mul_nm(mod,BDY(m),pi);        tp = nd_mul_nm(mod,BDY(m),pi);
Line 8809  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
Line 10867  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
   if ( mod == -2 )    if ( mod == -2 )
     return nd_btog_lf(f,v,ord,tlist,rp);      return nd_btog_lf(f,v,ord,tlist,rp);
   
   parse_nd_option(current_option);  
   get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
     parse_nd_option(vv,current_option);
   for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );    for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
   switch ( ord->id ) {    switch ( ord->id ) {
     case 1:      case 1:
Line 8832  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
Line 10890  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
   ind = BDY((LIST)ARG4(BDY(tlist)));    ind = BDY((LIST)ARG4(BDY(tlist)));
   perm = BDY((LIST)BDY(permtrace)); trace =NEXT(permtrace);    perm = BDY((LIST)BDY(permtrace)); trace =NEXT(permtrace);
   for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {    for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {
     j = ZTOS((Q)BDY(BDY((LIST)BDY(t))));      j = (int)ZTOS((Q)BDY(BDY((LIST)BDY(t))));
   if ( j > i ) i = j;    if ( j > i ) i = j;
   }    }
   n = i+1;    n = i+1;
Line 8840  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
Line 10898  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
   p = (ND **)MALLOC(n*sizeof(ND *));    p = (ND **)MALLOC(n*sizeof(ND *));
   for ( t = perm, i = 0; t; t = NEXT(t), i++ ) {    for ( t = perm, i = 0; t; t = NEXT(t), i++ ) {
     pi = BDY((LIST)BDY(t));      pi = BDY((LIST)BDY(t));
     pi0 = ZTOS((Q)ARG0(pi)); pi1 = ZTOS((Q)ARG1(pi));      pi0 = (int)ZTOS((Q)ARG0(pi)); pi1 = (int)ZTOS((Q)ARG1(pi));
     p[pi0] = c = (ND *)MALLOC(nb*sizeof(ND));      p[pi0] = c = (ND *)MALLOC(nb*sizeof(ND));
     ptomp(mod,(P)ARG2(pi),&inv);      ptomp(mod,(P)ARG2(pi),&inv);
     ((MQ)inv)->cont = invm(((MQ)inv)->cont,mod);      ((MQ)inv)->cont = invm(((MQ)inv)->cont,mod);
Line 8851  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
Line 10909  MAT nd_btog(LIST f,LIST v,int mod,struct order_spec *o
   for ( t = trace,i=0; t; t = NEXT(t), i++ ) {    for ( t = trace,i=0; t; t = NEXT(t), i++ ) {
     printf("%d ",i); fflush(stdout);      printf("%d ",i); fflush(stdout);
     ti = BDY((LIST)BDY(t));      ti = BDY((LIST)BDY(t));
     p[j=ZTOS((Q)ARG0(ti))] = btog(BDY((LIST)ARG1(ti)),p,nb,mod);      p[j=(int)ZTOS((Q)ARG0(ti))] = btog(BDY((LIST)ARG1(ti)),p,nb,mod);
   }    }
   for ( t = intred, i=0; t; t = NEXT(t), i++ ) {    for ( t = intred, i=0; t; t = NEXT(t), i++ ) {
     printf("%d ",i); fflush(stdout);      printf("%d ",i); fflush(stdout);
     ti = BDY((LIST)BDY(t));      ti = BDY((LIST)BDY(t));
     p[j=ZTOS((Q)ARG0(ti))] = btog(BDY((LIST)ARG1(ti)),p,nb,mod);      p[j=(int)ZTOS((Q)ARG0(ti))] = btog(BDY((LIST)ARG1(ti)),p,nb,mod);
   }    }
   m = length(ind);    m = length(ind);
   MKMAT(mat,nb,m);    MKMAT(mat,nb,m);
Line 8878  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
Line 10936  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
   LM lm;    LM lm;
   Z lf,inv;    Z lf,inv;
   
   parse_nd_option(current_option);  
   get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
     parse_nd_option(vv,current_option);
   for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );    for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
   switch ( ord->id ) {    switch ( ord->id ) {
     case 1:      case 1:
Line 8901  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
Line 10959  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
   ind = BDY((LIST)ARG4(BDY(tlist)));    ind = BDY((LIST)ARG4(BDY(tlist)));
   perm = BDY((LIST)BDY(permtrace)); trace =NEXT(permtrace);    perm = BDY((LIST)BDY(permtrace)); trace =NEXT(permtrace);
   for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {    for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {
     j = ZTOS((Q)BDY(BDY((LIST)BDY(t))));      j = (int)ZTOS((Q)BDY(BDY((LIST)BDY(t))));
   if ( j > i ) i = j;    if ( j > i ) i = j;
   }    }
   n = i+1;    n = i+1;
Line 8909  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
Line 10967  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
   p = (ND **)MALLOC(n*sizeof(ND *));    p = (ND **)MALLOC(n*sizeof(ND *));
   for ( t = perm, i = 0; t; t = NEXT(t), i++ ) {    for ( t = perm, i = 0; t; t = NEXT(t), i++ ) {
     pi = BDY((LIST)BDY(t));      pi = BDY((LIST)BDY(t));
     pi0 = ZTOS((Q)ARG0(pi)); pi1 = ZTOS((Q)ARG1(pi));      pi0 = (int)ZTOS((Q)ARG0(pi)); pi1 = (int)ZTOS((Q)ARG1(pi));
     p[pi0] = c = (ND *)MALLOC(nb*sizeof(ND));      p[pi0] = c = (ND *)MALLOC(nb*sizeof(ND));
     simp_ff((Obj)ARG2(pi),(Obj *)&lm); lmtolf(lm,&lf); invz(lf,current_mod_lf,&inv);      simp_ff((Obj)ARG2(pi),(Obj *)&lm); lmtolf(lm,&lf); invz(lf,current_mod_lf,&inv);
     u = ptond(CO,vv,(P)ONE);      u = ptond(CO,vv,(P)ONE);
Line 8919  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
Line 10977  MAT nd_btog_lf(LIST f,LIST v,struct order_spec *ord,LI
   for ( t = trace,i=0; t; t = NEXT(t), i++ ) {    for ( t = trace,i=0; t; t = NEXT(t), i++ ) {
     printf("%d ",i); fflush(stdout);      printf("%d ",i); fflush(stdout);
     ti = BDY((LIST)BDY(t));      ti = BDY((LIST)BDY(t));
     p[j=ZTOS((Q)ARG0(ti))] = btog_lf(BDY((LIST)ARG1(ti)),p,nb);      p[j=(int)ZTOS((Q)ARG0(ti))] = btog_lf(BDY((LIST)ARG1(ti)),p,nb);
   }    }
   for ( t = intred, i=0; t; t = NEXT(t), i++ ) {    for ( t = intred, i=0; t; t = NEXT(t), i++ ) {
     printf("%d ",i); fflush(stdout);      printf("%d ",i); fflush(stdout);
     ti = BDY((LIST)BDY(t));      ti = BDY((LIST)BDY(t));
     p[j=ZTOS((Q)ARG0(ti))] = btog_lf(BDY((LIST)ARG1(ti)),p,nb);      p[j=(int)ZTOS((Q)ARG0(ti))] = btog_lf(BDY((LIST)ARG1(ti)),p,nb);
   }    }
   m = length(ind);    m = length(ind);
   MKMAT(mat,nb,m);    MKMAT(mat,nb,m);
Line 8941  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 10999  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
   VL fv,tv,vv;    VL fv,tv,vv;
   NODE permtrace,perm,trace,intred,ind,t,pi,ti;    NODE permtrace,perm,trace,intred,ind,t,pi,ti;
   ND *p;    ND *p;
   ND *c;  
   ND u;    ND u;
   P inv;    P inv;
   VECT vect;    VECT vect;
Line 8949  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 11006  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
   if ( mod == -2 )    if ( mod == -2 )
     error("nd_btog_one : not implemented yet for a large finite field");      error("nd_btog_one : not implemented yet for a large finite field");
   
   parse_nd_option(current_option);  
   get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);    get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
     parse_nd_option(vv,current_option);
   for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );    for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
   switch ( ord->id ) {    switch ( ord->id ) {
     case 1:      case 1:
Line 8972  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 11029  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
   ind = BDY((LIST)ARG4(BDY(tlist)));    ind = BDY((LIST)ARG4(BDY(tlist)));
   perm = BDY((LIST)BDY(permtrace)); trace =NEXT(permtrace);    perm = BDY((LIST)BDY(permtrace)); trace =NEXT(permtrace);
   for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {    for ( i = length(perm)-1, t = trace; t; t = NEXT(t) ) {
     j = ZTOS((Q)BDY(BDY((LIST)BDY(t))));      j = (int)ZTOS((Q)BDY(BDY((LIST)BDY(t))));
   if ( j > i ) i = j;    if ( j > i ) i = j;
   }    }
   n = i+1;    n = i+1;
Line 8980  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 11037  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
   p = (ND *)MALLOC(n*sizeof(ND *));    p = (ND *)MALLOC(n*sizeof(ND *));
   for ( t = perm, i = 0; t; t = NEXT(t), i++ ) {    for ( t = perm, i = 0; t; t = NEXT(t), i++ ) {
     pi = BDY((LIST)BDY(t));      pi = BDY((LIST)BDY(t));
   pi0 = ZTOS((Q)ARG0(pi)); pi1 = ZTOS((Q)ARG1(pi));    pi0 = (int)ZTOS((Q)ARG0(pi)); pi1 = (int)ZTOS((Q)ARG1(pi));
   if ( pi1 == pos ) {    if ( pi1 == pos ) {
     ptomp(mod,(P)ARG2(pi),&inv);      ptomp(mod,(P)ARG2(pi),&inv);
     ((MQ)inv)->cont = invm(((MQ)inv)->cont,mod);      ((MQ)inv)->cont = invm(((MQ)inv)->cont,mod);
Line 8992  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 11049  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
   for ( t = trace,i=0; t; t = NEXT(t), i++ ) {    for ( t = trace,i=0; t; t = NEXT(t), i++ ) {
   printf("%d ",i); fflush(stdout);    printf("%d ",i); fflush(stdout);
     ti = BDY((LIST)BDY(t));      ti = BDY((LIST)BDY(t));
     p[j=ZTOS((Q)ARG0(ti))] = btog_one(BDY((LIST)ARG1(ti)),p,nb,mod);      p[j=(int)ZTOS((Q)ARG0(ti))] = btog_one(BDY((LIST)ARG1(ti)),p,nb,mod);
     if ( Demand ) {      if ( Demand ) {
         nd_save_mod(p[j],j); nd_free(p[j]); p[j] = 0;          nd_save_mod(p[j],j); nd_free(p[j]); p[j] = 0;
   }    }
Line 9000  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 11057  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
   for ( t = intred, i=0; t; t = NEXT(t), i++ ) {    for ( t = intred, i=0; t; t = NEXT(t), i++ ) {
   printf("%d ",i); fflush(stdout);    printf("%d ",i); fflush(stdout);
     ti = BDY((LIST)BDY(t));      ti = BDY((LIST)BDY(t));
     p[j=ZTOS((Q)ARG0(ti))] = btog_one(BDY((LIST)ARG1(ti)),p,nb,mod);      p[j=(int)ZTOS((Q)ARG0(ti))] = btog_one(BDY((LIST)ARG1(ti)),p,nb,mod);
     if ( Demand ) {      if ( Demand ) {
         nd_save_mod(p[j],j); nd_free(p[j]); p[j] = 0;          nd_save_mod(p[j],j); nd_free(p[j]); p[j] = 0;
   }    }
Line 9010  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
Line 11067  VECT nd_btog_one(LIST f,LIST v,int mod,struct order_sp
   for ( j = 0, t = ind; j < m; j++, t = NEXT(t) ) {    for ( j = 0, t = ind; j < m; j++, t = NEXT(t) ) {
   u = p[ZTOS((Q)BDY(t))];    u = p[ZTOS((Q)BDY(t))];
   if ( !u ) {    if ( !u ) {
     u = nd_load_mod(ZTOS((Q)BDY(t)));      u = nd_load_mod((int)ZTOS((Q)BDY(t)));
     BDY(vect)[j] = ndtodp(mod,u);      BDY(vect)[j] = ndtodp(mod,u);
     nd_free(u);      nd_free(u);
   } else    } else
Line 9039  void ndv_print_lf(NDV p)
Line 11096  void ndv_print_lf(NDV p)
   
 void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp)  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp)
 {  {
     VL tv,fv,vv,vc,av;      VL tv,fv,vv;
     NODE fd,fd0,in0,in,r,r0,t,s,cand,alist;      NODE fd,fd0,in0,in,r,t,s,cand;
     int m,nocheck,nvar,mindex,e,max;      int m,nvar,mindex,e,max;
     NDV c;      NDV c;
     NMV a;      NMV a;
     P p,zp;  
     Q dmy;  
     EPOS oepos;      EPOS oepos;
     int obpe,oadv,wmax,i,len,cbpe,ishomo,nalg,mrank,trank,ompos;      int obpe,oadv,wmax,i,len,cbpe,ishomo,mrank,trank,ompos;
     Alg alpha,dp;  
     P poly;  
     LIST f1,f2,zpl;  
     Obj obj;  
     NumberField nf;  
     struct order_spec *ord1;      struct order_spec *ord1;
     struct oEGT eg_check,eg0,eg1;      struct oEGT eg_check,eg0,eg1;
     NODE tr,tl1,tl2,tl3,tl4;  
     LIST l1,l2,l3,l4,l5;  
     int *perm;      int *perm;
     int j,ret;      int ret;
     NODE retn;      NODE retn;
     Q jq,bpe;  
   
     nd_module = 0;      nd_module = 0;
     parse_nd_option(current_option);  
     get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);      get_vars((Obj)f,&fv); pltovl(v,&vv); vlminus(fv,vv,&nd_vc);
       parse_nd_option(vv,current_option);
     if ( nd_vc )      if ( nd_vc )
       error("nd_f4_lf_trace : computation over a rational function field is not implemented");        error("nd_f4_lf_trace : computation over a rational function field is not implemented");
     for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );      for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ );
Line 9080  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
Line 11127  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
     nd_ntrans = nvar;      nd_ntrans = nvar;
     nd_nalg = 0;      nd_nalg = 0;
   
     nocheck = 0;  //    nocheck = 0;
     mindex = 0;      mindex = 0;
   
     /* do not use on-demand load/save */      /* do not use on-demand load/save */
Line 9127  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
Line 11174  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
         for ( t = in0, wmax = max; t; t = NEXT(t) ) {          for ( t = in0, wmax = max; t; t = NEXT(t) ) {
             c = (NDV)BDY(t); len = LEN(c);              c = (NDV)BDY(t); len = LEN(c);
             for ( a = BDY(c), i = 0; i < len; i++, NMV_ADV(a) )              for ( a = BDY(c), i = 0; i < len; i++, NMV_ADV(a) )
                 wmax = MAX(TD(DL(a)),wmax);                  wmax = MAX((int)TD(DL(a)),wmax);
         }          }
         homogenize_order(ord,nvar,&ord1);          homogenize_order(ord,nvar,&ord1);
         nd_init_ord(ord1);          nd_init_ord(ord1);
Line 9135  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
Line 11182  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
         for ( t = fd0; t; t = NEXT(t) )          for ( t = fd0; t; t = NEXT(t) )
             ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos,ompos);              ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos,ompos);
     }      }
     if ( MaxDeg > 0 ) nocheck = 1;  //    if ( MaxDeg > 0 ) nocheck = 1;
     ret = ndv_setup(-2,m,fd0,nd_gbblock?1:0,0);      ret = ndv_setup(-2,m,fd0,nd_gbblock?1:0,0,0);
     if ( ret )      if ( ret )
       cand = nd_f4_lf_trace_main(m,&perm);        cand = nd_f4_lf_trace_main(m,&perm);
     if ( !ret || !cand ) {      if ( !ret || !cand ) {
Line 9155  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
Line 11202  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
     if ( (ret = ndv_check_membership(-2,in0,obpe,oadv,oepos,cand)) != 0 ) {      if ( (ret = ndv_check_membership(-2,in0,obpe,oadv,oepos,cand)) != 0 ) {
       /* gbcheck : cand is a GB of Id(cand) ? */        /* gbcheck : cand is a GB of Id(cand) ? */
       retn = nd_f4(-2,0,0);        retn = nd_f4(-2,0,0);
         if ( !retn ) {
           /* failure */
           *rp = 0; return;
         }
     }      }
   if ( !retn ) {  
       /* failure */  
       *rp = 0; return;  
     }  
     get_eg(&eg1); init_eg(&eg_check); add_eg(&eg_check,&eg0,&eg1);      get_eg(&eg1); init_eg(&eg_check); add_eg(&eg_check,&eg0,&eg1);
     if ( DP_Print )      if ( DP_Print )
         fprintf(asir_out,"check=%.3fsec\n",eg_check.exectime);          fprintf(asir_out,"check=%.3fsec\n",eg_check.exectime);
Line 9175  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
Line 11222  void nd_f4_lf_trace(LIST f,LIST v,int trace,int homo,s
   
 NODE nd_f4_lf_trace_main(int m,int **indp)  NODE nd_f4_lf_trace_main(int m,int **indp)
 {  {
     int i,nh,stat,index;      int i,nh,stat;
     NODE r,rm,g;      NODE r,rm,g;
     ND_pairs d,l,l0,t;      ND_pairs d,l,l0,t;
     ND spol,red;      NDV nfqv,nfv;
     NDV nf,redv,nfqv,nfv;      NODE rp0,nflist,nflist_lf;
     NM s0,s;      UINT *s0vect;
     NODE rp0,srp0,nflist,nflist_lf;      int col;
     int nsp,nred,col,rank,len,k,j,a;  
     UINT c;  
     UINT **spmat;  
     UINT *s0vect,*svect,*p,*v;  
     int *colstat;  
     IndArray *imat;  
     int *rhead;  
     int spcol,sprow;  
     int sugar;      int sugar;
     PGeoBucket bucket;      PGeoBucket bucket;
     struct oEGT eg0,eg1,eg_f4;      struct oEGT eg0,eg1,eg_f4;
Line 9252  NODE nd_f4_lf_trace_main(int m,int **indp)
Line 11291  NODE nd_f4_lf_trace_main(int m,int **indp)
             if ( DL_COMPARE(HDL(nfv),HDL(nfqv)) ) return 0;              if ( DL_COMPARE(HDL(nfv),HDL(nfqv)) ) return 0;
             ndv_removecont(m,nfv);              ndv_removecont(m,nfv);
             ndv_removecont(-2,nfqv);              ndv_removecont(-2,nfqv);
             nh = ndv_newps(-2,nfv,nfqv,1);              nh = ndv_newps(-2,nfv,nfqv);
             d = update_pairs(d,g,nh,0);              d = update_pairs(d,g,nh,0);
             g = update_base(g,nh);              g = update_base(g,nh);
         }          }
Line 9292  NDV vect64_to_ndv(mp_limb_t *vect,int spcol,int col,in
Line 11331  NDV vect64_to_ndv(mp_limb_t *vect,int spcol,int col,in
     }      }
 }  }
   
   NDV vect64_to_ndv_s(mp_limb_t *vect,int col,UINT *s0vect)
   {
       int j,k,len;
       UINT *p;
       UINT c;
       NDV r;
       NMV mr0,mr;
   
       for ( j = 0, len = 0; j < col; j++ ) if ( vect[j] ) len++;
       if ( !len ) return 0;
       else {
           mr0 = (NMV)MALLOC_ATOMIC_IGNORE_OFF_PAGE(nmv_adv*len);
           mr = mr0;
           p = s0vect;
           for ( j = k = 0; j < col; j++, p += nd_wpd )
             if ( (c = (UINT)vect[k++]) != 0 ) {
               ndl_copy(p,DL(mr)); CM(mr) = c; NMV_ADV(mr);
             }
           MKNDV(nd_nvar,mr0,len,r);
           return r;
       }
   }
   
 int nd_to_vect64(int mod,UINT *s0,int n,ND d,mp_limb_t *r)  int nd_to_vect64(int mod,UINT *s0,int n,ND d,mp_limb_t *r)
 {  {
     NM m;      NM m;
Line 9321  int nd_to_vect64(int mod,UINT *s0,int n,ND d,mp_limb_t
Line 11383  int nd_to_vect64(int mod,UINT *s0,int n,ND d,mp_limb_t
   
 #define MOD128(a,c,m) ((a)=(((c)!=0||((a)>=(m)))?(((((U128)(c))<<64)+(a))%(m)):(a)))  #define MOD128(a,c,m) ((a)=(((c)!=0||((a)>=(m)))?(((((U128)(c))<<64)+(a))%(m)):(a)))
   
 int ndv_reduce_vect64(int m,mp_limb_t *svect,mp_limb_t *cvect,int col,IndArray *imat,NM_ind_pair *rp0,int nred)  int ndv_reduce_vect64(int m,mp_limb_t *svect,mp_limb_t *cvect,int col,IndArray *imat,NM_ind_pair *rp0,int nred,SIG sig)
 {  {
     int i,j,k,len,pos,prev;      int i,j,k,len,pos,prev;
     mp_limb_t a,c,c1,c2;      mp_limb_t a,c,c1,c2;
Line 9342  int ndv_reduce_vect64(int m,mp_limb_t *svect,mp_limb_t
Line 11404  int ndv_reduce_vect64(int m,mp_limb_t *svect,mp_limb_t
         a = svect[k]; c = cvect[k];          a = svect[k]; c = cvect[k];
         MOD128(a,c,m);          MOD128(a,c,m);
         svect[k] = a; cvect[k] = 0;          svect[k] = a; cvect[k] = 0;
         if ( (c = svect[k]) != 0 ) {          if ( (c = svect[k]) != 0 && (sig == 0 || comp_sig(sig,rp0[i]->sig) > 0 ) ) {
             Nf4_red++;              Nf4_red++;
             maxrs = MAX(maxrs,rp0[i]->sugar);              maxrs = MAX(maxrs,rp0[i]->sugar);
             c = m-c; redv = nd_ps[rp0[i]->index];              c = m-c; redv = nd_ps[rp0[i]->index];
Line 9416  NODE nd_f4_red_mod64_main(int m,ND_pairs sp0,int nsp,U
Line 11478  NODE nd_f4_red_mod64_main(int m,ND_pairs sp0,int nsp,U
         nd_sp(m,0,sp,&spol);          nd_sp(m,0,sp,&spol);
         if ( !spol ) continue;          if ( !spol ) continue;
         nd_to_vect64(m,s0vect,col,spol,svect);          nd_to_vect64(m,s0vect,col,spol,svect);
         maxrs = ndv_reduce_vect64(m,svect,cvect,col,imat,rvect,nred);          maxrs = ndv_reduce_vect64(m,svect,cvect,col,imat,rvect,nred,0);
         for ( i = 0; i < col; i++ ) if ( svect[i] ) break;          for ( i = 0; i < col; i++ ) if ( svect[i] ) break;
         if ( i < col ) {          if ( i < col ) {
             spmat[sprow] = v = (mp_limb_t *)MALLOC_ATOMIC(spcol*sizeof(mp_limb_t));              spmat[sprow] = v = (mp_limb_t *)MALLOC_ATOMIC(spcol*sizeof(mp_limb_t));
Line 9546  int nd_gauss_elim_mod64(mp_limb_t **mat,int *sugar,ND_
Line 11608  int nd_gauss_elim_mod64(mp_limb_t **mat,int *sugar,ND_
   GCFREE(cmat);    GCFREE(cmat);
   return rank;    return rank;
 }  }
   
   int nd_gauss_elim_mod64_s(mp_limb_t **mat,int *sugar,ND_pairs *spactive,int row,int col,int md,int *colstat,SIG *sig)
   {
     int i,j,k,l,rank,s,imin;
     mp_limb_t inv;
     mp_limb_t a;
     UINT c;
     mp_limb_t *t,*pivot,*pk;
     UINT *ck;
     UINT **cmat;
     UINT *ct;
     ND_pairs pair;
     SIG sg;
     int *used;
   
     used = (int *)MALLOC(row*sizeof(int));
     cmat = (UINT **)MALLOC(row*sizeof(UINT *));
     for ( i = 0; i < row; i++ ) {
       cmat[i] = MALLOC_ATOMIC(col*sizeof(UINT));
       bzero(cmat[i],col*sizeof(UINT));
     }
   
     for ( j = 0; j < col; j++ ) {
       for ( i = 0; i < row; i++ ) {
         a = mat[i][j]; c = cmat[i][j];
         MOD128(a,c,md);
         mat[i][j] = a; cmat[i][j] = 0;
       }
       for ( i = 0; i < row; i++ )
         if ( !used[i] && mat[i][j] ) break;
       if ( i == row ) {
         colstat[j] = 0;
         continue;
       } else {
         colstat[j] = 1;
         used[i] = 1;
       }
       /* column j is normalized */
       s = sugar[i];
       inv = invm((UINT)mat[i][j],md);
       /* normalize pivot row */
       for ( k = j, pk = mat[i]+j, ck = cmat[i]+j; k < col; k++, pk++, ck++ ) {
         a = *pk; c = *ck; MOD128(a,c,md); *pk = (a*inv)%md; *ck = 0;
       }
       for ( k = i+1; k < row; k++ ) {
         if ( (a = mat[k][j]) != 0 ) {
           sugar[k] = MAX(sugar[k],s);
           red_by_vect64(md,mat[k]+j,cmat[k]+j,mat[i]+j,(int)(md-a),col-j);
           Nf4_red++;
         }
       }
     }
     rank = 0;
     for ( i = 0; i < row; i++ ) {
       for ( j = 0; j < col; j++ )
         if ( mat[i][j] ) break;
       if ( j == col ) sugar[i] = -1;
       else rank++;
     }
     for ( i = 0; i < row; i++ ) GCFREE(cmat[i]);
     GCFREE(cmat);
     return rank;
   }
   
   NODE nd_f4_red_mod64_main_s(int m,ND_pairs sp0,int nsp,UINT *s0vect,int col,
           NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred,NODE *syzlistp)
   {
       int spcol,sprow,a;
       int i,j,k,l,rank;
       NODE r0,r;
       ND_pairs sp;
       ND spol;
       mp_limb_t **spmat;
       mp_limb_t *svect,*cvect;
       mp_limb_t *v;
       int *colstat;
       struct oEGT eg0,eg1,eg2,eg_f4,eg_f4_1,eg_f4_2;
       int maxrs;
       int *spsugar;
       ND_pairs *spactive;
       SIG *spsig;
   
       get_eg(&eg0);
       /* elimination (1st step) */
       spmat = (mp_limb_t **)MALLOC(nsp*sizeof(mp_limb_t *));
       cvect = (mp_limb_t *)MALLOC(col*sizeof(mp_limb_t));
       spsugar = (int *)MALLOC(nsp*sizeof(int));
       spsig = (SIG *)MALLOC(nsp*sizeof(SIG));
       for ( a = sprow = 0, sp = sp0; a < nsp; a++, sp = NEXT(sp) ) {
           nd_sp(m,0,sp,&spol);
           if ( !spol ) {
             syzlistp[sp->sig->pos] = insert_sig(syzlistp[sp->sig->pos],sp->sig);
             continue;
           }
           svect = (mp_limb_t *)MALLOC(col*sizeof(mp_limb_t));
           nd_to_vect64(m,s0vect,col,spol,svect);
           maxrs = ndv_reduce_vect64(m,svect,cvect,col,imat,rvect,nred,spol->sig);
           for ( i = 0; i < col; i++ ) if ( svect[i] ) break;
           if ( i < col ) {
               spmat[sprow] = svect;
               spsugar[sprow] = MAX(maxrs,SG(spol));
               spsig[sprow] = sp->sig;
               sprow++;
           } else {
             syzlistp[sp->sig->pos] = insert_sig(syzlistp[sp->sig->pos],sp->sig);
           }
           nd_free(spol);
       }
       get_eg(&eg1); init_eg(&eg_f4_1); add_eg(&eg_f4_1,&eg0,&eg1); add_eg(&f4_elim1,&eg0,&eg1);
       if ( DP_Print ) {
           fprintf(asir_out,"elim1=%.3fsec,",eg_f4_1.exectime);
           fflush(asir_out);
       }
       /* free index arrays */
       for ( i = 0; i < nred; i++ ) GCFREE(imat[i]->index.c);
   
       /* elimination (2nd step) */
       colstat = (int *)MALLOC(col*sizeof(int));
       rank = nd_gauss_elim_mod64_s(spmat,spsugar,0,sprow,col,m,colstat,spsig);
       r0 = 0;
       for ( i = 0; i < sprow; i++ ) {
           if ( spsugar[i] >= 0 ) {
             NEXTNODE(r0,r);
             BDY(r) = vect64_to_ndv_s(spmat[i],col,s0vect);
             SG((NDV)BDY(r)) = spsugar[i];
             ((NDV)BDY(r))->sig = spsig[i];
           } else
             syzlistp[spsig[i]->pos] = insert_sig(syzlistp[spsig[i]->pos],spsig[i]);
           GCFREE(spmat[i]);
       }
       if ( r0 ) NEXT(r) = 0;
       get_eg(&eg2); init_eg(&eg_f4_2); add_eg(&eg_f4_2,&eg1,&eg2); add_eg(&f4_elim2,&eg1,&eg2);
       init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg2);
       if ( DP_Print ) {
           fprintf(asir_out,"elim2=%.3fsec,",eg_f4_2.exectime);
           fprintf(asir_out,"nsp=%d,nred=%d,spmat=(%d,%d),rank=%d ",
               nsp,nred,sprow,col,rank);
           fprintf(asir_out,"%.3fsec,",eg_f4.exectime);
       }
       return r0;
   }
 #endif  #endif
   
   NODE nd_f4_red_s(int m,ND_pairs sp0,int trace,UINT *s0vect,int col,NODE rp0,NODE *syzlistp)
   {
     IndArray *imat;
     int nsp,nred,i,start;
     int *rhead;
     NODE r0,rp;
     ND_pairs sp;
     NM_ind_pair *rvect;
     struct oEGT eg0,eg1,eg_conv;
   
     for ( sp = sp0, nsp = 0; sp; sp = NEXT(sp), nsp++ );
     nred = length(rp0);
     imat = (IndArray *)MALLOC(nred*sizeof(IndArray));
     rhead = (int *)MALLOC(col*sizeof(int));
     for ( i = 0; i < col; i++ ) rhead[i] = 0;
   
     /* construction of index arrays */
     get_eg(&eg0);
     if ( DP_Print ) {
       fprintf(asir_out,"%dx%d,",nsp+nred,col);
       fflush(asir_out);
     }
     rvect = (NM_ind_pair *)MALLOC(nred*sizeof(NM_ind_pair));
     for ( start = 0, rp = rp0, i = 0; rp; i++, rp = NEXT(rp) ) {
       rvect[i] = (NM_ind_pair)BDY(rp);
       imat[i] = nm_ind_pair_to_vect_compress(trace,s0vect,col,rvect[i],start);
       rhead[imat[i]->head] = 1;
       start = imat[i]->head;
     }
     get_eg(&eg1); init_eg(&eg_conv); add_eg(&eg_conv,&eg0,&eg1); add_eg(&f4_conv,&eg0,&eg1);
     if ( DP_Print ) {
       fprintf(asir_out,"conv=%.3fsec,",eg_conv.exectime);
       fflush(asir_out);
     }
     if ( m > 0 )
   #if SIZEOF_LONG==8
       r0 = nd_f4_red_mod64_main_s(m,sp0,nsp,s0vect,col,rvect,rhead,imat,nred,syzlistp);
   #else
       r0 = nd_f4_red_main_s(m,sp0,nsp,s0vect,col,rvect,rhead,imat,nred,syzlistp);
   #endif
     else {
   //    r0 = nd_f4_red_q_main_s(sp0,nsp,trace,s0vect,col,rvect,rhead,imat,nred);
       r0 = 0;
       error("nd_f4_red_q_main_s : not implemented yet");
     }
     return r0;
   }
   
   INLINE int ndl_find_reducer_minsig(UINT *dg)
   {
     RHist r;
     int i,imin;
     SIG t;
     static int wpd,nvar;
     static SIG quo,quomin;
     static UINT *tmp;
   
     if ( !quo || nvar != nd_nvar ) { NEWSIG(quo); NEWSIG(quomin); }
     if ( wpd != nd_wpd ) {
       wpd = nd_wpd;
       tmp = (UINT *)MALLOC(wpd*sizeof(UINT));
     }
   #if 0
     d = ndl_hash_value(dg);
     for ( r = nd_red[d], k = 0; r; r = NEXT(r), k++ ) {
       if ( ndl_equal(dg,DL(r)) ) {
         return r->index;
       }
     }
   #endif
     imin = -1;
     for ( i = 0; i < nd_psn; i++ ) {
       r = nd_psh[i];
       if ( ndl_reducible(dg,DL(r)) ) {
         ndl_sub(dg,DL(r),tmp);
         _ndltodl(tmp,DL(quo));
         _addtodl(nd_nvar,DL(nd_psh[i]->sig),DL(quo));
         quo->pos = nd_psh[i]->sig->pos;
         _adddl(nd_nvar,DL(quo),nd_sba_hm[quo->pos],DL2(quo));
         if ( imin < 0 || comp_sig(quomin,quo) > 0 ) {
           t = quo; quo = quomin; quomin = t;
           imin = i;
         }
       }
     }
     if ( imin == -1 ) return nd_psn;
     else {
   #if 0
       nd_append_red(dg,i);
   #endif
       return imin;
     }
   }
   
   int nd_symbolic_preproc_s(PGeoBucket bucket,int trace,UINT **s0vect,NODE *r)
   {
     NODE rp0,rp;
     NM mul,head,s0,s;
     int index,col,i,sugar;
     RHist h;
     UINT *s0v,*p;
     NM_ind_pair pair;
     ND red;
     NDV *ps;
     SIG sig;
   
     s0 = 0; rp0 = 0; col = 0;
     if ( nd_demand )
       ps = trace?nd_ps_trace_sym:nd_ps_sym;
     else
       ps = trace?nd_ps_trace:nd_ps;
     while ( 1 ) {
       head = remove_head_pbucket_symbolic(bucket);
       if ( !head ) break;
       if ( !s0 ) s0 = head;
       else NEXT(s) = head;
       s = head;
       index = ndl_find_reducer_minsig(DL(head));
       if ( index >= 0 && index < nd_psn ) {
         h = nd_psh[index];
         NEWNM(mul);
         ndl_sub(DL(head),DL(h),DL(mul));
         if ( ndl_check_bound2(index,DL(mul)) )
           return 0;
         sugar = TD(DL(mul))+SG(ps[index]);
         NEWSIG(sig);
         _ndltodl(DL(mul),DL(sig));
         _addtodl(nd_nvar,DL(nd_psh[index]->sig),DL(sig));
         sig->pos = nd_psh[index]->sig->pos;
         _adddl(nd_nvar,DL(sig),nd_sba_hm[sig->pos],DL2(sig));
         MKNM_ind_pair(pair,mul,index,sugar,sig);
         red = ndv_mul_nm_symbolic(mul,ps[index]);
         add_pbucket_symbolic(bucket,nd_remove_head(red));
         NEXTNODE(rp0,rp); BDY(rp) = (pointer)pair;
       }
       col++;
     }
     if ( rp0 ) NEXT(rp) = 0;
     NEXT(s) = 0;
     s0v = (UINT *)MALLOC_ATOMIC(col*nd_wpd*sizeof(UINT));
     for ( i = 0, p = s0v, s = s0; i < col;
       i++, p += nd_wpd, s = NEXT(s) ) ndl_copy(DL(s),p);
     *s0vect = s0v;
     *r = rp0;
   
     return col;
   }
   
   NODE nd_sba_f4(int m,int **indp)
   {
     int i,nh,stat,f4red,f4step;
     int col,psugar,ms;
     NODE r,g,rp0,nflist;
     ND_pairs d,l,t,l1;
     ND h,nf;
     NDV nfv;
   //  union oNDC hc;
     UINT *s0vect;
     PGeoBucket bucket;
     NODE *syzlist;
     SIG sig;
     struct oEGT eg1;
     struct oEGT eg2,eg_remove,eg_nf,eg_nfzero;
   
     Nf4_red=0;
     d = 0;
     syzlist = (NODE *)MALLOC(nd_psn*sizeof(NODE));
     for ( i = 0; i < nd_psn; i++ ) {
       d = update_pairs_s(d,i,syzlist);
     }
     nd_nbase = nd_psn;
     f4red = 1;
     psugar = 0;
     f4step = 0;
     while ( d ) {
       for ( t = d, ms = SG(d); t; t = NEXT(t) )
         if ( SG(t) < ms ) ms = SG(t);
       if ( ms == psugar && f4step >= nd_sba_f4step ) {
   again:
         l = d; d = d->next;
   #if 0
         if ( small_lcm(l) ) {
           if ( DP_Print ) fprintf(asir_out,"M");
           continue;
         }
         sig = l->sig;
         stat = nd_sp(m,0,l,&h);
   #else
         l1 = find_smallest_lcm(l);
         if ( l1 == 0 ) {
           if ( DP_Print ) fprintf(asir_out,"M");
           continue;
         }
         sig = l1->sig;
         stat = nd_sp(m,0,l1,&h);
   #endif
         if ( !stat ) {
           NEXT(l) = d; d = l;
           d = nd_reconstruct(0,d);
           goto again;
         }
     get_eg(&eg1);
     #if USE_GEOBUCKET
         stat = m?nd_nf_pbucket_s(m,h,nd_ps,!nd_top&&!Top,&nf):nd_nf_s(m,0,h,nd_ps,!nd_top&&!Top,&nf);
     #else
         stat = nd_nf_s(m,0,h,nd_ps,!nd_top&&!Top,&nf);
     #endif
     get_eg(&eg2);
         if ( !stat ) {
           NEXT(l) = d; d = l;
           d = nd_reconstruct(0,d);
           goto again;
         } else if ( stat == -1 ) {
           if ( DP_Print ) { printf("S"); fflush(stdout); }
           FREENDP(l);
         } else if ( nf ) {
           if ( DP_Print ) { printf("+"); fflush(stdout); }
           add_eg(&eg_nf,&eg1,&eg2);
   //        hc = HCU(nf);
           nd_removecont(m,nf);
           nfv = ndtondv(m,nf); nd_free(nf);
           nh = ndv_newps(m,nfv,0);
   
           d = update_pairs_s(d,nh,syzlist);
           nd_sba_pos[sig->pos] = append_one(nd_sba_pos[sig->pos],nh);
           FREENDP(l);
         } else {
           add_eg(&eg_nfzero,&eg1,&eg2);
          // syzygy
     get_eg(&eg1);
           d = remove_spair_s(d,sig);
     get_eg(&eg2); add_eg(&eg_remove,&eg1,&eg2);
           syzlist[sig->pos] = insert_sig(syzlist[sig->pos],sig);
           if ( DP_Print ) { printf("."); fflush(stdout); }
           FREENDP(l);
         }
       } else {
         if ( ms != psugar ) f4step = 1;
         else f4step++;
   again2:
         psugar = ms;
         l = nd_minsugarp_s(d,&d);
   //      sugar = nd_sugarweight?d->sugar2:SG(d);
         bucket = create_pbucket();
         stat = nd_sp_f4(m,0,l,bucket);
         if ( !stat ) {
           for ( t = l; NEXT(t); t = NEXT(t) );
           NEXT(t) = d; d = l;
           d = nd_reconstruct(0,d);
           goto again2;
         }
         if ( bucket->m < 0 ) continue;
         col = nd_symbolic_preproc_s(bucket,0,&s0vect,&rp0);
         if ( !col ) {
           for ( t = l; NEXT(t); t = NEXT(t) )
             ;
           NEXT(t) = d; d = l;
           d = nd_reconstruct(0,d);
           goto again2;
         }
         if ( DP_Print ) fprintf(asir_out,"\nsugar=%d,",psugar);
         nflist = nd_f4_red_s(m,l,0,s0vect,col,rp0,syzlist);
         /* adding new bases */
         for ( r = nflist; r; r = NEXT(r) ) {
           nfv = (NDV)BDY(r);
           if ( nd_f4_td ) SG(nfv) = nd_tdeg(nfv);
           ndv_removecont(m,nfv);
           nh = ndv_newps(m,nfv,0);
           d = update_pairs_s(d,nh,syzlist);
           nd_sba_pos[nfv->sig->pos] = append_one(nd_sba_pos[nfv->sig->pos],nh);
         }
         for ( i = 0; i < nd_nbase; i++ )
           for ( r = syzlist[i]; r; r = NEXT(r) )
               d = remove_spair_s(d,(SIG)BDY(r));
         d = remove_large_lcm(d);
         if ( DP_Print ) {
           fprintf(asir_out,"f4red=%d,gblen=%d",f4red,nd_psn); fflush(asir_out);
         }
         f4red++;
       }
     }
     if ( DP_Print ) {
       fprintf(asir_out,"\nnumber of red=%d,",Nf4_red);
     }
     g = conv_ilist_s(nd_demand,0,indp);
     return g;
   }
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.66
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>