[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.146 and 1.151

version 1.146, 2006/11/27 07:31:25 version 1.151, 2006/12/14 02:49:59
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.145 2006/08/26 05:38:06 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.150 2006/12/12 11:50:37 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 8  int nd_dcomp;
Line 8  int nd_dcomp;
 NM _nm_free_list;  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;
   
 static NODE nd_subst;  static NODE nd_subst;
 static VL nd_vc;  static VL nd_vc;
Line 52  void nd_det_reconstruct(NDV **dm,int n,int j,NDV d);
Line 53  void nd_det_reconstruct(NDV **dm,int n,int j,NDV d);
 int nd_monic(int m,ND *p);  int nd_monic(int m,ND *p);
 NDV plain_vect_to_ndv_q(Q *mat,int col,UINT *s0vect);  NDV plain_vect_to_ndv_q(Q *mat,int col,UINT *s0vect);
   
   extern int Denominator,DP_Multiple;
   
 void nd_free_private_storage()  void nd_free_private_storage()
 {  {
         _nm_free_list = 0;          _nm_free_list = 0;
Line 1316  int ndv_check_candidate(NODE input,int obpe,int oadv,E
Line 1319  int ndv_check_candidate(NODE input,int obpe,int oadv,E
         NODE t,s;          NODE t,s;
         union oNDC dn;          union oNDC dn;
   
         ndv_setup(0,0,cand,0);          ndv_setup(0,0,cand,0,1);
         n = length(cand);          n = length(cand);
   
         /* membercheck : list is a subset of Id(cand) ? */          /* membercheck : list is a subset of Id(cand) ? */
Line 1569  ND normalize_pbucket(int mod,PGeoBucket g)
Line 1572  ND normalize_pbucket(int mod,PGeoBucket g)
         return r;          return r;
 }  }
   
   #if 0
   void register_hcf(NDV p)
   {
           DCP dc,t;
           P hc,h;
           int c;
           NODE l,l1,prev;
   
           hc = p->body->c.p;
           if ( !nd_vc || NUM(hc) ) return;
           fctrp(nd_vc,hc,&dc);
           for ( t = dc; t; t = NEXT(t) ) {
                   h = t->c;
                   if ( NUM(h) ) continue;
                   for ( prev = 0, l = nd_hcf; l; prev = l, l = NEXT(l) ) {
                           c = compp(nd_vc,h,(P)BDY(l));
                           if ( c >= 0 ) break;
                   }
                   if ( !l || c > 0  ) {
                           MKNODE(l1,h,l);
                           if ( !prev )
                                   nd_hcf = l1;
                           else
                                   NEXT(prev) = l1;
                   }
           }
   }
   #else
   void register_hcf(NDV p)
   {
           DCP dc,t;
           P hc,h,q;
           Q dmy;
           int c;
           NODE l,l1,prev;
   
           hc = p->body->c.p;
           if ( NUM(hc) ) return;
           ptozp(hc,1,&dmy,&h);
   #if 1
           for ( l = nd_hcf; l; l = NEXT(l) ) {
                   while ( 1 ) {
                           if ( divtpz(nd_vc,h,(P)BDY(l),&q) ) h = q;
                           else break;
                   }
           }
           if ( NUM(h) ) return;
   #endif
           for ( prev = 0, l = nd_hcf; l; prev = l, l = NEXT(l) ) {
                   c = compp(nd_vc,h,(P)BDY(l));
                   if ( c >= 0 ) break;
           }
           if ( !l || c > 0  ) {
                   MKNODE(l1,h,l);
                   if ( !prev )
                           nd_hcf = l1;
                   else
                           NEXT(prev) = l1;
           }
   }
   #endif
   
 int do_diagonalize(int sugar,int m)  int do_diagonalize(int sugar,int m)
 {  {
         int i,nh,stat;          int i,nh,stat;
Line 1592  int do_diagonalize(int sugar,int m)
Line 1657  int do_diagonalize(int sugar,int m)
                 nfv = ndtondv(m,nf);                  nfv = ndtondv(m,nf);
                 nd_free(nf);                  nd_free(nf);
                 nd_bound[i] = ndv_compute_bound(nfv);                  nd_bound[i] = ndv_compute_bound(nfv);
                   if ( !m ) register_hcf(nfv);
                 if ( nd_demand ) {                  if ( nd_demand ) {
                         ndv_save(nfv,i);                          ndv_save(nfv,i);
                         ndv_free(nfv);                          ndv_free(nfv);
Line 1685  again:
Line 1751  again:
         else          else
                 for ( t = g; t; t = NEXT(t) )                  for ( t = g; t; t = NEXT(t) )
                         BDY(t) = (pointer)nd_ps[(int)BDY(t)];                          BDY(t) = (pointer)nd_ps[(int)BDY(t)];
           if ( !checkonly && DP_Print ) { printf("nd_gb done.\n"); fflush(stdout); }
         return g;          return g;
 }  }
   
Line 1721  int do_diagonalize_trace(int sugar,int m)
Line 1788  int do_diagonalize_trace(int sugar,int m)
                 nfv = ndtondv(0,nf);                  nfv = ndtondv(0,nf);
                 nd_free(nf);                  nd_free(nf);
                 nd_bound[i] = ndv_compute_bound(nfv);                  nd_bound[i] = ndv_compute_bound(nfv);
                   register_hcf(nfv);
                 if ( nd_demand ) {                  if ( nd_demand ) {
                         ndv_save(nfv,i);                          ndv_save(nfv,i);
                         ndv_free(nfv);                          ndv_free(nfv);
Line 1733  int do_diagonalize_trace(int sugar,int m)
Line 1801  int do_diagonalize_trace(int sugar,int m)
 static struct oEGT eg_invdalg;  static struct oEGT eg_invdalg;
 struct oEGT eg_le;  struct oEGT eg_le;
   
   void nd_subst_vector(VL vl,P p,NODE subst,P *r)
   {
           NODE tn;
           P p1;
   
           for ( tn = subst; tn; tn = NEXT(NEXT(tn)) ) {
                   substp(vl,p,BDY(tn),BDY(NEXT(tn)),&p1); p = p1;
           }
           *r = p;
   }
   
 NODE nd_gb_trace(int m,int ishomo)  NODE nd_gb_trace(int m,int ishomo)
 {  {
         int i,nh,sugar,stat;          int i,nh,sugar,stat;
Line 1742  NODE nd_gb_trace(int m,int ishomo)
Line 1821  NODE nd_gb_trace(int m,int ishomo)
         ND h,nf,nfq,s,head;          ND h,nf,nfq,s,head;
         NDV nfv,nfqv;          NDV nfv,nfqv;
         Q q,den,num;          Q q,den,num;
           P hc;
         union oNDC dn;          union oNDC dn;
         struct oEGT eg_monic,egm0,egm1;          struct oEGT eg_monic,egm0,egm1;
         int diag_count = 0;          int diag_count = 0;
Line 1761  again:
Line 1841  again:
                 if ( SG(l) != sugar ) {                  if ( SG(l) != sugar ) {
 #if 1  #if 1
                         if ( ishomo ) {                          if ( ishomo ) {
                                   if ( DP_Print > 2 ) fprintf(asir_out,"|");
                                 stat = do_diagonalize_trace(sugar,m);                                  stat = do_diagonalize_trace(sugar,m);
                                   if ( DP_Print > 2 ) fprintf(asir_out,"|");
                                 diag_count = 0;                                  diag_count = 0;
                                 if ( !stat ) {                                  if ( !stat ) {
                                         NEXT(l) = d; d = l;                                          NEXT(l) = d; d = l;
Line 1803  again:
Line 1885  again:
                         }                          }
                         if ( nfq ) {                          if ( nfq ) {
                                 /* m|HC(nfq) => failure */                                  /* m|HC(nfq) => failure */
                                 if ( !rem(NM(HCQ(nfq)),m) ) return 0;                                  if ( nd_vc ) {
                                           nd_subst_vector(nd_vc,HCP(nfq),nd_subst,&hc); q = (Q)hc;
                                   } else
                                           q = HCQ(nfq);
                                   if ( !rem(NM(q),m) ) return 0;
   
                                 if ( DP_Print ) { printf("+"); fflush(stdout); }                                  if ( DP_Print ) { printf("+"); fflush(stdout); }
                                 if ( nd_nalg ) {                                  if ( nd_nalg ) {
Line 1820  again:
Line 1906  again:
                                 nh = ndv_newps(0,nfv,nfqv);                                  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,"|");
                                         stat = do_diagonalize_trace(sugar,m);                                          stat = do_diagonalize_trace(sugar,m);
                                           if ( DP_Print > 2 ) fprintf(asir_out,"|");
                                         if ( !stat ) {                                          if ( !stat ) {
                                                 NEXT(l) = d; d = l;                                                  NEXT(l) = d; d = l;
                                                 d = nd_reconstruct(1,d);                                                  d = nd_reconstruct(1,d);
Line 1848  again:
Line 1936  again:
                 print_eg("invdalg",&eg_invdalg);                  print_eg("invdalg",&eg_invdalg);
                 print_eg("le",&eg_le);                  print_eg("le",&eg_le);
         }          }
           if ( DP_Print ) { printf("nd_gb_trace done.\n"); fflush(stdout); }
         return g;          return g;
 }  }
   
Line 1878  NODE ndv_reduceall(int m,NODE f)
Line 1967  NODE ndv_reduceall(int m,NODE f)
                 (int (*)(const void *,const void *))ndv_compare);                  (int (*)(const void *,const void *))ndv_compare);
         for ( t = f, i = 0; t; i++, t = NEXT(t) ) BDY(t) = (pointer)w[i];          for ( t = f, i = 0; t; i++, t = NEXT(t) ) BDY(t) = (pointer)w[i];
 #endif  #endif
         ndv_setup(m,0,f,0);          ndv_setup(m,0,f,0,1);
         for ( i = 0; i < n; ) {          for ( i = 0; i < n; ) {
                 g = ndvtond(m,nd_ps[i]);                  g = ndvtond(m,nd_ps[i]);
                 g = nd_separate_head(g,&head);                  g = nd_separate_head(g,&head);
Line 2191  int ndv_newps(int m,NDV a,NDV aq)
Line 2280  int ndv_newps(int m,NDV a,NDV aq)
         nd_ps[nd_psn] = a;          nd_ps[nd_psn] = a;
         if ( aq ) {          if ( aq ) {
                 nd_ps_trace[nd_psn] = aq;                  nd_ps_trace[nd_psn] = aq;
                   register_hcf(aq);
                 nd_bound[nd_psn] = ndv_compute_bound(aq);                  nd_bound[nd_psn] = ndv_compute_bound(aq);
                 SG(r) = SG(aq); ndl_copy(HDL(aq),DL(r));                  SG(r) = SG(aq); ndl_copy(HDL(aq),DL(r));
         } else {          } else {
                   if ( !m ) register_hcf(a);
                 nd_bound[nd_psn] = ndv_compute_bound(a);                  nd_bound[nd_psn] = ndv_compute_bound(a);
                 SG(r) = SG(a); ndl_copy(HDL(a),DL(r));                  SG(r) = SG(a); ndl_copy(HDL(a),DL(r));
         }          }
Line 2209  int ndv_newps(int m,NDV a,NDV aq)
Line 2300  int ndv_newps(int m,NDV a,NDV aq)
         return nd_psn++;          return nd_psn++;
 }  }
   
 void ndv_setup(int mod,int trace,NODE f,int dont_sort)  void ndv_setup(int mod,int trace,NODE f,int dont_sort,int dont_removecont)
 {  {
         int i,j,td,len,max;          int i,j,td,len,max;
         NODE s,s0,f0;          NODE s,s0,f0;
Line 2237  void ndv_setup(int mod,int trace,NODE f,int dont_sort)
Line 2328  void ndv_setup(int mod,int trace,NODE f,int dont_sort)
         nd_ps_trace = (NDV *)MALLOC(nd_pslen*sizeof(NDV));          nd_ps_trace = (NDV *)MALLOC(nd_pslen*sizeof(NDV));
         nd_psh = (RHist *)MALLOC(nd_pslen*sizeof(RHist));          nd_psh = (RHist *)MALLOC(nd_pslen*sizeof(RHist));
         nd_bound = (UINT **)MALLOC(nd_pslen*sizeof(UINT *));          nd_bound = (UINT **)MALLOC(nd_pslen*sizeof(UINT *));
           nd_hcf = 0;
   
         if ( trace && nd_vc )          if ( trace && nd_vc )
                 makesubst(nd_vc,&nd_subst);                  makesubst(nd_vc,&nd_subst);
Line 2249  void ndv_setup(int mod,int trace,NODE f,int dont_sort)
Line 2341  void ndv_setup(int mod,int trace,NODE f,int dont_sort)
         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]);
                         ndv_removecont(0,a);                          if ( !dont_removecont) ndv_removecont(0,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);
                         ndv_removecont(mod,am);                          ndv_removecont(mod,am);
                 } else {                  } else {
                         a = nd_ps[i] = ndv_dup(mod,w[i]);                          a = nd_ps[i] = ndv_dup(mod,w[i]);
                         ndv_removecont(mod,a);                          if ( mod || !dont_removecont ) ndv_removecont(mod,a);
                           if ( !mod ) register_hcf(a);
                 }                  }
                 NEWRHist(r); SG(r) = HTD(a); ndl_copy(HDL(a),DL(r));                  NEWRHist(r); SG(r) = HTD(a); ndl_copy(HDL(a),DL(r));
                 nd_bound[i] = ndv_compute_bound(a);                  nd_bound[i] = ndv_compute_bound(a);
Line 2395  void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
Line 2489  void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
         if ( !m && Demand ) nd_demand = 1;          if ( !m && Demand ) nd_demand = 1;
         else nd_demand = 0;          else nd_demand = 0;
   
           if ( DP_Multiple )
                   nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1);
 #if 0  #if 0
         ndv_alloc = 0;          ndv_alloc = 0;
 #endif  #endif
Line 2440  void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
Line 2536  void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe
                 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);          ndv_setup(m,0,fd0,0,0);
         x = f4?nd_f4(m):nd_gb(m,ishomo,0);          x = f4?nd_f4(m):nd_gb(m,ishomo,0);
         nd_demand = 0;          nd_demand = 0;
         x = ndv_reducebase(x);          x = ndv_reducebase(x);
Line 2515  void nd_gr_postproc(LIST f,LIST v,int m,struct order_s
Line 2611  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);          ndv_setup(m,0,fd0,0,1);
         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 2559  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 2655  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
         struct order_spec *ord1;          struct order_spec *ord1;
         struct oEGT eg_check,eg0,eg1;          struct oEGT eg_check,eg0,eg1;
   
           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);          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 ) {
Line 2630  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
Line 2729  void nd_gr_trace(LIST f,LIST v,int trace,int homo,int 
         while ( 1 ) {          while ( 1 ) {
                 if ( Demand )                  if ( Demand )
                         nd_demand = 1;                          nd_demand = 1;
                 ndv_setup(m,1,fd0,0);                  ndv_setup(m,1,fd0,0,0);
                 cand = f4?nd_f4_trace(m):nd_gb_trace(m,ishomo || homo);                  cand = f4?nd_f4_trace(m):nd_gb_trace(m,ishomo || homo);
                 if ( !cand ) {                  if ( !cand ) {
                         /* failure */                          /* failure */
Line 2809  void nd_removecont(int mod,ND p)
Line 2908  void nd_removecont(int mod,ND p)
                 v.len = n;                  v.len = n;
                 v.body = (pointer *)w;                  v.body = (pointer *)w;
                 for ( m = BDY(p), i = 0; i < n; m = NEXT(m), i++ ) w[i] = CQ(m);                  for ( m = BDY(p), i = 0; i < n; m = NEXT(m), i++ ) w[i] = CQ(m);
                 removecont_array((P *)w,n);                  removecont_array((P *)w,n,1);
                 for ( m = BDY(p), i = 0; i < n; m = NEXT(m), i++ ) CQ(m) = w[i];                  for ( m = BDY(p), i = 0; i < n; m = NEXT(m), i++ ) CQ(m) = w[i];
         }          }
 }  }
Line 2823  void nd_removecont2(ND p1,ND p2)
Line 2922  void nd_removecont2(ND p1,ND p2)
         struct oVECT v;          struct oVECT v;
         N q,r;          N q,r;
   
         if ( !p1 ) {  
                 nd_removecont(0,p2); return;  
         } else if ( !p2 ) {  
                 nd_removecont(0,p1); return;  
         }  
         n1 = nd_length(p1);          n1 = nd_length(p1);
         n2 = nd_length(p2);          n2 = nd_length(p2);
         n = n1+n2;          n = n1+n2;
         w = (Q *)ALLOCA(n*sizeof(Q));          w = (Q *)ALLOCA(n*sizeof(Q));
         v.len = n;          v.len = n;
         v.body = (pointer *)w;          v.body = (pointer *)w;
         for ( m = BDY(p1), i = 0; i < n1; m = NEXT(m), i++ ) w[i] = CQ(m);          i = 0;
         for ( m = BDY(p2); i < n; m = NEXT(m), i++ ) w[i] = CQ(m);          if ( p1 )
         removecont_array((P *)w,n);                  for ( m = BDY(p1); i < n1; m = NEXT(m), i++ ) w[i] = CQ(m);
         for ( m = BDY(p1), i = 0; i < n1; m = NEXT(m), i++ ) CQ(m) = w[i];          if ( p2 )
         for ( m = BDY(p2); i < n; m = NEXT(m), i++ ) CQ(m) = w[i];                  for ( m = BDY(p2); i < n; m = NEXT(m), i++ ) w[i] = CQ(m);
           removecont_array((P *)w,n,1);
           i = 0;
           if ( p1 )
                   for ( m = BDY(p1); i < n1; m = NEXT(m), i++ ) CQ(m) = w[i];
           if ( p2 )
                   for ( m = BDY(p2); i < n; m = NEXT(m), i++ ) CQ(m) = w[i];
 }  }
   
 void ndv_removecont(int mod,NDV p)  void ndv_removecont(int mod,NDV p)
Line 2863  void ndv_removecont(int mod,NDV p)
Line 2963  void ndv_removecont(int mod,NDV p)
                         all_p = all_p && !NUM(w[i]);                          all_p = all_p && !NUM(w[i]);
                 }                  }
                 if ( all_p ) {                  if ( all_p ) {
                         qltozl(c,len,&dvr); heu_nezgcdnpz(nd_vc,w,len,&g);                          qltozl(c,len,&dvr); nd_heu_nezgcdnpz(nd_vc,w,len,1,&g);
                         mulp(nd_vc,(P)dvr,g,&cont);                          mulp(nd_vc,(P)dvr,g,&cont);
                         for ( m = BDY(p), i = 0; i < len; NMV_ADV(m), i++ ) {                          for ( m = BDY(p), i = 0; i < len; NMV_ADV(m), i++ ) {
                                 divsp(nd_vc,CP(m),cont,&tp); CP(m) = tp;                                  divsp(nd_vc,CP(m),cont,&tp); CP(m) = tp;
Line 2924  void ndv_dehomogenize(NDV p,struct order_spec *ord)
Line 3024  void ndv_dehomogenize(NDV p,struct order_spec *ord)
         NV(p)--;          NV(p)--;
 }  }
   
 void removecont_array(P *p,int n)  void nd_heu_nezgcdnpz(VL vl,P *pl,int m,int full,P *pr)
 {  {
           int i;
           P *tpl,*tpl1;
           NODE l;
           P h,gcd,t;
   
           tpl = (P *)ALLOCA(m*sizeof(P));
           tpl1 = (P *)ALLOCA(m*sizeof(P));
           bcopy(pl,tpl,m*sizeof(P));
           gcd = (P)ONE;
           for ( l = nd_hcf; l; l = NEXT(l) ) {
                   h = (P)BDY(l);
                   while ( 1 ) {
                           for ( i = 0; i < m; i++ )
                                   if ( !divtpz(vl,tpl[i],h,&tpl1[i]) )
                                           break;
                           if ( i == m ) {
                                   bcopy(tpl1,tpl,m*sizeof(P));
                                   mulp(vl,gcd,h,&t); gcd = t;
                           } else
                                   break;
                   }
           }
           if ( DP_Print > 2 ){fprintf(asir_out,"[%d]",nmonop(gcd)); fflush(asir_out);}
           if ( full ) {
                   heu_nezgcdnpz(vl,tpl,m,&t);
                   mulp(vl,gcd,t,pr);
           } else
                   *pr = gcd;
   }
   
   void removecont_array(P *p,int n,int full)
   {
         int all_p,all_q,i;          int all_p,all_q,i;
         Q *c;          Q *c;
         P *w;          P *w;
Line 2942  void removecont_array(P *p,int n)
Line 3074  void removecont_array(P *p,int n)
                         ptozp(p[i],1,&c[i],&w[i]);                          ptozp(p[i],1,&c[i],&w[i]);
                 }                  }
                 removecont_array_q(c,n);                  removecont_array_q(c,n);
                 heu_nezgcdnpz(nd_vc,w,n,&t);                  nd_heu_nezgcdnpz(nd_vc,w,n,full,&t);
                 for ( i = 0; i < n; i++ ) {                  for ( i = 0; i < n; i++ ) {
                         divsp(nd_vc,w[i],t,&s); mulp(nd_vc,s,(P)c[i],&p[i]);                          divsp(nd_vc,w[i],t,&s); mulp(nd_vc,s,(P)c[i],&p[i]);
                 }                  }
Line 3298  ND nd_copy(ND p)
Line 3430  ND nd_copy(ND p)
   
 int nd_sp(int mod,int trace,ND_pairs p,ND *rp)  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
 {  {
         NM m;          NM m1,m2;
         NDV p1,p2;          NDV p1,p2;
         ND t1,t2;          ND t1,t2;
         UINT *lcm;          UINT *lcm;
           P gp,tp;
           Q g,t;
         int td;          int td;
   
         if ( !mod && nd_demand ) {          if ( !mod && nd_demand ) {
Line 3314  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
Line 3448  int nd_sp(int mod,int trace,ND_pairs p,ND *rp)
                 }                  }
         }          }
         lcm = LCM(p);          lcm = LCM(p);
         NEWNM(m);          NEWNM(m1); ndl_sub(lcm,HDL(p1),DL(m1));
         CQ(m) = HCQ(p2);          if ( ndl_check_bound2(p->i1,DL(m1)) ) {
         ndl_sub(lcm,HDL(p1),DL(m));                  FREENM(m1); return 0;
         if ( ndl_check_bound2(p->i1,DL(m)) )  
                 return 0;  
         t1 = ndv_mul_nm(mod,m,p1);  
         if ( mod == -1 ) CM(m) = _chsgnsf(HCM(p1));  
         else if ( mod ) CM(m) = mod-HCM(p1);  
         else chsgnp(HCP(p1),&CP(m));  
         ndl_sub(lcm,HDL(p2),DL(m));  
         if ( ndl_check_bound2(p->i2,DL(m)) ) {  
                 nd_free(t1);  
                 return 0;  
         }          }
         t2 = ndv_mul_nm(mod,m,p2);          NEWNM(m2); ndl_sub(lcm,HDL(p2),DL(m2));
           if ( ndl_check_bound2(p->i2,DL(m2)) ) {
                   FREENM(m1); FREENM(m2); return 0;
           }
   
           if ( mod == -1 ) {
                   CM(m1) = HCM(p2); CM(m2) = _chsgnsf(HCM(p1));
           } else if ( mod ) {
                   CM(m1) = HCM(p2); CM(m2) = mod-HCM(p1);
           } else if ( nd_vc ) {
                   ezgcdpz(nd_vc,HCP(p1),HCP(p2),&gp);
                   divsp(nd_vc,HCP(p2),gp,&CP(m1));
                   divsp(nd_vc,HCP(p1),gp,&tp); chsgnp(tp,&CP(m2));
           } else {
                   igcd_cofactor(HCQ(p1),HCQ(p2),&g,&t,&CQ(m1)); chsgnq(t,&CQ(m2));
           }
           t1 = ndv_mul_nm(mod,m1,p1); t2 = ndv_mul_nm(mod,m2,p2);
         *rp = nd_add(mod,t1,t2);          *rp = nd_add(mod,t1,t2);
         FREENM(m);          FREENM(m1); FREENM(m2);
         return 1;          return 1;
 }  }
   
Line 3710  void ndv_mod(int mod,NDV p)
Line 3850  void ndv_mod(int mod,NDV p)
         NMV t,d;          NMV t,d;
         int r,s,u;          int r,s,u;
         int i,len,dlen;          int i,len,dlen;
           P cp;
         Q c;          Q c;
         P cp,cp1;  
         NODE tn;  
         Obj gfs;          Obj gfs;
   
         if ( !p ) return;          if ( !p ) return;
Line 3730  void ndv_mod(int mod,NDV p)
Line 3869  void ndv_mod(int mod,NDV p)
         else          else
                 for ( t = d = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) {                  for ( t = d = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) {
                         if ( nd_vc ) {                          if ( nd_vc ) {
                                 for ( tn = nd_subst, cp = CP(t); tn; tn = NEXT(NEXT(tn)) ) {                                  nd_subst_vector(nd_vc,CP(t),nd_subst,&cp);
                                         substp(nd_vc,cp,BDY(tn),BDY(NEXT(tn)),&cp1); cp = cp1;  
                                 }  
                                 c = (Q)cp;                                  c = (Q)cp;
                         } else                          } else
                                 c = CQ(t);                                  c = CQ(t);
Line 4137  void nd_nf_p(P f,LIST g,LIST v,int m,struct order_spec
Line 4274  void nd_nf_p(P f,LIST g,LIST v,int m,struct order_spec
         if ( m ) ndv_mod(m,(NDV)BDY(in));          if ( m ) ndv_mod(m,(NDV)BDY(in));
         NEXT(in) = 0;          NEXT(in) = 0;
   
         /* dont sort */          /* dont sort, dont removecont */
         ndv_setup(m,0,in0,1);          ndv_setup(m,0,in0,1,1);
         nd_psn--;          nd_psn--;
         nd_scale=2;          nd_scale=2;
         while ( 1 ) {          while ( 1 ) {
Line 4336  int ndv_reduce_vect_q(Q *svect,int trace,int col,IndAr
Line 4473  int ndv_reduce_vect_q(Q *svect,int trace,int col,IndAr
                         if ( j == col ) break;                          if ( j == col ) break;
                         if ( hmag && ((double)p_mag((P)svect[j]) > hmag) ) {                          if ( hmag && ((double)p_mag((P)svect[j]) > hmag) ) {
                                 nz = compress_array(svect,cvect,col);                                  nz = compress_array(svect,cvect,col);
                                 removecont_array((P *)cvect,nz);                                  removecont_array((P *)cvect,nz,1);
                                 expand_array(svect,cvect,nz);                                  expand_array(svect,cvect,nz);
                                 hmag = ((double)p_mag((P)svect[j]))*nd_scale;                                  hmag = ((double)p_mag((P)svect[j]))*nd_scale;
                         }                          }
                 }                  }
         }          }
         nz = compress_array(svect,cvect,col);          nz = compress_array(svect,cvect,col);
         removecont_array((P *)cvect,nz);          removecont_array((P *)cvect,nz,1);
         expand_array(svect,cvect,nz);          expand_array(svect,cvect,nz);
         if ( DP_Print ) {          if ( DP_Print ) {
                 fprintf(asir_out,"-"); fflush(asir_out);                  fprintf(asir_out,"-"); fflush(asir_out);
Legend:
Removed from v.1.146  
changed lines
  Added in v.1.151
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>