version 1.38, 2005/11/27 00:07:05 |
version 1.43, 2007/09/17 12:47:45 |
|
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* |
* |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.37 2005/11/26 01:28:12 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/dist.c,v 1.42 2006/08/27 22:17:27 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
|
|
|
|
#define ORD_HOMO_WW_DRL_ZIGZAG 13 |
#define ORD_HOMO_WW_DRL_ZIGZAG 13 |
|
|
int cmpdl_drl_zigzag(), cmpdl_homo_ww_drl_zigzag(); |
int cmpdl_drl_zigzag(), cmpdl_homo_ww_drl_zigzag(); |
|
int cmpdl_top_weight(); |
|
|
int (*cmpdl)()=cmpdl_revgradlex; |
int (*cmpdl)()=cmpdl_revgradlex; |
|
int (*cmpdl_tie_breaker)(); |
int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex}; |
int (*primitive_cmpdl[3])() = {cmpdl_revgradlex,cmpdl_gradlex,cmpdl_lex}; |
|
|
int do_weyl; |
int do_weyl; |
Line 118 int has_sfcoef_p(P f) |
|
Line 120 int has_sfcoef_p(P f) |
|
} |
} |
} |
} |
|
|
|
extern N *current_top_weight_vector; |
|
static int current_top_weight_len; |
|
|
void initd(struct order_spec *spec) |
void initd(struct order_spec *spec) |
{ |
{ |
|
int len,i; |
|
|
switch ( spec->id ) { |
switch ( spec->id ) { |
case 3: |
case 3: |
cmpdl = cmpdl_composite; |
cmpdl = cmpdl_composite; |
Line 165 void initd(struct order_spec *spec) |
|
Line 172 void initd(struct order_spec *spec) |
|
} |
} |
break; |
break; |
} |
} |
|
if ( current_top_weight_vector ) { |
|
cmpdl_tie_breaker = cmpdl; |
|
cmpdl = cmpdl_top_weight; |
|
for ( len = 0, i = 0; i < spec->nv; i++ ) |
|
if ( current_top_weight_vector[i] ) |
|
len = MAX(PL(current_top_weight_vector[i]),len); |
|
current_top_weight_len = len; |
|
} |
dp_current_spec = spec; |
dp_current_spec = spec; |
} |
} |
|
|
Line 1121 int compd(VL vl,DP p1,DP p2) |
|
Line 1136 int compd(VL vl,DP p1,DP p2) |
|
return p2 ? -1 : 0; |
return p2 ? -1 : 0; |
else if ( !p2 ) |
else if ( !p2 ) |
return 1; |
return 1; |
|
else if ( NV(p1) != NV(p2) ) |
|
error("compd : size mismatch"); |
else { |
else { |
for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2); |
for ( n = NV(p1), m1 = BDY(p1), m2 = BDY(p2); |
m1 && m2; m1 = NEXT(m1), m2 = NEXT(m2) ) |
m1 && m2; m1 = NEXT(m1), m2 = NEXT(m2) ) |
Line 1496 int cmpdl_order_pair(int n,DL d1,DL d2) |
|
Line 1513 int cmpdl_order_pair(int n,DL d1,DL d2) |
|
struct order_pair *pair; |
struct order_pair *pair; |
|
|
len = dp_current_spec->ord.block.length; |
len = dp_current_spec->ord.block.length; |
|
if ( n != dp_current_spec->nv ) |
|
error("cmpdl_order_pair : incompatible order specification"); |
pair = dp_current_spec->ord.block.order_pair; |
pair = dp_current_spec->ord.block.order_pair; |
|
|
head = 0; |
head = 0; |
Line 1636 int cmpdl_matrix(int n,DL d1,DL d2) |
|
Line 1655 int cmpdl_matrix(int n,DL d1,DL d2) |
|
return 0; |
return 0; |
} |
} |
|
|
|
int cmpdl_top_weight(int n,DL d1,DL d2) |
|
{ |
|
int *w; |
|
N sum,wm,wma,t; |
|
N *a; |
|
struct oN tn; |
|
int len,i,sgn,tsgn; |
|
int *t1,*t2; |
|
|
|
w = (int *)ALLOCA(n*sizeof(int)); |
|
len = current_top_weight_len+3; |
|
t1 = d1->d; t2 = d2->d; |
|
for ( i = 0; i < n; i++ ) w[i] = t1[i]-t2[i]; |
|
sum = (N)W_ALLOC(len); sgn = 0; |
|
wm = (N)W_ALLOC(len); |
|
wma = (N)W_ALLOC(len); |
|
a = current_top_weight_vector; |
|
for ( i = 0; i < n; i++ ) { |
|
if ( !a[i] || !w[i] ) continue; |
|
tn.p = 1; |
|
if ( w[i] > 0 ) { |
|
tn.b[0] = w[i]; tsgn = 1; |
|
} else { |
|
tn.b[0] = -w[i]; tsgn = -1; |
|
} |
|
_muln(a[i],&tn,wm); |
|
if ( !sgn ) { |
|
sgn = tsgn; |
|
t = wm; wm = sum; sum = t; |
|
} else if ( sgn == tsgn ) { |
|
_addn(sum,wm,wma); |
|
if ( !PL(wma) ) |
|
sgn = 0; |
|
t = wma; wma = sum; sum = t; |
|
} else { |
|
sgn *= _subn(sum,wm,wma); |
|
t = wma; wma = sum; sum = t; |
|
} |
|
} |
|
if ( sgn > 0 ) return 1; |
|
else if ( sgn < 0 ) return -1; |
|
else return (*cmpdl_tie_breaker)(n,d1,d2); |
|
} |
|
|
GeoBucket create_bucket() |
GeoBucket create_bucket() |
{ |
{ |
GeoBucket g; |
GeoBucket g; |
Line 1799 int ni_next(int *a,int n) |
|
Line 1862 int ni_next(int *a,int n) |
|
|
|
int comp_nbm(NBM a,NBM b) |
int comp_nbm(NBM a,NBM b) |
{ |
{ |
int d,i,w; |
int d,i,w,ai,bi; |
int *ab,*bb; |
int *ab,*bb; |
|
|
if ( a->d > b->d ) return 1; |
if ( a->d > b->d ) return 1; |
else if ( a->d < b->d ) return -1; |
else if ( a->d < b->d ) return -1; |
else { |
else { |
d = a->d; ab = a->b; bb = b->b; |
d = a->d; ab = a->b; bb = b->b; |
|
#if 0 |
w = (d+31)/32; |
w = (d+31)/32; |
for ( i = 0; i < w; i++ ) |
for ( i = 0; i < w; i++ ) |
if ( ab[i] > bb[i] ) return 1; |
if ( ab[i] > bb[i] ) return 1; |
else if ( ab[i] < bb[i] ) return -1; |
else if ( ab[i] < bb[i] ) return -1; |
|
#else |
|
for ( i = 0; i < d; i++ ) { |
|
ai = NBM_GET(ab,i); |
|
bi = NBM_GET(bb,i); |
|
if ( ai > bi ) return 1; |
|
else if ( ai < bi ) return -1; |
|
} |
|
#endif |
return 0; |
return 0; |
} |
} |
} |
} |
Line 1819 NBM mul_nbm(NBM a,NBM b) |
|
Line 1891 NBM mul_nbm(NBM a,NBM b) |
|
int ad,bd,d,i,j; |
int ad,bd,d,i,j; |
int *ab,*bb,*mb; |
int *ab,*bb,*mb; |
NBM m; |
NBM m; |
Q c,c1; |
|
NODE r; |
NODE r; |
NBP u; |
NBP u; |
|
|
ad = a->d; bd = b->d; ab = a->b; bb = b->b; |
ad = a->d; bd = b->d; ab = a->b; bb = b->b; |
d = ad + bd; |
d = ad + bd; |
NEWNBM(m); NEWNBMBDY(m,d); |
NEWNBM(m); NEWNBMBDY(m,d); |
m->d = d; mulq(a->c,b->c,&m->c); mb = m->b; |
m->d = d; mulp(CO,a->c,b->c,&m->c); mb = m->b; |
j = 0; |
j = 0; |
for ( i = 0; i < ad; i++, j++ ) |
for ( i = 0; i < ad; i++, j++ ) |
if ( NBM_GET(ab,i) ) NBM_SET(mb,j); |
if ( NBM_GET(ab,i) ) NBM_SET(mb,j); |
Line 1849 NBP nbmtonbp(NBM m) |
|
Line 1920 NBP nbmtonbp(NBM m) |
|
|
|
/* a=c*x*rest -> a0= x*rest, ah=x, ar=rest */ |
/* a=c*x*rest -> a0= x*rest, ah=x, ar=rest */ |
|
|
Q separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
P separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
{ |
{ |
int i,d1; |
int i,d1; |
NBM t; |
NBM t; |
Line 1857 Q separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
Line 1928 Q separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
if ( !a->d ) error("separate_nbm : invalid argument"); |
if ( !a->d ) error("separate_nbm : invalid argument"); |
|
|
if ( a0 ) { |
if ( a0 ) { |
NEWNBM(t); t->d = a->d; t->b = a->b; t->c = ONE; |
NEWNBM(t); t->d = a->d; t->b = a->b; t->c = (P)ONE; |
*a0 = nbmtonbp(t); |
*a0 = nbmtonbp(t); |
} |
} |
|
|
if ( ah ) { |
if ( ah ) { |
NEWNBM(t); NEWNBMBDY(t,1); t->d = 1; t->c = ONE; |
NEWNBM(t); NEWNBMBDY(t,1); t->d = 1; t->c = (P)ONE; |
if ( NBM_GET(a->b,0) ) NBM_SET(t->b,0); |
if ( NBM_GET(a->b,0) ) NBM_SET(t->b,0); |
else NBM_CLR(t->b,0); |
else NBM_CLR(t->b,0); |
*ah = nbmtonbp(t); |
*ah = nbmtonbp(t); |
Line 1870 Q separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
Line 1941 Q separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
|
|
if ( ar ) { |
if ( ar ) { |
d1 = a->d-1; |
d1 = a->d-1; |
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = ONE; |
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = (P)ONE; |
for ( i = 0; i < d1; i++ ) { |
for ( i = 0; i < d1; i++ ) { |
if ( NBM_GET(a->b,i+1) ) NBM_SET(t->b,i); |
if ( NBM_GET(a->b,i+1) ) NBM_SET(t->b,i); |
else NBM_CLR(t->b,i); |
else NBM_CLR(t->b,i); |
Line 1881 Q separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
Line 1952 Q separate_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
return a->c; |
return a->c; |
} |
} |
|
|
|
/* a=c*rest*x -> a0= rest*x, ar=rest, at=x */ |
|
|
|
P separate_tail_nbm(NBM a,NBP *a0,NBP *ar,NBP *at) |
|
{ |
|
int i,d,d1; |
|
NBM t; |
|
|
|
if ( !(d=a->d) ) error("separate_tail_nbm : invalid argument"); |
|
|
|
if ( a0 ) { |
|
NEWNBM(t); t->d = a->d; t->b = a->b; t->c = (P)ONE; |
|
*a0 = nbmtonbp(t); |
|
} |
|
|
|
d1 = a->d-1; |
|
if ( at ) { |
|
NEWNBM(t); NEWNBMBDY(t,1); t->d = 1; t->c = (P)ONE; |
|
if ( NBM_GET(a->b,d1) ) NBM_SET(t->b,0); |
|
else NBM_CLR(t->b,0); |
|
*at = nbmtonbp(t); |
|
} |
|
|
|
if ( ar ) { |
|
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = (P)ONE; |
|
for ( i = 0; i < d1; i++ ) { |
|
if ( NBM_GET(a->b,i) ) NBM_SET(t->b,i); |
|
else NBM_CLR(t->b,i); |
|
} |
|
*ar = nbmtonbp(t); |
|
} |
|
|
|
return a->c; |
|
} |
|
|
NBP make_xky(int k) |
NBP make_xky(int k) |
{ |
{ |
int k1,i; |
int k1,i; |
NBM t; |
NBM t; |
|
|
NEWNBM(t); NEWNBMBDY(t,k); t->d = k; t->c = ONE; |
NEWNBM(t); NEWNBMBDY(t,k); t->d = k; t->c = (P)ONE; |
k1 = k-1; |
k1 = k-1; |
for ( i = 0; i < k1; i++ ) NBM_SET(t->b,i); |
for ( i = 0; i < k1; i++ ) NBM_SET(t->b,i); |
NBM_CLR(t->b,i); |
NBM_CLR(t->b,i); |
Line 1895 NBP make_xky(int k) |
|
Line 2000 NBP make_xky(int k) |
|
|
|
/* a=c*x^(k-1)*y*rest -> a0= x^(k-1)*y*rest, ah=x^(k-1)*y, ar=rest */ |
/* a=c*x^(k-1)*y*rest -> a0= x^(k-1)*y*rest, ah=x^(k-1)*y, ar=rest */ |
|
|
Q separate_xky_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
P separate_xky_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
{ |
{ |
int i,d1,k,k1; |
int i,d1,k,k1; |
NBM t; |
NBM t; |
Line 1909 Q separate_xky_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
Line 2014 Q separate_xky_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
k = i+1; |
k = i+1; |
|
|
if ( a0 ) { |
if ( a0 ) { |
NEWNBM(t); t->d = a->d; t->b = a->b; t->c = ONE; |
NEWNBM(t); t->d = a->d; t->b = a->b; t->c = (P)ONE; |
*a0 = nbmtonbp(t); |
*a0 = nbmtonbp(t); |
} |
} |
|
|
if ( ah ) { |
if ( ah ) { |
NEWNBM(t); NEWNBMBDY(t,k); t->d = k; t->c = ONE; |
NEWNBM(t); NEWNBMBDY(t,k); t->d = k; t->c = (P)ONE; |
for ( i = 0; i < k1; i++ ) NBM_SET(t->b,i); |
for ( i = 0; i < k1; i++ ) NBM_SET(t->b,i); |
NBM_CLR(t->b,i); |
NBM_CLR(t->b,i); |
*ah = nbmtonbp(t); |
*ah = nbmtonbp(t); |
Line 1922 Q separate_xky_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
Line 2027 Q separate_xky_nbm(NBM a,NBP *a0,NBP *ah,NBP *ar) |
|
|
|
if ( ar ) { |
if ( ar ) { |
d1 = a->d-k; |
d1 = a->d-k; |
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = ONE; |
NEWNBM(t); NEWNBMBDY(t,d1); t->d = d1; t->c = (P)ONE; |
for ( i = 0; i < d1; i++ ) { |
for ( i = 0; i < d1; i++ ) { |
if ( NBM_GET(a->b,i+k) ) NBM_SET(t->b,i); |
if ( NBM_GET(a->b,i+k) ) NBM_SET(t->b,i); |
else NBM_CLR(t->b,i); |
else NBM_CLR(t->b,i); |
Line 1941 void mulnbpnbm(VL vl,NBP p,NBM m, NBP *rp); |
|
Line 2046 void mulnbpnbm(VL vl,NBP p,NBM m, NBP *rp); |
|
NBP shuffle_mul_nbm(NBM a,NBM b) |
NBP shuffle_mul_nbm(NBM a,NBM b) |
{ |
{ |
NBP u,a0,ah,ar,b0,bh,br,a1,b1,t; |
NBP u,a0,ah,ar,b0,bh,br,a1,b1,t; |
Q ac,bc,c; |
P ac,bc,c; |
|
|
if ( !a->d || !b->d ) |
if ( !a->d || !b->d ) |
u = nbmtonbp(mul_nbm(a,b)); |
u = nbmtonbp(mul_nbm(a,b)); |
else { |
else { |
ac = separate_nbm(a,&a0,&ah,&ar); |
ac = separate_nbm(a,&a0,&ah,&ar); |
bc = separate_nbm(b,&b0,&bh,&br); |
bc = separate_nbm(b,&b0,&bh,&br); |
mulq(ac,bc,&c); |
mulp(CO,ac,bc,&c); |
shuffle_mulnbp(CO,ar,b0,&t); mulnbp(CO,ah,t,&a1); |
shuffle_mulnbp(CO,ar,b0,&t); mulnbp(CO,ah,t,&a1); |
shuffle_mulnbp(CO,a0,br,&t); mulnbp(CO,bh,t,&b1); |
shuffle_mulnbp(CO,a0,br,&t); mulnbp(CO,bh,t,&b1); |
addnbp(CO,a1,b1,&t); mulnbp(CO,(NBP)c,t,&u); |
addnbp(CO,a1,b1,&t); mulnbp(CO,(NBP)c,t,&u); |
Line 1959 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
Line 2064 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
NBP harmonic_mul_nbm(NBM a,NBM b) |
NBP harmonic_mul_nbm(NBM a,NBM b) |
{ |
{ |
NBP u,a0,ah,ar,b0,bh,br,a1,b1,t,s,abk,ab1; |
NBP u,a0,ah,ar,b0,bh,br,a1,b1,t,s,abk,ab1; |
Q ac,bc,c; |
P ac,bc,c; |
|
|
if ( !a->d || !b->d ) |
if ( !a->d || !b->d ) |
u = nbmtonbp(mul_nbm(a,b)); |
u = nbmtonbp(mul_nbm(a,b)); |
else { |
else { |
mulq(a->c,b->c,&c); |
mulp(CO,a->c,b->c,&c); |
ac = separate_xky_nbm(a,&a0,&ah,&ar); |
ac = separate_xky_nbm(a,&a0,&ah,&ar); |
bc = separate_xky_nbm(b,&b0,&bh,&br); |
bc = separate_xky_nbm(b,&b0,&bh,&br); |
mulq(ac,bc,&c); |
mulp(CO,ac,bc,&c); |
harmonic_mulnbp(CO,ar,b0,&t); mulnbp(CO,ah,t,&a1); |
harmonic_mulnbp(CO,ar,b0,&t); mulnbp(CO,ah,t,&a1); |
harmonic_mulnbp(CO,a0,br,&t); mulnbp(CO,bh,t,&b1); |
harmonic_mulnbp(CO,a0,br,&t); mulnbp(CO,bh,t,&b1); |
abk = make_xky(((NBM)BDY(BDY(ah)))->d+((NBM)BDY(BDY(bh)))->d); |
abk = make_xky(((NBM)BDY(BDY(ah)))->d+((NBM)BDY(BDY(bh)))->d); |
Line 1982 void addnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
Line 2087 void addnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
{ |
{ |
NODE b1,b2,br,br0; |
NODE b1,b2,br,br0; |
NBM m1,m2,m; |
NBM m1,m2,m; |
Q c; |
P c; |
|
|
if ( !p1 ) |
if ( !p1 ) |
*rp = p2; |
*rp = p2; |
Line 1993 void addnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
Line 2098 void addnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
m1 = (NBM)BDY(b1); m2 = (NBM)BDY(b2); |
m1 = (NBM)BDY(b1); m2 = (NBM)BDY(b2); |
switch ( comp_nbm(m1,m2) ) { |
switch ( comp_nbm(m1,m2) ) { |
case 0: |
case 0: |
addq(m1->c,m2->c,&c); |
addp(CO,m1->c,m2->c,&c); |
if ( c ) { |
if ( c ) { |
NEXTNODE(br0,br); |
NEXTNODE(br0,br); |
NEWNBM(m); m->d = m1->d; m->c = c; m->b = m1->b; |
NEWNBM(m); m->d = m1->d; m->c = c; m->b = m1->b; |
Line 2043 void chsgnnbp(NBP p,NBP *rp) |
|
Line 2148 void chsgnnbp(NBP p,NBP *rp) |
|
for ( r0 = 0, b = BDY(p); b; b = NEXT(b) ) { |
for ( r0 = 0, b = BDY(p); b; b = NEXT(b) ) { |
NEXTNODE(r0,r); |
NEXTNODE(r0,r); |
m = (NBM)BDY(b); |
m = (NBM)BDY(b); |
NEWNBM(m1); m1->d = m->d; m1->b = m->b; chsgnq(m->c,&m1->c); |
NEWNBM(m1); m1->d = m->d; m1->b = m->b; chsgnp(m->c,&m1->c); |
BDY(r) = m1; |
BDY(r) = m1; |
} |
} |
if ( r0 ) NEXT(r) = 0; |
if ( r0 ) NEXT(r) = 0; |
Line 2060 void mulnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
Line 2165 void mulnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
*rp = 0; return; |
*rp = 0; return; |
} |
} |
if ( OID(p1) != O_NBP ) { |
if ( OID(p1) != O_NBP ) { |
if ( !NUM(p1) || !RATN(p1) ) |
if ( !POLY(p1) ) |
error("mulnbp : invalid argument"); |
error("mulnbp : invalid argument"); |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (Q)p1; |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (P)p1; |
MKNODE(n,m,0); MKNBP(p1,n); |
MKNODE(n,m,0); MKNBP(p1,n); |
} |
} |
if ( OID(p2) != O_NBP ) { |
if ( OID(p2) != O_NBP ) { |
if ( !NUM(p2) || !RATN(p2) ) |
if ( !POLY(p2) ) |
error("mulnbp : invalid argument"); |
error("mulnbp : invalid argument"); |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (Q)p2; |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (P)p2; |
MKNODE(n,m,0); MKNBP(p2,n); |
MKNODE(n,m,0); MKNBP(p2,n); |
} |
} |
if ( length(BDY(p1)) < length(BDY(p2)) ) { |
if ( length(BDY(p1)) < length(BDY(p2)) ) { |
Line 2126 void pwrnbp(VL vl,NBP a,Q q,NBP *c) |
|
Line 2231 void pwrnbp(VL vl,NBP a,Q q,NBP *c) |
|
NODE r; |
NODE r; |
|
|
if ( !q ) { |
if ( !q ) { |
NEWNBM(m); m->d = 0; m->c = ONE; m->b = 0; |
NEWNBM(m); m->d = 0; m->c = (P)ONE; m->b = 0; |
MKNODE(r,m,0); MKNBP(*c,r); |
MKNODE(r,m,0); MKNBP(*c,r); |
} else if ( !a ) |
} else if ( !a ) |
*c = 0; |
*c = 0; |
Line 2157 int compnbp(VL vl,NBP p1,NBP p2) |
|
Line 2262 int compnbp(VL vl,NBP p1,NBP p2) |
|
for ( n1 = BDY(p1), n2 = BDY(p2); |
for ( n1 = BDY(p1), n2 = BDY(p2); |
n1 && n2; n1 = NEXT(n1), n2 = NEXT(n2) ) { |
n1 && n2; n1 = NEXT(n1), n2 = NEXT(n2) ) { |
m1 = (NBM)BDY(n1); m2 = (NBM)BDY(n2); |
m1 = (NBM)BDY(n1); m2 = (NBM)BDY(n2); |
if ( (t = comp_nbm(m1,m2)) || (t = cmpq(m1->c,m2->c) ) ) |
if ( (t = comp_nbm(m1,m2)) || (t = compp(CO,m1->c,m2->c) ) ) |
return t; |
return t; |
} |
} |
if ( n1 ) |
if ( n1 ) |
Line 2179 void shuffle_mulnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
Line 2284 void shuffle_mulnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
*rp = 0; return; |
*rp = 0; return; |
} |
} |
if ( OID(p1) != O_NBP ) { |
if ( OID(p1) != O_NBP ) { |
if ( !NUM(p1) || !RATN(p1) ) |
if ( !POLY(p1) ) |
error("shuffle_mulnbp : invalid argument"); |
error("shuffle_mulnbp : invalid argument"); |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (Q)p1; |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (P)p1; |
MKNODE(n,m,0); MKNBP(p1,n); |
MKNODE(n,m,0); MKNBP(p1,n); |
} |
} |
if ( OID(p2) != O_NBP ) { |
if ( OID(p2) != O_NBP ) { |
if ( !NUM(p2) || !RATN(p2) ) |
if ( !POLY(p2) ) |
error("shuffle_mulnbp : invalid argument"); |
error("shuffle_mulnbp : invalid argument"); |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (Q)p2; |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (P)p2; |
MKNODE(n,m,0); MKNBP(p2,n); |
MKNODE(n,m,0); MKNBP(p2,n); |
} |
} |
for ( r = 0, b1 = BDY(p1); b1; b1 = NEXT(b1) ) |
for ( r = 0, b1 = BDY(p1); b1; b1 = NEXT(b1) ) |
Line 2208 void harmonic_mulnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
Line 2313 void harmonic_mulnbp(VL vl,NBP p1,NBP p2, NBP *rp) |
|
*rp = 0; return; |
*rp = 0; return; |
} |
} |
if ( OID(p1) != O_NBP ) { |
if ( OID(p1) != O_NBP ) { |
if ( !NUM(p1) || !RATN(p1) ) |
if ( !POLY(p1) ) |
error("harmonic_mulnbp : invalid argument"); |
error("harmonic_mulnbp : invalid argument"); |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (Q)p1; |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (P)p1; |
MKNODE(n,m,0); MKNBP(p1,n); |
MKNODE(n,m,0); MKNBP(p1,n); |
} |
} |
if ( OID(p2) != O_NBP ) { |
if ( OID(p2) != O_NBP ) { |
if ( !NUM(p2) || !RATN(p2) ) |
if ( !POLY(p2) ) |
error("harmonic_mulnbp : invalid argument"); |
error("harmonic_mulnbp : invalid argument"); |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (Q)p2; |
NEWNBM(m); m->d = 0; m->b = 0; m->c = (P)p2; |
MKNODE(n,m,0); MKNBP(p2,n); |
MKNODE(n,m,0); MKNBP(p2,n); |
} |
} |
for ( r = 0, b1 = BDY(p1); b1; b1 = NEXT(b1) ) |
for ( r = 0, b1 = BDY(p1); b1; b1 = NEXT(b1) ) |
Line 2233 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
Line 2338 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
int ad,bd,d,i,ai,bi,bit,s; |
int ad,bd,d,i,ai,bi,bit,s; |
int *ab,*bb,*wmb,*w; |
int *ab,*bb,*wmb,*w; |
NBM wm,tm; |
NBM wm,tm; |
Q c,c1; |
P c,c1; |
NODE r,t,t1,p; |
NODE r,t,t1,p; |
NBP u; |
NBP u; |
|
|
Line 2243 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
Line 2348 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
NEWNBM(wm); NEWNBMBDY(wm,d); wmb = wm->b; |
NEWNBM(wm); NEWNBMBDY(wm,d); wmb = wm->b; |
for ( i = 0; i < ad; i++ ) w[i] = 1; |
for ( i = 0; i < ad; i++ ) w[i] = 1; |
for ( ; i < d; i++ ) w[i] = 0; |
for ( ; i < d; i++ ) w[i] = 0; |
mulq(a->c,b->c,&c); |
mulp(CO,a->c,b->c,&c); |
r = 0; |
r = 0; |
do { |
do { |
wm->d = d; wm->c = c; |
wm->d = d; wm->c = c; |
Line 2266 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
Line 2371 NBP shuffle_mul_nbm(NBM a,NBM b) |
|
break; |
break; |
} else if ( s == 0 ) { |
} else if ( s == 0 ) { |
/* add coefs */ |
/* add coefs */ |
addq(tm->c,c,&c1); |
addp(CO,tm->c,c,&c1); |
if ( c1 ) tm->c = c1; |
if ( c1 ) tm->c = c1; |
else NEXT(p) = NEXT(t); |
else NEXT(p) = NEXT(t); |
break; |
break; |
Line 2304 NBP harmonic_mul_nbm(NBM a,NBM b) |
|
Line 2409 NBP harmonic_mul_nbm(NBM a,NBM b) |
|
int da,db,d,la,lb,lmax,lmin,l,lab,la1,lb1,lab1; |
int da,db,d,la,lb,lmax,lmin,l,lab,la1,lb1,lab1; |
int i,j,k,ia,ib,s; |
int i,j,k,ia,ib,s; |
int *wa,*wb,*w,*wab,*wa1,*wmb; |
int *wa,*wb,*w,*wab,*wa1,*wmb; |
Q c,c1; |
P c,c1; |
NBM wm,tm; |
NBM wm,tm; |
NODE r,t1,t,p; |
NODE r,t1,t,p; |
NBP u; |
NBP u; |
Line 2314 NBP harmonic_mul_nbm(NBM a,NBM b) |
|
Line 2419 NBP harmonic_mul_nbm(NBM a,NBM b) |
|
wb = (int *)ALLOCA(db*sizeof(int)); |
wb = (int *)ALLOCA(db*sizeof(int)); |
la = nbmtoxky(a,wa); |
la = nbmtoxky(a,wa); |
lb = nbmtoxky(b,wb); |
lb = nbmtoxky(b,wb); |
mulq(a->c,b->c,&c); |
mulp(CO,a->c,b->c,&c); |
/* wa[0],..,wa[la-1] <-> x^wa[0]y x^wa[1]y .. */ |
/* wa[0],..,wa[la-1] <-> x^wa[0]y x^wa[1]y .. */ |
/* lmax : total length */ |
/* lmax : total length */ |
lmax = la+lb; |
lmax = la+lb; |
Line 2361 NBP harmonic_mul_nbm(NBM a,NBM b) |
|
Line 2466 NBP harmonic_mul_nbm(NBM a,NBM b) |
|
break; |
break; |
} else if ( s == 0 ) { |
} else if ( s == 0 ) { |
/* add coefs */ |
/* add coefs */ |
addq(tm->c,c,&c1); |
addp(CO,tm->c,c,&c1); |
if ( c1 ) tm->c = c1; |
if ( c1 ) tm->c = c1; |
else NEXT(p) = NEXT(t); |
else NEXT(p) = NEXT(t); |
break; |
break; |