[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.209 and 1.210

version 1.209, 2013/09/25 02:36:24 version 1.210, 2013/09/26 00:38:47
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.208 2013/09/15 04:30:28 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.209 2013/09/25 02:36:24 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 5363  Q *nm_ind_pair_to_vect(int mod,UINT *s0,int n,NM_ind_p
Line 5363  Q *nm_ind_pair_to_vect(int mod,UINT *s0,int n,NM_ind_p
     return r;      return r;
 }  }
   
 IndArray nm_ind_pair_to_vect_compress(int mod,UINT *s0,int n,NM_ind_pair pair)  IndArray nm_ind_pair_to_vect_compress(int mod,UINT *s0,int n,int *s0hash,NM_ind_pair pair)
 {  {
     NM m;      NM m;
     NMV mr;      NMV mr;
Line 5372  IndArray nm_ind_pair_to_vect_compress(int mod,UINT *s0
Line 5372  IndArray nm_ind_pair_to_vect_compress(int mod,UINT *s0
     unsigned char *ivc;      unsigned char *ivc;
     unsigned short *ivs;      unsigned short *ivs;
     UINT *v,*ivi,*s0v;      UINT *v,*ivi,*s0v;
     int i,j,len,prev,diff,cdiff;      int i,j,len,prev,diff,cdiff,h;
     IndArray r;      IndArray r;
 struct oEGT eg0,eg1;  struct oEGT eg0,eg1;
   
Line 5385  struct oEGT eg0,eg1;
Line 5385  struct oEGT eg0,eg1;
 get_eg(&eg0);  get_eg(&eg0);
     for ( i = j = 0, s = s0, mr = BDY(p); j < len; j++, NMV_ADV(mr) ) {      for ( i = j = 0, s = s0, mr = BDY(p); j < len; j++, NMV_ADV(mr) ) {
         ndl_add(d,DL(mr),t);          ndl_add(d,DL(mr),t);
         for ( ; !ndl_equal(t,s); s += nd_wpd, i++ );                  h = ndl_hash_value(t);
           for ( ; h != s0hash[i] || !ndl_equal(t,s); s += nd_wpd, i++ );
         v[j] = i;          v[j] = i;
     }      }
 get_eg(&eg1); add_eg(&eg_search,&eg0,&eg1);  get_eg(&eg1); add_eg(&eg_search,&eg0,&eg1);
Line 5845  NODE nd_f4(int m,int **indp)
Line 5846  NODE nd_f4(int m,int **indp)
         if ( DP_Print )          if ( DP_Print )
             fprintf(asir_out,"sugar=%d,symb=%fsec,",              fprintf(asir_out,"sugar=%d,symb=%fsec,",
                 sugar,eg_f4.exectime+eg_f4.gctime);                  sugar,eg_f4.exectime+eg_f4.gctime);
         if ( 1 )          nflist = nd_f4_red(m,l,0,s0vect,col,rp0,nd_gentrace?&ll:0);
             nflist = nd_f4_red(m,l,0,s0vect,col,rp0,nd_gentrace?&ll:0);  
         else  
             nflist = nd_f4_red_dist(m,l,s0vect,col,rp0,nd_gentrace?&ll:0);  
         /* adding new bases */          /* adding new bases */
         for ( r = nflist; r; r = NEXT(r) ) {          for ( r = nflist; r; r = NEXT(r) ) {
             nf = (NDV)BDY(r);              nf = (NDV)BDY(r);
Line 6099  NODE nd_f4_red(int m,ND_pairs sp0,int trace,UINT *s0ve
Line 6097  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;
   
 init_eg(&eg_search);  init_eg(&eg_search);
     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);
Line 6108  init_eg(&eg_search);
Line 6109  init_eg(&eg_search);
   
     /* construction of index arrays */      /* construction of index arrays */
     rvect = (NM_ind_pair *)ALLOCA(nred*sizeof(NM_ind_pair));      rvect = (NM_ind_pair *)ALLOCA(nred*sizeof(NM_ind_pair));
       s0hash = (int *)ALLOCA(col*sizeof(int));
       for ( i = 0, s = s0vect; i < col; i++, s += nd_wpd )
           s0hash[i] = ndl_hash_value(s);
     for ( rp = rp0, i = 0; rp; i++, rp = NEXT(rp) ) {      for ( rp = rp0, i = 0; rp; i++, rp = NEXT(rp) ) {
         rvect[i] = (NM_ind_pair)BDY(rp);          rvect[i] = (NM_ind_pair)BDY(rp);
         imat[i] = nm_ind_pair_to_vect_compress(m,s0vect,col,rvect[i]);          imat[i] = nm_ind_pair_to_vect_compress(m,s0vect,col,s0hash,rvect[i]);
         rhead[imat[i]->head] = 1;          rhead[imat[i]->head] = 1;
     }      }
     if ( m )      if ( m )
Line 6435  int ox_exec_f4_red(Q proc)
Line 6439  int ox_exec_f4_red(Q proc)
     return s;      return s;
 }  }
   
   #if 0
 NODE nd_f4_red_dist(int m,ND_pairs sp0,UINT *s0vect,int col,NODE rp0,ND_pairs *nz)  NODE nd_f4_red_dist(int m,ND_pairs sp0,UINT *s0vect,int col,NODE rp0,ND_pairs *nz)
 {  {
     int nsp,nred;      int nsp,nred;
Line 6597  void nd_exec_f4_red_dist()
Line 6602  void nd_exec_f4_red_dist()
     }      }
     fflush(nd_write);      fflush(nd_write);
 }  }
   #endif
   
 int nd_gauss_elim_q(Q **mat0,int *sugar,int row,int col,int *colstat)  int nd_gauss_elim_q(Q **mat0,int *sugar,int row,int col,int *colstat)
 {  {
Legend:
Removed from v.1.209  
changed lines
  Added in v.1.210
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>