OpenXM_contrib2/asir2000/builtin/dp.c - diff

Return to dp.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / builtin

Diff for /OpenXM_contrib2/asir2000/builtin/dp.c between version 1.9 and 1.12

version 1.9, 2000/12/08 02:39:05

version 1.12, 2000/12/11 02:00:40

Line 45

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

* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.

* $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.8 2000/12/05 06:59:15 noro Exp $

* $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.11 2000/12/08 08:26:08 noro Exp $

#include "ca.h"

#include "base.h"

Line 57 extern int dp_order_pair_length;

extern struct order_pair *dp_order_pair;

extern struct order_spec dp_current_spec;

int do_weyl;

void Pdp_ord(), Pdp_ptod(), Pdp_dtop();

void Pdp_ptozp(), Pdp_ptozp2(), Pdp_red(), Pdp_red2(), Pdp_lcm(), Pdp_redble();

Line 76 void Pdp_mbase(),Pdp_lnf_mod(),Pdp_nf_tab_mod(),Pdp_md

Line 77 void Pdp_mbase(),Pdp_lnf_mod(),Pdp_nf_tab_mod(),Pdp_md

void Pdp_vtoe(), Pdp_etov(), Pdp_dtov(), Pdp_idiv(), Pdp_sep();

void Pdp_cont();

void Pdp_weyl_red(),Pdp_weyl_sp(),Pdp_weyl_nf();

void Pdp_weyl_gr_main(),Pdp_weyl_gr_mod_main();

void Pdp_weyl_f4_main(),Pdp_weyl_f4_mod_main();

void Pdp_weyl_mul();

struct ftab dp_tab[] = {

/* content reduction */

{"dp_ptozp",Pdp_ptozp,1},

Line 84 struct ftab dp_tab[] = {

Line 90 struct ftab dp_tab[] = {

{"dp_red_coef",Pdp_red_coef,2},

{"dp_cont",Pdp_cont,1},

/* polynomial ring */

/* s-poly */

{"dp_sp",Pdp_sp,2},

{"dp_sp_mod",Pdp_sp_mod,3},

Line 108 struct ftab dp_tab[] = {

Line 115 struct ftab dp_tab[] = {

{"dp_f4_main",Pdp_f4_main,3},

{"dp_f4_mod_main",Pdp_f4_mod_main,4},

/* weyl algebra */

/* multiplication */

{"dp_weyl_mul",Pdp_weyl_mul,2},

/* s-poly */

{"dp_weyl_sp",Pdp_weyl_sp,2},

/* m-reduction */

{"dp_weyl_red",Pdp_weyl_red,3},

/* normal form */

{"dp_weyl_nf",Pdp_weyl_nf,4},

/* Buchberger algorithm */

{"dp_weyl_gr_main",Pdp_weyl_gr_main,5},

{"dp_weyl_gr_mod_main",Pdp_weyl_gr_mod_main,5},

/* F4 algorithm */

{"dp_weyl_f4_main",Pdp_weyl_f4_main,3},

{"dp_weyl_f4_mod_main",Pdp_weyl_f4_mod_main,4},

{0,0,0},

};

Line 380 NODE arg;

Line 408 NODE arg;

DP *rp;

{

asir_assert(ARG0(arg),O_DP,"dp_ptozp");

if ( Dist )

dp_ptozp((DP)ARG0(arg),rp);

dp_ptozp_d(BDY(Dist),length(BDY(Dist)),(DP)ARG0(arg),rp);

else

dp_ptozp((DP)ARG0(arg),rp);

}

void Pdp_ptozp2(arg,rp)

Line 396 LIST *rp;

Line 421 LIST *rp;

p0 = (DP)ARG0(arg); p1 = (DP)ARG1(arg);

asir_assert(p0,O_DP,"dp_ptozp2");

asir_assert(p1,O_DP,"dp_ptozp2");

if ( Dist )

dp_ptozp2(p0,p1,&h,&r);

dp_ptozp2_d(BDY(Dist),length(BDY(Dist)),p0,p1,&h,&r);

else

dp_ptozp2(p0,p1,&h,&r);

NEWNODE(n0); BDY(n0) = (pointer)h;

NEWNODE(NEXT(n0)); BDY(NEXT(n0)) = (pointer)r;

NEXT(NEXT(n0)) = 0;

Line 451 DP *rp;

Line 473 DP *rp;

DP g;

int full;

do_weyl = 0;

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

asir_assert(ARG1(arg),O_DP,"dp_nf");

asir_assert(ARG2(arg),O_VECT,"dp_nf");

Line 463 DP *rp;

Line 486 DP *rp;

dp_nf_ptozp(b,g,ps,full,DP_Multiple,rp);

}

void Pdp_weyl_nf(arg,rp)

NODE arg;

DP *rp;

{

NODE b;

DP *ps;

DP g;

int full;

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

asir_assert(ARG1(arg),O_DP,"dp_weyl_nf");

asir_assert(ARG2(arg),O_VECT,"dp_weyl_nf");

asir_assert(ARG3(arg),O_N,"dp_weyl_nf");

if ( !(g = (DP)ARG1(arg)) ) {

*rp = 0; return;

}

b = BDY((LIST)ARG0(arg)); ps = (DP *)BDY((VECT)ARG2(arg));

full = (Q)ARG3(arg) ? 1 : 0;

do_weyl = 1;

dp_nf_ptozp(b,g,ps,full,DP_Multiple,rp);

do_weyl = 0;

}

void Pdp_true_nf(arg,rp)

NODE arg;

LIST *rp;

Line 474 LIST *rp;

Line 520 LIST *rp;

P dn;

int full;

do_weyl = 0;

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

asir_assert(ARG1(arg),O_DP,"dp_true_nf");

asir_assert(ARG2(arg),O_VECT,"dp_true_nf");

Line 500 DP *rp;

Line 547 DP *rp;

int mod,full,ac;

NODE n,n0;

do_weyl = 0;

ac = argc(arg);

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

asir_assert(ARG1(arg),O_DP,"dp_nf_mod");

Line 531 LIST *rp;

Line 579 LIST *rp;

int mod,full;

NODE n;

do_weyl = 0;

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

asir_assert(ARG1(arg),O_DP,"dp_nf_mod");

asir_assert(ARG2(arg),O_VECT,"dp_nf_mod");

Line 625 LIST *rp;

Line 674 LIST *rp;

P dmy;

NODE n;

do_weyl = 0;

asir_assert(ARG0(arg),O_DP,"dp_red_mod");

asir_assert(ARG1(arg),O_DP,"dp_red_mod");

asir_assert(ARG2(arg),O_DP,"dp_red_mod");

Line 647 DP *rp;

Line 697 DP *rp;

dp_subd(p1,p2,rp);

}

void Pdp_weyl_mul(arg,rp)

NODE arg;

DP *rp;

{

DP p1,p2;

p1 = (DP)ARG0(arg); p2 = (DP)ARG1(arg);

asir_assert(p1,O_DP,"dp_weyl_mul"); asir_assert(p2,O_DP,"dp_mul");

do_weyl = 1;

muld(CO,p1,p2,rp);

do_weyl = 0;

}

void Pdp_red(arg,rp)

NODE arg;

LIST *rp;

Line 655 LIST *rp;

Line 718 LIST *rp;

DP head,rest,dmy1;

P dmy;

do_weyl = 0;

asir_assert(ARG0(arg),O_DP,"dp_red");

asir_assert(ARG1(arg),O_DP,"dp_red");

asir_assert(ARG2(arg),O_DP,"dp_red");

Line 664 LIST *rp;

Line 728 LIST *rp;

NEXT(NEXT(n)) = 0; MKLIST(*rp,n);

}

void Pdp_weyl_red(arg,rp)

NODE arg;

LIST *rp;

{

NODE n;

DP head,rest,dmy1;

P dmy;

asir_assert(ARG0(arg),O_DP,"dp_weyl_red");

asir_assert(ARG1(arg),O_DP,"dp_weyl_red");

asir_assert(ARG2(arg),O_DP,"dp_weyl_red");

do_weyl = 1;

dp_red((DP)ARG0(arg),(DP)ARG1(arg),(DP)ARG2(arg),&head,&rest,&dmy,&dmy1);

do_weyl = 0;

NEWNODE(n); BDY(n) = (pointer)head;

NEWNODE(NEXT(n)); BDY(NEXT(n)) = (pointer)rest;

NEXT(NEXT(n)) = 0; MKLIST(*rp,n);

}

void Pdp_sp(arg,rp)

NODE arg;

DP *rp;

{

DP p1,p2;

do_weyl = 0;

p1 = (DP)ARG0(arg); p2 = (DP)ARG1(arg);

asir_assert(p1,O_DP,"dp_sp"); asir_assert(p2,O_DP,"dp_sp");

dp_sp(p1,p2,rp);

}

void Pdp_weyl_sp(arg,rp)

NODE arg;

DP *rp;

{

DP p1,p2;

p1 = (DP)ARG0(arg); p2 = (DP)ARG1(arg);

asir_assert(p1,O_DP,"dp_weyl_sp"); asir_assert(p2,O_DP,"dp_sp");

do_weyl = 1;

dp_sp(p1,p2,rp);

do_weyl = 0;

}

void Pdp_sp_mod(arg,rp)

NODE arg;

DP *rp;

Line 682 DP *rp;

Line 779 DP *rp;

DP p1,p2;

int mod;

do_weyl = 0;

p1 = (DP)ARG0(arg); p2 = (DP)ARG1(arg);

asir_assert(p1,O_DP,"dp_sp_mod"); asir_assert(p2,O_DP,"dp_sp_mod");

asir_assert(ARG2(arg),O_N,"dp_sp_mod");

Line 1012 LIST *rp;

Line 1110 LIST *rp;

int modular;

struct order_spec ord;

do_weyl = 0;

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

asir_assert(ARG1(arg),O_LIST,"dp_gr_main");

asir_assert(ARG2(arg),O_N,"dp_gr_main");

Line 1036 LIST *rp;

Line 1135 LIST *rp;

LIST f,v;

struct order_spec ord;

do_weyl = 0;

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

asir_assert(ARG1(arg),O_LIST,"dp_f4_main");

f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);

Line 1051 LIST *rp;

Line 1151 LIST *rp;

int m;

struct order_spec ord;

do_weyl = 0;

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

asir_assert(ARG1(arg),O_LIST,"dp_f4_main");

asir_assert(ARG2(arg),O_N,"dp_f4_main");

Line 1068 LIST *rp;

Line 1169 LIST *rp;

int m;

struct order_spec ord;

do_weyl = 0;

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

asir_assert(ARG1(arg),O_LIST,"dp_gr_mod_main");

asir_assert(ARG2(arg),O_N,"dp_gr_mod_main");

Line 1076 LIST *rp;

Line 1178 LIST *rp;

homo = (Num)ARG2(arg); m = QTOS((Q)ARG3(arg));

create_order_spec(ARG4(arg),&ord);

dp_gr_mod_main(f,v,homo,m,&ord,rp);

}

/* for Weyl algebra */

void Pdp_weyl_gr_main(arg,rp)

NODE arg;

LIST *rp;

{

LIST f,v;

Num homo;

Q m;

int modular;

struct order_spec ord;

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

asir_assert(ARG1(arg),O_LIST,"dp_weyl_gr_main");

asir_assert(ARG2(arg),O_N,"dp_weyl_gr_main");

asir_assert(ARG3(arg),O_N,"dp_weyl_gr_main");

f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);

homo = (Num)ARG2(arg);

m = (Q)ARG3(arg);

if ( !m )

modular = 0;

else if ( PL(NM(m))>1 || (PL(NM(m)) == 1 && BD(NM(m))[0] >= 0x80000000) )

error("dp_gr_main : too large modulus");

else

modular = QTOS(m);

create_order_spec(ARG4(arg),&ord);

do_weyl = 1;

dp_gr_main(f,v,homo,modular,&ord,rp);

do_weyl = 0;

}

void Pdp_weyl_f4_main(arg,rp)

NODE arg;

LIST *rp;

{

LIST f,v;

struct order_spec ord;

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

asir_assert(ARG1(arg),O_LIST,"dp_weyl_f4_main");

f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);

create_order_spec(ARG2(arg),&ord);

do_weyl = 1;

dp_f4_main(f,v,&ord,rp);

do_weyl = 0;

}

void Pdp_weyl_f4_mod_main(arg,rp)

NODE arg;

LIST *rp;

{

LIST f,v;

int m;

struct order_spec ord;

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

asir_assert(ARG1(arg),O_LIST,"dp_weyl_f4_main");

asir_assert(ARG2(arg),O_N,"dp_f4_main");

f = (LIST)ARG0(arg); v = (LIST)ARG1(arg); m = QTOS((Q)ARG2(arg));

create_order_spec(ARG3(arg),&ord);

do_weyl = 1;

dp_f4_mod_main(f,v,m,&ord,rp);

do_weyl = 0;

}

void Pdp_weyl_gr_mod_main(arg,rp)

NODE arg;

LIST *rp;

{

LIST f,v;

Num homo;

int m;

struct order_spec ord;

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

asir_assert(ARG1(arg),O_LIST,"dp_weyl_gr_mod_main");

asir_assert(ARG2(arg),O_N,"dp_weyl_gr_mod_main");

asir_assert(ARG3(arg),O_N,"dp_weyl_gr_mod_main");

f = (LIST)ARG0(arg); v = (LIST)ARG1(arg);

homo = (Num)ARG2(arg); m = QTOS((Q)ARG3(arg));

create_order_spec(ARG4(arg),&ord);

do_weyl = 1;

dp_gr_mod_main(f,v,homo,m,&ord,rp);

do_weyl = 0;

}

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