version 1.198, 2012/04/10 07:15:07 |
version 1.201, 2013/01/30 07:44:48 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.197 2011/08/15 09:21:00 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.200 2012/12/17 07:20:44 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); |
|
|
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 2807 NODE postprocess_algcoef(VL av,NODE alist,NODE r) |
|
Line 2806 NODE postprocess_algcoef(VL av,NODE alist,NODE r) |
|
return u0; |
return u0; |
} |
} |
|
|
void nd_gr(LIST f,LIST v,int m,int homo,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,vc,av; |
NODE fd,fd0,r,r0,t,x,s,xx,alist; |
NODE fd,fd0,r,r0,t,x,s,xx,alist; |
|
|
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 ) BDY(r) = ndvtopl(m,CO,vv,BDY(t),mrank); |
if ( nd_module ) BDY(r) = ndvtopl(m,CO,vv,BDY(t),mrank); |
else BDY(r) = ndvtop(m,CO,vv,BDY(t)); |
else if ( retdp ) BDY(r) = ndvtodp(m,BDY(t)); |
|
else BDY(r) = ndvtop(m,CO,vv,BDY(t)); |
} |
} |
if ( r0 ) NEXT(r) = 0; |
if ( r0 ) NEXT(r) = 0; |
if ( nalg ) |
if ( nalg ) |
Line 3841 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 3952 void nd_setup_parameters(int nvar,int max) { |
|
Line 3952 void nd_setup_parameters(int nvar,int max) { |
|
else if ( max < 65536 ) nd_bpe = 16; |
else if ( max < 65536 ) nd_bpe = 16; |
else nd_bpe = 32; |
else nd_bpe = 32; |
} |
} |
|
if ( !do_weyl && current_dl_weight_vector ) { |
|
UINT t; |
|
/* t = max(weights) */ |
|
for ( i = 0, t = 0; i < nd_nvar; i++ ) |
|
if ( t < current_dl_weight_vector[i] ) |
|
t = current_dl_weight_vector[i]; |
|
/* i = bitsize of t */ |
|
for ( i = 0; t; t >>=1, i++ ); |
|
/* i += bitsize of nd_nvar */ |
|
for ( t = nd_nvar; t; t >>=1, i++); |
|
/* nd_bpe+i = bitsize of max(weights)*max(exp)*nd_nvar */ |
|
if ( (nd_bpe+i) >= 31 ) |
|
error("nd_setup_parameters : too large weight"); |
|
} |
nd_epw = (sizeof(UINT)*8)/nd_bpe; |
nd_epw = (sizeof(UINT)*8)/nd_bpe; |
elen = nd_nvar/nd_epw+(nd_nvar%nd_epw?1:0); |
elen = nd_nvar/nd_epw+(nd_nvar%nd_epw?1:0); |
nd_exporigin = nd_get_exporigin(nd_ord); |
nd_exporigin = nd_get_exporigin(nd_ord); |
Line 4806 NDV ndtondv(int mod,ND p) |
|
Line 4820 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 5501 NDV vect_to_ndv(UINT *vect,int spcol,int col,int *rhea |
|
Line 5515 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 5531 NDV vect_to_ndv_q(Q *vect,int spcol,int col,int *rhead |
|
Line 5545 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 5563 NDV plain_vect_to_ndv_q(Q *vect,int col,UINT *s0vect) |
|
Line 5577 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 6004 NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s |
|
Line 6018 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 6017 NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s |
|
Line 6031 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 6085 NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U |
|
Line 6099 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 6094 NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U |
|
Line 6108 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 6106 NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,U |
|
Line 6120 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 6238 NDV nd_recv_ndv() |
|
Line 6252 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 |