version 1.3, 2004/03/25 01:31:03 |
version 1.19, 2009/01/04 05:44:51 |
|
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.h,v 1.18 2006/12/14 02:49:59 noro Exp $ */ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
#include "ox.h" |
#include "ox.h" |
Line 47 typedef union oNDC { |
|
Line 48 typedef union oNDC { |
|
int m; |
int m; |
Q z; |
Q z; |
P p; |
P p; |
|
R r; |
|
DAlg a; |
} *NDC; |
} *NDC; |
|
|
/* monomial; linked list rep. */ |
/* monomial; linked list rep. */ |
Line 133 extern int *current_weyl_weight_vector; |
|
Line 136 extern int *current_weyl_weight_vector; |
|
#define HCM(d) ((d)->body->c.m) |
#define HCM(d) ((d)->body->c.m) |
#define HCQ(d) ((d)->body->c.z) |
#define HCQ(d) ((d)->body->c.z) |
#define HCP(d) ((d)->body->c.p) |
#define HCP(d) ((d)->body->c.p) |
#define CM(a) ((a)->c.m) |
#define HCA(d) ((d)->body->c.a) |
#define CQ(a) ((a)->c.z) |
#define CM(x) ((x)->c.m) |
#define CP(a) ((a)->c.p) |
#define CQ(x) ((x)->c.z) |
#define DL(a) ((a)->dl) |
#define CP(x) ((x)->c.p) |
#define SG(a) ((a)->sugar) |
#define CA(x) ((x)->c.a) |
#define LEN(a) ((a)->len) |
#define DL(x) ((x)->dl) |
#define LCM(a) ((a)->lcm) |
#define SG(x) ((x)->sugar) |
|
#define LEN(x) ((x)->len) |
|
#define LCM(x) ((x)->lcm) |
#define GET_EXP(d,a) (((d)[nd_epos[a].i]>>nd_epos[a].s)&nd_mask0) |
#define GET_EXP(d,a) (((d)[nd_epos[a].i]>>nd_epos[a].s)&nd_mask0) |
#define GET_EXP_MASK(d,a,m) ((((d)[nd_epos[a].i]&(m)[nd_epos[a].i])>>nd_epos[a].s)&nd_mask0) |
#define GET_EXP_MASK(d,a,m) ((((d)[nd_epos[a].i]&(m)[nd_epos[a].i])>>nd_epos[a].s)&nd_mask0) |
#define PUT_EXP(r,a,e) ((r)[nd_epos[a].i] |= ((e)<<nd_epos[a].s)) |
#define PUT_EXP(r,a,e) ((r)[nd_epos[a].i] |= ((e)<<nd_epos[a].s)) |
Line 147 extern int *current_weyl_weight_vector; |
|
Line 152 extern int *current_weyl_weight_vector; |
|
|
|
#define GET_EXP_OLD(d,a) (((d)[oepos[a].i]>>oepos[a].s)&omask0) |
#define GET_EXP_OLD(d,a) (((d)[oepos[a].i]>>oepos[a].s)&omask0) |
#define PUT_EXP_OLD(r,a,e) ((r)[oepos[a].i] |= ((e)<<oepos[a].s)) |
#define PUT_EXP_OLD(r,a,e) ((r)[oepos[a].i] |= ((e)<<oepos[a].s)) |
|
#define MPOS(d) (d[nd_mpos]) |
|
|
#define ROUND_FOR_ALIGN(s) ((((s)+sizeof(void *)-1)/sizeof(void *))*sizeof(void *)) |
#define ROUND_FOR_ALIGN(s) ((((s)+sizeof(void *)-1)/sizeof(void *))*sizeof(void *)) |
|
|
Line 214 NODE append_one(NODE,int); |
|
Line 220 NODE append_one(NODE,int); |
|
/* manipulation of coefficients */ |
/* manipulation of coefficients */ |
void nd_removecont(int mod,ND p); |
void nd_removecont(int mod,ND p); |
void nd_removecont2(ND p1,ND p2); |
void nd_removecont2(ND p1,ND p2); |
void removecont_array(Q *c,int n); |
void removecont_array(P *c,int n,int full); |
|
void removecont_array_q(Q *c,int n); |
|
|
/* GeoBucket functions */ |
/* GeoBucket functions */ |
ND normalize_pbucket(int mod,PGeoBucket g); |
ND normalize_pbucket(int mod,PGeoBucket g); |
Line 243 int ndv_newps(int m,NDV a,NDV aq); |
|
Line 250 int ndv_newps(int m,NDV a,NDV aq); |
|
|
|
/* top level functions */ |
/* top level functions */ |
void nd_gr(LIST f,LIST v,int m,int f4,struct order_spec *ord,LIST *rp); |
void nd_gr(LIST f,LIST v,int m,int f4,struct order_spec *ord,LIST *rp); |
void nd_gr_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp); |
void nd_gr_trace(LIST f,LIST v,int trace,int homo,int f4,struct order_spec *ord,LIST *rp); |
NODE nd_f4(int m); |
NODE nd_f4(int m); |
NODE nd_gb(int m,int ishomo,int checkonly); |
NODE nd_gb(int m,int ishomo,int checkonly); |
NODE nd_gb_trace(int m,int ishomo); |
NODE nd_gb_trace(int m,int ishomo); |
|
NODE nd_f4_trace(int m); |
|
|
/* ndl functions */ |
/* ndl functions */ |
int ndl_weight(UINT *d); |
int ndl_weight(UINT *d); |
Line 270 INLINE int ndl_hash_value(UINT *d); |
|
Line 278 INLINE int ndl_hash_value(UINT *d); |
|
/* normal forms */ |
/* normal forms */ |
INLINE int ndl_find_reducer(UINT *g); |
INLINE int ndl_find_reducer(UINT *g); |
int nd_sp(int mod,int trace,ND_pairs p,ND *nf); |
int nd_sp(int mod,int trace,ND_pairs p,ND *nf); |
int nd_nf(int mod,ND g,NDV *ps,int full,NDC dn,ND *nf); |
int nd_nf(int mod,ND d,ND g,NDV *ps,int full,NDC dn,ND *nf); |
int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *nf); |
int nd_nf_pbucket(int mod,ND g,NDV *ps,int full,ND *nf); |
|
|
/* finalizers */ |
/* finalizers */ |
Line 293 void ndp_print(ND_pairs d); |
|
Line 301 void ndp_print(ND_pairs d); |
|
|
|
/* setup, reconstruct */ |
/* setup, reconstruct */ |
void nd_init_ord(struct order_spec *spec); |
void nd_init_ord(struct order_spec *spec); |
ND_pairs nd_reconstruct(int mod,int trace,ND_pairs ndp); |
ND_pairs nd_reconstruct(int trace,ND_pairs ndp); |
void ndv_setup(int mod,int trace,NODE f); |
void ndv_setup(int mod,int trace,NODE f,int dont_sort,int dont_removecont); |
void nd_setup_parameters(int nvar,int max); |
void nd_setup_parameters(int nvar,int max); |
BlockMask nd_create_blockmask(struct order_spec *ord); |
BlockMask nd_create_blockmask(struct order_spec *ord); |
EPOS nd_create_epos(struct order_spec *ord); |
EPOS nd_create_epos(struct order_spec *ord); |
Line 306 ND nd_dup(ND p); |
|
Line 314 ND nd_dup(ND p); |
|
/* ND functions */ |
/* ND functions */ |
int ndv_check_candidate(NODE input,int obpe,int oadv,EPOS oepos,NODE cand); |
int ndv_check_candidate(NODE input,int obpe,int oadv,EPOS oepos,NODE cand); |
void nd_mul_c(int mod,ND p,int mul); |
void nd_mul_c(int mod,ND p,int mul); |
void nd_mul_c_q(ND p,Q mul); |
void nd_mul_c_q(ND p,P mul); |
void nd_mul_c_p(VL vl,ND p,P mul); |
void nd_mul_c_p(VL vl,ND p,P mul); |
ND nd_remove_head(ND p); |
ND nd_remove_head(ND p); |
ND nd_separate_head(ND p,ND *head); |
ND nd_separate_head(ND p,ND *head); |
int nd_length(ND p); |
int nd_length(ND p); |
void nd_append_red(UINT *d,int i); |
void nd_append_red(UINT *d,int i); |
UINT *ndv_compute_bound(NDV p); |
UINT *ndv_compute_bound(NDV p); |
|
UINT *nd_compute_bound(ND p); |
ND nd_copy(ND p); |
ND nd_copy(ND p); |
ND nd_merge(ND p1,ND p2); |
ND nd_merge(ND p1,ND p2); |
ND nd_add(int mod,ND p1,ND p2); |
ND nd_add(int mod,ND p1,ND p2); |
ND nd_add_q(ND p1,ND p2); |
ND nd_add_q(ND p1,ND p2); |
ND nd_add_sf(ND p1,ND p2); |
ND nd_add_sf(ND p1,ND p2); |
|
ND nd_quo(int mod,PGeoBucket p,NDV d); |
INLINE int nd_length(ND p); |
INLINE int nd_length(ND p); |
NODE nd_f4_red(int m,ND_pairs sp0,UINT *s0vect,int col,NODE rp0); |
NODE nd_f4_red(int m,ND_pairs sp0,int trace,UINT *s0vect,int col,NODE rp0,ND_pairs *nz); |
NODE nd_f4_red_dist(int m,ND_pairs sp0,UINT *s0vect,int col,NODE rp0); |
NODE nd_f4_red_dist(int m,ND_pairs sp0,UINT *s0vect,int col,NODE rp0, ND_pairs *nz); |
|
NODE nd_f4_red_main(int m,ND_pairs sp0,int nsp,UINT *s0vect,int col, |
|
NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred,ND_pairs *nz); |
|
NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,int trace,UINT *s0vect,int col, |
|
NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred); |
|
|
/* NDV functions */ |
/* NDV functions */ |
ND weyl_ndv_mul_nm(int mod,NM m0,NDV p); |
ND weyl_ndv_mul_nm(int mod,NM m0,NDV p); |
Line 329 void ndv_mul_c(int mod,NDV p,int mul); |
|
Line 343 void ndv_mul_c(int mod,NDV p,int mul); |
|
void ndv_mul_c_q(NDV p,Q mul); |
void ndv_mul_c_q(NDV p,Q mul); |
ND ndv_mul_nm_symbolic(NM m0,NDV p); |
ND ndv_mul_nm_symbolic(NM m0,NDV p); |
ND ndv_mul_nm(int mod,NM m0,NDV p); |
ND ndv_mul_nm(int mod,NM m0,NDV p); |
|
ND ndv_mul_nmv_trunc(int mod,NMV m0,NDV p,UINT *d); |
void ndv_realloc(NDV p,int obpe,int oadv,EPOS oepos); |
void ndv_realloc(NDV p,int obpe,int oadv,EPOS oepos); |
NDV ndv_dup_realloc(NDV p,int obpe,int oadv,EPOS oepos); |
NDV ndv_dup_realloc(NDV p,int obpe,int oadv,EPOS oepos); |
void ndv_homogenize(NDV p,int obpe,int oadv,EPOS eops); |
void ndv_homogenize(NDV p,int obpe,int oadv,EPOS eops); |
Line 346 NDV ptondv(VL vl,VL dvl,P p); |
|
Line 361 NDV ptondv(VL vl,VL dvl,P p); |
|
P ndvtop(int mod,VL vl,VL dvl,NDV p); |
P ndvtop(int mod,VL vl,VL dvl,NDV p); |
NDV ndtondv(int mod,ND p); |
NDV ndtondv(int mod,ND p); |
ND ndvtond(int mod,NDV p); |
ND ndvtond(int mod,NDV p); |
int nm_ind_pair_to_vect(int m,UINT *s0,int n,NM_ind_pair pair,UINT *r); |
Q *nm_ind_pair_to_vect(int m,UINT *s0,int n,NM_ind_pair pair); |
IndArray nm_ind_pair_to_vect_compress(int m,UINT *s0,int n,NM_ind_pair pair); |
IndArray nm_ind_pair_to_vect_compress(int m,UINT *s0,int n,NM_ind_pair pair); |
int nd_to_vect(int mod,UINT *s0,int n,ND d,UINT *r); |
int nd_to_vect(int mod,UINT *s0,int n,ND d,UINT *r); |
|
int nd_to_vect_q(UINT *s0,int n,ND d,Q *r); |
|
NDV vect_to_ndv_q(Q *vect,int spcol,int col,int *rhead,UINT *s0vect); |
|
|
/* elimination */ |
/* elimination */ |
int nd_gauss_elim_mod(int **mat0,int *sugar,int row,int col,int md,int *colstat); |
int nd_gauss_elim_mod(int **mat0,int *sugar,ND_pairs *spactive,int row,int col,int md,int *colstat); |
int nd_gauss_elim_sf(int **mat0,int *sugar,int row,int col,int md,int *colstat); |
int nd_gauss_elim_sf(int **mat0,int *sugar,int row,int col,int md,int *colstat); |
|
int nd_gauss_elim_q(Q **mat0,int *sugar,int row,int col,int *colstat); |
|
|
int ndl_ww_lex_compare(UINT *a1,UINT *a2); |
int ndl_ww_lex_compare(UINT *a1,UINT *a2); |
|
|
|
|