version 1.16, 2019/10/11 03:45:56 |
version 1.19, 2019/11/21 04:03:16 |
|
|
* 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.15 2019/09/19 06:29:47 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2018/builtin/dp.c,v 1.18 2019/11/19 10:50:31 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 90 void Pdp_vtoe(), Pdp_etov(), Pdp_dtov(), Pdp_idiv(), P |
|
Line 90 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_set_schreyer(),Pdpm_nf(),Pdpm_weyl_nf(),Pdpm_sp(),Pdpm_weyl_sp(),Pdpm_nf_and_quotient(); |
void Pdpm_ltod(),Pdpm_dtol(),Pdpm_set_schreyer(),Pdpm_nf(),Pdpm_weyl_nf(),Pdpm_sp(),Pdpm_weyl_sp(),Pdpm_nf_and_quotient(),Pdpm_nf_and_quotient2(); |
|
void Pdpm_schreyer_frame(),Pdpm_set_schreyer_level(); |
|
void Pdpm_list_to_array(),Pdpm_sp_nf(),Pdpm_insert_to_zlist(); |
void Pdpm_hm(),Pdpm_ht(),Pdpm_hc(),Pdpm_hp(),Pdpm_rest(),Pdpm_shift(),Pdpm_split(),Pdpm_sort(),Pdpm_dptodpm(),Pdpm_redble(); |
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(),Pdpm_td(); |
void Pdpm_schreyer_base(),Pdpm_simplify_syz(),Pdpm_td(); |
|
|
Line 160 struct ftab dp_tab[] = { |
|
Line 162 struct ftab dp_tab[] = { |
|
{"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_and_quotient",Pdpm_nf_and_quotient,-3}, |
{"dpm_nf_and_quotient",Pdpm_nf_and_quotient,-3}, |
|
{"dpm_nf_and_quotient2",Pdpm_nf_and_quotient2,-3}, |
{"dpm_nf_f",Pdpm_nf_f,-4}, |
{"dpm_nf_f",Pdpm_nf_f,-4}, |
{"dpm_weyl_nf_f",Pdpm_weyl_nf_f,-4}, |
{"dpm_weyl_nf_f",Pdpm_weyl_nf_f,-4}, |
{"dpm_nf",Pdpm_nf,-4}, |
{"dpm_nf",Pdpm_nf,-4}, |
Line 256 struct ftab dp_supp_tab[] = { |
|
Line 259 struct ftab dp_supp_tab[] = { |
|
{"dp_sort",Pdp_sort,1}, |
{"dp_sort",Pdp_sort,1}, |
{"dp_ord",Pdp_ord,-1}, |
{"dp_ord",Pdp_ord,-1}, |
{"dpm_set_schreyer",Pdpm_set_schreyer,-1}, |
{"dpm_set_schreyer",Pdpm_set_schreyer,-1}, |
|
{"dpm_set_schreyer_level",Pdpm_set_schreyer_level,1}, |
|
{"dpm_schreyer_frame",Pdpm_schreyer_frame,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 330 struct ftab dp_supp_tab[] = { |
|
Line 335 struct ftab dp_supp_tab[] = { |
|
{"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_schreyer_base",Pdpm_schreyer_base,1}, |
|
{"dpm_list_to_array",Pdpm_list_to_array,1}, |
|
{"dpm_sp_nf",Pdpm_sp_nf,4}, |
|
{"dpm_insert_to_zlist",Pdpm_insert_to_zlist,3}, |
{"dpm_simplify_syz",Pdpm_simplify_syz,2}, |
{"dpm_simplify_syz",Pdpm_simplify_syz,2}, |
|
|
{"dp_rref2",Pdp_rref2,2}, |
{"dp_rref2",Pdp_rref2,2}, |
Line 1344 void Pdpm_nf(NODE arg,DPM *rp) |
|
Line 1352 void Pdpm_nf(NODE arg,DPM *rp) |
|
} |
} |
|
|
DP *dpm_nf_and_quotient(NODE b,DPM g,VECT ps,DPM *rp,P *dnp); |
DP *dpm_nf_and_quotient(NODE b,DPM g,VECT ps,DPM *rp,P *dnp); |
|
DPM dpm_nf_and_quotient2(NODE b,DPM g,VECT ps,DPM *rp,P *dnp); |
|
|
void Pdpm_nf_and_quotient(NODE arg,LIST *rp) |
void Pdpm_nf_and_quotient(NODE arg,LIST *rp) |
{ |
{ |
Line 1378 void Pdpm_nf_and_quotient(NODE arg,LIST *rp) |
|
Line 1387 void Pdpm_nf_and_quotient(NODE arg,LIST *rp) |
|
MKLIST(*rp,n); |
MKLIST(*rp,n); |
} |
} |
|
|
|
void Pdpm_nf_and_quotient2(NODE arg,LIST *rp) |
|
{ |
|
NODE b; |
|
VECT ps; |
|
DPM g,nm,q; |
|
P dn; |
|
NODE n; |
|
int ac; |
|
|
|
do_weyl = 0; dp_fcoeffs = 0; |
|
ac = argc(arg); |
|
if ( ac < 2 ) |
|
error("dpm_nf_and_quotient2 : invalid arguments"); |
|
else if ( ac == 2 ) { |
|
asir_assert(ARG1(arg),O_VECT,"dpm_nf_and_quotient2"); |
|
b = 0; g = (DPM)ARG0(arg); ps = (VECT)ARG1(arg); |
|
} else if ( ac == 3 ) { |
|
asir_assert(ARG0(arg),O_LIST,"dpm_nf_and_quotient2"); |
|
asir_assert(ARG2(arg),O_VECT,"dpm_nf_and_quotient2"); |
|
b = BDY((LIST)ARG0(arg)); g = (DPM)ARG1(arg); ps = (VECT)ARG2(arg); |
|
} |
|
if ( g ) { |
|
q = dpm_nf_and_quotient2(b,g,ps,&nm,&dn); |
|
} else { |
|
q = 0; nm = 0; dn = (P)ONE; |
|
} |
|
n = mknode(3,nm,dn,q); |
|
MKLIST(*rp,n); |
|
} |
|
|
void Pdpm_weyl_nf(NODE arg,DPM *rp) |
void Pdpm_weyl_nf(NODE arg,DPM *rp) |
{ |
{ |
NODE b; |
NODE b; |
Line 1889 void Pdpm_redble(NODE arg,Z *rp) |
|
Line 1928 void Pdpm_redble(NODE arg,Z *rp) |
|
} |
} |
|
|
void dpm_schreyer_base(LIST g,LIST *s); |
void dpm_schreyer_base(LIST g,LIST *s); |
|
void dpm_schreyer_base_zlist(LIST g,LIST *s); |
|
|
void Pdpm_schreyer_base(NODE arg,LIST *rp) |
void Pdpm_schreyer_base(NODE arg,LIST *rp) |
{ |
{ |
asir_assert(ARG0(arg),O_LIST,"dpm_schreyer_base"); |
asir_assert(ARG0(arg),O_LIST,"dpm_schreyer_base"); |
dpm_schreyer_base((LIST)ARG0(arg),rp); |
dpm_schreyer_base_zlist((LIST)ARG0(arg),rp); |
} |
} |
|
|
|
void dpm_list_to_array(LIST g,VECT *psv,VECT *psiv); |
|
|
|
void Pdpm_list_to_array(NODE arg,LIST *rp) |
|
{ |
|
VECT psv,psiv; |
|
NODE nd; |
|
|
|
asir_assert(ARG0(arg),O_LIST,"dpm_list_to_array"); |
|
dpm_list_to_array((LIST)ARG0(arg),&psv,&psiv); |
|
nd = mknode(2,psv,psiv); |
|
MKLIST(*rp,nd); |
|
} |
|
|
|
/* [quo,nf] = dpm_sp_nf(psv,psiv,i,j) */ |
|
DPM dpm_sp_nf_zlist(VECT psv,VECT psiv,int i,int j,DPM *nf); |
|
|
|
void Pdpm_sp_nf(NODE arg,LIST *rp) |
|
{ |
|
VECT psv,psiv; |
|
DPM quo,nf; |
|
int i,j; |
|
NODE nd; |
|
|
|
asir_assert(ARG0(arg),O_VECT,"dpm_sp_nf"); psv = (VECT)ARG0(arg); |
|
asir_assert(ARG1(arg),O_VECT,"dpm_sp_nf"); psiv = (VECT)ARG1(arg); |
|
asir_assert(ARG2(arg),O_N,"dpm_sp_nf"); i = ZTOS((Q)ARG2(arg)); |
|
asir_assert(ARG3(arg),O_N,"dpm_sp_nf"); j = ZTOS((Q)ARG3(arg)); |
|
quo = dpm_sp_nf_zlist(psv,psiv,i,j,&nf); |
|
nd = mknode(2,quo,nf); |
|
MKLIST(*rp,nd); |
|
} |
|
|
|
void dpm_insert_to_zlist(VECT psiv,int pos,int i); |
|
|
|
/* insert_to_zlist(indarray,dpm_hp(f),i) */ |
|
void Pdpm_insert_to_zlist(NODE arg,VECT *rp) |
|
{ |
|
VECT psiv; |
|
int i,pos; |
|
|
|
asir_assert(ARG0(arg),O_VECT,"dpm_insert_to_zlist"); psiv = (VECT)ARG0(arg); |
|
asir_assert(ARG1(arg),O_N,"dpm_insert_to_zlist"); pos = ZTOS((Q)ARG1(arg)); |
|
asir_assert(ARG2(arg),O_N,"dpm_insert_to_zlist"); i = ZTOS((Q)ARG2(arg)); |
|
dpm_insert_to_zlist(psiv,pos,i); |
|
*rp = psiv; |
|
} |
|
|
|
|
void dpm_simplify_syz(LIST m,LIST s,LIST *m1,LIST *s1,LIST *w1); |
void dpm_simplify_syz(LIST m,LIST s,LIST *m1,LIST *s1,LIST *w1); |
|
|
void Pdpm_simplify_syz(NODE arg,LIST *rp) |
void Pdpm_simplify_syz(NODE arg,LIST *rp) |
Line 3966 void set_schreyer_order(LIST n); |
|
Line 4054 void set_schreyer_order(LIST n); |
|
void Pdpm_set_schreyer(NODE arg,LIST *rp) |
void Pdpm_set_schreyer(NODE arg,LIST *rp) |
{ |
{ |
if ( argc(arg) ) { |
if ( argc(arg) ) { |
set_schreyer_order((LIST)ARG0(arg)); |
set_schreyer_order(ARG0(arg)?(LIST)ARG0(arg):0); |
} |
} |
*rp = dmm_stack->obj; |
if ( dmm_stack ) |
|
*rp = dmm_stack->obj; |
|
else |
|
*rp = 0; |
} |
} |
|
|
|
DMMstack_array Schreyer_Frame; |
|
DMMstack_array dpm_schreyer_frame(NODE n); |
|
void set_schreyer_level(DMMstack_array array,int level); |
|
|
|
void Pdpm_set_schreyer_level(NODE arg,Q *rp) |
|
{ |
|
set_schreyer_level(Schreyer_Frame,ZTOS((Q)ARG0(arg))); |
|
*rp = (Q)ARG0(arg); |
|
} |
|
|
|
void Pdpm_schreyer_frame(NODE arg,LIST *rp) |
|
{ |
|
DMMstack_array a; |
|
DMMstack *body; |
|
NODE b,b1; |
|
LIST l; |
|
int len,i; |
|
|
|
Schreyer_Frame = a = dpm_schreyer_frame(BDY((LIST)ARG0(arg))); |
|
len = a->len; |
|
body = a->body; |
|
b = 0; |
|
for ( i = 0; i < len; i++ ) { |
|
MKNODE(b1,(pointer)body[i]->obj,b); |
|
b = b1; |
|
} |
|
MKLIST(l,b); |
|
*rp = l; |
|
} |
|
|
|
|
void Pdpm_hm(NODE arg,DPM *rp) |
void Pdpm_hm(NODE arg,DPM *rp) |
{ |
{ |