version 1.1, 2004/08/18 00:17:02 |
version 1.2, 2018/03/29 01:32:52 |
|
|
P p1,p2; |
P p1,p2; |
P *q,*r; |
P *q,*r; |
{ |
{ |
register int i,j; |
register int i,j; |
register DCP dc1,dc2,dc; |
register DCP dc1,dc2,dc; |
P m,s,dvr; |
P m,s,dvr; |
P *pq,*pr,*pd; |
P *pq,*pr,*pd; |
V v1,v2; |
V v1,v2; |
Q deg1,deg2; |
Q deg1,deg2; |
int d1,d2,sgn; |
int d1,d2,sgn; |
|
|
if ( !p1 ) { |
if ( !p1 ) { |
*q = 0; *r = 0; |
*q = 0; *r = 0; |
} else if ( NUM(p2) ) |
} else if ( NUM(p2) ) |
if ( NUM(p1) ) { |
if ( NUM(p1) ) { |
DIVNUM(p1,p2,q); *r = 0; |
DIVNUM(p1,p2,q); *r = 0; |
} else { |
} else { |
DIVSDCP(vl,p1,p2,q); *r = 0; |
DIVSDCP(vl,p1,p2,q); *r = 0; |
} |
} |
else if ( NUM(p1) ) { |
else if ( NUM(p1) ) { |
*q = 0; *r = p1; |
*q = 0; *r = p1; |
} else if ( ( v1 = VR(p1) ) == ( v2 = VR(p2) ) ) { |
} else if ( ( v1 = VR(p1) ) == ( v2 = VR(p2) ) ) { |
dc1 = DC(p1); dc2 = DC(p2); |
dc1 = DC(p1); dc2 = DC(p2); |
deg1 = DEG(dc1); deg2 = DEG(dc2); |
deg1 = DEG(dc1); deg2 = DEG(dc2); |
sgn = cmpq(deg1,deg2); |
sgn = cmpq(deg1,deg2); |
if ( sgn == 0 ) { |
if ( sgn == 0 ) { |
DIVSP(vl,COEF(dc1),COEF(dc2),q); |
DIVSP(vl,COEF(dc1),COEF(dc2),q); |
MULP(vl,p2,*q,&m); SUBP(vl,p1,m,r); |
MULP(vl,p2,*q,&m); SUBP(vl,p1,m,r); |
} else if ( sgn < 0 ) { |
} else if ( sgn < 0 ) { |
*q = 0; *r = p1; |
*q = 0; *r = p1; |
} else { |
} else { |
if ( (PL(NM(deg1)) > 1) ) |
if ( (PL(NM(deg1)) > 1) ) |
error("divsrp : invalid input"); |
error("divsrp : invalid input"); |
d1 = QTOS(deg1); d2 = QTOS(deg2); |
d1 = QTOS(deg1); d2 = QTOS(deg2); |
W_CALLOC(d1-d2,P,pq); W_CALLOC(d1,P,pr); W_CALLOC(d2,P,pd); |
W_CALLOC(d1-d2,P,pq); W_CALLOC(d1,P,pr); W_CALLOC(d2,P,pd); |
for ( dc = dc1; dc; dc = NEXT(dc) ) |
for ( dc = dc1; dc; dc = NEXT(dc) ) |
pr[QTOS(DEG(dc))] = COEF(dc); |
pr[QTOS(DEG(dc))] = COEF(dc); |
for ( dc = dc2; dc; dc = NEXT(dc) ) |
for ( dc = dc2; dc; dc = NEXT(dc) ) |
pd[QTOS(DEG(dc))] = COEF(dc); |
pd[QTOS(DEG(dc))] = COEF(dc); |
for ( dvr = COEF(dc2), i = d1 - d2; i >= 0; i-- ) { |
for ( dvr = COEF(dc2), i = d1 - d2; i >= 0; i-- ) { |
if ( !pr[i+d2] ) |
if ( !pr[i+d2] ) |
continue; |
continue; |
DIVSP(vl,pr[i+d2],dvr,&pq[i]); |
DIVSP(vl,pr[i+d2],dvr,&pq[i]); |
for ( j = d2; j >= 0; j-- ) { |
for ( j = d2; j >= 0; j-- ) { |
MULP(vl,pq[i],pd[j],&m); |
MULP(vl,pq[i],pd[j],&m); |
SUBP(vl,pr[i + j],m,&s); pr[i + j] = s; |
SUBP(vl,pr[i + j],m,&s); pr[i + j] = s; |
} |
} |
} |
} |
plisttop(pq,v1,d1 - d2,q); plisttop(pr,v1,d1 - 1,r); |
plisttop(pq,v1,d1 - d2,q); plisttop(pr,v1,d1 - 1,r); |
} |
} |
} else { |
} else { |
for ( ; (v1 != vl->v) && (v2 != vl->v); vl = NEXT(vl) ); |
for ( ; (v1 != vl->v) && (v2 != vl->v); vl = NEXT(vl) ); |
if ( v2 == vl->v ) { |
if ( v2 == vl->v ) { |
*q = 0; *r = p1; |
*q = 0; *r = p1; |
} else |
} else |
DIVSRDCP(vl,p1,p2,q,r); |
DIVSRDCP(vl,p1,p2,q,r); |
} |
} |
} |
} |
|
|
void D_DIVSRDCP(vl,p1,p2,q,r) |
void D_DIVSRDCP(vl,p1,p2,q,r) |
|
|
P p1,p2,*q,*r; |
P p1,p2,*q,*r; |
{ |
{ |
|
|
P qc,rc; |
P qc,rc; |
DCP dc,dcq,dcq0,dcr,dcr0; |
DCP dc,dcq,dcq0,dcr,dcr0; |
|
|
for ( dc = DC(p1), dcq0 = 0, dcr0 = 0; dc; dc = NEXT(dc) ) { |
for ( dc = DC(p1), dcq0 = 0, dcr0 = 0; dc; dc = NEXT(dc) ) { |
DIVSRP(vl,COEF(dc),p2,&qc,&rc); |
DIVSRP(vl,COEF(dc),p2,&qc,&rc); |
if ( qc ) { |
if ( qc ) { |
NEXTDC(dcq0,dcq); DEG(dcq) = DEG(dc); COEF(dcq) = qc; |
NEXTDC(dcq0,dcq); DEG(dcq) = DEG(dc); COEF(dcq) = qc; |
} |
} |
if ( rc ) { |
if ( rc ) { |
NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = rc; |
NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = rc; |
} |
} |
} |
} |
if ( dcq0 ) { |
if ( dcq0 ) { |
NEXT(dcq) = 0; MKP(VR(p1),dcq0,*q); |
NEXT(dcq) = 0; MKP(VR(p1),dcq0,*q); |
} else |
} else |
*q = 0; |
*q = 0; |
if ( dcr0 ) { |
if ( dcr0 ) { |
NEXT(dcr) = 0; MKP(VR(p1),dcr0,*r); |
NEXT(dcr) = 0; MKP(VR(p1),dcr0,*r); |
} else |
} else |
*r = 0; |
*r = 0; |
} |
} |
|
|
void D_DIVSP(vl,p1,p2,q) |
void D_DIVSP(vl,p1,p2,q) |
VL vl; |
VL vl; |
P p1,p2,*q; |
P p1,p2,*q; |
{ |
{ |
P t; |
P t; |
|
|
DIVSRP(vl,p1,p2,q,&t); |
DIVSRP(vl,p1,p2,q,&t); |
if ( t ) |
if ( t ) |
error("divsp: cannot happen"); |
error("divsp: cannot happen"); |
} |
} |
|
|
void D_DIVSDCP(vl,p1,p2,q) |
void D_DIVSDCP(vl,p1,p2,q) |
|
|
P p1,p2,*q; |
P p1,p2,*q; |
{ |
{ |
|
|
P m; |
P m; |
register DCP dc,dcr,dcr0; |
register DCP dc,dcr,dcr0; |
|
|
for ( dc = DC(p1), dcr0 = 0; dc; dc = NEXT(dc) ) { |
for ( dc = DC(p1), dcr0 = 0; dc; dc = NEXT(dc) ) { |
DIVSP(vl,COEF(dc),p2,&m); |
DIVSP(vl,COEF(dc),p2,&m); |
NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = m; NEXT(dcr) = 0; |
NEXTDC(dcr0,dcr); DEG(dcr) = DEG(dc); COEF(dcr) = m; NEXT(dcr) = 0; |
} |
} |
MKP(VR(p1),dcr0,*q); |
MKP(VR(p1),dcr0,*q); |
} |
} |
|
|