version 1.8, 2019/03/18 10:30:41 |
version 1.12, 2019/09/04 05:32:10 |
|
|
* 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/asir2018/builtin/dp.c,v 1.7 2019/03/18 07:09:58 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2018/builtin/dp.c,v 1.11 2019/09/04 01:12:02 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 89 void Pdp_vtoe(), Pdp_etov(), Pdp_dtov(), Pdp_idiv(), P |
|
Line 89 void Pdp_vtoe(), Pdp_etov(), Pdp_dtov(), Pdp_idiv(), P |
|
void Pdp_cont(); |
void Pdp_cont(); |
void Pdp_gr_checklist(); |
void Pdp_gr_checklist(); |
void Pdp_ltod(),Pdpv_ord(),Pdpv_ht(),Pdpv_hm(),Pdpv_hc(); |
void Pdp_ltod(),Pdpv_ord(),Pdpv_ht(),Pdpv_hm(),Pdpv_hc(); |
void Pdpm_ltod(),Pdpm_dtol(),Pdpm_ord(),Pdpm_nf(),Pdpm_weyl_nf(),Pdpm_sp(),Pdpm_weyl_sp(); |
void Pdpm_ltod(),Pdpm_dtol(),Pdpm_set_schreyer(),Pdpm_nf(),Pdpm_weyl_nf(),Pdpm_sp(),Pdpm_weyl_sp(),Pdpm_nf_and_quotient(); |
void Pdpm_hm(),Pdpm_ht(),Pdpm_hc(); |
void Pdpm_hm(),Pdpm_ht(),Pdpm_hc(),Pdpm_hp(),Pdpm_rest(),Pdpm_shift(),Pdpm_split(),Pdpm_sort(),Pdpm_dptodpm(),Pdpm_redble(); |
|
void Pdpm_schreyer_base(),Pdpm_simplify_syz(); |
|
|
void Pdp_weyl_red(); |
void Pdp_weyl_red(); |
void Pdp_weyl_sp(); |
void Pdp_weyl_sp(); |
Line 157 struct ftab dp_tab[] = { |
|
Line 158 struct ftab dp_tab[] = { |
|
{"dp_nf",Pdp_nf,4}, |
{"dp_nf",Pdp_nf,4}, |
{"dp_nf_mod",Pdp_nf_mod,5}, |
{"dp_nf_mod",Pdp_nf_mod,5}, |
{"dp_nf_f",Pdp_nf_f,4}, |
{"dp_nf_f",Pdp_nf_f,4}, |
{"dpm_nf_f",Pdpm_nf_f,4}, |
{"dpm_nf_and_quotient",Pdpm_nf_and_quotient,-3}, |
{"dpm_weyl_nf_f",Pdpm_weyl_nf_f,4}, |
{"dpm_nf_f",Pdpm_nf_f,-4}, |
{"dpm_nf",Pdpm_nf,4}, |
{"dpm_weyl_nf_f",Pdpm_weyl_nf_f,-4}, |
|
{"dpm_nf",Pdpm_nf,-4}, |
{"dpm_sp",Pdpm_sp,2}, |
{"dpm_sp",Pdpm_sp,2}, |
{"dpm_weyl_sp",Pdpm_weyl_sp,2}, |
{"dpm_weyl_sp",Pdpm_weyl_sp,2}, |
|
|
Line 215 struct ftab dp_tab[] = { |
|
Line 217 struct ftab dp_tab[] = { |
|
|
|
/* normal form */ |
/* normal form */ |
{"dp_weyl_nf",Pdp_weyl_nf,4}, |
{"dp_weyl_nf",Pdp_weyl_nf,4}, |
{"dpm_weyl_nf",Pdpm_weyl_nf,4}, |
{"dpm_weyl_nf",Pdpm_weyl_nf,-4}, |
{"dp_weyl_nf_mod",Pdp_weyl_nf_mod,5}, |
{"dp_weyl_nf_mod",Pdp_weyl_nf_mod,5}, |
{"dp_weyl_nf_f",Pdp_weyl_nf_f,4}, |
{"dp_weyl_nf_f",Pdp_weyl_nf_f,4}, |
|
|
Line 252 struct ftab dp_supp_tab[] = { |
|
Line 254 struct ftab dp_supp_tab[] = { |
|
/* setting flags */ |
/* setting flags */ |
{"dp_sort",Pdp_sort,1}, |
{"dp_sort",Pdp_sort,1}, |
{"dp_ord",Pdp_ord,-1}, |
{"dp_ord",Pdp_ord,-1}, |
{"dpm_ord",Pdpm_ord,-1}, |
{"dpm_set_schreyer",Pdpm_set_schreyer,-1}, |
{"dpv_ord",Pdpv_ord,-2}, |
{"dpv_ord",Pdpv_ord,-2}, |
{"dp_set_kara",Pdp_set_kara,-1}, |
{"dp_set_kara",Pdp_set_kara,-1}, |
{"dp_nelim",Pdp_nelim,-1}, |
{"dp_nelim",Pdp_nelim,-1}, |
Line 274 struct ftab dp_supp_tab[] = { |
|
Line 276 struct ftab dp_supp_tab[] = { |
|
{"dp_ltod",Pdp_ltod,-2}, |
{"dp_ltod",Pdp_ltod,-2}, |
|
|
{"dpm_ltod",Pdpm_ltod,2}, |
{"dpm_ltod",Pdpm_ltod,2}, |
|
{"dpm_dptodpm",Pdpm_dptodpm,2}, |
{"dpm_dtol",Pdpm_dtol,3}, |
{"dpm_dtol",Pdpm_dtol,3}, |
|
|
/* criteria */ |
/* criteria */ |
Line 293 struct ftab dp_supp_tab[] = { |
|
Line 296 struct ftab dp_supp_tab[] = { |
|
{"dpm_hm",Pdpm_hm,1}, |
{"dpm_hm",Pdpm_hm,1}, |
{"dpm_ht",Pdpm_ht,1}, |
{"dpm_ht",Pdpm_ht,1}, |
{"dpm_hc",Pdpm_hc,1}, |
{"dpm_hc",Pdpm_hc,1}, |
|
{"dpm_hp",Pdpm_hp,1}, |
|
{"dpm_rest",Pdpm_rest,1}, |
|
{"dpm_shift",Pdpm_shift,2}, |
|
{"dpm_split",Pdpm_split,2}, |
|
{"dpm_sort",Pdpm_sort,1}, |
{"dp_rest",Pdp_rest,1}, |
{"dp_rest",Pdp_rest,1}, |
{"dp_initial_term",Pdp_initial_term,1}, |
{"dp_initial_term",Pdp_initial_term,1}, |
{"dp_order",Pdp_order,1}, |
{"dp_order",Pdp_order,1}, |
Line 307 struct ftab dp_supp_tab[] = { |
|
Line 315 struct ftab dp_supp_tab[] = { |
|
/* misc */ |
/* misc */ |
{"dp_mbase",Pdp_mbase,1}, |
{"dp_mbase",Pdp_mbase,1}, |
{"dp_redble",Pdp_redble,2}, |
{"dp_redble",Pdp_redble,2}, |
|
{"dpm_redble",Pdpm_redble,2}, |
{"dp_sep",Pdp_sep,2}, |
{"dp_sep",Pdp_sep,2}, |
{"dp_idiv",Pdp_idiv,2}, |
{"dp_idiv",Pdp_idiv,2}, |
{"dp_tdiv",Pdp_tdiv,2}, |
{"dp_tdiv",Pdp_tdiv,2}, |
Line 316 struct ftab dp_supp_tab[] = { |
|
Line 325 struct ftab dp_supp_tab[] = { |
|
{"dp_compute_essential_df",Pdp_compute_essential_df,2}, |
{"dp_compute_essential_df",Pdp_compute_essential_df,2}, |
{"dp_mono_raddec",Pdp_mono_raddec,2}, |
{"dp_mono_raddec",Pdp_mono_raddec,2}, |
{"dp_mono_reduce",Pdp_mono_reduce,2}, |
{"dp_mono_reduce",Pdp_mono_reduce,2}, |
|
{"dpm_schreyer_base",Pdpm_schreyer_base,1}, |
|
{"dpm_simplify_syz",Pdpm_simplify_syz,2}, |
|
|
{"dp_rref2",Pdp_rref2,2}, |
{"dp_rref2",Pdp_rref2,2}, |
{"sumi_updatepairs",Psumi_updatepairs,3}, |
{"sumi_updatepairs",Psumi_updatepairs,3}, |
Line 503 P binpoly(P n,int a,int b) |
|
Line 514 P binpoly(P n,int a,int b) |
|
return r; |
return r; |
} |
} |
|
|
|
void ibin(unsigned long int n,unsigned long int k,Z *r); |
|
|
void mhp_to_hf(VL vl,P hp,int n,P *plist,VECT *head,P *hf) |
void mhp_to_hf(VL vl,P hp,int n,P *plist,VECT *head,P *hf) |
{ |
{ |
P tv,gcd,q,h,hphead,tt,ai,hpoly,nv,bp,w; |
P tv,gcd,q,h,hphead,tt,ai,hpoly,nv,bp,w; |
Line 538 void mhp_to_hf(VL vl,P hp,int n,P *plist,VECT *head,P |
|
Line 551 void mhp_to_hf(VL vl,P hp,int n,P *plist,VECT *head,P |
|
topdc = 0; |
topdc = 0; |
for ( i = 0; i < qd; i++ ) { |
for ( i = 0; i < qd; i++ ) { |
NEWDC(dc); NEXT(dc) = topdc; |
NEWDC(dc); NEXT(dc) = topdc; |
ibin(i+s-1,s-1,&COEF(dc)); |
ibin(i+s-1,s-1,(Z *)&COEF(dc)); |
STOZ(i,d); DEG(dc) = d; |
STOZ(i,d); DEG(dc) = d; |
topdc = dc; |
topdc = dc; |
} |
} |
Line 897 void Pdp_nf_tab_f(NODE arg,DP *rp) |
|
Line 910 void Pdp_nf_tab_f(NODE arg,DP *rp) |
|
dp_nf_tab_f((DP)ARG0(arg),(LIST *)BDY((VECT)ARG1(arg)),rp); |
dp_nf_tab_f((DP)ARG0(arg),(LIST *)BDY((VECT)ARG1(arg)),rp); |
} |
} |
|
|
|
extern int dpm_ordtype; |
|
|
void Pdp_ord(NODE arg,Obj *rp) |
void Pdp_ord(NODE arg,Obj *rp) |
{ |
{ |
struct order_spec *spec; |
struct order_spec *spec; |
Line 914 void Pdp_ord(NODE arg,Obj *rp) |
|
Line 929 void Pdp_ord(NODE arg,Obj *rp) |
|
else if ( !create_order_spec(0,(Obj)ARG0(arg),&spec) ) |
else if ( !create_order_spec(0,(Obj)ARG0(arg),&spec) ) |
error("dp_ord : invalid order specification"); |
error("dp_ord : invalid order specification"); |
initd(spec); *rp = spec->obj; |
initd(spec); *rp = spec->obj; |
|
if ( spec->id >= 256 ) dpm_ordtype = spec->ispot; |
} |
} |
} |
} |
|
|
Line 1064 void Pdpm_ltod(NODE arg,DPM *rp) |
|
Line 1080 void Pdpm_ltod(NODE arg,DPM *rp) |
|
|
|
nd = BDY(f); |
nd = BDY(f); |
len = length(nd); |
len = length(nd); |
for ( i = 0, t = nd, s = 0; i < len; i++, t = NEXT(t) ) { |
for ( i = 1, t = nd, s = 0; i <= len; i++, t = NEXT(t) ) { |
ptod(CO,vl,(P)BDY(t),&d); |
ptod(CO,vl,(P)BDY(t),&d); |
dtodpm(d,i,&u); |
dtodpm(d,i,&u); |
adddpm(CO,s,u,&w); s = w; |
adddpm(CO,s,u,&w); s = w; |
Line 1072 void Pdpm_ltod(NODE arg,DPM *rp) |
|
Line 1088 void Pdpm_ltod(NODE arg,DPM *rp) |
|
*rp = s; |
*rp = s; |
} |
} |
|
|
|
// c*[monomial,i]+... -> c*<<monomial:i>>+... |
|
|
|
void Pdpm_dptodpm(NODE arg,DPM *rp) |
|
{ |
|
DP p; |
|
MP mp; |
|
int pos; |
|
DMM m0,m; |
|
|
|
p = (DP)ARG0(arg); |
|
pos = ZTOS((Z)ARG1(arg)); |
|
if ( pos <= 0 ) |
|
error("dpm_mtod : position must be positive"); |
|
if ( !p ) *rp = 0; |
|
else { |
|
for ( m0 = 0, mp = BDY(p); mp; mp = NEXT(mp) ) { |
|
NEXTDMM(m0,m); m->dl = mp->dl; m->c = mp->c; m->pos = pos; |
|
} |
|
MKDPM(p->nv,m0,*rp); (*rp)->sugar = p->sugar; |
|
} |
|
} |
|
|
void Pdpm_dtol(NODE arg,LIST *rp) |
void Pdpm_dtol(NODE arg,LIST *rp) |
{ |
{ |
DPM a; |
DPM a; |
Line 1265 void Pdp_weyl_nf(NODE arg,DP *rp) |
|
Line 1303 void Pdp_weyl_nf(NODE arg,DP *rp) |
|
void Pdpm_nf(NODE arg,DPM *rp) |
void Pdpm_nf(NODE arg,DPM *rp) |
{ |
{ |
NODE b; |
NODE b; |
DPM *ps; |
VECT ps; |
DPM g; |
DPM g; |
int full; |
int ac,full; |
|
|
if ( !(g = (DPM)ARG1(arg)) ) { |
if ( !(g = (DPM)ARG1(arg)) ) { |
*rp = 0; return; |
*rp = 0; return; |
} |
} |
do_weyl = 0; dp_fcoeffs = 0; |
do_weyl = 0; dp_fcoeffs = 0; |
asir_assert(ARG0(arg),O_LIST,"dpm_nf"); |
ac = argc(arg); |
asir_assert(ARG1(arg),O_DPM,"dpm_nf"); |
if ( ac < 3 ) |
asir_assert(ARG2(arg),O_VECT,"dpm_nf"); |
error("dpm_nf: invalid arguments"); |
asir_assert(ARG3(arg),O_N,"dpm_nf"); |
else if ( ac == 3 ) { |
b = BDY((LIST)ARG0(arg)); ps = (DPM *)BDY((VECT)ARG2(arg)); |
asir_assert(ARG1(arg),O_VECT,"dpm_nf"); |
full = (Q)ARG3(arg) ? 1 : 0; |
b = 0; g = (DPM)ARG0(arg); ps = (VECT)ARG1(arg); |
|
} else if ( ac == 4 ) { |
|
asir_assert(ARG0(arg),O_LIST,"dpm_nf"); |
|
asir_assert(ARG2(arg),O_VECT,"dpm_nf"); |
|
b = BDY((LIST)ARG0(arg)); g = (DPM)ARG1(arg); ps = (VECT)ARG2(arg); |
|
full = (Q)ARG3(arg) ? 1 : 0; |
|
} |
dpm_nf_z(b,g,ps,full,DP_Multiple,rp); |
dpm_nf_z(b,g,ps,full,DP_Multiple,rp); |
} |
} |
|
|
|
DP *dpm_nf_and_quotient(NODE b,DPM g,VECT ps,DPM *rp,P *dnp); |
|
|
|
void Pdpm_nf_and_quotient(NODE arg,LIST *rp) |
|
{ |
|
NODE b; |
|
VECT ps; |
|
DPM g,nm; |
|
P dn; |
|
VECT quo; |
|
NODE n; |
|
int ac; |
|
|
|
do_weyl = 0; dp_fcoeffs = 0; |
|
ac = argc(arg); |
|
if ( ac < 2 ) |
|
error("dpm_nf_and_quotient : invalid arguments"); |
|
else if ( ac == 2 ) { |
|
asir_assert(ARG1(arg),O_VECT,"dpm_nf_and_quotient"); |
|
b = 0; g = (DPM)ARG0(arg); ps = (VECT)ARG1(arg); |
|
} else if ( ac == 3 ) { |
|
asir_assert(ARG0(arg),O_LIST,"dpm_nf_and_quotient"); |
|
asir_assert(ARG2(arg),O_VECT,"dpm_nf_and_quotient"); |
|
b = BDY((LIST)ARG0(arg)); g = (DPM)ARG1(arg); ps = (VECT)ARG2(arg); |
|
} |
|
NEWVECT(quo); quo->len = ps->len; |
|
if ( g ) { |
|
quo->body = (pointer *)dpm_nf_and_quotient(b,g,ps,&nm,&dn); |
|
} else { |
|
quo->body = (pointer *)MALLOC(quo->len*sizeof(pointer)); |
|
nm = 0; dn = (P)ONE; |
|
} |
|
n = mknode(3,nm,dn,quo); |
|
MKLIST(*rp,n); |
|
} |
|
|
void Pdpm_weyl_nf(NODE arg,DPM *rp) |
void Pdpm_weyl_nf(NODE arg,DPM *rp) |
{ |
{ |
NODE b; |
NODE b; |
DPM *ps; |
VECT ps; |
DPM g; |
DPM g; |
int full; |
int ac,full; |
|
|
if ( !(g = (DPM)ARG1(arg)) ) { |
if ( !(g = (DPM)ARG1(arg)) ) { |
*rp = 0; return; |
*rp = 0; return; |
} |
} |
asir_assert(ARG0(arg),O_LIST,"dpm_weyl_nf"); |
do_weyl = 1; dp_fcoeffs = 0; |
asir_assert(ARG1(arg),O_DPM,"dpm_weyl_nf"); |
ac = argc(arg); |
asir_assert(ARG2(arg),O_VECT,"dpm_weyl_nf"); |
if ( ac < 3 ) |
asir_assert(ARG3(arg),O_N,"dpm_weyl_nf"); |
error("dpm_weyl_nf: invalid arguments"); |
b = BDY((LIST)ARG0(arg)); ps = (DPM *)BDY((VECT)ARG2(arg)); |
else if ( ac == 3 ) { |
full = (Q)ARG3(arg) ? 1 : 0; |
asir_assert(ARG1(arg),O_VECT,"dpm_nf"); |
do_weyl = 1; |
b = 0; g = (DPM)ARG0(arg); ps = (VECT)ARG1(arg); |
|
} else if ( ac == 4 ) { |
|
asir_assert(ARG0(arg),O_LIST,"dpm_weyl_nf"); |
|
asir_assert(ARG2(arg),O_VECT,"dpm_weyl_nf"); |
|
b = BDY((LIST)ARG0(arg)); g = (DPM)ARG1(arg); ps = (VECT)ARG2(arg); |
|
full = (Q)ARG3(arg) ? 1 : 0; |
|
} |
dpm_nf_z(b,g,ps,full,DP_Multiple,rp); |
dpm_nf_z(b,g,ps,full,DP_Multiple,rp); |
do_weyl = 0; |
do_weyl = 0; |
} |
} |
Line 1349 void Pdp_weyl_nf_f(NODE arg,DP *rp) |
|
Line 1434 void Pdp_weyl_nf_f(NODE arg,DP *rp) |
|
void Pdpm_nf_f(NODE arg,DPM *rp) |
void Pdpm_nf_f(NODE arg,DPM *rp) |
{ |
{ |
NODE b; |
NODE b; |
DPM *ps; |
VECT ps; |
DPM g; |
DPM g; |
int full; |
int ac,full; |
|
|
if ( !(g = (DPM)ARG1(arg)) ) { |
if ( !(g = (DPM)ARG1(arg)) ) { |
*rp = 0; return; |
*rp = 0; return; |
} |
} |
asir_assert(ARG0(arg),O_LIST,"dpm_nf_f"); |
ac = argc(arg); |
asir_assert(ARG1(arg),O_DPM,"dpm_nf_f"); |
if ( ac < 3 ) |
asir_assert(ARG2(arg),O_VECT,"dpm_nf_f"); |
error("dpm_nf_f: invalid arguments"); |
asir_assert(ARG3(arg),O_N,"dpm_nf_f"); |
else if ( ac == 3 ) { |
b = BDY((LIST)ARG0(arg)); ps = (DPM *)BDY((VECT)ARG2(arg)); |
asir_assert(ARG1(arg),O_VECT,"dpm_nf_f"); |
full = (Q)ARG3(arg) ? 1 : 0; |
b = 0; g = (DPM)ARG0(arg); ps = (VECT)ARG1(arg); |
|
} else if ( ac == 4 ) { |
|
asir_assert(ARG0(arg),O_LIST,"dpm_nf_f"); |
|
asir_assert(ARG2(arg),O_VECT,"dpm_nf_f"); |
|
b = BDY((LIST)ARG0(arg)); g = (DPM)ARG1(arg); ps = (VECT)ARG2(arg); |
|
full = (Q)ARG3(arg) ? 1 : 0; |
|
} |
|
do_weyl = 0; |
dpm_nf_f(b,g,ps,full,rp); |
dpm_nf_f(b,g,ps,full,rp); |
} |
} |
|
|
void Pdpm_weyl_nf_f(NODE arg,DPM *rp) |
void Pdpm_weyl_nf_f(NODE arg,DPM *rp) |
{ |
{ |
NODE b; |
NODE b; |
DPM *ps; |
VECT ps; |
DPM g; |
DPM g; |
int full; |
int ac,full; |
|
|
if ( !(g = (DPM)ARG1(arg)) ) { |
if ( !(g = (DPM)ARG1(arg)) ) { |
*rp = 0; return; |
*rp = 0; return; |
} |
} |
asir_assert(ARG0(arg),O_LIST,"dpm_weyl_nf_f"); |
ac = argc(arg); |
asir_assert(ARG1(arg),O_DP,"dpm_weyl_nf_f"); |
if ( ac < 3 ) |
asir_assert(ARG2(arg),O_VECT,"dpm_weyl_nf_f"); |
error("dpm_weyl_nf_f: invalid arguments"); |
asir_assert(ARG3(arg),O_N,"dpm_weyl_nf_f"); |
else if ( ac == 3 ) { |
b = BDY((LIST)ARG0(arg)); ps = (DPM *)BDY((VECT)ARG2(arg)); |
asir_assert(ARG1(arg),O_VECT,"dpm_weyl_nf_f"); |
full = (Q)ARG3(arg) ? 1 : 0; |
b = 0; g = (DPM)ARG0(arg); ps = (VECT)ARG1(arg); |
|
} else if ( ac == 4 ) { |
|
asir_assert(ARG0(arg),O_LIST,"dpm_weyl_nf_f"); |
|
asir_assert(ARG2(arg),O_VECT,"dpm_weyl_nf_f"); |
|
b = BDY((LIST)ARG0(arg)); g = (DPM)ARG1(arg); ps = (VECT)ARG2(arg); |
|
full = (Q)ARG3(arg) ? 1 : 0; |
|
} |
do_weyl = 1; |
do_weyl = 1; |
dpm_nf_f(b,g,ps,full,rp); |
dpm_nf_f(b,g,ps,full,rp); |
do_weyl = 0; |
do_weyl = 0; |
Line 1763 void Pdp_redble(NODE arg,Z *rp) |
|
Line 1861 void Pdp_redble(NODE arg,Z *rp) |
|
*rp = 0; |
*rp = 0; |
} |
} |
|
|
|
void Pdpm_redble(NODE arg,Z *rp) |
|
{ |
|
asir_assert(ARG0(arg),O_DPM,"dpm_redble"); |
|
asir_assert(ARG1(arg),O_DPM,"dpm_redble"); |
|
if ( dpm_redble((DPM)ARG0(arg),(DPM)ARG1(arg)) ) |
|
*rp = ONE; |
|
else |
|
*rp = 0; |
|
} |
|
|
|
void dpm_schreyer_base(LIST g,LIST *s); |
|
|
|
void Pdpm_schreyer_base(NODE arg,LIST *rp) |
|
{ |
|
asir_assert(ARG0(arg),O_LIST,"dpm_schreyer_base"); |
|
dpm_schreyer_base((LIST)ARG0(arg),rp); |
|
} |
|
|
|
void dpm_simplify_syz(LIST m,LIST s,LIST *m1,LIST *s1); |
|
|
|
void Pdpm_simplify_syz(NODE arg,LIST *rp) |
|
{ |
|
LIST s1,m1; |
|
NODE t; |
|
|
|
asir_assert(ARG0(arg),O_LIST,"dpm_simplify_syz"); |
|
asir_assert(ARG1(arg),O_LIST,"dpm_simplify_syz"); |
|
dpm_simplify_syz((LIST)ARG0(arg),(LIST)ARG1(arg),&s1,&m1); |
|
t = mknode(2,s1,m1); |
|
MKLIST(*rp,t); |
|
} |
|
|
|
|
void Pdp_red_mod(NODE arg,LIST *rp) |
void Pdp_red_mod(NODE arg,LIST *rp) |
{ |
{ |
DP h,r; |
DP h,r; |
Line 1925 void Pdp_weyl_sp(NODE arg,DP *rp) |
|
Line 2056 void Pdp_weyl_sp(NODE arg,DP *rp) |
|
do_weyl = 0; |
do_weyl = 0; |
} |
} |
|
|
void Pdpm_sp(NODE arg,DPM *rp) |
void Pdpm_sp(NODE arg,Obj *rp) |
{ |
{ |
DPM p1,p2; |
DPM p1,p2,sp; |
|
DP mul1,mul2; |
|
Obj val; |
|
NODE nd; |
|
LIST l; |
|
|
do_weyl = 0; |
do_weyl = 0; |
p1 = (DPM)ARG0(arg); p2 = (DPM)ARG1(arg); |
p1 = (DPM)ARG0(arg); p2 = (DPM)ARG1(arg); |
asir_assert(p1,O_DPM,"dpm_sp"); asir_assert(p2,O_DPM,"dpm_sp"); |
asir_assert(p1,O_DPM,"dpm_sp"); asir_assert(p2,O_DPM,"dpm_sp"); |
dpm_sp(p1,p2,rp); |
dpm_sp(p1,p2,&sp,&mul1,&mul2); |
|
if ( get_opt("coef",&val) && val ) { |
|
nd = mknode(3,sp,mul1,mul2); |
|
MKLIST(l,nd); |
|
*rp = (Obj)l; |
|
} else { |
|
*rp = (Obj)sp; |
|
} |
} |
} |
|
|
void Pdpm_weyl_sp(NODE arg,DPM *rp) |
void Pdpm_weyl_sp(NODE arg,Obj *rp) |
{ |
{ |
DPM p1,p2; |
DPM p1,p2,sp; |
|
DP mul1,mul2; |
|
Obj val; |
|
NODE nd; |
|
LIST l; |
|
|
p1 = (DPM)ARG0(arg); p2 = (DPM)ARG1(arg); |
p1 = (DPM)ARG0(arg); p2 = (DPM)ARG1(arg); |
asir_assert(p1,O_DPM,"dpm_weyl_sp"); asir_assert(p2,O_DPM,"dpm_weyl_sp"); |
asir_assert(p1,O_DPM,"dpm_weyl_sp"); asir_assert(p2,O_DPM,"dpm_weyl_sp"); |
do_weyl = 1; |
do_weyl = 1; |
dpm_sp(p1,p2,rp); |
dpm_sp(p1,p2,&sp,&mul1,&mul2); |
do_weyl = 0; |
do_weyl = 0; |
|
if ( get_opt("coef",&val) && val ) { |
|
nd = mknode(3,sp,mul1,mul2); |
|
MKLIST(l,nd); |
|
*rp = (Obj)l; |
|
} else { |
|
*rp = (Obj)sp; |
|
} |
} |
} |
|
|
void Pdp_sp_mod(NODE arg,DP *rp) |
void Pdp_sp_mod(NODE arg,DP *rp) |
Line 2813 void Pnd_gr_trace(NODE arg,LIST *rp) |
|
Line 2966 void Pnd_gr_trace(NODE arg,LIST *rp) |
|
{ |
{ |
LIST f,v; |
LIST f,v; |
int m,homo,ac; |
int m,homo,ac; |
|
Obj val; |
|
int retdp; |
Num nhomo; |
Num nhomo; |
struct order_spec *ord; |
struct order_spec *ord; |
|
|
Line 2836 void Pnd_gr_trace(NODE arg,LIST *rp) |
|
Line 2991 void Pnd_gr_trace(NODE arg,LIST *rp) |
|
homo = ZTOS((Q)nhomo); |
homo = ZTOS((Q)nhomo); |
} else |
} else |
error("nd_gr_trace : invalid argument"); |
error("nd_gr_trace : invalid argument"); |
nd_gr_trace(f,v,m,homo,0,ord,rp); |
retdp = 0; |
|
if ( get_opt("dp",&val) && val ) retdp = 1; |
|
nd_gr_trace(f,v,m,homo,retdp,0,ord,rp); |
} |
} |
|
|
void Pnd_f4_trace(NODE arg,LIST *rp) |
void Pnd_f4_trace(NODE arg,LIST *rp) |
{ |
{ |
LIST f,v; |
LIST f,v; |
int m,homo,ac; |
int m,homo,ac; |
|
int retdp; |
|
Obj val; |
Num nhomo; |
Num nhomo; |
struct order_spec *ord; |
struct order_spec *ord; |
|
|
Line 2866 void Pnd_f4_trace(NODE arg,LIST *rp) |
|
Line 3025 void Pnd_f4_trace(NODE arg,LIST *rp) |
|
homo = ZTOS((Q)nhomo); |
homo = ZTOS((Q)nhomo); |
} else |
} else |
error("nd_gr_trace : invalid argument"); |
error("nd_gr_trace : invalid argument"); |
nd_gr_trace(f,v,m,homo,1,ord,rp); |
retdp = 0; |
|
if ( get_opt("dp",&val) && val ) retdp = 1; |
|
nd_gr_trace(f,v,m,homo,retdp,1,ord,rp); |
} |
} |
|
|
void Pnd_weyl_gr(NODE arg,LIST *rp) |
void Pnd_weyl_gr(NODE arg,LIST *rp) |
Line 2907 void Pnd_weyl_gr(NODE arg,LIST *rp) |
|
Line 3068 void Pnd_weyl_gr(NODE arg,LIST *rp) |
|
void Pnd_weyl_gr_trace(NODE arg,LIST *rp) |
void Pnd_weyl_gr_trace(NODE arg,LIST *rp) |
{ |
{ |
LIST f,v; |
LIST f,v; |
int m,homo,ac; |
int m,homo,ac,retdp; |
|
Obj val; |
Num nhomo; |
Num nhomo; |
struct order_spec *ord; |
struct order_spec *ord; |
|
|
Line 2931 void Pnd_weyl_gr_trace(NODE arg,LIST *rp) |
|
Line 3093 void Pnd_weyl_gr_trace(NODE arg,LIST *rp) |
|
homo = ZTOS((Q)nhomo); |
homo = ZTOS((Q)nhomo); |
} else |
} else |
error("nd_weyl_gr_trace : invalid argument"); |
error("nd_weyl_gr_trace : invalid argument"); |
nd_gr_trace(f,v,m,homo,0,ord,rp); |
retdp = 0; |
|
if ( get_opt("dp",&val) && val ) retdp = 1; |
|
nd_gr_trace(f,v,m,homo,retdp,0,ord,rp); |
do_weyl = 0; |
do_weyl = 0; |
} |
} |
|
|
Line 3750 void Pdpv_ord(NODE arg,Obj *rp) |
|
Line 3914 void Pdpv_ord(NODE arg,Obj *rp) |
|
*rp = dp_current_modspec->obj; |
*rp = dp_current_modspec->obj; |
} |
} |
|
|
extern int dpm_ispot; |
extern int dpm_ordtype; |
|
|
void Pdpm_ord(NODE arg,LIST *rp) |
void set_schreyer_order(LIST n); |
|
|
|
LIST schreyer_obj; |
|
|
|
void Pdpm_set_schreyer(NODE arg,LIST *rp) |
{ |
{ |
Z q; |
if ( argc(arg) ) { |
NODE nd; |
schreyer_obj = (LIST)ARG0(arg); |
struct order_spec *spec; |
set_schreyer_order(schreyer_obj); |
|
|
if ( arg ) { |
|
nd = BDY((LIST)ARG0(arg)); |
|
if ( !create_order_spec(0,(Obj)ARG1(nd),&spec) ) |
|
error("dpm_ord : invalid order specification"); |
|
initdpm(spec,ZTOS((Q)ARG0(nd))); |
|
} |
} |
STOZ(dpm_ispot,q); |
*rp = schreyer_obj; |
nd = mknode(2,q,dp_current_spec->obj); |
|
MKLIST(*rp,nd); |
|
} |
} |
|
|
void Pdpm_hm(NODE arg,DPM *rp) |
void Pdpm_hm(NODE arg,DPM *rp) |
Line 3781 void Pdpm_ht(NODE arg,DPM *rp) |
|
Line 3941 void Pdpm_ht(NODE arg,DPM *rp) |
|
{ |
{ |
DPM p; |
DPM p; |
|
|
p = (DPM)ARG0(arg); asir_assert(p,O_DPM,"dp_ht"); |
p = (DPM)ARG0(arg); asir_assert(p,O_DPM,"dpm_ht"); |
dpm_ht(p,rp); |
dpm_ht(p,rp); |
} |
} |
|
|
|
void dpm_rest(DPM p,DPM *r); |
|
|
|
void Pdpm_rest(NODE arg,DPM *rp) |
|
{ |
|
DPM p; |
|
|
|
p = (DPM)ARG0(arg); asir_assert(p,O_DPM,"dpm_ht"); |
|
dpm_rest(p,rp); |
|
} |
|
|
|
|
|
void Pdpm_hp(NODE arg,Z *rp) |
|
{ |
|
DPM p; |
|
int pos; |
|
|
|
p = (DPM)ARG0(arg); asir_assert(p,O_DPM,"dpm_ht"); |
|
pos = BDY(p)->pos; |
|
STOZ(pos,*rp); |
|
} |
|
|
|
void dpm_shift(DPM p,int s,DPM *rp); |
|
|
|
void Pdpm_shift(NODE arg,DPM *rp) |
|
{ |
|
DPM p; |
|
int s; |
|
|
|
p = (DPM)ARG0(arg); asir_assert(p,O_DPM,"dpm_shift"); |
|
s = ZTOS((Z)ARG1(arg)); |
|
dpm_shift(p,s,rp); |
|
} |
|
|
|
void dpm_sort(DPM p,DPM *rp); |
|
|
|
void Pdpm_sort(NODE arg,DPM *rp) |
|
{ |
|
DPM p; |
|
int s; |
|
|
|
p = (DPM)ARG0(arg); |
|
if ( !p ) *rp = 0; |
|
dpm_sort(p,rp); |
|
} |
|
|
|
void dpm_split(DPM p,int s,DPM *up,DPM *lo); |
|
|
|
void Pdpm_split(NODE arg,LIST *rp) |
|
{ |
|
DPM p,up,lo; |
|
int s; |
|
NODE nd; |
|
|
|
p = (DPM)ARG0(arg); |
|
s = ZTOS((Z)ARG1(arg)); |
|
dpm_split(p,s,&up,&lo); |
|
nd = mknode(2,up,lo); |
|
MKLIST(*rp,nd); |
|
} |
|
|
|
|
void Pdpm_hc(NODE arg,Obj *rp) |
void Pdpm_hc(NODE arg,Obj *rp) |
{ |
{ |
asir_assert(ARG0(arg),O_DPM,"dpm_hc"); |
asir_assert(ARG0(arg),O_DPM,"dpm_hc"); |
Line 3894 int dpv_hp(DPV p) |
|
Line 4115 int dpv_hp(DPV p) |
|
case ORD_LEX: |
case ORD_LEX: |
for ( i = 0; i < len; i++ ) |
for ( i = 0; i < len; i++ ) |
if ( e[i] ) return i; |
if ( e[i] ) return i; |
|
return -1; |
|
break; |
|
default: |
|
error("dpv_hp : unsupported term ordering"); |
return -1; |
return -1; |
break; |
break; |
} |
} |