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

Diff for /OpenXM_contrib2/asir2000/engine/nd.c between version 1.70 and 1.71

version 1.70, 2003/09/15 10:51:45 version 1.71, 2003/09/17 07:16:53
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.69 2003/09/15 09:49:44 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.70 2003/09/15 10:51:45 noro Exp $ */
   
 #include "ca.h"  #include "ca.h"
 #include "inline.h"  #include "inline.h"
Line 143  static int nmv_adv;
Line 143  static int nmv_adv;
 static int nd_dcomp;  static int nd_dcomp;
   
 extern VL CO;  extern VL CO;
 extern int Top,Reverse,dp_nelim,do_weyl;  extern int Top,Reverse,DP_Print,dp_nelim,do_weyl;
 extern int *current_weyl_weight_vector;  extern int *current_weyl_weight_vector;
   
 /* fundamental macros */  /* fundamental macros */
Line 224  if(!(r)){NEWND_pairs(r);(c)=(r);}else{NEWND_pairs(NEXT
Line 224  if(!(r)){NEWND_pairs(r);(c)=(r);}else{NEWND_pairs(NEXT
 #define NMV_OPREV(m) (m = (NMV)(((char *)m)-oadv))  #define NMV_OPREV(m) (m = (NMV)(((char *)m)-oadv))
   
 /* external functions */  /* external functions */
   #if 1
 void GC_gcollect();  void GC_gcollect();
   #endif
 NODE append_one(NODE,int);  NODE append_one(NODE,int);
   
 /* manipulation of coefficients */  /* manipulation of coefficients */
Line 331  ND nd_copy(ND p);
Line 333  ND nd_copy(ND p);
 ND nd_merge(ND p1,ND p2);  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);
 ND nd_add_q(ND p1,ND p2);  ND nd_add_q(ND p1,ND p2);
   ND nd_add_sf(ND p1,ND p2);
 INLINE int nd_length(ND p);  INLINE int nd_length(ND p);
   
 /* NDV functions */  /* NDV functions */
Line 361  int nd_to_vect(int mod,UINT *s0,int n,ND d,UINT *r);
Line 364  int nd_to_vect(int mod,UINT *s0,int n,ND d,UINT *r);
   
 void nd_free_private_storage()  void nd_free_private_storage()
 {  {
         _nd_free_list = 0;  
         _nm_free_list = 0;          _nm_free_list = 0;
         _ndp_free_list = 0;          _ndp_free_list = 0;
         bzero(nd_red,sizeof(REDTAB_LEN*sizeof(RHist)));  #if 0
         GC_gcollect();          GC_gcollect();
   #endif
 }  }
   
 void _NM_alloc()  void _NM_alloc()
Line 1141  ND nd_add(int mod,ND p1,ND p2)
Line 1144  ND nd_add(int mod,ND p1,ND p2)
   
         if ( !p1 ) return p2;          if ( !p1 ) return p2;
         else if ( !p2 ) return p1;          else if ( !p2 ) return p1;
           else if ( mod == -1 ) return nd_add_sf(p1,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;
Line 1228  ND nd_add_q(ND p1,ND p2)
Line 1232  ND nd_add_q(ND p1,ND p2)
         }          }
 }  }
   
   ND nd_add_sf(ND p1,ND p2)
   {
           int n,c,can;
           ND r;
           NM m1,m2,mr0,mr,s;
           int t;
   
           if ( !p1 ) return p2;
           else if ( !p2 ) return p1;
           else {
                   can = 0;
                   for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2), mr0 = 0; m1 && m2; ) {
                           c = DL_COMPARE(DL(m1),DL(m2));
                           switch ( c ) {
                                   case 0:
                                           t = _addsf(CM(m1),CM(m2));
                                           s = m1; m1 = NEXT(m1);
                                           if ( t ) {
                                                   can++; NEXTNM2(mr0,mr,s); CM(mr) = (t);
                                           } else {
                                                   can += 2; FREENM(s);
                                           }
                                           s = m2; m2 = NEXT(m2); FREENM(s);
                                           break;
                                   case 1:
                                           s = m1; m1 = NEXT(m1); NEXTNM2(mr0,mr,s);
                                           break;
                                   case -1:
                                           s = m2; m2 = NEXT(m2); NEXTNM2(mr0,mr,s);
                                           break;
                           }
                   }
                   if ( !mr0 )
                           if ( m1 ) mr0 = m1;
                           else if ( m2 ) mr0 = m2;
                           else return 0;
                   else if ( m1 ) NEXT(mr) = m1;
                   else if ( m2 ) NEXT(mr) = m2;
                   else NEXT(mr) = 0;
                   BDY(p1) = mr0;
                   SG(p1) = MAX(SG(p1),SG(p2));
                   LEN(p1) = LEN(p1)+LEN(p2)-can;
                   FREEND(p2);
                   return p1;
           }
   }
   
 /* ret=1 : success, ret=0 : overflow */  /* ret=1 : success, ret=0 : overflow */
 int nd_nf(int mod,ND g,NDV *ps,int full,NDC dn,ND *rp)  int nd_nf(int mod,ND g,NDV *ps,int full,NDC dn,ND *rp)
 {  {
Line 1266  int nd_nf(int mod,ND g,NDV *ps,int full,NDC dn,ND *rp)
Line 1317  int nd_nf(int mod,ND g,NDV *ps,int full,NDC dn,ND *rp)
                                 return 0;                                  return 0;
                         }                          }
                         p = ps[index];                          p = ps[index];
                         if ( mod ) {                          if ( mod == -1 )
                                   CM(mul) = _mulsf(_invsf(HCM(p)),_chsgnsf(HCM(g)));
                           else if ( mod ) {
                                 c1 = invm(HCM(p),mod); c2 = mod-HCM(g);                                  c1 = invm(HCM(p),mod); c2 = mod-HCM(g);
                                 DMAR(c1,c2,0,mod,c); CM(mul) = c;                                  DMAR(c1,c2,0,mod,c); CM(mul) = c;
                         } else {                          } else {
Line 1353  int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *rp
Line 1406  int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *rp
                                 return 0;                                  return 0;
                         }                          }
                         p = ps[index];                          p = ps[index];
                         if ( mod ) {                          if ( mod == -1 )
                                   CM(mul) = _mulsf(_invsf(HCM(p)),_chsgnsf(HCM(g)));
                           else if ( mod ) {
                                 c1 = invm(HCM(p),mod); c2 = mod-HCM(g);                                  c1 = invm(HCM(p),mod); c2 = mod-HCM(g);
                                 DMAR(c1,c2,0,mod,c); CM(mul) = c;                                  DMAR(c1,c2,0,mod,c); CM(mul) = c;
                         } else {                          } else {
Line 1430  again:
Line 1485  again:
                         nd_reconstruct(0,0,0);                          nd_reconstruct(0,0,0);
                         goto again;                          goto again;
                 } else if ( nf ) return 0;                  } else if ( nf ) return 0;
                 printf("."); fflush(stdout);                  if ( DP_Print ) { printf("."); fflush(stdout); }
         }          }
         printf("\n");          if ( DP_Print ) { printf("\n"); }
         /* gbcheck : cand is a GB of Id(cand) ? */          /* gbcheck : cand is a GB of Id(cand) ? */
         if ( !nd_gb(0,1) ) return 0;          if ( !nd_gb(0,1) ) return 0;
         /* XXX */          /* XXX */
Line 1594  int head_pbucket(int mod,PGeoBucket g)
Line 1649  int head_pbucket(int mod,PGeoBucket g)
                                         dj = HDL(gj);                                          dj = HDL(gj);
                                         sum = HCM(gj);                                          sum = HCM(gj);
                                 } else if ( c == 0 ) {                                  } else if ( c == 0 ) {
                                         sum = sum+HCM(gi)-mod;                                          if ( mod == -1 )
                                         if ( sum < 0 ) sum += mod;                                                  sum = _addsf(sum,HCM(gi));
                                           else {
                                                   sum = sum+HCM(gi)-mod;
                                                   if ( sum < 0 ) sum += mod;
                                           }
                                         g->body[i] = nd_remove_head(gi);                                          g->body[i] = nd_remove_head(gi);
                                 }                                  }
                         }                          }
Line 1686  again:
Line 1745  again:
                 l = nd_minp(d,&d);                  l = nd_minp(d,&d);
                 if ( SG(l) != sugar ) {                  if ( SG(l) != sugar ) {
                         sugar = SG(l);                          sugar = SG(l);
                         fprintf(asir_out,"%d",sugar);                          if ( DP_Print ) fprintf(asir_out,"%d",sugar);
                 }                  }
                 stat = nd_sp(m,0,l,&h);                  stat = nd_sp(m,0,l,&h);
                 if ( !stat ) {                  if ( !stat ) {
Line 1705  again:
Line 1764  again:
                         goto again;                          goto again;
                 } else if ( nf ) {                  } else if ( nf ) {
                         if ( checkonly ) return 0;                          if ( checkonly ) return 0;
                         printf("+"); fflush(stdout);                          if ( DP_Print ) { printf("+"); fflush(stdout); }
                         nd_removecont(m,nf);                          nd_removecont(m,nf);
                         nfv = ndtondv(m,nf); nd_free(nf);                          nfv = ndtondv(m,nf); nd_free(nf);
                         nh = ndv_newps(nfv,0);                          nh = ndv_newps(nfv,0);
Line 1713  again:
Line 1772  again:
                         g = update_base(g,nh);                          g = update_base(g,nh);
                         FREENDP(l);                          FREENDP(l);
                 } else {                  } else {
                         printf("."); fflush(stdout);                          if ( DP_Print ) { printf("."); fflush(stdout); }
                         FREENDP(l);                          FREENDP(l);
                 }                  }
         }          }
Line 1742  again:
Line 1801  again:
                 l = nd_minp(d,&d);                  l = nd_minp(d,&d);
                 if ( SG(l) != sugar ) {                  if ( SG(l) != sugar ) {
                         sugar = SG(l);                          sugar = SG(l);
                         fprintf(asir_out,"%d",sugar);                          if ( DP_Print ) fprintf(asir_out,"%d",sugar);
                 }                  }
                 stat = nd_sp(m,0,l,&h);                  stat = nd_sp(m,0,l,&h);
                 if ( !stat ) {                  if ( !stat ) {
Line 1769  again:
Line 1828  again:
                                 /* m|HC(nfq) => failure */                                  /* m|HC(nfq) => failure */
                                 if ( !rem(NM(HCQ(nfq)),m) ) return 0;                                  if ( !rem(NM(HCQ(nfq)),m) ) return 0;
   
                                 printf("+"); fflush(stdout);                                  if ( DP_Print ) { printf("+"); fflush(stdout); }
                                 nd_removecont(m,nf); nfv = ndtondv(m,nf); nd_free(nf);                                  nd_removecont(m,nf); nfv = ndtondv(m,nf); nd_free(nf);
                                 nd_removecont(0,nfq); nfqv = ndtondv(0,nfq); nd_free(nfq);                                  nd_removecont(0,nfq); nfqv = ndtondv(0,nfq); nd_free(nfq);
                                 nh = ndv_newps(nfv,nfqv);                                  nh = ndv_newps(nfv,nfqv);
                                 d = update_pairs(d,g,nh);                                  d = update_pairs(d,g,nh);
                                 g = update_base(g,nh);                                  g = update_base(g,nh);
                         } else {                          } else {
                                 printf("*"); fflush(stdout);                                  if ( DP_Print ) { printf("*"); fflush(stdout); }
                         }                          }
                 } else {                  } else {
                         printf("."); fflush(stdout);                          if ( DP_Print ) { printf("."); fflush(stdout); }
                 }                  }
                 FREENDP(l);                  FREENDP(l);
         }          }
Line 1823  NODE ndv_reduceall(int m,NODE f)
Line 1882  NODE ndv_reduceall(int m,NODE f)
                 if ( !stat )                  if ( !stat )
                         nd_reconstruct(m,0,0);                          nd_reconstruct(m,0,0);
                 else {                  else {
                         printf("."); fflush(stdout);                          if ( DP_Print ) { printf("."); fflush(stdout); }
                         if ( !m ) { mulq(HCQ(head),dn.z,&q); HCQ(head) = q; }                          if ( !m ) { mulq(HCQ(head),dn.z,&q); HCQ(head) = q; }
                         nf = nd_add(m,head,nf);                          nf = nd_add(m,head,nf);
                         ndv_free(nd_ps[i]);                          ndv_free(nd_ps[i]);
Line 1833  NODE ndv_reduceall(int m,NODE f)
Line 1892  NODE ndv_reduceall(int m,NODE f)
                         i++;                          i++;
                 }                  }
         }          }
         printf("\n");          if ( DP_Print ) { printf("\n"); }
         for ( a0 = 0, i = 0; i < n; i++ ) {          for ( a0 = 0, i = 0; i < n; i++ ) {
                 NEXTNODE(a0,a);                  NEXTNODE(a0,a);
                 BDY(a) = (pointer)nd_ps[i];                  BDY(a) = (pointer)nd_ps[i];
Line 2156  void ndv_setup(int mod,int trace,NODE f)
Line 2215  void ndv_setup(int mod,int trace,NODE f)
   
         if ( !nd_red )          if ( !nd_red )
                 nd_red = (RHist *)MALLOC(REDTAB_LEN*sizeof(RHist));                  nd_red = (RHist *)MALLOC(REDTAB_LEN*sizeof(RHist));
         bzero(nd_red,REDTAB_LEN*sizeof(RHist));          for ( i = 0; i < REDTAB_LEN; i++ ) nd_red[i] = 0;
         nd_free_private_storage();  
         for ( i = 0; i < nd_psn; i++ ) {          for ( i = 0; i < nd_psn; i++ ) {
                 if ( trace ) {                  if ( trace ) {
                         a = nd_ps_trace[i] = ndv_dup(0,w[i]);                          a = nd_ps_trace[i] = ndv_dup(0,w[i]);
Line 2182  void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
Line 2240  void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
         int e,max,nvar;          int e,max,nvar;
         NDV b;          NDV b;
   
           nd_free_private_storage();
         get_vars((Obj)f,&fv); pltovl(v,&vv);          get_vars((Obj)f,&fv); pltovl(v,&vv);
         nvar = length(vv);          nvar = length(vv);
         nd_init_ord(ord);          nd_init_ord(ord);
Line 2226  void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru
Line 2285  void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru
         nocheck = 0;          nocheck = 0;
         mindex = 0;          mindex = 0;
   
           nd_free_private_storage();
         /* setup modulus */          /* setup modulus */
         if ( trace < 0 ) {          if ( trace < 0 ) {
                 trace = -trace;                  trace = -trace;
Line 2383  void nd_print(ND p)
Line 2443  void nd_print(ND p)
                 printf("0\n");                  printf("0\n");
         else {          else {
                 for ( m = BDY(p); m; m = NEXT(m) ) {                  for ( m = BDY(p); m; m = NEXT(m) ) {
                         printf("+%d*",CM(m));                          if ( CM(m) & 0x80000000 ) printf("+@_%d*",IFTOF(CM(m)));
                           else printf("+%d*",CM(m));
                         ndl_print(DL(m));                          ndl_print(DL(m));
                 }                  }
                 printf("\n");                  printf("\n");
Line 2424  void nd_removecont(int mod,ND p)
Line 2485  void nd_removecont(int mod,ND p)
         struct oVECT v;          struct oVECT v;
         N q,r;          N q,r;
   
         if ( mod ) nd_mul_c(mod,p,invm(HCM(p),mod));          if ( mod == -1 ) nd_mul_c(mod,p,_invsf(HCM(p)));
           else if ( mod ) nd_mul_c(mod,p,invm(HCM(p),mod));
         else {          else {
                 for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ );                  for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ );
                 w = (Q *)ALLOCA(n*sizeof(Q));                  w = (Q *)ALLOCA(n*sizeof(Q));
Line 2470  void ndv_removecont(int mod,NDV p)
Line 2532  void ndv_removecont(int mod,NDV p)
         Q dvr,t;          Q dvr,t;
         NMV m;          NMV m;
   
         if ( mod )          if ( mod == -1 )
                   ndv_mul_c(mod,p,_invsf(HCM(p)));
           else if ( mod )
                 ndv_mul_c(mod,p,invm(HCM(p),mod));                  ndv_mul_c(mod,p,invm(HCM(p),mod));
         else {          else {
                 len = p->len;                  len = p->len;
Line 2571  void nd_mul_c(int mod,ND p,int mul)
Line 2635  void nd_mul_c(int mod,ND p,int mul)
         int c,c1;          int c,c1;
   
         if ( !p ) return;          if ( !p ) return;
         for ( m = BDY(p); m; m = NEXT(m) ) {          if ( mod == -1 )
                 c1 = CM(m);                  for ( m = BDY(p); m; m = NEXT(m) )
                 DMAR(c1,mul,0,mod,c);                          CM(m) = _mulsf(CM(m),mul);
                 CM(m) = c;          else
         }                  for ( m = BDY(p); m; m = NEXT(m) ) {
                           c1 = CM(m); DMAR(c1,mul,0,mod,c); CM(m) = c;
                   }
 }  }
   
 void nd_mul_c_q(ND p,Q mul)  void nd_mul_c_q(ND p,Q mul)
Line 2721  void nd_setup_parameters(int nvar,int max) {
Line 2787  void nd_setup_parameters(int nvar,int max) {
 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
 {  {
         int i,obpe,oadv,h;          int i,obpe,oadv,h;
         NM prev_nm_free_list;          static NM prev_nm_free_list;
           static ND_pairs prev_ndp_free_list;
         RHist mr0,mr;          RHist mr0,mr;
         RHist r;          RHist r;
         RHist *old_red;          RHist *old_red;
         ND_pairs s0,s,t,prev_ndp_free_list;          ND_pairs s0,s,t;
         EPOS oepos;          EPOS oepos;
   
         obpe = nd_bpe;          obpe = nd_bpe;
Line 2785  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
Line 2852  ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d)
         if ( s0 ) NEXT(s) = 0;          if ( s0 ) NEXT(s) = 0;
         prev_nm_free_list = 0;          prev_nm_free_list = 0;
         prev_ndp_free_list = 0;          prev_ndp_free_list = 0;
   #if 0
         GC_gcollect();          GC_gcollect();
   #endif
         return s0;          return s0;
 }  }
   
Line 2860  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
Line 2929  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
         if ( ndl_check_bound2(p->i1,DL(m)) )          if ( ndl_check_bound2(p->i1,DL(m)) )
                 return 0;                  return 0;
         t1 = ndv_mul_nm(mod,m,p1);          t1 = ndv_mul_nm(mod,m,p1);
         if ( mod ) CM(m) = mod-HCM(p1);          if ( mod == -1 ) CM(m) = _chsgnsf(HCM(p1));
           else if ( mod ) CM(m) = mod-HCM(p1);
         else chsgnq(HCQ(p1),&CQ(m));          else chsgnq(HCQ(p1),&CQ(m));
         ndl_sub(lcm,HDL(p2),DL(m));          ndl_sub(lcm,HDL(p2),DL(m));
         if ( ndl_check_bound2(p->i2,DL(m)) ) {          if ( ndl_check_bound2(p->i2,DL(m)) ) {
Line 2880  void ndv_mul_c(int mod,NDV p,int mul)
Line 2950  void ndv_mul_c(int mod,NDV p,int mul)
   
         if ( !p ) return;          if ( !p ) return;
         len = LEN(p);          len = LEN(p);
         for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) {          if ( mod == -1 )
                 c1 = CM(m); DMAR(c1,mul,0,mod,c); CM(m) = c;                  for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) )
         }                          CM(m) = _mulsf(CM(m),mul);
           else
                   for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) {
                           c1 = CM(m); DMAR(c1,mul,0,mod,c); CM(m) = c;
                   }
 }  }
   
 void ndv_mul_c_q(NDV p,Q mul)  void ndv_mul_c_q(NDV p,Q mul)
Line 3072  ND ndv_mul_nm(int mod,NM m0,NDV p)
Line 3146  ND ndv_mul_nm(int mod,NM m0,NDV p)
   
         if ( !p ) return 0;          if ( !p ) return 0;
         else if ( do_weyl )          else if ( do_weyl )
                 return weyl_ndv_mul_nm(mod,m0,p);                  if ( mod == -1 )
                           error("ndv_mul_nm : not implemented (weyl)");
                   else
                           return weyl_ndv_mul_nm(mod,m0,p);
         else {          else {
                 n = NV(p); m = BDY(p);                  n = NV(p); m = BDY(p);
                 d = DL(m0);                  d = DL(m0);
                 len = LEN(p);                  len = LEN(p);
                 mr0 = 0;                  mr0 = 0;
                 td = TD(d);                  td = TD(d);
                 if ( mod ) {                  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) ) {
                                 NEXTNM(mr0,mr);                                  NEXTNM(mr0,mr);
                                   CM(mr) = _mulsf(CM(m),c);
                                   ndl_add(DL(m),d,DL(mr));
                           }
                   } else if ( mod ) {
                           c = CM(m0);
                           for ( i = 0; i < len; i++, NMV_ADV(m) ) {
                                   NEXTNM(mr0,mr);
                                 c1 = CM(m);                                  c1 = CM(m);
                                 DMAR(c1,c,0,mod,c2);                                  DMAR(c1,c,0,mod,c2);
                                 CM(mr) = c2;                                  CM(mr) = c2;
Line 3187  void ndv_mod(int mod,NDV p)
Line 3271  void ndv_mod(int mod,NDV p)
         NMV t,d;          NMV t,d;
         int r;          int r;
         int i,len,dlen;          int i,len,dlen;
           Obj gfs;
   
         if ( !p ) return;          if ( !p ) return;
         len = LEN(p);          len = LEN(p);
         dlen = 0;          dlen = 0;
         for ( t = d = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) {          if ( mod == -1 )
                 r = rem(NM(CQ(t)),mod);                  for ( t = d = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) {
                 if ( r ) {                          simp_ff((Obj)CP(t),&gfs);
                         if ( SGN(CQ(t)) < 0 )                          r = FTOIF(CONT((GFS)gfs));
                                 r = mod-r;  
                         CM(d) = r;                          CM(d) = r;
                         ndl_copy(DL(t),DL(d));                          ndl_copy(DL(t),DL(d));
                         NMV_ADV(d);                          NMV_ADV(d);
                         dlen++;                          dlen++;
                 }                  }
         }          else
                   for ( t = d = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) {
                           r = rem(NM(CQ(t)),mod);
                           if ( r ) {
                                   if ( SGN(CQ(t)) < 0 )
                                           r = mod-r;
                                   CM(d) = r;
                                   ndl_copy(DL(t),DL(d));
                                   NMV_ADV(d);
                                   dlen++;
                           }
                   }
         LEN(p) = dlen;          LEN(p) = dlen;
 }  }
   
Line 3276  P ndvtop(int mod,VL vl,VL dvl,NDV p)
Line 3371  P ndvtop(int mod,VL vl,VL dvl,NDV p)
         P c;          P c;
         UINT *d;          UINT *d;
         P s,r,u,t,w;          P s,r,u,t,w;
           GFS gfs;
   
         if ( !p ) return 0;          if ( !p ) return 0;
         else {          else {
Line 3283  P ndvtop(int mod,VL vl,VL dvl,NDV p)
Line 3379  P ndvtop(int mod,VL vl,VL dvl,NDV p)
                 n = NV(p);                  n = NV(p);
                 m = (NMV)(((char *)BDY(p))+nmv_adv*(len-1));                  m = (NMV)(((char *)BDY(p))+nmv_adv*(len-1));
                 for ( j = len-1, s = 0; j >= 0; j--, NMV_PREV(m) ) {                  for ( j = len-1, s = 0; j >= 0; j--, NMV_PREV(m) ) {
                         if ( mod ) {                          if ( mod == -1 ) {
                                   e = IFTOF(CM(m)); MKGFS(e,gfs); c = (P)gfs;
                           } else if ( mod ) {
                                 STOQ(CM(m),q); c = (P)q;                                  STOQ(CM(m),q); c = (P)q;
                         } else                          } else
                                 c = CP(m);                                  c = CP(m);
Line 3347  void ndv_print(NDV p)
Line 3445  void ndv_print(NDV p)
         else {          else {
                 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) ) {
                         printf("+%d*",CM(m));                          if ( CM(m) & 0x80000000 ) printf("+@_%d*",IFTOF(CM(m)));
                           else printf("+%d*",CM(m));
                         ndl_print(DL(m));                          ndl_print(DL(m));
                 }                  }
                 printf("\n");                  printf("\n");
Line 3888  NODE nd_f4(int m)
Line 3987  NODE nd_f4(int m)
                 rank = generic_gauss_elim_mod(spmat,sprow,spcol,m,colstat);                  rank = generic_gauss_elim_mod(spmat,sprow,spcol,m,colstat);
   
                 get_eg(&eg1); init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg1);                  get_eg(&eg1); init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg1);
                 fprintf(asir_out,"sugar=%d,nsp=%d,nred=%d,spmat=(%d,%d),rank=%d  ",                  if ( DP_Print ) {
                         sugar,nsp,nred,sprow,spcol,rank);                          fprintf(asir_out,"sugar=%d,nsp=%d,nred=%d,spmat=(%d,%d),rank=%d  ",
                 fprintf(asir_out,"%fsec\n",eg_f4.exectime+eg_f4.gctime);                                  sugar,nsp,nred,sprow,spcol,rank);
                           fprintf(asir_out,"%fsec\n",eg_f4.exectime+eg_f4.gctime);
                   }
   
                 /* adding new bases */                  /* adding new bases */
                 for ( i = 0; i < rank; i++ ) {                  for ( i = 0; i < rank; i++ ) {
Line 3989  NODE nd_f4(int m)
Line 4090  NODE nd_f4(int m)
                 rank = generic_gauss_elim_mod(spmat,sprow,spcol,m,colstat);                  rank = generic_gauss_elim_mod(spmat,sprow,spcol,m,colstat);
   
                 get_eg(&eg1); init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg1);                  get_eg(&eg1); init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg1);
                 fprintf(asir_out,"sugar=%d,nsp=%d,nred=%d,spmat=(%d,%d),rank=%d  ",                  if ( DP_Print ) {
                         sugar,nsp,nred,sprow,spcol,rank);                          fprintf(asir_out,"sugar=%d,nsp=%d,nred=%d,spmat=(%d,%d),rank=%d  ",
                 fprintf(asir_out,"%fsec\n",eg_f4.exectime+eg_f4.gctime);                                  sugar,nsp,nred,sprow,spcol,rank);
                           fprintf(asir_out,"%fsec\n",eg_f4.exectime+eg_f4.gctime);
                   }
   
                 /* adding new bases */                  /* adding new bases */
                 for ( i = 0; i < rank; i++ ) {                  for ( i = 0; i < rank; i++ ) {
Legend:
Removed from v.1.70  
changed lines
  Added in v.1.71
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>