version 1.38, 2005/11/27 00:07:05 |
version 1.42, 2006/08/27 22:17:27 |
|
|
* 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.41 2006/08/08 02:39:17 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
|
|
Line 1121 int compd(VL vl,DP p1,DP p2) |
|
Line 1121 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 1498 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 1799 int ni_next(int *a,int n) |
|
Line 1803 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 1832 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 1861 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 1869 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 1882 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 1893 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 1941 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 1955 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 1968 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 1987 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 2005 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 2028 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 2039 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 2089 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 2106 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 2172 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 2203 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 2225 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 2254 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 2279 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 2289 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 2312 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 2350 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 2360 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 2407 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; |