[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.199 and 1.200

version 1.199, 2012/08/27 05:38:00 version 1.200, 2012/12/17 07:20:44
Line 1 
Line 1 
 /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.198 2012/04/10 07:15:07 noro Exp $ */  /* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.199 2012/08/27 05:38:00 noro Exp $ */
   
 #include "nd.h"  #include "nd.h"
   
Line 64  LIST ndvtopl(int mod,VL vl,VL dvl,NDV p,int rank);
Line 64  LIST ndvtopl(int mod,VL vl,VL dvl,NDV p,int rank);
 NDV pltondv(VL vl,VL dvl,LIST p);  NDV pltondv(VL vl,VL dvl,LIST p);
 void pltozpl(LIST l,Q *cont,LIST *pp);  void pltozpl(LIST l,Q *cont,LIST *pp);
 void ndl_max(UINT *d1,unsigned *d2,UINT *d);  void ndl_max(UINT *d1,unsigned *d2,UINT *d);
 pointer GC_malloc_atomic_ignore_off_page(int);  
 void nmtodp(int mod,NM m,DP *r);  void nmtodp(int mod,NM m,DP *r);
 NODE reverse_node(NODE n);  NODE reverse_node(NODE n);
 P ndc_div(int mod,union oNDC a,union oNDC b);  P ndc_div(int mod,union oNDC a,union oNDC b);
Line 92  void _NM_alloc()
Line 91  void _NM_alloc()
     int i;      int i;
   
     for ( i = 0; i < 1024; i++ ) {      for ( i = 0; i < 1024; i++ ) {
         p = (NM)GC_malloc(sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT));          p = (NM)MALLOC(sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT));
         p->next = _nm_free_list; _nm_free_list = p;          p->next = _nm_free_list; _nm_free_list = p;
     }      }
 }  }
Line 103  void _ND_alloc()
Line 102  void _ND_alloc()
     int i;      int i;
   
     for ( i = 0; i < 1024; i++ ) {      for ( i = 0; i < 1024; i++ ) {
         p = (ND)GC_malloc(sizeof(struct oND));          p = (ND)MALLOC(sizeof(struct oND));
         p->body = (NM)_nd_free_list; _nd_free_list = p;          p->body = (NM)_nd_free_list; _nd_free_list = p;
     }      }
 }  }
Line 114  void _NDP_alloc()
Line 113  void _NDP_alloc()
     int i;      int i;
   
     for ( i = 0; i < 1024; i++ ) {      for ( i = 0; i < 1024; i++ ) {
         p = (ND_pairs)GC_malloc(sizeof(struct oND_pairs)          p = (ND_pairs)MALLOC(sizeof(struct oND_pairs)
             +(nd_wpd-1)*sizeof(UINT));              +(nd_wpd-1)*sizeof(UINT));
         p->next = _ndp_free_list; _ndp_free_list = p;          p->next = _ndp_free_list; _ndp_free_list = p;
     }      }
Line 1581  void free_pbucket(PGeoBucket b) {
Line 1580  void free_pbucket(PGeoBucket b) {
             nd_free(b->body[i]);              nd_free(b->body[i]);
             b->body[i] = 0;              b->body[i] = 0;
         }          }
     GC_free(b);      GCFREE(b);
 }  }
   
 void add_pbucket_symbolic(PGeoBucket g,ND d)  void add_pbucket_symbolic(PGeoBucket g,ND d)
Line 3842  void nd_free(ND p)
Line 3841  void nd_free(ND p)
   
 void ndv_free(NDV p)  void ndv_free(NDV p)
 {  {
     GC_free(BDY(p));      GCFREE(BDY(p));
 }  }
   
 void nd_append_red(UINT *d,int i)  void nd_append_red(UINT *d,int i)
Line 4807  NDV ndtondv(int mod,ND p)
Line 4806  NDV ndtondv(int mod,ND p)
     if ( !p ) return 0;      if ( !p ) return 0;
     len = LEN(p);      len = LEN(p);
     if ( mod )      if ( mod )
         m0 = m = (NMV)GC_malloc_atomic_ignore_off_page(len*nmv_adv);          m0 = m = (NMV)MALLOC_ATOMIC_IGNORE_OFF_PAGE(len*nmv_adv);
     else      else
         m0 = m = MALLOC(len*nmv_adv);          m0 = m = MALLOC(len*nmv_adv);
 #if 0  #if 0
Line 5502  NDV vect_to_ndv(UINT *vect,int spcol,int col,int *rhea
Line 5501  NDV vect_to_ndv(UINT *vect,int spcol,int col,int *rhea
     for ( j = 0, len = 0; j < spcol; j++ ) if ( vect[j] ) len++;      for ( j = 0, len = 0; j < spcol; j++ ) if ( vect[j] ) len++;
     if ( !len ) return 0;      if ( !len ) return 0;
     else {      else {
         mr0 = (NMV)GC_malloc_atomic_ignore_off_page(nmv_adv*len);          mr0 = (NMV)MALLOC_ATOMIC_IGNORE_OFF_PAGE(nmv_adv*len);
 #if 0  #if 0
         ndv_alloc += nmv_adv*len;          ndv_alloc += nmv_adv*len;
 #endif  #endif
Line 5532  NDV vect_to_ndv_q(Q *vect,int spcol,int col,int *rhead
Line 5531  NDV vect_to_ndv_q(Q *vect,int spcol,int col,int *rhead
     for ( j = 0, len = 0; j < spcol; j++ ) if ( vect[j] ) len++;      for ( j = 0, len = 0; j < spcol; j++ ) if ( vect[j] ) len++;
     if ( !len ) return 0;      if ( !len ) return 0;
     else {      else {
         mr0 = (NMV)GC_malloc(nmv_adv*len);          mr0 = (NMV)MALLOC(nmv_adv*len);
 #if 0  #if 0
         ndv_alloc += nmv_adv*len;          ndv_alloc += nmv_adv*len;
 #endif  #endif
Line 5564  NDV plain_vect_to_ndv_q(Q *vect,int col,UINT *s0vect)
Line 5563  NDV plain_vect_to_ndv_q(Q *vect,int col,UINT *s0vect)
     for ( j = 0, len = 0; j < col; j++ ) if ( vect[j] ) len++;      for ( j = 0, len = 0; j < col; j++ ) if ( vect[j] ) len++;
     if ( !len ) return 0;      if ( !len ) return 0;
     else {      else {
         mr0 = (NMV)GC_malloc(nmv_adv*len);          mr0 = (NMV)MALLOC(nmv_adv*len);
 #if 0  #if 0
         ndv_alloc += nmv_adv*len;          ndv_alloc += nmv_adv*len;
 #endif  #endif
Line 6005  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
Line 6004  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
         fflush(asir_out);          fflush(asir_out);
     }      }
     /* free index arrays */      /* free index arrays */
     for ( i = 0; i < nred; i++ ) GC_free(imat[i]->index.c);      for ( i = 0; i < nred; i++ ) GCFREE(imat[i]->index.c);
   
     /* elimination (2nd step) */      /* elimination (2nd step) */
     colstat = (int *)ALLOCA(spcol*sizeof(int));      colstat = (int *)ALLOCA(spcol*sizeof(int));
Line 6018  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
Line 6017  NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s
         NEXTNODE(r0,r); BDY(r) =          NEXTNODE(r0,r); BDY(r) =
             (pointer)vect_to_ndv(spmat[i],spcol,col,rhead,s0vect);              (pointer)vect_to_ndv(spmat[i],spcol,col,rhead,s0vect);
         SG((NDV)BDY(r)) = spsugar[i];          SG((NDV)BDY(r)) = spsugar[i];
         GC_free(spmat[i]);          GCFREE(spmat[i]);
     }      }
     if ( r0 ) NEXT(r) = 0;      if ( r0 ) NEXT(r) = 0;
   
     for ( ; i < sprow; i++ ) GC_free(spmat[i]);      for ( ; i < sprow; i++ ) GCFREE(spmat[i]);
     get_eg(&eg2); init_eg(&eg_f4_2); add_eg(&eg_f4_2,&eg1,&eg2);      get_eg(&eg2); init_eg(&eg_f4_2); add_eg(&eg_f4_2,&eg1,&eg2);
     init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg2);      init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg2);
     if ( DP_Print ) {      if ( DP_Print ) {
Line 6086  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
Line 6085  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
         fflush(asir_out);          fflush(asir_out);
     }      }
     /* free index arrays */      /* free index arrays */
 /*    for ( i = 0; i < nred; i++ ) GC_free(imat[i]->index.c); */  /*    for ( i = 0; i < nred; i++ ) GCFREE(imat[i]->index.c); */
   
     /* elimination (2nd step) */      /* elimination (2nd step) */
     colstat = (int *)ALLOCA(spcol*sizeof(int));      colstat = (int *)ALLOCA(spcol*sizeof(int));
Line 6095  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
Line 6094  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
     for ( i = 0; i < rank; i++ ) {      for ( i = 0; i < rank; i++ ) {
         w[rank-i-1] = (pointer)vect_to_ndv_q(spmat[i],spcol,col,rhead,s0vect);          w[rank-i-1] = (pointer)vect_to_ndv_q(spmat[i],spcol,col,rhead,s0vect);
         SG((NDV)w[rank-i-1]) = spsugar[i];          SG((NDV)w[rank-i-1]) = spsugar[i];
 /*        GC_free(spmat[i]); */  /*        GCFREE(spmat[i]); */
     }      }
 #if 0  #if 0
     qsort(w,rank,sizeof(NDV),      qsort(w,rank,sizeof(NDV),
Line 6107  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
Line 6106  NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U
     }      }
     if ( r0 ) NEXT(r) = 0;      if ( r0 ) NEXT(r) = 0;
   
 /*    for ( ; i < sprow; i++ ) GC_free(spmat[i]); */  /*    for ( ; i < sprow; i++ ) GCFREE(spmat[i]); */
     get_eg(&eg2); init_eg(&eg_f4_2); add_eg(&eg_f4_2,&eg1,&eg2);      get_eg(&eg2); init_eg(&eg_f4_2); add_eg(&eg_f4_2,&eg1,&eg2);
     init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg2);      init_eg(&eg_f4); add_eg(&eg_f4,&eg0,&eg2);
     if ( DP_Print ) {      if ( DP_Print ) {
Line 6239  NDV nd_recv_ndv()
Line 6238  NDV nd_recv_ndv()
     len = nd_recv_int();      len = nd_recv_int();
     if ( !len ) return 0;      if ( !len ) return 0;
     else {      else {
         m0 = m = (NMV)GC_malloc_atomic_ignore_off_page(nmv_adv*len);          m0 = m = (NMV)MALLOC_ATOMIC_IGNORE_OFF_PAGE(nmv_adv*len);
 #if 0  #if 0
         ndv_alloc += len*nmv_adv;          ndv_alloc += len*nmv_adv;
 #endif  #endif
Legend:
Removed from v.1.199  
changed lines
  Added in v.1.200
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>