version 1.160, 2009/01/05 01:47:30 |
version 1.163, 2009/01/05 06:29:46 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.159 2009/01/05 00:52:20 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.162 2009/01/05 02:47:39 noro Exp $ */ |
|
|
#include "nd.h" |
#include "nd.h" |
|
|
Line 230 void ndl_homogenize(UINT *d,UINT *r,int obpe,EPOS oepo |
|
Line 230 void ndl_homogenize(UINT *d,UINT *r,int obpe,EPOS oepo |
|
} |
} |
w = TD(d); |
w = TD(d); |
PUT_EXP(r,nd_nvar-1,weight-w); |
PUT_EXP(r,nd_nvar-1,weight-w); |
|
if ( nd_module ) MPOS(r) = MPOS(d); |
TD(r) = weight; |
TD(r) = weight; |
if ( nd_blockmask ) ndl_weight_mask(r); |
if ( nd_blockmask ) ndl_weight_mask(r); |
if ( nd_module ) MPOS(r) = MPOS(d); |
|
} |
} |
|
|
void ndl_dehomogenize(UINT *d) |
void ndl_dehomogenize(UINT *d) |
Line 364 void ndl_lcm(UINT *d1,unsigned *d2,UINT *d) |
|
Line 364 void ndl_lcm(UINT *d1,unsigned *d2,UINT *d) |
|
d[i] = u; |
d[i] = u; |
} |
} |
#endif |
#endif |
|
if ( nd_module ) MPOS(d) = MPOS(d1); |
TD(d) = ndl_weight(d); |
TD(d) = ndl_weight(d); |
if ( nd_blockmask ) ndl_weight_mask(d); |
if ( nd_blockmask ) ndl_weight_mask(d); |
if ( nd_module ) MPOS(d) = MPOS(d1); |
|
} |
} |
|
|
void ndl_max(UINT *d1,unsigned *d2,UINT *d) |
void ndl_max(UINT *d1,unsigned *d2,UINT *d) |
Line 399 int ndl_weight(UINT *d) |
|
Line 399 int ndl_weight(UINT *d) |
|
for ( j = 0; j < nd_epw; j++, u>>=nd_bpe ) |
for ( j = 0; j < nd_epw; j++, u>>=nd_bpe ) |
t += (u&nd_mask0); |
t += (u&nd_mask0); |
} |
} |
|
if ( nd_module && current_module_weight_vector && MPOS(d) ) |
|
t += current_module_weight_vector[MPOS(d)]; |
return t; |
return t; |
} |
} |
|
|
Line 725 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
Line 727 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
{ |
{ |
int i; |
int i; |
|
|
|
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_add : invalid operation"); |
|
} |
#if 1 |
#if 1 |
switch ( nd_wpd ) { |
switch ( nd_wpd ) { |
case 2: |
case 2: |
Line 743 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
Line 749 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
#else |
#else |
for ( i = 0; i < nd_wpd; i++ ) d[i] = d1[i]+d2[i]; |
for ( i = 0; i < nd_wpd; i++ ) d[i] = d1[i]+d2[i]; |
#endif |
#endif |
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_add : invalid operation"); |
|
MPOS(d) = MPOS(d1); |
|
} |
|
} |
} |
|
|
/* d1 += d2 */ |
/* d1 += d2 */ |
Line 755 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
Line 756 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
{ |
{ |
int i; |
int i; |
|
|
|
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_addto : invalid operation"); |
|
} |
#if 1 |
#if 1 |
switch ( nd_wpd ) { |
switch ( nd_wpd ) { |
case 2: |
case 2: |
Line 773 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
Line 778 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
#else |
#else |
for ( i = 0; i < nd_wpd; i++ ) d1[i] += d2[i]; |
for ( i = 0; i < nd_wpd; i++ ) d1[i] += d2[i]; |
#endif |
#endif |
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_addto : invalid operation"); |
|
MPOS(d1) = MPOS(d2); |
|
} |
|
} |
} |
|
|
INLINE void ndl_sub(UINT *d1,UINT *d2,UINT *d) |
INLINE void ndl_sub(UINT *d1,UINT *d2,UINT *d) |
Line 2209 ND_pairs crit_B( ND_pairs d, int s ) |
|
Line 2209 ND_pairs crit_B( ND_pairs d, int s ) |
|
lcm = (UINT *)ALLOCA(nd_wpd*sizeof(UINT)); |
lcm = (UINT *)ALLOCA(nd_wpd*sizeof(UINT)); |
while ( cur ) { |
while ( cur ) { |
tl = cur->lcm; |
tl = cur->lcm; |
if ( ndl_reducible(tl,t) |
if ( ndl_reducible(tl,t) ) { |
&& (ndl_lcm(DL(nd_psh[cur->i1]),t,lcm),!ndl_equal(lcm,tl)) |
ndl_lcm(DL(nd_psh[cur->i1]),t,lcm); |
&& (ndl_lcm(DL(nd_psh[cur->i2]),t,lcm),!ndl_equal(lcm,tl)) ) { |
if ( !ndl_equal(lcm,tl) ) { |
remove = cur; |
ndl_lcm(DL(nd_psh[cur->i2]),t,lcm); |
if ( !prev ) { |
if (!ndl_equal(lcm,tl)) { |
head = cur = NEXT(cur); |
remove = cur; |
} else { |
if ( !prev ) { |
cur = NEXT(prev) = NEXT(cur); |
head = cur = NEXT(cur); |
} |
} else { |
FREENDP(remove); |
cur = NEXT(prev) = NEXT(cur); |
|
} |
|
FREENDP(remove); |
|
} else { |
|
prev = cur; cur = NEXT(cur); |
|
} |
|
} else { |
|
prev = cur; cur = NEXT(cur); |
|
} |
} else { |
} else { |
prev = cur; cur = NEXT(cur); |
prev = cur; cur = NEXT(cur); |
} |
} |
Line 2637 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe |
|
Line 2645 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe |
|
NumberField nf; |
NumberField nf; |
struct order_spec *ord1; |
struct order_spec *ord1; |
|
|
|
nd_module = 0; |
if ( !m && Demand ) nd_demand = 1; |
if ( !m && Demand ) nd_demand = 1; |
else nd_demand = 0; |
else nd_demand = 0; |
|
|
Line 2823 void nd_gr_trace(LIST f,LIST v,int trace,int homo,int |
|
Line 2832 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; |
|
|
|
nd_module = 0; |
if ( DP_Multiple ) |
if ( DP_Multiple ) |
nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1); |
nd_scale = ((double)DP_Multiple)/(double)(Denominator?Denominator:1); |
|
|
Line 3773 void weyl_mul_nm_nmv(int n,int mod,NM m0,NMV m1,NM *ta |
|
Line 3783 void weyl_mul_nm_nmv(int n,int mod,NM m0,NMV m1,NM *ta |
|
} |
} |
if ( nd_module ) { |
if ( nd_module ) { |
mpos = MPOS(d1); |
mpos = MPOS(d1); |
|
TD(d1) = ndl_weight(d1); |
if ( MPOS(d0) ) error("weyl_mul_nm_nmv : invalid operation"); |
if ( MPOS(d0) ) error("weyl_mul_nm_nmv : invalid operation"); |
} |
} |
tab[0] = m; |
tab[0] = m; |
Line 3809 void weyl_mul_nm_nmv(int n,int mod,NM m0,NMV m1,NM *ta |
|
Line 3820 void weyl_mul_nm_nmv(int n,int mod,NM m0,NMV m1,NM *ta |
|
TD(d) = s; |
TD(d) = s; |
PUT_EXP(d,n-1,s-h); |
PUT_EXP(d,n-1,s-h); |
} else TD(d) = h; |
} else TD(d) = h; |
|
if ( nd_module ) { |
|
MPOS(d) = mpos; |
|
TD(d) = ndl_weight(d); |
|
} |
if ( nd_blockmask ) ndl_weight_mask(d); |
if ( nd_blockmask ) ndl_weight_mask(d); |
if ( nd_module ) MPOS(d) = mpos; |
|
if ( mod ) c = ctab[j]; |
if ( mod ) c = ctab[j]; |
else q = ctab_q[j]; |
else q = ctab_q[j]; |
p = tab+curlen*j; |
p = tab+curlen*j; |
Line 4150 NDV pltondv(VL vl,VL dvl,LIST p) |
|
Line 4164 NDV pltondv(VL vl,VL dvl,LIST p) |
|
r = 0; |
r = 0; |
for ( i = 1, t = BDY(p); t; t = NEXT(t), i++ ) { |
for ( i = 1, t = BDY(p); t; t = NEXT(t), i++ ) { |
ri = ptond(vl,dvl,(P)BDY(t)); |
ri = ptond(vl,dvl,(P)BDY(t)); |
if ( ri ) |
if ( ri ) |
for ( m = BDY(ri); m; m = NEXT(m) ) |
for ( m = BDY(ri); m; m = NEXT(m) ) { |
MPOS(DL(m)) = i; |
MPOS(DL(m)) = i; |
|
TD(DL(m)) = ndl_weight(DL(m)); |
|
if ( nd_blockmask ) ndl_weight_mask(DL(m)); |
|
} |
r = nd_add(0,r,ri); |
r = nd_add(0,r,ri); |
} |
} |
return ndtondv(0,r); |
return ndtondv(0,r); |
Line 6416 int nd_monic(int mod,ND *p) |
|
Line 6433 int nd_monic(int mod,ND *p) |
|
PUT_EXP(DL(mr),i+nd_ntrans,e); |
PUT_EXP(DL(mr),i+nd_ntrans,e); |
td += MUL_WEIGHT(e,i+nd_ntrans); |
td += MUL_WEIGHT(e,i+nd_ntrans); |
} |
} |
|
if ( nd_module ) MPOS(DL(mr)) = MPOS(DL(m)); |
TD(DL(mr)) = td; |
TD(DL(mr)) = td; |
if ( nd_blockmask) ndl_weight_mask(DL(mr)); |
if ( nd_blockmask) ndl_weight_mask(DL(mr)); |
if ( nd_module ) MPOS(DL(mr)) = MPOS(DL(m)); |
|
} |
} |
} |
} |
NEXT(mr) = 0; |
NEXT(mr) = 0; |