version 1.15, 2003/07/31 01:02:39 |
version 1.16, 2003/07/31 02:56:13 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.14 2003/07/30 07:57:45 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.15 2003/07/31 01:02:39 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
#include "inline.h" |
#include "inline.h" |
Line 92 NODE append_one(NODE,int); |
|
Line 92 NODE append_one(NODE,int); |
|
#define HTD(d) ((d)->body->td) |
#define HTD(d) ((d)->body->td) |
#define HDL(d) ((d)->body->dl) |
#define HDL(d) ((d)->body->dl) |
#define HCM(d) ((d)->body->c.m) |
#define HCM(d) ((d)->body->c.m) |
|
#define HCQ(d) ((d)->body->c.z) |
#define CM(a) ((a)->c.m) |
#define CM(a) ((a)->c.m) |
|
#define CQ(a) ((a)->c.z) |
#define DL(a) ((a)->dl) |
#define DL(a) ((a)->dl) |
#define TD(a) ((a)->td) |
#define TD(a) ((a)->td) |
#define SG(a) ((a)->sugar) |
#define SG(a) ((a)->sugar) |
Line 117 if(!(r)){(c)=(r)=(s);}else{NEXT(c)=(s);(c)=(s);} |
|
Line 119 if(!(r)){(c)=(r)=(s);}else{NEXT(c)=(s);(c)=(s);} |
|
#define NEXTND_pairs(r,c) \ |
#define NEXTND_pairs(r,c) \ |
if(!(r)){NEWND_pairs(r);(c)=(r);}else{NEWND_pairs(NEXT(c));(c)=NEXT(c);} |
if(!(r)){NEWND_pairs(r);(c)=(r);}else{NEWND_pairs(NEXT(c));(c)=NEXT(c);} |
|
|
|
void nd_removecont(ND p); |
|
void ndv_removecont(NDV p); |
|
void ndv_mul_c_q(NDV p,Q mul); |
|
void nd_mul_c_q(ND p,Q mul); |
|
|
ND_pairs crit_B( ND_pairs d, int s ); |
ND_pairs crit_B( ND_pairs d, int s ); |
void nd_gr(LIST f,LIST v,int m,struct order_spec *ord,LIST *rp); |
void nd_gr(LIST f,LIST v,int m,struct order_spec *ord,LIST *rp); |
NODE nd_setup(NODE f); |
NODE nd_setup(int mod,NODE f); |
int nd_newps(ND a); |
int nd_newps(int mod,ND a); |
ND_pairs nd_minp( ND_pairs d, ND_pairs *prest ); |
ND_pairs nd_minp( ND_pairs d, ND_pairs *prest ); |
NODE update_base(NODE nd,int ndp); |
NODE update_base(NODE nd,int ndp); |
static ND_pairs equivalent_pairs( ND_pairs d1, ND_pairs *prest ); |
static ND_pairs equivalent_pairs( ND_pairs d1, ND_pairs *prest ); |
Line 129 ND_pairs crit_F( ND_pairs d1 ); |
|
Line 136 ND_pairs crit_F( ND_pairs d1 ); |
|
ND_pairs crit_M( ND_pairs d1 ); |
ND_pairs crit_M( ND_pairs d1 ); |
ND_pairs nd_newpairs( NODE g, int t ); |
ND_pairs nd_newpairs( NODE g, int t ); |
ND_pairs update_pairs( ND_pairs d, NODE /* of index */ g, int t); |
ND_pairs update_pairs( ND_pairs d, NODE /* of index */ g, int t); |
NODE nd_gb(NODE f); |
NODE nd_gb(int m,NODE f); |
void nd_free_private_storage(); |
void nd_free_private_storage(); |
void _NM_alloc(); |
void _NM_alloc(); |
void _ND_alloc(); |
void _ND_alloc(); |
Line 138 ND nd_add(ND p1,ND p2); |
|
Line 145 ND nd_add(ND p1,ND p2); |
|
ND nd_mul_nm(ND p,NM m0); |
ND nd_mul_nm(ND p,NM m0); |
ND nd_mul_ind_nm(int index,NM m0); |
ND nd_mul_ind_nm(int index,NM m0); |
ND nd_mul_term(ND p,int td,unsigned int *d); |
ND nd_mul_term(ND p,int td,unsigned int *d); |
int nd_sp(ND_pairs p,ND *nf); |
int nd_sp(int mod,ND_pairs p,ND *nf); |
int nd_find_reducer(ND g); |
int nd_find_reducer(ND g); |
int nd_nf(ND g,int full,ND *nf); |
int nd_nf(int mod,ND g,int full,ND *nf); |
ND nd_reduce(ND p1,ND p2); |
ND nd_reduce(ND p1,ND p2); |
ND nd_reduce_special(ND p1,ND p2); |
ND nd_reduce_special(ND p1,ND p2); |
void nd_free(ND p); |
void nd_free(ND p); |
void ndl_print(unsigned int *dl); |
void ndl_print(unsigned int *dl); |
void nd_print(ND p); |
void nd_print(ND p); |
|
void nd_print_q(ND p); |
|
void ndv_print(NDV p); |
|
void ndv_print_q(NDV p); |
void ndp_print(ND_pairs d); |
void ndp_print(ND_pairs d); |
int nd_length(ND p); |
int nd_length(ND p); |
void nd_monic(ND p); |
void nd_monic(ND p); |
|
|
|
|
void ndv_mul_c(NDV p,int mul); |
void ndv_mul_c(NDV p,int mul); |
ND ndv_add(ND p1,NDV p2); |
ND ndv_add(ND p1,NDV p2); |
NDV ndtondv(ND p); |
ND ndv_add_q(ND p1,NDV p2); |
void ndv_mul_nm(NDV pv,NM m,NDV r); |
NDV ndtondv(int mod,ND p); |
ND ndv_mul_nm_create(NDV p,NM m0); |
void ndv_mul_nm(int mod,NDV pv,NM m,NDV r); |
|
ND ndv_mul_nm_create(int mod,NDV p,NM m0); |
void ndv_realloc(NDV p,int obpe,int oadv); |
void ndv_realloc(NDV p,int obpe,int oadv); |
NDV dptondv(DP); |
NDV dptondv(int,DP); |
DP ndvtodp(NDV); |
DP ndvtodp(int,NDV); |
#else |
#else |
static ND *nd_ps; |
static ND *nd_ps; |
ND dptond(DP); |
ND dptond(DP); |
Line 915 ND nd_add(ND p1,ND p2) |
|
Line 926 ND nd_add(ND p1,ND p2) |
|
|
|
#if 1 |
#if 1 |
/* ret=1 : success, ret=0 : overflow */ |
/* ret=1 : success, ret=0 : overflow */ |
int nd_nf(ND g,int full,ND *rp) |
int nd_nf(int mod,ND g,int full,ND *rp) |
{ |
{ |
ND d; |
ND d; |
NM m,mrd,tail; |
NM m,mrd,tail; |
Line 927 int nd_nf(ND g,int full,ND *rp) |
|
Line 938 int nd_nf(ND g,int full,ND *rp) |
|
#else |
#else |
ND p,red; |
ND p,red; |
#endif |
#endif |
|
Q cg,cred,gcd; |
|
|
if ( !g ) { |
if ( !g ) { |
*rp = 0; |
*rp = 0; |
Line 950 int nd_nf(ND g,int full,ND *rp) |
|
Line 962 int nd_nf(ND g,int full,ND *rp) |
|
nd_free(g); nd_free(d); |
nd_free(g); nd_free(d); |
return 0; |
return 0; |
} |
} |
c1 = invm(HCM(p),nd_mod); c2 = nd_mod-HCM(g); |
|
DMAR(c1,c2,0,nd_mod,c); CM(mul) = c; |
|
#if USE_NDV |
#if USE_NDV |
ndv_mul_nm(nd_ps[index],mul,ndv_red); |
if ( mod ) { |
g = ndv_add(g,ndv_red); |
c1 = invm(HCM(p),nd_mod); c2 = nd_mod-HCM(g); |
|
DMAR(c1,c2,0,nd_mod,c); CM(mul) = c; |
|
ndv_mul_nm(mod,nd_ps[index],mul,ndv_red); |
|
g = ndv_add(g,ndv_red); |
|
} else { |
|
igcd_cofactor(HCQ(g),HCQ(nd_ps[index]),&gcd,&cg,&cred); |
|
chsgnq(cg,&CQ(mul)); |
|
nd_mul_c_q(d,cred); |
|
nd_mul_c_q(g,cred); |
|
ndv_mul_nm(mod,nd_ps[index],mul,ndv_red); |
|
g = ndv_add_q(g,ndv_red); |
|
} |
sugar = MAX(sugar,SG(ndv_red)); |
sugar = MAX(sugar,SG(ndv_red)); |
#else |
#else |
|
c1 = invm(HCM(p),nd_mod); c2 = nd_mod-HCM(g); |
|
DMAR(c1,c2,0,nd_mod,c); CM(mul) = c; |
red = nd_mul_ind_nm(index,mul); |
red = nd_mul_ind_nm(index,mul); |
g = nd_add(g,red); |
g = nd_add(g,red); |
sugar = MAX(sugar,SG(red)); |
sugar = MAX(sugar,SG(red)); |
Line 1146 ND nd_nf(ND g,int full) |
|
Line 1169 ND nd_nf(ND g,int full) |
|
} |
} |
#endif |
#endif |
|
|
NODE nd_gb(NODE f) |
NODE nd_gb(int m,NODE f) |
{ |
{ |
int i,nh,sugar,stat; |
int i,nh,sugar,stat; |
NODE r,g,gall; |
NODE r,g,gall; |
|
|
sugar = SG(l); |
sugar = SG(l); |
fprintf(asir_out,"%d",sugar); |
fprintf(asir_out,"%d",sugar); |
} |
} |
stat = nd_sp(l,&h); |
stat = nd_sp(m,l,&h); |
if ( !stat ) { |
if ( !stat ) { |
NEXT(l) = d; d = l; |
NEXT(l) = d; d = l; |
d = nd_reconstruct(d); |
d = nd_reconstruct(d); |
goto again; |
goto again; |
} |
} |
stat = nd_nf(h,!Top,&nf); |
stat = nd_nf(m,h,!Top,&nf); |
if ( !stat ) { |
if ( !stat ) { |
NEXT(l) = d; d = l; |
NEXT(l) = d; d = l; |
d = nd_reconstruct(d); |
d = nd_reconstruct(d); |
goto again; |
goto again; |
} else if ( nf ) { |
} else if ( nf ) { |
printf("+"); fflush(stdout); |
printf("+"); fflush(stdout); |
nh = nd_newps(nf); |
nh = nd_newps(m,nf); |
d = update_pairs(d,g,nh); |
d = update_pairs(d,g,nh); |
g = update_base(g,nh); |
g = update_base(g,nh); |
gall = append_one(gall,nh); |
gall = append_one(gall,nh); |
|
|
return m; |
return m; |
} |
} |
|
|
int nd_newps(ND a) |
int nd_newps(int mod,ND a) |
{ |
{ |
int len; |
int len; |
RHist r; |
RHist r; |
Line 1492 int nd_newps(ND a) |
|
Line 1515 int nd_newps(ND a) |
|
nd_bound = (unsigned int **) |
nd_bound = (unsigned int **) |
REALLOC((char *)nd_bound,nd_pslen*sizeof(unsigned int *)); |
REALLOC((char *)nd_bound,nd_pslen*sizeof(unsigned int *)); |
} |
} |
nd_monic(a); |
if ( mod ) |
|
nd_monic(a); |
|
else |
|
nd_removecont(a); |
nd_bound[nd_psn] = nd_compute_bound(a); |
nd_bound[nd_psn] = nd_compute_bound(a); |
NEWRHist(r); TD(r) = HTD(a); ndl_copy(HDL(a),DL(r)); nd_psh[nd_psn] = r; |
NEWRHist(r); TD(r) = HTD(a); ndl_copy(HDL(a),DL(r)); nd_psh[nd_psn] = r; |
#if USE_NDV |
#if USE_NDV |
nd_ps[nd_psn]= ndtondv(a); |
nd_ps[nd_psn]= ndtondv(mod,a); |
nd_free(a); |
nd_free(a); |
nd_psl[nd_psn] = len = LEN(nd_ps[nd_psn]); |
nd_psl[nd_psn] = len = LEN(nd_ps[nd_psn]); |
if ( len > nmv_len ) { |
if ( len > nmv_len ) { |
Line 1510 int nd_newps(ND a) |
|
Line 1536 int nd_newps(ND a) |
|
return nd_psn++; |
return nd_psn++; |
} |
} |
|
|
NODE nd_setup(NODE f) |
NODE nd_setup(int mod,NODE f) |
{ |
{ |
int i,j,td,len,max; |
int i,j,td,len,max; |
NODE s,s0,f0; |
NODE s,s0,f0; |
Line 1555 NODE nd_setup(NODE f) |
|
Line 1581 NODE nd_setup(NODE f) |
|
len = 0; |
len = 0; |
for ( i = 0; i < nd_psn; i++, f = NEXT(f) ) { |
for ( i = 0; i < nd_psn; i++, f = NEXT(f) ) { |
#if USE_NDV |
#if USE_NDV |
nd_ps[i] = dptondv((DP)BDY(f)); |
nd_ps[i] = dptondv(mod,(DP)BDY(f)); |
ndv_mul_c(nd_ps[i],invm(HCM(nd_ps[i]),nd_mod)); |
if ( mod ) |
|
ndv_mul_c(nd_ps[i],invm(HCM(nd_ps[i]),nd_mod)); |
|
else |
|
ndv_removecont(nd_ps[i]); |
len = MAX(len,LEN(nd_ps[i])); |
len = MAX(len,LEN(nd_ps[i])); |
#else |
#else |
nd_ps[i] = dptond((DP)BDY(f)); |
nd_ps[i] = dptond((DP)BDY(f)); |
Line 1569 NODE nd_setup(NODE f) |
|
Line 1598 NODE nd_setup(NODE f) |
|
#if USE_NDV |
#if USE_NDV |
nmv_len = 16*len; |
nmv_len = 16*len; |
NEWNDV(ndv_red); |
NEWNDV(ndv_red); |
BDY(ndv_red) = (NMV)MALLOC_ATOMIC(nmv_len*nmv_adv); |
if ( mod ) |
|
BDY(ndv_red) = (NMV)MALLOC_ATOMIC(nmv_len*nmv_adv); |
|
else |
|
BDY(ndv_red) = (NMV)MALLOC(nmv_len*nmv_adv); |
#endif |
#endif |
for ( s0 = 0, i = 0; i < nd_psn; i++ ) { |
for ( s0 = 0, i = 0; i < nd_psn; i++ ) { |
NEXTNODE(s0,s); BDY(s) = (pointer)i; |
NEXTNODE(s0,s); BDY(s) = (pointer)i; |
Line 1603 void nd_gr(LIST f,LIST v,int m,struct order_spec *ord, |
|
Line 1635 void nd_gr(LIST f,LIST v,int m,struct order_spec *ord, |
|
nd_mod = m; |
nd_mod = m; |
for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) { |
for ( fd0 = 0, t = BDY(f); t; t = NEXT(t) ) { |
ptod(CO,vv,(P)BDY(t),&b); |
ptod(CO,vv,(P)BDY(t),&b); |
_dp_mod(b,m,0,&c); |
if ( m ) |
|
_dp_mod(b,m,0,&c); |
|
else |
|
c = b; |
if ( c ) { |
if ( c ) { |
NEXTNODE(fd0,fd); BDY(fd) = (pointer)c; |
NEXTNODE(fd0,fd); BDY(fd) = (pointer)c; |
} |
} |
} |
} |
if ( fd0 ) NEXT(fd) = 0; |
if ( fd0 ) NEXT(fd) = 0; |
s = nd_setup(fd0); |
s = nd_setup(m,fd0); |
x = nd_gb(s); |
x = nd_gb(m,s); |
#if 0 |
#if 0 |
x = nd_reduceall(x,m); |
x = nd_reduceall(x,m); |
#endif |
#endif |
for ( r0 = 0; x; x = NEXT(x) ) { |
for ( r0 = 0; x; x = NEXT(x) ) { |
NEXTNODE(r0,r); |
NEXTNODE(r0,r); |
#if USE_NDV |
#if USE_NDV |
a = ndvtodp(nd_ps[(int)BDY(x)]); |
a = ndvtodp(m,nd_ps[(int)BDY(x)]); |
#else |
#else |
a = ndtodp(nd_ps[(int)BDY(x)]); |
a = ndtodp(m,nd_ps[(int)BDY(x)]); |
#endif |
#endif |
_dtop_mod(CO,vv,a,(P *)&BDY(r)); |
if ( m ) |
|
_dtop_mod(CO,vv,a,(P *)&BDY(r)); |
|
else |
|
dtop(CO,vv,a,(P *)&BDY(r)); |
} |
} |
if ( r0 ) NEXT(r) = 0; |
if ( r0 ) NEXT(r) = 0; |
MKLIST(*rp,r0); |
MKLIST(*rp,r0); |
Line 1745 void nd_print(ND p) |
|
Line 1783 void nd_print(ND p) |
|
} |
} |
} |
} |
|
|
|
void nd_print_q(ND p) |
|
{ |
|
NM m; |
|
|
|
if ( !p ) |
|
printf("0\n"); |
|
else { |
|
for ( m = BDY(p); m; m = NEXT(m) ) { |
|
printf("+"); |
|
printexpr(CO,CQ(m)); |
|
printf("*"); |
|
ndl_print(DL(m)); |
|
} |
|
printf("\n"); |
|
} |
|
} |
|
|
void ndp_print(ND_pairs d) |
void ndp_print(ND_pairs d) |
{ |
{ |
ND_pairs t; |
ND_pairs t; |
Line 1755 void ndp_print(ND_pairs d) |
|
Line 1810 void ndp_print(ND_pairs d) |
|
printf("\n"); |
printf("\n"); |
} |
} |
|
|
|
void nd_removecont(ND p) |
|
{ |
|
int i,n; |
|
Q *w; |
|
Q dvr,t; |
|
NM m; |
|
|
|
for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ ); |
|
w = (Q *)ALLOCA(n*sizeof(Q)); |
|
for ( m = BDY(p), i = 0; i < n; m = NEXT(m), i++ ) |
|
w[i] = CQ(m); |
|
sortbynm(w,n); |
|
qltozl(w,n,&dvr); |
|
for ( m = BDY(p); m; m = NEXT(m) ) { |
|
divq(CQ(m),dvr,&t); CQ(m) = t; |
|
} |
|
} |
|
|
|
void ndv_removecont(NDV p) |
|
{ |
|
int i,len; |
|
Q *w; |
|
Q dvr,t; |
|
NMV m; |
|
|
|
len = p->len; |
|
w = (Q *)ALLOCA(len*sizeof(Q)); |
|
for ( m = BDY(p), i = 0; i < len; NMV_ADV(m), i++ ) |
|
w[i] = CQ(m); |
|
sortbynm(w,len); |
|
qltozl(w,len,&dvr); |
|
for ( m = BDY(p), i = 0; i < len; NMV_ADV(m), i++ ) { |
|
divq(CQ(m),dvr,&t); CQ(m) = t; |
|
} |
|
} |
|
|
void nd_monic(ND p) |
void nd_monic(ND p) |
{ |
{ |
if ( !p ) |
if ( !p ) |
Line 1777 void nd_mul_c(ND p,int mul) |
|
Line 1868 void nd_mul_c(ND p,int mul) |
|
} |
} |
} |
} |
|
|
|
void nd_mul_c_q(ND p,Q mul) |
|
{ |
|
NM m; |
|
Q c; |
|
|
|
if ( !p ) |
|
return; |
|
for ( m = BDY(p); m; m = NEXT(m) ) { |
|
mulq(CQ(m),mul,&c); CQ(m) = c; |
|
} |
|
} |
|
|
void nd_free(ND p) |
void nd_free(ND p) |
{ |
{ |
NM t,s; |
NM t,s; |
Line 1924 ND_pairs nd_reconstruct(ND_pairs d) |
|
Line 2027 ND_pairs nd_reconstruct(ND_pairs d) |
|
} |
} |
for ( i = 0; i < REDTAB_LEN; i++ ) { |
for ( i = 0; i < REDTAB_LEN; i++ ) { |
for ( mr0 = 0, r = nd_red[i]; r; r = NEXT(r) ) { |
for ( mr0 = 0, r = nd_red[i]; r; r = NEXT(r) ) { |
NEXTRHist(mr0,r); |
NEXTRHist(mr0,mr); |
mr->index = r->index; |
mr->index = r->index; |
TD(mr) = TD(r); |
TD(mr) = TD(r); |
ndl_dup(obpe,DL(r),DL(mr)); |
ndl_dup(obpe,DL(r),DL(mr)); |
Line 2011 ND nd_copy(ND p) |
|
Line 2114 ND nd_copy(ND p) |
|
} |
} |
|
|
#if USE_NDV |
#if USE_NDV |
int nd_sp(ND_pairs p,ND *rp) |
int nd_sp(int mod,ND_pairs p,ND *rp) |
{ |
{ |
NM m; |
NM m; |
NDV p1,p2; |
NDV p1,p2; |
Line 2024 int nd_sp(ND_pairs p,ND *rp) |
|
Line 2127 int nd_sp(ND_pairs p,ND *rp) |
|
lcm = p->lcm; |
lcm = p->lcm; |
td = TD(p); |
td = TD(p); |
NEWNM(m); |
NEWNM(m); |
CM(m) = HCM(p2); TD(m) = td-HTD(p1); ndl_sub(lcm,HDL(p1),DL(m)); |
if ( mod ) CM(m) = HCM(p2); |
|
else CQ(m) = HCQ(p2); |
|
TD(m) = td-HTD(p1); ndl_sub(lcm,HDL(p1),DL(m)); |
if ( ndl_check_bound2(p->i1,DL(m)) ) |
if ( ndl_check_bound2(p->i1,DL(m)) ) |
return 0; |
return 0; |
t1 = ndv_mul_nm_create(p1,m); |
t1 = ndv_mul_nm_create(mod,p1,m); |
CM(m) = nd_mod-HCM(p1); TD(m) = td-HTD(p2); ndl_sub(lcm,HDL(p2),DL(m)); |
if ( mod ) |
|
CM(m) = mod-HCM(p1); |
|
else |
|
chsgnq(HCQ(p1),&CQ(m)); |
|
TD(m) = td-HTD(p2); ndl_sub(lcm,HDL(p2),DL(m)); |
if ( ndl_check_bound2(p->i2,DL(m)) ) { |
if ( ndl_check_bound2(p->i2,DL(m)) ) { |
nd_free(t1); |
nd_free(t1); |
return 0; |
return 0; |
} |
} |
ndv_mul_nm(p2,m,ndv_red); |
ndv_mul_nm(mod,p2,m,ndv_red); |
FREENM(m); |
FREENM(m); |
*rp = ndv_add(t1,ndv_red); |
if ( mod ) |
|
*rp = ndv_add(t1,ndv_red); |
|
else |
|
*rp = ndv_add_q(t1,ndv_red); |
return 1; |
return 1; |
} |
} |
|
|
Line 2054 void ndv_mul_c(NDV p,int mul) |
|
Line 2166 void ndv_mul_c(NDV p,int mul) |
|
} |
} |
} |
} |
|
|
void ndv_mul_nm(NDV p,NM m0,NDV r) |
void ndv_mul_c_q(NDV p,Q mul) |
{ |
{ |
|
NMV m; |
|
Q c; |
|
int len,i; |
|
|
|
if ( !p ) |
|
return; |
|
len = LEN(p); |
|
for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) { |
|
mulq(CQ(m),mul,&c); CQ(m) = c; |
|
} |
|
} |
|
|
|
void ndv_mul_nm(int mod,NDV p,NM m0,NDV r) |
|
{ |
NMV m,mr,mr0; |
NMV m,mr,mr0; |
unsigned int *d,*dt,*dm; |
unsigned int *d,*dt,*dm; |
int c,n,td,i,c1,c2,len; |
int c,n,td,i,c1,c2,len; |
|
Q q; |
|
|
if ( !p ) |
if ( !p ) |
/* XXX */ |
/* XXX */ |
LEN(r) = 0; |
LEN(r) = 0; |
else { |
else { |
n = NV(p); m = BDY(p); len = LEN(p); |
n = NV(p); m = BDY(p); len = LEN(p); |
d = DL(m0); td = TD(m0); c = CM(m0); |
d = DL(m0); td = TD(m0); |
mr = BDY(r); |
mr = BDY(r); |
for ( ; len > 0; len--, NMV_ADV(m), NMV_ADV(mr) ) { |
if ( mod ) { |
c1 = CM(m); DMAR(c1,c,0,nd_mod,c2); CM(mr) = c2; |
c = CM(m0); |
TD(mr) = TD(m)+td; ndl_add(DL(m),d,DL(mr)); |
for ( ; len > 0; len--, NMV_ADV(m), NMV_ADV(mr) ) { |
|
c1 = CM(m); DMAR(c1,c,0,nd_mod,c2); CM(mr) = c2; |
|
TD(mr) = TD(m)+td; ndl_add(DL(m),d,DL(mr)); |
|
} |
|
} else { |
|
q = CQ(m0); |
|
for ( ; len > 0; len--, NMV_ADV(m), NMV_ADV(mr) ) { |
|
mulq(CQ(m),q,&CQ(mr)); |
|
TD(mr) = TD(m)+td; ndl_add(DL(m),d,DL(mr)); |
|
} |
} |
} |
NV(r) = NV(p); |
NV(r) = NV(p); |
LEN(r) = LEN(p); |
LEN(r) = LEN(p); |
Line 2077 void ndv_mul_nm(NDV p,NM m0,NDV r) |
|
Line 2213 void ndv_mul_nm(NDV p,NM m0,NDV r) |
|
} |
} |
} |
} |
|
|
ND ndv_mul_nm_create(NDV p,NM m0) |
ND ndv_mul_nm_create(int mod,NDV p,NM m0) |
{ |
{ |
NM mr,mr0; |
NM mr,mr0; |
NMV m; |
NMV m; |
unsigned int *d,*dt,*dm; |
unsigned int *d,*dt,*dm; |
int c,n,td,i,c1,c2,len; |
int c,n,td,i,c1,c2,len; |
|
Q q; |
ND r; |
ND r; |
|
|
if ( !p ) |
if ( !p ) |
return 0; |
return 0; |
else { |
else { |
n = NV(p); m = BDY(p); |
n = NV(p); m = BDY(p); |
d = DL(m0); td = TD(m0); c = CM(m0); |
d = DL(m0); td = TD(m0); |
len = LEN(p); |
len = LEN(p); |
mr0 = 0; |
mr0 = 0; |
for ( i = 0; i < len; i++, NMV_ADV(m) ) { |
if ( mod ) { |
NEXTNM(mr0,mr); |
c = CM(m0); |
c1 = CM(m); |
for ( i = 0; i < len; i++, NMV_ADV(m) ) { |
DMAR(c1,c,0,nd_mod,c2); |
NEXTNM(mr0,mr); |
CM(mr) = c2; |
c1 = CM(m); |
TD(mr) = TD(m)+td; |
DMAR(c1,c,0,nd_mod,c2); |
ndl_add(DL(m),d,DL(mr)); |
CM(mr) = c2; |
|
TD(mr) = TD(m)+td; |
|
ndl_add(DL(m),d,DL(mr)); |
|
} |
|
} else { |
|
q = CQ(m0); |
|
for ( i = 0; i < len; i++, NMV_ADV(m) ) { |
|
NEXTNM(mr0,mr); |
|
mulq(CQ(m),q,&CQ(mr)); |
|
TD(mr) = TD(m)+td; |
|
ndl_add(DL(m),d,DL(mr)); |
|
} |
} |
} |
NEXT(mr) = 0; |
NEXT(mr) = 0; |
MKND(NV(p),mr0,r); |
MKND(NV(p),mr0,r); |
Line 2114 ND ndv_add(ND p1,NDV p2) |
|
Line 2262 ND ndv_add(ND p1,NDV p2) |
|
NM head; |
NM head; |
unsigned int *d; |
unsigned int *d; |
NMV m2; |
NMV m2; |
|
Q q; |
|
|
if ( !p1 ) |
if ( !p1 ) |
return 0; |
return 0; |
Line 2149 ND ndv_add(ND p1,NDV p2) |
|
Line 2298 ND ndv_add(ND p1,NDV p2) |
|
break; |
break; |
case -1: |
case -1: |
ndl_copy(DL(m2),DL(new)); |
ndl_copy(DL(m2),DL(new)); |
CM(new) = CM(m2); |
CQ(new) = CQ(m2); |
if ( !prev ) { |
if ( !prev ) { |
/* cur = head */ |
/* cur = head */ |
prev = new; NEXT(prev) = head; head = prev; |
prev = new; NEXT(prev) = head; head = prev; |
Line 2161 ND ndv_add(ND p1,NDV p2) |
|
Line 2310 ND ndv_add(ND p1,NDV p2) |
|
} |
} |
} |
} |
for ( ; i < len; i++, NMV_ADV(m2) ) { |
for ( ; i < len; i++, NMV_ADV(m2) ) { |
td2 = TD(new) = TD(m2); CM(new) = CM(m2); ndl_copy(DL(m2),DL(new)); |
td2 = TD(new) = TD(m2); CQ(new) = CQ(m2); ndl_copy(DL(m2),DL(new)); |
if ( !prev ) { |
if ( !prev ) { |
prev = new; NEXT(prev) = 0; head = prev; |
prev = new; NEXT(prev) = 0; head = prev; |
} else { |
} else { |
Line 2181 ND ndv_add(ND p1,NDV p2) |
|
Line 2330 ND ndv_add(ND p1,NDV p2) |
|
} |
} |
} |
} |
|
|
|
ND ndv_add_q(ND p1,NDV p2) |
|
{ |
|
register NM prev,cur,new; |
|
int c,c1,c2,t,td,td2,mul,len,i; |
|
NM head; |
|
unsigned int *d; |
|
NMV m2; |
|
Q q; |
|
|
|
if ( !p1 ) |
|
return 0; |
|
else { |
|
prev = 0; head = cur = BDY(p1); |
|
NEWNM(new); len = LEN(p2); |
|
for ( m2 = BDY(p2), i = 0; cur && i < len; ) { |
|
td2 = TD(new) = TD(m2); |
|
if ( TD(cur) > td2 ) { |
|
prev = cur; cur = NEXT(cur); |
|
continue; |
|
} else if ( TD(cur) < td2 ) c = -1; |
|
else if ( nd_wpd == 1 ) { |
|
if ( DL(cur)[0] > DL(m2)[0] ) c = is_rlex ? -1 : 1; |
|
else if ( DL(cur)[0] < DL(m2)[0] ) c = is_rlex ? 1 : -1; |
|
else c = 0; |
|
} |
|
else c = ndl_compare(DL(cur),DL(m2)); |
|
switch ( c ) { |
|
case 0: |
|
addq(CQ(cur),CQ(m2),&q); |
|
if ( q ) |
|
CQ(cur) = q; |
|
else if ( !prev ) { |
|
head = NEXT(cur); FREENM(cur); cur = head; |
|
} else { |
|
NEXT(prev) = NEXT(cur); FREENM(cur); cur = NEXT(prev); |
|
} |
|
NMV_ADV(m2); i++; |
|
break; |
|
case 1: |
|
prev = cur; cur = NEXT(cur); |
|
break; |
|
case -1: |
|
ndl_copy(DL(m2),DL(new)); |
|
CQ(new) = CQ(m2); |
|
if ( !prev ) { |
|
/* cur = head */ |
|
prev = new; NEXT(prev) = head; head = prev; |
|
} else { |
|
NEXT(prev) = new; NEXT(new) = cur; prev = new; |
|
} |
|
NEWNM(new); NMV_ADV(m2); i++; |
|
break; |
|
} |
|
} |
|
for ( ; i < len; i++, NMV_ADV(m2) ) { |
|
td2 = TD(new) = TD(m2); CQ(new) = CQ(m2); ndl_copy(DL(m2),DL(new)); |
|
if ( !prev ) { |
|
prev = new; NEXT(prev) = 0; head = prev; |
|
} else { |
|
NEXT(prev) = new; NEXT(new) = 0; prev = new; |
|
} |
|
NEWNM(new); |
|
} |
|
FREENM(new); |
|
if ( head ) { |
|
BDY(p1) = head; SG(p1) = MAX(SG(p1),SG(p2)); |
|
return p1; |
|
} else { |
|
FREEND(p1); |
|
return 0; |
|
} |
|
|
|
} |
|
} |
|
|
void ndv_realloc(NDV p,int obpe,int oadv) |
void ndv_realloc(NDV p,int obpe,int oadv) |
{ |
{ |
NMV m,mr,mr0,t; |
NMV m,mr,mr0,t; |
Line 2199 void ndv_realloc(NDV p,int obpe,int oadv) |
|
Line 2423 void ndv_realloc(NDV p,int obpe,int oadv) |
|
mr = (NMV)((char *)mr0+(len-1)*nmv_adv); |
mr = (NMV)((char *)mr0+(len-1)*nmv_adv); |
t = (NMV)ALLOCA(nmv_adv); |
t = (NMV)ALLOCA(nmv_adv); |
for ( i = 0; i < len; i++, NMV_OPREV(m), NMV_PREV(mr) ) { |
for ( i = 0; i < len; i++, NMV_OPREV(m), NMV_PREV(mr) ) { |
CM(t) = CM(m); |
CQ(t) = CQ(m); |
TD(t) = TD(m); |
TD(t) = TD(m); |
for ( k = 0; k < nd_wpd; k++ ) DL(t)[k] = 0; |
for ( k = 0; k < nd_wpd; k++ ) DL(t)[k] = 0; |
ndl_dup(obpe,DL(m),DL(t)); |
ndl_dup(obpe,DL(m),DL(t)); |
CM(mr) = CM(t); |
CQ(mr) = CQ(t); |
TD(mr) = TD(t); |
TD(mr) = TD(t); |
ndl_copy(DL(t),DL(mr)); |
ndl_copy(DL(t),DL(mr)); |
} |
} |
Line 2211 void ndv_realloc(NDV p,int obpe,int oadv) |
|
Line 2435 void ndv_realloc(NDV p,int obpe,int oadv) |
|
} |
} |
} |
} |
|
|
NDV ndtondv(ND p) |
NDV ndtondv(int mod,ND p) |
{ |
{ |
NDV d; |
NDV d; |
NMV m,m0; |
NMV m,m0; |
Line 2221 NDV ndtondv(ND p) |
|
Line 2445 NDV ndtondv(ND p) |
|
if ( !p ) |
if ( !p ) |
return 0; |
return 0; |
len = nd_length(p); |
len = nd_length(p); |
m0 = m = (NMV)MALLOC_ATOMIC(len*nmv_adv); |
if ( mod ) |
|
m0 = m = (NMV)MALLOC_ATOMIC(len*nmv_adv); |
|
else |
|
m0 = m = (NMV)MALLOC(len*nmv_adv); |
for ( t = BDY(p), i = 0; t; t = NEXT(t), i++, NMV_ADV(m) ) { |
for ( t = BDY(p), i = 0; t; t = NEXT(t), i++, NMV_ADV(m) ) { |
TD(m) = TD(t); |
TD(m) = TD(t); |
ndl_copy(DL(t),DL(m)); |
ndl_copy(DL(t),DL(m)); |
CM(m) = CM(t); |
CQ(m) = CQ(t); |
} |
} |
MKNDV(NV(p),m0,len,d); |
MKNDV(NV(p),m0,len,d); |
SG(d) = SG(p); |
SG(d) = SG(p); |
return d; |
return d; |
} |
} |
|
|
NDV dptondv(DP p) |
NDV dptondv(int mod,DP p) |
{ |
{ |
NDV d; |
NDV d; |
NMV m0,m; |
NMV m0,m; |
Line 2242 NDV dptondv(DP p) |
|
Line 2469 NDV dptondv(DP p) |
|
if ( !p ) |
if ( !p ) |
return 0; |
return 0; |
for ( t = BDY(p), l = 0; t; t = NEXT(t), l++ ); |
for ( t = BDY(p), l = 0; t; t = NEXT(t), l++ ); |
m0 = m = (NMV)MALLOC_ATOMIC(l*nmv_adv); |
if ( mod ) |
|
m0 = m = (NMV)MALLOC_ATOMIC(l*nmv_adv); |
|
else |
|
m0 = m = (NMV)MALLOC(l*nmv_adv); |
n = NV(p); |
n = NV(p); |
for ( t = BDY(p), i = 0; i < l; i++, t = NEXT(t), NMV_ADV(m) ) { |
if ( mod ) { |
CM(m) = ITOS(C(t)); |
for ( t = BDY(p), i = 0; i < l; i++, t = NEXT(t), NMV_ADV(m) ) { |
TD(m) = TD(DL(t)); |
CM(m) = ITOS(C(t)); |
dltondl(n,DL(t),DL(m)); |
TD(m) = TD(DL(t)); |
|
dltondl(n,DL(t),DL(m)); |
|
} |
|
} else { |
|
for ( t = BDY(p), i = 0; i < l; i++, t = NEXT(t), NMV_ADV(m) ) { |
|
CQ(m) = (Q)C(t); |
|
TD(m) = TD(DL(t)); |
|
dltondl(n,DL(t),DL(m)); |
|
} |
} |
} |
MKNDV(n,m0,l,d); |
MKNDV(n,m0,l,d); |
SG(d) = SG(p); |
SG(d) = SG(p); |
return d; |
return d; |
} |
} |
|
|
DP ndvtodp(NDV p) |
DP ndvtodp(int mod,NDV p) |
{ |
{ |
DP d; |
DP d; |
MP m0,m; |
MP m0,m; |
Line 2266 DP ndvtodp(NDV p) |
|
Line 2504 DP ndvtodp(NDV p) |
|
m0 = 0; |
m0 = 0; |
len = LEN(p); |
len = LEN(p); |
n = NV(p); |
n = NV(p); |
for ( t = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) { |
if ( mod ) { |
NEXTMP(m0,m); |
for ( t = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) { |
C(m) = STOI(CM(t)); |
NEXTMP(m0,m); |
DL(m) = ndltodl(n,TD(t),DL(t)); |
C(m) = STOI(CM(t)); |
|
DL(m) = ndltodl(n,TD(t),DL(t)); |
|
} |
|
} else { |
|
for ( t = BDY(p), i = 0; i < len; i++, NMV_ADV(t) ) { |
|
NEXTMP(m0,m); |
|
C(m) = (P)CQ(t); |
|
DL(m) = ndltodl(n,TD(t),DL(t)); |
|
} |
} |
} |
NEXT(m) = 0; |
NEXT(m) = 0; |
MKDP(NV(p),m0,d); |
MKDP(NV(p),m0,d); |
Line 2288 void ndv_print(NDV p) |
|
Line 2534 void ndv_print(NDV p) |
|
len = LEN(p); |
len = LEN(p); |
for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) { |
for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) { |
printf("+%d*",CM(m)); |
printf("+%d*",CM(m)); |
|
ndl_print(DL(m)); |
|
} |
|
printf("\n"); |
|
} |
|
} |
|
|
|
void ndv_print_q(NDV p) |
|
{ |
|
NMV m; |
|
int i,len; |
|
|
|
if ( !p ) |
|
printf("0\n"); |
|
else { |
|
len = LEN(p); |
|
for ( m = BDY(p), i = 0; i < len; i++, NMV_ADV(m) ) { |
|
printf("+"); |
|
printexpr(CO,CQ(m)); |
|
printf("*"); |
ndl_print(DL(m)); |
ndl_print(DL(m)); |
} |
} |
printf("\n"); |
printf("\n"); |