version 1.7, 2004/09/17 01:27:48 |
version 1.28, 2012/08/27 05:38:00 |
|
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.h,v 1.27 2010/04/16 07:13:42 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 113 typedef struct oIndArray |
|
Line 116 typedef struct oIndArray |
|
} index; |
} index; |
} *IndArray; |
} *IndArray; |
|
|
|
typedef struct oNDVI { |
|
NDV p; |
|
int i; |
|
} *NDVI; |
|
|
extern int (*ndl_compare_function)(UINT *a1,UINT *a2); |
extern int (*ndl_compare_function)(UINT *a1,UINT *a2); |
extern int nd_dcomp; |
extern int nd_dcomp; |
|
|
Line 125 extern char *Demand; |
|
Line 133 extern char *Demand; |
|
extern VL CO; |
extern VL CO; |
extern int Top,Reverse,DP_Print,dp_nelim,do_weyl,NoSugar; |
extern int Top,Reverse,DP_Print,dp_nelim,do_weyl,NoSugar; |
extern int *current_weyl_weight_vector; |
extern int *current_weyl_weight_vector; |
|
extern int *current_module_weight_vector; |
|
|
/* fundamental macros */ |
/* fundamental macros */ |
#define TD(d) (d[0]) |
#define TD(d) (d[0]) |
#define HDL(d) ((d)->body->dl) |
#define HDL(d) ((d)->body->dl) |
#define HTD(d) (TD(HDL(d))) |
#define HTD(d) (TD(HDL(d))) |
|
#define HCU(d) ((d)->body->c) |
#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 159 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 227 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 233 ND_pairs nd_newpairs( NODE g, int t ); |
|
Line 247 ND_pairs nd_newpairs( NODE g, int t ); |
|
ND_pairs nd_minp( ND_pairs d, ND_pairs *prest ); |
ND_pairs nd_minp( ND_pairs d, ND_pairs *prest ); |
ND_pairs nd_minsugarp( ND_pairs d, ND_pairs *prest ); |
ND_pairs nd_minsugarp( ND_pairs d, ND_pairs *prest ); |
NODE update_base(NODE nd,int ndp); |
NODE update_base(NODE nd,int ndp); |
ND_pairs update_pairs( ND_pairs d, NODE /* of index */ g, int t); |
ND_pairs update_pairs( ND_pairs d, NODE /* of index */ g, int t, int gensyz); |
ND_pairs equivalent_pairs( ND_pairs d1, ND_pairs *prest ); |
ND_pairs equivalent_pairs( ND_pairs d1, ND_pairs *prest ); |
ND_pairs crit_B( ND_pairs d, int s ); |
ND_pairs crit_B( ND_pairs d, int s ); |
ND_pairs crit_M( ND_pairs d1 ); |
ND_pairs crit_M( ND_pairs d1 ); |
Line 242 int crit_2( int dp1, int dp2 ); |
|
Line 256 int crit_2( int dp1, int dp2 ); |
|
int ndv_newps(int m,NDV a,NDV aq); |
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 homo,int retdp,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,int **indp); |
NODE nd_gb(int m,int ishomo,int checkonly); |
NODE nd_gb(int m,int ishomo,int checkonly,int gensyz,int **indp); |
NODE nd_gb_trace(int m,int ishomo); |
NODE nd_gb_trace(int m,int ishomo,int **indp); |
|
NODE nd_f4_trace(int m,int **indp); |
|
|
/* ndl functions */ |
/* ndl functions */ |
int ndl_weight(UINT *d); |
int ndl_weight(UINT *d); |
void ndl_weight_mask(UINT *d); |
void ndl_weight_mask(UINT *d); |
void ndl_homogenize(UINT *d,UINT *r,int obpe,EPOS oepos,int weight); |
void ndl_homogenize(UINT *d,UINT *r,int obpe,EPOS oepos,int ompos,int weight); |
void ndl_dehomogenize(UINT *p); |
void ndl_dehomogenize(UINT *p); |
void ndl_reconstruct(UINT *d,UINT *r,int obpe,EPOS oepos); |
void ndl_reconstruct(UINT *d,UINT *r,int obpe,EPOS oepos); |
INLINE int ndl_reducible(UINT *d1,UINT *d2); |
INLINE int ndl_reducible(UINT *d1,UINT *d2); |
Line 270 INLINE int ndl_hash_value(UINT *d); |
|
Line 285 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 */ |
NODE ndv_reducebase(NODE x); |
NODE ndv_reducebase(NODE x,int *perm); |
NODE ndv_reduceall(int m,NODE f); |
NODE ndv_reduceall(int m,NODE f); |
|
|
/* allocators */ |
/* allocators */ |
Line 294 void ndp_print(ND_pairs d); |
|
Line 309 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 trace,ND_pairs ndp); |
ND_pairs nd_reconstruct(int trace,ND_pairs ndp); |
void ndv_setup(int mod,int trace,NODE f); |
int 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 304 NDV ndv_dup(int mod,NDV p); |
|
Line 319 NDV ndv_dup(int mod,NDV p); |
|
ND nd_dup(ND p); |
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_membership(int m,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); |
Line 321 ND nd_add_q(ND p1,ND p2); |
|
Line 336 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); |
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 334 ND ndv_mul_nm(int mod,NM m0,NDV p); |
|
Line 353 ND ndv_mul_nm(int mod,NM m0,NDV p); |
|
ND ndv_mul_nmv_trunc(int mod,NMV m0,NDV p,UINT *d); |
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,int ompos); |
void ndv_dehomogenize(NDV p,struct order_spec *spec); |
void ndv_dehomogenize(NDV p,struct order_spec *spec); |
void ndv_removecont(int mod,NDV p); |
void ndv_removecont(int mod,NDV p); |
void ndv_print(NDV p); |
void ndv_print(NDV p); |
Line 349 NDV ptondv(VL vl,VL dvl,P p); |
|
Line 368 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); |
|
|
|
|