| version 1.21, 2002/01/28 00:54:41 |
version 1.23, 2003/01/04 09:06:15 |
|
|
| * 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/builtin/dp-supp.c,v 1.20 2001/10/09 01:36:05 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/dp-supp.c,v 1.22 2002/12/27 07:37:57 noro Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #include "base.h" |
#include "base.h" |
| Line 371 void dp_ptozp2_d(DP p0,DP p1,DP *hp,DP *rp) |
|
| Line 371 void dp_ptozp2_d(DP p0,DP p1,DP *hp,DP *rp) |
|
| *hp = h; *rp = r; |
*hp = h; *rp = r; |
| } |
} |
| |
|
| |
int have_sf_coef(P p) |
| |
{ |
| |
DCP dc; |
| |
|
| |
if ( !p ) |
| |
return 0; |
| |
else if ( NUM(p) ) |
| |
return NID((Num)p) == N_GFS ? 1 : 0; |
| |
else { |
| |
for ( dc = DC(p); dc; dc = NEXT(dc) ) |
| |
if ( have_sf_coef(COEF(dc)) ) |
| |
return 1; |
| |
return 0; |
| |
} |
| |
} |
| |
|
| void dp_prim(DP p,DP *rp) |
void dp_prim(DP p,DP *rp) |
| { |
{ |
| P t,g; |
P t,g; |
| Line 383 void dp_prim(DP p,DP *rp) |
|
| Line 399 void dp_prim(DP p,DP *rp) |
|
| |
|
| if ( !p ) |
if ( !p ) |
| *rp = 0; |
*rp = 0; |
| else if ( dp_fcoeffs ) |
else if ( dp_fcoeffs == N_GFS ) { |
| |
for ( m = BDY(p); m; m = NEXT(m) ) |
| |
if ( OID(m->c) == O_N ) { |
| |
/* GCD of coeffs = 1 */ |
| |
*rp = p; |
| |
return; |
| |
} else break; |
| |
/* compute GCD over the finite fieid */ |
| |
for ( m = BDY(p), n = 0; m; m = NEXT(m), n++ ); |
| |
w = (P *)ALLOCA(n*sizeof(P)); |
| |
for ( m = BDY(p), i = 0; i < n; m = NEXT(m), i++ ) |
| |
w[i] = m->c; |
| |
gcdsf(CO,w,n,&g); |
| |
if ( NUM(g) ) |
| |
*rp = p; |
| |
else { |
| |
for ( mr0 = 0, m = BDY(p); m; m = NEXT(m) ) { |
| |
NEXTMP(mr0,mr); divsp(CO,m->c,g,&mr->c); mr->dl = m->dl; |
| |
} |
| |
NEXT(mr) = 0; MKDP(p->nv,mr0,*rp); (*rp)->sugar = p->sugar; |
| |
} |
| |
return; |
| |
} else if ( dp_fcoeffs ) |
| *rp = p; |
*rp = p; |
| else if ( NoGCD ) |
else if ( NoGCD ) |
| dp_ptozp(p,rp); |
dp_ptozp(p,rp); |
| Line 688 void dp_red(DP p0,DP p1,DP p2,DP *head,DP *rest,P *dnp |
|
| Line 726 void dp_red(DP p0,DP p1,DP p2,DP *head,DP *rest,P *dnp |
|
| Q c,c1,c2; |
Q c,c1,c2; |
| N gn,tn; |
N gn,tn; |
| P g,a; |
P g,a; |
| |
P p[2]; |
| |
|
| n = p1->nv; d1 = BDY(p1)->dl; d2 = BDY(p2)->dl; |
n = p1->nv; d1 = BDY(p1)->dl; d2 = BDY(p2)->dl; |
| NEWDL(d,n); d->td = d1->td - d2->td; |
NEWDL(d,n); d->td = d1->td - d2->td; |
| for ( i = 0; i < n; i++ ) |
for ( i = 0; i < n; i++ ) |
| d->d[i] = d1->d[i]-d2->d[i]; |
d->d[i] = d1->d[i]-d2->d[i]; |
| c1 = (Q)BDY(p1)->c; c2 = (Q)BDY(p2)->c; |
c1 = (Q)BDY(p1)->c; c2 = (Q)BDY(p2)->c; |
| if ( dp_fcoeffs ) { |
if ( dp_fcoeffs == N_GFS ) { |
| |
p[0] = (P)c1; p[1] = (P)c2; |
| |
gcdsf(CO,p,2,&g); |
| |
divsp(CO,(P)c1,g,&a); c1 = (Q)a; divsp(CO,(P)c2,g,&a); c2 = (Q)a; |
| |
} else if ( dp_fcoeffs ) { |
| /* do nothing */ |
/* do nothing */ |
| } else if ( INT(c1) && INT(c2) ) { |
} else if ( INT(c1) && INT(c2) ) { |
| gcdn(NM(c1),NM(c2),&gn); |
gcdn(NM(c1),NM(c2),&gn); |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to dp-supp.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / builtin