version 1.42, 2003/08/21 08:20:06 |
version 1.43, 2003/08/22 07:12:49 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.41 2003/08/21 07:39:25 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.42 2003/08/21 08:20:06 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
#include "inline.h" |
#include "inline.h" |
Line 79 typedef struct oEPOS { |
|
Line 79 typedef struct oEPOS { |
|
int s; /* shift */ |
int s; /* shift */ |
} *EPOS; |
} *EPOS; |
|
|
|
typedef struct oBlockMask { |
|
int n; |
|
struct order_pair *order_pair; |
|
unsigned int **mask; |
|
} *BlockMask; |
|
|
int (*nd_compare_function)(unsigned int *a1,unsigned int *a2); |
int (*nd_compare_function)(unsigned int *a1,unsigned int *a2); |
|
|
static double nd_scale=2; |
static double nd_scale=2; |
static unsigned int **nd_bound; |
static unsigned int **nd_bound; |
static struct order_spec *nd_ord; |
static struct order_spec *nd_ord; |
static EPOS nd_epos; |
static EPOS nd_epos; |
|
static BlockMask nd_blockmask; |
static int nd_nvar; |
static int nd_nvar; |
static int nd_isrlex; |
static int nd_isrlex; |
static int nd_epw,nd_bpe,nd_wpd,nd_exporigin; |
static int nd_epw,nd_bpe,nd_wpd,nd_exporigin; |
Line 127 extern int Top,Reverse,dp_nelim; |
|
Line 134 extern int Top,Reverse,dp_nelim; |
|
/* macros for term comparison */ |
/* macros for term comparison */ |
#define TD_DL_COMPARE(d1,d2)\ |
#define TD_DL_COMPARE(d1,d2)\ |
(TD(d1)>TD(d2)?1:(TD(d1)<TD(d2)?-1:ndl_lex_compare(d1,d2))) |
(TD(d1)>TD(d2)?1:(TD(d1)<TD(d2)?-1:ndl_lex_compare(d1,d2))) |
|
#if 0 |
#define DL_COMPARE(d1,d2)\ |
#define DL_COMPARE(d1,d2)\ |
(nd_dcomp?TD_DL_COMPARE(d1,d2):(*nd_compare_function)(d1,d2)) |
(nd_dcomp>0?TD_DL_COMPARE(d1,d2)\ |
|
:(nd_dcomp==0?ndl_lex_compare(d1,d2)\ |
|
:(nd_blockmask?ndl_block_compare(d1,d2)\ |
|
:(*nd_compare_function)(d1,d2)))) |
|
#else |
|
#define DL_COMPARE(d1,d2)\ |
|
(nd_dcomp>0?TD_DL_COMPARE(d1,d2):(*nd_compare_function)(d1,d2)) |
|
#endif |
|
|
/* allocators */ |
/* allocators */ |
#define NEWRHist(r) \ |
#define NEWRHist(r) \ |
Line 202 NODE nd_gb_trace(int m); |
|
Line 217 NODE nd_gb_trace(int m); |
|
|
|
/* ndl functions */ |
/* ndl functions */ |
int ndl_weight(unsigned int *d); |
int ndl_weight(unsigned int *d); |
|
int ndl_weight_mask(unsigned int *d,int i); |
void ndl_dehomogenize(unsigned int *p); |
void ndl_dehomogenize(unsigned int *p); |
void ndl_reconstruct(int obpe,unsigned int *d,unsigned int *r); |
void ndl_reconstruct(int obpe,EPOS oepos,unsigned int *d,unsigned int *r); |
INLINE int ndl_reducible(unsigned int *d1,unsigned int *d2); |
INLINE int ndl_reducible(unsigned int *d1,unsigned int *d2); |
INLINE int ndl_lex_compare(unsigned int *d1,unsigned int *d2); |
INLINE int ndl_lex_compare(unsigned int *d1,unsigned int *d2); |
|
INLINE int ndl_block_compare(unsigned int *d1,unsigned int *d2); |
INLINE int ndl_equal(unsigned int *d1,unsigned int *d2); |
INLINE int ndl_equal(unsigned int *d1,unsigned int *d2); |
INLINE void ndl_copy(unsigned int *d1,unsigned int *d2); |
INLINE void ndl_copy(unsigned int *d1,unsigned int *d2); |
INLINE void ndl_add(unsigned int *d1,unsigned int *d2,unsigned int *d); |
INLINE void ndl_add(unsigned int *d1,unsigned int *d2,unsigned int *d); |
Line 234 void nd_free_private_storage(); |
|
Line 251 void nd_free_private_storage(); |
|
void _NM_alloc(); |
void _NM_alloc(); |
void _ND_alloc(); |
void _ND_alloc(); |
void nd_free(ND p); |
void nd_free(ND p); |
void nd_realloc(ND p,int obpe); |
|
void nd_free_redlist(); |
void nd_free_redlist(); |
|
|
/* printing */ |
/* printing */ |
Line 250 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs ndp |
|
Line 266 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs ndp |
|
void nd_reconstruct_direct(int mod,NDV *ps,int len); |
void nd_reconstruct_direct(int mod,NDV *ps,int len); |
void nd_setup(int mod,int trace,NODE f); |
void nd_setup(int mod,int trace,NODE f); |
void nd_setup_parameters(); |
void nd_setup_parameters(); |
|
BlockMask nd_create_blockmask(struct order_spec *ord); |
|
|
/* ND functions */ |
/* ND functions */ |
int nd_check_candidate(NODE input,NODE cand); |
int nd_check_candidate(NODE input,NODE cand); |
Line 268 INLINE int nd_length(ND p); |
|
Line 285 INLINE int nd_length(ND p); |
|
/* NDV functions */ |
/* NDV functions */ |
void ndv_mul_c(int mod,NDV p,int mul); |
void ndv_mul_c(int mod,NDV p,int mul); |
void ndv_mul_c_q(NDV p,Q mul); |
void ndv_mul_c_q(NDV p,Q mul); |
void ndv_realloc(NDV p,int obpe,int oadv); |
void ndv_realloc(NDV p,int obpe,int oadv,EPOS oepos); |
ND ndv_mul_nm(int mod,NDV p,NM m0); |
ND ndv_mul_nm(int mod,NDV p,NM m0); |
void ndv_dehomogenize(NDV p); |
void ndv_dehomogenize(NDV p); |
void ndv_removecont(int mod,NDV p); |
void ndv_removecont(int mod,NDV p); |
Line 405 INLINE int ndl_reducible(unsigned int *d1,unsigned int |
|
Line 422 INLINE int ndl_reducible(unsigned int *d1,unsigned int |
|
} |
} |
} |
} |
|
|
|
/* XXX : block order not supported */ |
|
|
void ndl_dehomogenize(unsigned int *d) |
void ndl_dehomogenize(unsigned int *d) |
{ |
{ |
unsigned int mask; |
unsigned int mask; |
Line 435 void ndl_dehomogenize(unsigned int *d) |
|
Line 454 void ndl_dehomogenize(unsigned int *d) |
|
void ndl_lcm(unsigned int *d1,unsigned *d2,unsigned int *d) |
void ndl_lcm(unsigned int *d1,unsigned *d2,unsigned int *d) |
{ |
{ |
unsigned int t1,t2,u,u1,u2; |
unsigned int t1,t2,u,u1,u2; |
int i,j; |
int i,j,l; |
|
|
switch ( nd_bpe ) { |
switch ( nd_bpe ) { |
case 4: |
case 4: |
Line 498 void ndl_lcm(unsigned int *d1,unsigned *d2,unsigned in |
|
Line 517 void ndl_lcm(unsigned int *d1,unsigned *d2,unsigned in |
|
break; |
break; |
} |
} |
TD(d) = ndl_weight(d); |
TD(d) = ndl_weight(d); |
|
if ( nd_blockmask ) { |
|
l = nd_blockmask->n; |
|
for ( j = 0; j < l; j++ ) |
|
d[j+1] = ndl_weight_mask(d,j); |
|
} |
} |
} |
|
|
int ndl_weight(unsigned int *d) |
int ndl_weight(unsigned int *d) |
Line 513 int ndl_weight(unsigned int *d) |
|
Line 537 int ndl_weight(unsigned int *d) |
|
return t; |
return t; |
} |
} |
|
|
INLINE int ndl_lex_compare(unsigned int *d1,unsigned int *d2) |
int ndl_weight_mask(unsigned int *d,int index) |
{ |
{ |
|
unsigned int t,u; |
|
unsigned int *mask; |
|
int i,j; |
|
|
|
mask = nd_blockmask->mask[index]; |
|
for ( t = 0, i = nd_exporigin; i < nd_wpd; i++ ) { |
|
u = d[i]&mask[i]; |
|
for ( j = 0; j < nd_epw; j++, u>>=nd_bpe ) |
|
t += (u&nd_mask0); |
|
} |
|
return t; |
|
} |
|
|
|
int ndl_lex_compare(unsigned int *d1,unsigned int *d2) |
|
{ |
int i; |
int i; |
|
|
d1 += nd_exporigin; |
d1 += nd_exporigin; |
Line 527 INLINE int ndl_lex_compare(unsigned int *d1,unsigned i |
|
Line 566 INLINE int ndl_lex_compare(unsigned int *d1,unsigned i |
|
return 0; |
return 0; |
} |
} |
|
|
|
int ndl_block_compare(unsigned int *d1,unsigned int *d2) |
|
{ |
|
int i,l,j,ord_o,ord_l; |
|
struct order_pair *op; |
|
unsigned int t1,t2; |
|
unsigned int *mask; |
|
|
|
l = nd_blockmask->n; |
|
op = nd_blockmask->order_pair; |
|
for ( j = 0; j < l; j++ ) { |
|
mask = nd_blockmask->mask[j]; |
|
ord_o = op[j].order; |
|
if ( ord_o < 2 ) |
|
if ( d1[j+1] > d2[j+1] ) return 1; |
|
else if ( d1[j+1] < d2[j+1] ) return -1; |
|
for ( i = nd_exporigin; i < nd_wpd; i++ ) { |
|
t1 = d1[i]&mask[i]; |
|
t2 = d2[i]&mask[i]; |
|
if ( t1 > t2 ) |
|
return !ord_o ? -1 : 1; |
|
else if ( t1 < t2 ) |
|
return !ord_o ? 1 : -1; |
|
} |
|
} |
|
return 0; |
|
} |
|
|
INLINE int ndl_equal(unsigned int *d1,unsigned int *d2) |
INLINE int ndl_equal(unsigned int *d1,unsigned int *d2) |
{ |
{ |
int i; |
int i; |
Line 562 INLINE void ndl_add(unsigned int *d1,unsigned int *d2, |
|
Line 628 INLINE void ndl_add(unsigned int *d1,unsigned int *d2, |
|
{ |
{ |
int i; |
int i; |
|
|
|
#if 1 |
switch ( nd_wpd ) { |
switch ( nd_wpd ) { |
case 2: |
case 2: |
TD(d) = TD(d1)+TD(d2); |
TD(d) = TD(d1)+TD(d2); |
Line 573 INLINE void ndl_add(unsigned int *d1,unsigned int *d2, |
|
Line 640 INLINE void ndl_add(unsigned int *d1,unsigned int *d2, |
|
d[2] = d1[2]+d2[2]; |
d[2] = d1[2]+d2[2]; |
break; |
break; |
default: |
default: |
for ( i = 0; i < nd_wpd; i++ ) |
for ( i = 0; i < nd_wpd; i++ ) d[i] = d1[i]+d2[i]; |
d[i] = d1[i]+d2[i]; |
|
break; |
break; |
} |
} |
|
#else |
|
for ( i = 0; i < nd_wpd; i++ ) d[i] = d1[i]+d2[i]; |
|
#endif |
} |
} |
|
|
INLINE void ndl_sub(unsigned int *d1,unsigned int *d2,unsigned int *d) |
INLINE void ndl_sub(unsigned int *d1,unsigned int *d2,unsigned int *d) |
{ |
{ |
int i; |
int i; |
|
|
for ( i = 0; i < nd_wpd; i++ ) |
for ( i = 0; i < nd_wpd; i++ ) d[i] = d1[i]-d2[i]; |
d[i] = d1[i]-d2[i]; |
|
} |
} |
|
|
int ndl_disjoint(unsigned int *d1,unsigned int *d2) |
int ndl_disjoint(unsigned int *d1,unsigned int *d2) |
Line 2075 void nd_gr_trace(LIST f,LIST v,int m,int homo,struct o |
|
Line 2143 void nd_gr_trace(LIST f,LIST v,int m,int homo,struct o |
|
void dltondl(int n,DL dl,unsigned int *r) |
void dltondl(int n,DL dl,unsigned int *r) |
{ |
{ |
unsigned int *d; |
unsigned int *d; |
int i; |
int i,j,l,s,ord_l,ord_o; |
|
struct order_pair *op; |
|
|
d = dl->d; |
d = dl->d; |
for ( i = 0; i < nd_wpd; i++ ) r[i] = 0; |
for ( i = 0; i < nd_wpd; i++ ) r[i] = 0; |
if ( nd_isrlex ) |
if ( nd_blockmask ) { |
for ( i = 0; i < n; i++ ) PUT_EXP(r,n-1-i,d[i]); |
l = nd_blockmask->n; |
else |
op = nd_blockmask->order_pair; |
for ( i = 0; i < n; i++ ) PUT_EXP(r,i,d[i]); |
for ( j = 0, s = 0; j < l; j++ ) { |
TD(r) = ndl_weight(r); |
ord_o = op[j].order; |
|
ord_l = op[j].length; |
|
if ( !ord_o ) |
|
for ( i = 0; i < ord_l; i++ ) |
|
PUT_EXP(r,s+ord_l-i-1,d[s+i]); |
|
else |
|
for ( i = 0; i < ord_l; i++ ) |
|
PUT_EXP(r,s+i,d[s+i]); |
|
s += ord_l; |
|
} |
|
TD(r) = ndl_weight(r); |
|
for ( j = 0; j < l; j++ ) |
|
r[j+1] = ndl_weight_mask(r,j); |
|
} else { |
|
if ( nd_isrlex ) |
|
for ( i = 0; i < n; i++ ) PUT_EXP(r,n-1-i,d[i]); |
|
else |
|
for ( i = 0; i < n; i++ ) PUT_EXP(r,i,d[i]); |
|
TD(r) = ndl_weight(r); |
|
} |
} |
} |
|
|
DL ndltodl(int n,unsigned int *ndl) |
DL ndltodl(int n,unsigned int *ndl) |
{ |
{ |
DL dl; |
DL dl; |
int *d; |
int *d; |
int i; |
int i,j,l,s,ord_l,ord_o; |
|
struct order_pair *op; |
|
|
NEWDL(dl,n); |
NEWDL(dl,n); |
dl->td = TD(ndl); |
dl->td = TD(ndl); |
d = dl->d; |
d = dl->d; |
if ( nd_isrlex ) |
if ( nd_blockmask ) { |
for ( i = 0; i < n; i++ ) |
l = nd_blockmask->n; |
d[i] = GET_EXP(ndl,n-1-i); |
op = nd_blockmask->order_pair; |
else |
for ( j = 0, s = 0; j < l; j++ ) { |
for ( i = 0; i < n; i++ ) |
ord_o = op[j].order; |
d[i] = GET_EXP(ndl,i); |
ord_l = op[j].length; |
|
if ( !ord_o ) |
|
for ( i = 0; i < ord_l; i++ ) |
|
d[s+i] = GET_EXP(ndl,s+ord_l-i-1); |
|
else |
|
for ( i = 0; i < ord_l; i++ ) |
|
d[s+i] = GET_EXP(ndl,s+i); |
|
s += ord_l; |
|
} |
|
} else { |
|
if ( nd_isrlex ) |
|
for ( i = 0; i < n; i++ ) |
|
d[i] = GET_EXP(ndl,n-1-i); |
|
else |
|
for ( i = 0; i < n; i++ ) |
|
d[i] = GET_EXP(ndl,i); |
|
} |
return dl; |
return dl; |
} |
} |
|
|
Line 2153 DP ndtodp(int mod,ND p) |
|
Line 2258 DP ndtodp(int mod,ND p) |
|
void ndl_print(unsigned int *dl) |
void ndl_print(unsigned int *dl) |
{ |
{ |
int n; |
int n; |
int i; |
int i,j,l,ord_o,ord_l,s,s0; |
|
struct order_pair *op; |
|
|
n = nd_nvar; |
n = nd_nvar; |
printf("<<"); |
printf("<<"); |
if ( nd_isrlex ) |
if ( nd_blockmask ) { |
for ( i = 0; i < n; i++ ) printf(i==n-1?"%d":"%d,",GET_EXP(dl,n-1-i)); |
l = nd_blockmask->n; |
else |
op = nd_blockmask->order_pair; |
for ( i = 0; i < n; i++ ) printf(i==n-1?"%d":"%d,",GET_EXP(dl,i)); |
for ( j = 0, s = s0 = 0; j < l; j++ ) { |
|
ord_o = op[j].order; |
|
ord_l = op[j].length; |
|
if ( !ord_o ) |
|
for ( i = 0, s0 += ord_l; i < ord_l; i++, s++ ) |
|
printf(s==n-1?"%d":"%d,",GET_EXP(dl,s0-i-1)); |
|
else |
|
for ( i = 0; i < ord_l; i++, s++ ) |
|
printf(s==n-1?"%d":"%d,",GET_EXP(dl,s)); |
|
} |
|
} else { |
|
if ( nd_isrlex ) |
|
for ( i = 0; i < n; i++ ) printf(i==n-1?"%d":"%d,",GET_EXP(dl,n-1-i)); |
|
else |
|
for ( i = 0; i < n; i++ ) printf(i==n-1?"%d":"%d,",GET_EXP(dl,i)); |
|
} |
printf(">>"); |
printf(">>"); |
} |
} |
|
|
Line 2437 unsigned int *nd_compute_bound(ND p) |
|
Line 2558 unsigned int *nd_compute_bound(ND p) |
|
} |
} |
|
|
void nd_setup_parameters() { |
void nd_setup_parameters() { |
int i; |
int i,n,elen; |
|
|
nd_epw = (sizeof(unsigned int)*8)/nd_bpe; |
nd_epw = (sizeof(unsigned int)*8)/nd_bpe; |
nd_wpd = nd_nvar/nd_epw+(nd_nvar%nd_epw?1:0); |
elen = nd_nvar/nd_epw+(nd_nvar%nd_epw?1:0); |
|
|
switch ( nd_ord->id ) { |
switch ( nd_ord->id ) { |
case 0: |
case 0: |
nd_exporigin = 1; |
nd_exporigin = 1; |
break; |
break; |
case 1: |
case 1: |
/* block order */ |
/* block order */ |
nd_exporigin = nd_ord->ord.block.length; |
/* d[0]:weight d[1]:w0,...,d[nd_exporigin-1]:w(n-1) */ |
|
nd_exporigin = nd_ord->ord.block.length+1; |
break; |
break; |
case 2: |
case 2: |
error("nd_setup_parameters : matrix order is not supported yet."); |
error("nd_setup_parameters : matrix order is not supported yet."); |
break; |
break; |
} |
} |
nd_wpd += nd_exporigin; |
nd_wpd = nd_exporigin+elen; |
|
nd_epos = (EPOS)MALLOC_ATOMIC(nd_nvar*sizeof(struct oEPOS)); |
|
for ( i = 0; i < nd_nvar; i++ ) { |
|
nd_epos[i].i = nd_exporigin + i/nd_epw; |
|
nd_epos[i].s = (nd_epw-(i%nd_epw)-1)*nd_bpe; |
|
} |
if ( nd_bpe < 32 ) { |
if ( nd_bpe < 32 ) { |
nd_mask0 = (1<<nd_bpe)-1; |
nd_mask0 = (1<<nd_bpe)-1; |
} else { |
} else { |
Line 2467 void nd_setup_parameters() { |
|
Line 2595 void nd_setup_parameters() { |
|
} |
} |
nm_adv = sizeof(struct oNM)+(nd_wpd-1)*sizeof(unsigned int); |
nm_adv = sizeof(struct oNM)+(nd_wpd-1)*sizeof(unsigned int); |
nmv_adv = sizeof(struct oNMV)+(nd_wpd-1)*sizeof(unsigned int); |
nmv_adv = sizeof(struct oNMV)+(nd_wpd-1)*sizeof(unsigned int); |
|
nd_blockmask = nd_create_blockmask(nd_ord); |
nd_epos = (EPOS)MALLOC_ATOMIC(nd_nvar*sizeof(struct oEPOS)); |
|
for ( i = 0; i < nd_nvar; i++ ) { |
|
nd_epos[i].i = nd_exporigin + i/nd_epw; |
|
nd_epos[i].s = (nd_epw-(i%nd_epw)-1)*nd_bpe; |
|
} |
|
} |
} |
|
|
ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
Line 2483 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
Line 2606 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
RHist r; |
RHist r; |
RHist *old_red; |
RHist *old_red; |
ND_pairs s0,s,t,prev_ndp_free_list; |
ND_pairs s0,s,t,prev_ndp_free_list; |
|
EPOS oepos; |
|
|
obpe = nd_bpe; |
obpe = nd_bpe; |
oadv = nmv_adv; |
oadv = nmv_adv; |
|
oepos = nd_epos; |
if ( obpe < 4 ) nd_bpe = 4; |
if ( obpe < 4 ) nd_bpe = 4; |
else if ( obpe < 6 ) nd_bpe = 6; |
else if ( obpe < 6 ) nd_bpe = 6; |
else if ( obpe < 8 ) nd_bpe = 8; |
else if ( obpe < 8 ) nd_bpe = 8; |
Line 2499 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
Line 2624 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
_nm_free_list = 0; |
_nm_free_list = 0; |
_ndp_free_list = 0; |
_ndp_free_list = 0; |
if ( mod != 0 ) |
if ( mod != 0 ) |
for ( i = nd_psn-1; i >= 0; i-- ) ndv_realloc(nd_ps[i],obpe,oadv); |
for ( i = nd_psn-1; i >= 0; i-- ) ndv_realloc(nd_ps[i],obpe,oadv,oepos); |
if ( !mod || trace ) |
if ( !mod || trace ) |
for ( i = nd_psn-1; i >= 0; i-- ) ndv_realloc(nd_psq[i],obpe,oadv); |
for ( i = nd_psn-1; i >= 0; i-- ) ndv_realloc(nd_psq[i],obpe,oadv,oepos); |
s0 = 0; |
s0 = 0; |
for ( t = d; t; t = NEXT(t) ) { |
for ( t = d; t; t = NEXT(t) ) { |
NEXTND_pairs(s0,s); |
NEXTND_pairs(s0,s); |
s->i1 = t->i1; |
s->i1 = t->i1; |
s->i2 = t->i2; |
s->i2 = t->i2; |
SG(s) = SG(t); |
SG(s) = SG(t); |
ndl_reconstruct(obpe,LCM(t),LCM(s)); |
ndl_reconstruct(obpe,oepos,LCM(t),LCM(s)); |
} |
} |
|
|
old_red = (RHist *)ALLOCA(REDTAB_LEN*sizeof(RHist)); |
old_red = (RHist *)ALLOCA(REDTAB_LEN*sizeof(RHist)); |
Line 2521 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
Line 2646 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
NEWRHist(mr); |
NEWRHist(mr); |
mr->index = r->index; |
mr->index = r->index; |
SG(mr) = SG(r); |
SG(mr) = SG(r); |
ndl_reconstruct(obpe,DL(r),DL(mr)); |
ndl_reconstruct(obpe,oepos,DL(r),DL(mr)); |
h = ndl_hash_value(DL(mr)); |
h = ndl_hash_value(DL(mr)); |
NEXT(mr) = nd_red[h]; |
NEXT(mr) = nd_red[h]; |
nd_red[h] = mr; |
nd_red[h] = mr; |
Line 2530 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
Line 2655 ND_pairs nd_reconstruct(int mod,int trace,ND_pairs d) |
|
old_red = 0; |
old_red = 0; |
for ( i = 0; i < nd_psn; i++ ) { |
for ( i = 0; i < nd_psn; i++ ) { |
NEWRHist(r); SG(r) = SG(nd_psh[i]); |
NEWRHist(r); SG(r) = SG(nd_psh[i]); |
ndl_reconstruct(obpe,DL(nd_psh[i]),DL(r)); |
ndl_reconstruct(obpe,oepos,DL(nd_psh[i]),DL(r)); |
nd_psh[i] = r; |
nd_psh[i] = r; |
} |
} |
if ( s0 ) NEXT(s) = 0; |
if ( s0 ) NEXT(s) = 0; |
Line 2548 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
Line 2673 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
RHist r; |
RHist r; |
RHist *old_red; |
RHist *old_red; |
ND_pairs s0,s,t,prev_ndp_free_list; |
ND_pairs s0,s,t,prev_ndp_free_list; |
|
EPOS oepos; |
|
|
obpe = nd_bpe; |
obpe = nd_bpe; |
oadv = nmv_adv; |
oadv = nmv_adv; |
|
oepos = nd_epos; |
if ( obpe < 4 ) nd_bpe = 4; |
if ( obpe < 4 ) nd_bpe = 4; |
else if ( obpe < 6 ) nd_bpe = 6; |
else if ( obpe < 6 ) nd_bpe = 6; |
else if ( obpe < 8 ) nd_bpe = 8; |
else if ( obpe < 8 ) nd_bpe = 8; |
Line 2562 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
Line 2689 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
prev_nm_free_list = _nm_free_list; |
prev_nm_free_list = _nm_free_list; |
prev_ndp_free_list = _ndp_free_list; |
prev_ndp_free_list = _ndp_free_list; |
_nm_free_list = 0; _ndp_free_list = 0; |
_nm_free_list = 0; _ndp_free_list = 0; |
for ( i = len-1; i >= 0; i-- ) ndv_realloc(ps[i],obpe,oadv); |
for ( i = len-1; i >= 0; i-- ) ndv_realloc(ps[i],obpe,oadv,oepos); |
old_red = (RHist *)ALLOCA(REDTAB_LEN*sizeof(RHist)); |
old_red = (RHist *)ALLOCA(REDTAB_LEN*sizeof(RHist)); |
for ( i = 0; i < REDTAB_LEN; i++ ) { |
for ( i = 0; i < REDTAB_LEN; i++ ) { |
old_red[i] = nd_red[i]; |
old_red[i] = nd_red[i]; |
Line 2573 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
Line 2700 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
NEWRHist(mr); |
NEWRHist(mr); |
mr->index = r->index; |
mr->index = r->index; |
SG(mr) = SG(r); |
SG(mr) = SG(r); |
ndl_reconstruct(obpe,DL(r),DL(mr)); |
ndl_reconstruct(obpe,oepos,DL(r),DL(mr)); |
h = ndl_hash_value(DL(mr)); |
h = ndl_hash_value(DL(mr)); |
NEXT(mr) = nd_red[h]; |
NEXT(mr) = nd_red[h]; |
nd_red[h] = mr; |
nd_red[h] = mr; |
Line 2585 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
Line 2712 void nd_reconstruct_direct(int mod,NDV *ps,int len) |
|
GC_gcollect(); |
GC_gcollect(); |
} |
} |
|
|
void ndl_reconstruct(int obpe,unsigned int *d,unsigned int *r) |
void ndl_reconstruct(int obpe,EPOS oepos,unsigned int *d,unsigned int *r) |
{ |
{ |
int n,i,ei,oepw,cepw,cbpe; |
int n,i,ei,oepw,omask0,j,s,ord_l,ord_o,l; |
|
struct order_pair *op; |
|
#define GET_EXP_OLD(d,a) (((d)[oepos[a].i]>>oepos[a].s)&omask0) |
|
#define PUT_EXP_OLD(r,a,e) ((r)[oepos[a].i] |= ((e)<<oepos[a].s)) |
|
|
n = nd_nvar; |
n = nd_nvar; |
oepw = (sizeof(unsigned int)*8)/obpe; |
oepw = (sizeof(unsigned int)*8)/obpe; |
cepw = nd_epw; |
omask0 = (1<<obpe)-1; |
cbpe = nd_bpe; |
|
TD(r) = TD(d); |
TD(r) = TD(d); |
for ( i = nd_exporigin; i < nd_wpd; i++ ) r[i] = 0; |
for ( i = nd_exporigin; i < nd_wpd; i++ ) r[i] = 0; |
if ( nd_isrlex ) |
if ( nd_blockmask ) { |
for ( i = 0; i < n; i++ ) { |
l = nd_blockmask->n; |
ei = (d[(n-1-i)/oepw+nd_exporigin]>>((oepw-((n-1-i)%oepw)-1)*obpe)) |
op = nd_blockmask->order_pair; |
&((1<<obpe)-1); |
for ( i = 1; i < nd_exporigin; i++ ) |
r[(n-1-i)/cepw+nd_exporigin] |= (ei<<((cepw-((n-1-i)%cepw)-1)*cbpe)); |
r[i] = d[i]; |
|
for ( j = 0, s = 0; j < l; j++ ) { |
|
ord_o = op[j].order; |
|
ord_l = op[j].length; |
|
if ( !ord_o ) |
|
for ( i = 0; i < ord_l; i++ ) { |
|
ei = GET_EXP_OLD(d,s+ord_l-i-1); |
|
PUT_EXP(r,s+ord_l-i-1,ei); |
|
} |
|
else |
|
for ( i = 0; i < ord_l; i++ ) { |
|
ei = GET_EXP_OLD(d,s+i); |
|
PUT_EXP(r,s+i,ei); |
|
} |
|
s += ord_l; |
} |
} |
else |
} else { |
for ( i = 0; i < n; i++ ) { |
if ( nd_isrlex ) |
ei = (d[i/oepw+nd_exporigin]>>((oepw-(i%oepw)-1)*obpe)) |
for ( i = 0; i < n; i++ ) { |
&((1<<obpe)-1); |
ei = GET_EXP_OLD(d,n-1-i); |
r[i/cepw+nd_exporigin] |= (ei<<((cepw-(i%cepw)-1)*cbpe)); |
PUT_EXP(r,n-1-i,ei); |
} |
} |
} |
else |
|
for ( i = 0; i < n; i++ ) { |
void nd_realloc(ND p,int obpe) |
ei = GET_EXP_OLD(d,i); |
{ |
PUT_EXP(r,i,ei); |
NM m,mr,mr0; |
} |
|
|
if ( p ) { |
|
m = BDY(p); |
|
for ( mr0 = 0; m; m = NEXT(m) ) { |
|
NEXTNM(mr0,mr); |
|
CM(mr) = CM(m); |
|
ndl_reconstruct(obpe,DL(m),DL(mr)); |
|
} |
|
NEXT(mr) = 0; |
|
BDY(p) = mr0; |
|
} |
} |
} |
} |
|
|
Line 2743 ND ndv_mul_nm(int mod,NDV p,NM m0) |
|
Line 2876 ND ndv_mul_nm(int mod,NDV p,NM m0) |
|
} |
} |
} |
} |
|
|
void ndv_realloc(NDV p,int obpe,int oadv) |
void ndv_realloc(NDV p,int obpe,int oadv,EPOS oepos) |
{ |
{ |
NMV m,mr,mr0,t; |
NMV m,mr,mr0,t; |
int len,i,k; |
int len,i,k; |
Line 2760 void ndv_realloc(NDV p,int obpe,int oadv) |
|
Line 2893 void ndv_realloc(NDV p,int obpe,int oadv) |
|
for ( i = 0; i < len; i++, NMV_OPREV(m), NMV_PREV(mr) ) { |
for ( i = 0; i < len; i++, NMV_OPREV(m), NMV_PREV(mr) ) { |
CQ(t) = CQ(m); |
CQ(t) = CQ(m); |
for ( k = 0; k < nd_wpd; k++ ) DL(t)[k] = 0; |
for ( k = 0; k < nd_wpd; k++ ) DL(t)[k] = 0; |
ndl_reconstruct(obpe,DL(m),DL(t)); |
ndl_reconstruct(obpe,oepos,DL(m),DL(t)); |
CQ(mr) = CQ(t); |
CQ(mr) = CQ(t); |
ndl_copy(DL(t),DL(mr)); |
ndl_copy(DL(t),DL(mr)); |
} |
} |
Line 2914 NODE nd_reducebase(NODE x) |
|
Line 3047 NODE nd_reducebase(NODE x) |
|
return x; |
return x; |
} |
} |
|
|
|
/* XXX incomplete */ |
|
|
void nd_init_ord(struct order_spec *ord) |
void nd_init_ord(struct order_spec *ord) |
{ |
{ |
if ( ord->id ) |
switch ( ord->id ) { |
error("nd_gr : unsupported order"); |
|
nd_dcomp = 0; |
|
switch ( ord->ord.simple ) { |
|
case 0: |
case 0: |
nd_dcomp = 1; |
switch ( ord->ord.simple ) { |
nd_isrlex = 1; |
case 0: |
|
nd_dcomp = 1; |
|
nd_isrlex = 1; |
|
break; |
|
case 1: |
|
nd_dcomp = 1; |
|
nd_isrlex = 0; |
|
break; |
|
case 2: |
|
nd_dcomp = 0; |
|
nd_isrlex = 0; |
|
break; |
|
default: |
|
error("nd_gr : unsupported order"); |
|
} |
break; |
break; |
case 1: |
case 1: |
nd_dcomp = 1; |
/* XXX */ |
|
nd_dcomp = -1; |
nd_isrlex = 0; |
nd_isrlex = 0; |
|
nd_compare_function = ndl_block_compare; |
break; |
break; |
case 9: |
case 2: |
/* td1->td->rlex */ |
error("nd_init_ord : matrix order is not supported yet."); |
nd_dcomp = -dp_nelim; |
|
nd_isrlex = 1; |
|
break; |
break; |
default: |
|
error("nd_gr : unsupported order"); |
|
} |
} |
nd_ord = ord; |
nd_ord = ord; |
} |
} |
|
|
|
BlockMask nd_create_blockmask(struct order_spec *ord) |
|
{ |
|
int n,i,j,s,l; |
|
unsigned int *t; |
|
BlockMask bm; |
|
|
|
if ( !ord->id ) |
|
return 0; |
|
n = ord->ord.block.length; |
|
bm = (BlockMask)MALLOC(sizeof(struct oBlockMask)); |
|
bm->n = n; |
|
bm->order_pair = ord->ord.block.order_pair; |
|
bm->mask = (unsigned int **)MALLOC(n*sizeof(unsigned int *)); |
|
for ( i = 0, s = 0; i < n; i++ ) { |
|
bm->mask[i] = t |
|
= (unsigned int *)MALLOC_ATOMIC(nd_wpd*sizeof(unsigned int)); |
|
for ( j = 0; j < nd_wpd; j++ ) t[j] = 0; |
|
l = bm->order_pair[i].length; |
|
for ( j = 0; j < l; j++, s++ ) PUT_EXP(t,s,nd_mask0); |
|
} |
|
return bm; |
|
} |