OpenXM_contrib2/asir2000/builtin/poly.c - diff

Return to poly.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / builtin

Diff for /OpenXM_contrib2/asir2000/builtin/poly.c between version 1.20 and 1.22

version 1.20, 2003/06/21 02:09:16

version 1.22, 2011/03/30 02:43:18

Line 45

* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,

* 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.21 2003/06/24 09:49:35 noro Exp $

#include "ca.h"

#include "parse.h"

Line 54

void Pranp();

void Pheadsgn();

void Pmul_trunc();

void Pmul_trunc(),Pquo_trunc();

void Pumul(),Pumul_ff(),Pusquare(),Pusquare_ff(),Putmul(),Putmul_ff();

void Pkmul(),Pksquare(),Pktmul();

void Pord(), Pcoef0(), Pcoef(), Pdeg(), Pmindeg(), Psetmod();

void Pord(), Premove_vars(), Pcoef0(), Pcoef(), Pdeg(), Pmindeg(), Psetmod();

void Pcoef_gf2n();

void getcoef(), getdeglist(), mergedeglist(), change_mvar(), restore_mvar();

Line 114 int current_ff;

struct ftab poly_tab[] = {

{"headsgn",Pheadsgn,1},

{"quo_trunc",Pquo_trunc,2},

{"mul_trunc",Pmul_trunc,3},

{"homogeneous_deg",Phomogeneous_deg,-2},

{"homogeneous_part",Phomogeneous_part,-3},

Line 123 struct ftab poly_tab[] = {

Line 124 struct ftab poly_tab[] = {

{"p_mag",Pp_mag,1},

{"maxblen",Pmaxblen,1},

{"ord",Pord,-1},

{"remove_vars",Premove_vars,1},

{"coef0",Pcoef0,-3},

{"coef",Pcoef,-3},

{"coef_gf2n",Pcoef_gf2n,2},

Line 262 void Pmul_trunc(NODE arg,P *rp)

Line 264 void Pmul_trunc(NODE arg,P *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)

{

if ( argc(arg) == 2 )

Line 458 void Pord(NODE arg,LIST *listp)

Line 491 void Pord(NODE arg,LIST *listp)

NEXTNODE(n,tn); MKV(VR(vl),t); BDY(tn) = (pointer)t;

}

NEXT(tn) = 0; MKLIST(l,n); *listp = l;

}

void Premove_vars(NODE arg,LIST *listp)

{

NODE l,nd,tnd;

V *v,*va;

int n,na,i,j;

VL vl,vl1;

P t;

LIST list;

asir_assert(ARG0(arg),O_LIST,"remove_vars");

l = BDY((LIST)ARG0(arg)); n = length(l);

v = (V *)ALLOCA(n*sizeof(V));

for ( i = 0; i < n; i++, l = NEXT(l) )

if ( !(t = (P)BDY(l)) || (OID(t) != O_P) )

error("ord : invalid argument");

else v[i] = VR(t);

for ( na = 0, vl = CO; vl; vl = NEXT(vl), na++ );

va = (V *)ALLOCA(na*sizeof(V));

for ( i = 0, vl = CO; i < na; i++, vl = NEXT(vl) ) va[i] = VR(vl);

for ( i = 0; i < na; i++ )

for ( j = 0; j < n; j++ ) if ( va[i] == v[j] ) va[i] = 0;

for ( vl = 0, i = na-1; i >= 0; i-- )

if ( va[i] ) {

NEWVL(vl1); VR(vl1) = va[i]; NEXT(vl1) = vl; vl = vl1;

}

CO = vl;

for ( nd = 0, vl = CO; vl; vl = NEXT(vl) ) {

NEXTNODE(nd,tnd); MKV(VR(vl),t); BDY(tnd) = (pointer)t;

}

if ( nd ) NEXT(tnd) = 0;

MKLIST(list,nd); *listp = list;

}

void Pcoef0(NODE arg,Obj *rp)

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>