version 1.20, 2003/06/21 02:09:16 |
version 1.21, 2003/06/24 09:49:35 |
|
|
* 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/poly.c,v 1.19 2003/06/19 07:08:18 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/poly.c,v 1.20 2003/06/21 02:09:16 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
|
|
void Pranp(); |
void Pranp(); |
|
|
void Pheadsgn(); |
void Pheadsgn(); |
void Pmul_trunc(); |
void Pmul_trunc(),Pquo_trunc(); |
void Pumul(),Pumul_ff(),Pusquare(),Pusquare_ff(),Putmul(),Putmul_ff(); |
void Pumul(),Pumul_ff(),Pusquare(),Pusquare_ff(),Putmul(),Putmul_ff(); |
void Pkmul(),Pksquare(),Pktmul(); |
void Pkmul(),Pksquare(),Pktmul(); |
void Pord(), Pcoef0(), Pcoef(), Pdeg(), Pmindeg(), Psetmod(); |
void Pord(), Pcoef0(), Pcoef(), Pdeg(), Pmindeg(), Psetmod(); |
|
|
|
|
struct ftab poly_tab[] = { |
struct ftab poly_tab[] = { |
{"headsgn",Pheadsgn,1}, |
{"headsgn",Pheadsgn,1}, |
|
{"quo_trunc",Pquo_trunc,2}, |
{"mul_trunc",Pmul_trunc,3}, |
{"mul_trunc",Pmul_trunc,3}, |
{"homogeneous_deg",Phomogeneous_deg,-2}, |
{"homogeneous_deg",Phomogeneous_deg,-2}, |
{"homogeneous_part",Phomogeneous_part,-3}, |
{"homogeneous_part",Phomogeneous_part,-3}, |
Line 260 void Pmul_trunc(NODE arg,P *rp) |
|
Line 261 void Pmul_trunc(NODE arg,P *rp) |
|
vn[i].n = QTOS(DEG(DC(h))); |
vn[i].n = QTOS(DEG(DC(h))); |
} |
} |
mulp_trunc(vl,p1,p2,vn,rp); |
mulp_trunc(vl,p1,p2,vn,rp); |
|
} |
|
|
|
void Pquo_trunc(NODE arg,P *rp) |
|
{ |
|
P p1,p2,p,h; |
|
VL vl0,vl1,vl2,tvl,vl; |
|
VN vn; |
|
int i,n; |
|
|
|
p1 = (P)ARG0(arg); |
|
p2 = (P)ARG1(arg); |
|
if ( !p1 ) |
|
*rp = 0; |
|
else if ( NUM(p2) ) |
|
divsp(CO,p1,p2,rp); |
|
else { |
|
get_vars((Obj)p1,&vl1); get_vars((Obj)p2,&vl2); mergev(CO,vl1,vl2,&vl); |
|
for ( tvl = vl, n = 0; tvl; tvl = NEXT(tvl), n++ ); |
|
vn = (VN) ALLOCA((n+1)*sizeof(struct oVN)); |
|
for ( i = 0, tvl = vl; i < n; tvl = NEXT(tvl), i++ ) { |
|
vn[i].v = tvl->v; |
|
vn[i].n = 0; |
|
} |
|
vn[i].v = 0; |
|
vn[i].n = 0; |
|
for ( h = p2, i = 0; OID(h) == O_P; h = COEF(DC(h)) ) { |
|
for ( ; vn[i].v != VR(h); i++ ); |
|
vn[i].n = QTOS(DEG(DC(h))); |
|
} |
|
quop_trunc(vl,p1,p2,vn,rp); |
|
} |
} |
} |
|
|
void Phomogeneous_part(NODE arg,P *rp) |
void Phomogeneous_part(NODE arg,P *rp) |