version 1.12, 2004/12/04 09:39:27 |
version 1.37, 2017/09/14 01:34:54 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.h,v 1.11 2004/10/06 12:09:19 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.h,v 1.36 2017/02/21 09:20:23 noro Exp $ */ |
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
#include "ox.h" |
#include "ox.h" |
|
|
#include <time.h> |
#include <time.h> |
|
|
#if defined(__GNUC__) |
#if defined(__GNUC__) |
#define INLINE inline |
#define INLINE static inline |
#elif defined(VISUAL) |
#elif defined(VISUAL) || defined(__MINGW32__) |
#define INLINE __inline |
#define INLINE __inline |
#else |
#else |
#define INLINE |
#define INLINE |
Line 47 typedef struct oNDV { |
|
Line 47 typedef struct oNDV { |
|
typedef union oNDC { |
typedef union oNDC { |
int m; |
int m; |
Q z; |
Q z; |
|
LM lm; |
|
GZ gz; |
P p; |
P p; |
|
R r; |
DAlg a; |
DAlg a; |
} *NDC; |
} *NDC; |
|
|
Line 77 typedef struct oND_pairs { |
|
Line 80 typedef struct oND_pairs { |
|
struct oND_pairs *next; |
struct oND_pairs *next; |
int i1,i2; |
int i1,i2; |
int sugar; |
int sugar; |
|
int sugar2; |
UINT lcm[1]; |
UINT lcm[1]; |
} *ND_pairs; |
} *ND_pairs; |
|
|
Line 115 typedef struct oIndArray |
|
Line 119 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 127 extern char *Demand; |
|
Line 136 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 HCLM(d) ((d)->body->c.lm) |
#define HCQ(d) ((d)->body->c.z) |
#define HCQ(d) ((d)->body->c.z) |
|
#define HCZ(d) ((d)->body->c.gz) |
#define HCP(d) ((d)->body->c.p) |
#define HCP(d) ((d)->body->c.p) |
#define HCA(d) ((d)->body->c.a) |
#define HCA(d) ((d)->body->c.a) |
#define CM(x) ((x)->c.m) |
#define CM(x) ((x)->c.m) |
|
#define CLM(x) ((x)->c.lm) |
#define CQ(x) ((x)->c.z) |
#define CQ(x) ((x)->c.z) |
|
#define CZ(x) ((x)->c.gz) |
#define CP(x) ((x)->c.p) |
#define CP(x) ((x)->c.p) |
#define CA(x) ((x)->c.a) |
#define CA(x) ((x)->c.a) |
#define DL(x) ((x)->dl) |
#define DL(x) ((x)->dl) |
Line 151 extern int *current_weyl_weight_vector; |
|
Line 166 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 218 NODE append_one(NODE,int); |
|
Line 234 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 237 ND_pairs nd_newpairs( NODE g, int t ); |
|
Line 254 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 ); |
ND_pairs crit_F( ND_pairs d1 ); |
ND_pairs crit_F( ND_pairs d1 ); |
int crit_2( int dp1, int dp2 ); |
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,int f4); |
|
|
/* 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 checkonly,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 274 INLINE int ndl_hash_value(UINT *d); |
|
Line 292 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 298 void ndp_print(ND_pairs d); |
|
Line 316 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); |
int nd_get_exporigin(struct order_spec *ord); |
int nd_get_exporigin(struct order_spec *ord); |
void ndv_mod(int mod,NDV p); |
void ndv_mod(int mod,NDV p); |
NDV ndv_dup(int mod,NDV p); |
NDV ndv_dup(int mod,NDV p); |
|
NDV ndv_symbolic(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); |
int nd_length(ND p); |
INLINE 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); |
UINT *nd_compute_bound(ND p); |
Line 325 ND nd_add_q(ND p1,ND p2); |
|
Line 344 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, |
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_mod_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_sf_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); |
NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred); |
NODE nd_f4_red_q_main(ND_pairs sp0,int nsp,UINT *s0vect,int col, |
NODE nd_f4_red_gz_main(ND_pairs sp0,int nsp,int trace,UINT *s0vect,int col, |
NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred); |
NM_ind_pair *rvect,int *rhead,IndArray *imat,int nred); |
|
|
/* NDV functions */ |
/* NDV functions */ |
Line 342 ND ndv_mul_nm(int mod,NM m0,NDV p); |
|
Line 367 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 357 NDV ptondv(VL vl,VL dvl,P p); |
|
Line 382 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,int *s0hash,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); |
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); |
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 nd_gauss_elim_q(Q **mat0,int *sugar,int row,int col,int *colstat); |
|
|