OpenXM_contrib2/asir2000/include/ca.h - diff

Return to ca.h CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / include

Diff for /OpenXM_contrib2/asir2000/include/ca.h between version 1.85 and 1.95

version 1.85, 2013/12/20 02:02:24

version 1.95, 2016/03/31 01:40:11

Line 45

* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,

* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.

* $OpenXM: OpenXM_contrib2/asir2000/include/ca.h,v 1.84 2013/07/31 00:38:12 noro Exp $

* $OpenXM: OpenXM_contrib2/asir2000/include/ca.h,v 1.94 2015/08/20 08:42:07 noro Exp $

#include <stdio.h>

#include <stdlib.h>

#include <gmp.h>

#include <mpc.h>

/* from mpfr-impl.h */

#define MPFR_PREC(x) ((x)->_mpfr_prec)

#define MPFR_EXP(x) ((x)->_mpfr_exp)

#define MPFR_MANT(x) ((x)->_mpfr_d)

#define MPFR_LAST_LIMB(x) ((MPFR_PREC (x) - 1) / GMP_NUMB_BITS)

#define MPFR_LIMB_SIZE(x) (MPFR_LAST_LIMB (x) + 1)

#if defined(hpux)

#include <netinet/in.h>

# define setbuffer(FP,buf,siz) setvbuf(FP,buf,_IOFBF,siz)

#endif

#if !defined(VISUAL)

#if !defined(VISUAL) && !defined(__MINGW32__)

#include <unistd.h>

#include <sys/param.h>

#endif

Line 69

Line 77

#define alloca(x) __builtin_alloca(x)

#endif

#if defined(VISUAL)

#if defined(VISUAL) || defined(__MINGW32__)

#include <limits.h>

#include <malloc.h>

#endif

Line 86 typedef void * pointer;

Line 94 typedef void * pointer;

#include <strings.h>

#else

#include <string.h>

#if defined(VISUAL)

#if defined(VISUAL) || defined(__MINGW32__)

#define index(s,c) strchr(s,c)

#define bzero(s,len) memset(s,0,len)

#define bcopy(x,y,len) memcpy(y,x,len)

Line 138 typedef void * pointer;

Line 146 typedef void * pointer;

#define N_R 1

#define N_A 2

#define N_B 3

#define N_C 4

#define N_NEXT_B (N_B)

#define N_M 5

#if defined(INTERVAL)

#define N_LM 6

#define N_Quad (N_NEXT_B+1)

#define N_GF2N 7

#define N_IP (N_NEXT_B+2)

#define N_GFPN 8

#define N_IntervalDouble (N_NEXT_B+3)

#define N_GFS 9

#define N_IntervalQuad (N_NEXT_B+4)

#define N_GFSN 10

#define N_IntervalBigFloat (N_NEXT_B+5)

#define N_DA 11

#define N_PRE_C N_IntervalBigFloat

#define N_GZ 12

#define N_BASE (N_NEXT_B+5)

#define N_GQ 13

#else

#define N_BASE N_NEXT_B

#endif

#define N_C (N_BASE+1)

#define N_M (N_BASE+2)

#define N_LM (N_BASE+3)

#define N_GF2N (N_BASE+4)

#define N_GFPN (N_BASE+5)

#define N_GFS (N_BASE+6)

#define N_GFSN (N_BASE+7)

#define N_DA (N_BASE+8)

#define N_GZ (N_BASE+9)

#define N_GQ (N_BASE+10)

#define N_PARIB (N_BASE+11)

#define ORD_REVGRADLEX 0

#define ORD_GRADLEX 1

Line 229 typedef struct oBF {

Line 250 typedef struct oBF {

short id;

char nid;

char pad;

long body[1];

mpfr_t body;

} *BF;

typedef struct oC {

Line 642 struct order_spec {

Line 663 struct order_spec {

int nv;

int ispot; /* 1 means Position over Term (Pos then Term) */

int pot_nelim; /* size of positions for pot-elimination order */

int *top_weight;

int module_rank;

int *module_top_weight;

union {

int simple;

struct {

Line 751 typedef unsigned int ModNum;

Line 775 typedef unsigned int ModNum;

#define GCFREE(p) Risa_GC_free(p)

#define CALLOC(d,e) MALLOC((d)*(e))

#if !defined(__CYGWIN__) && (defined(__GNUC__) || defined(vax) || defined(apollo) || defined(alloca) || defined(VISUAL))

#if !defined(__CYGWIN__) && (defined(__GNUC__) || defined(vax) || defined(apollo) || defined(alloca) || defined(VISUAL) || defined(__MINGW32__))

#define ALLOCA(d) alloca(d)

#else

#define ALLOCA(d) MALLOC(d)

#endif

/* for setjmp/longjmp compatibility */

#if defined(__CYGWIN__) || defined(__x86_64)

#if defined(__CYGWIN__) || (defined(__x86_64) && !defined(__MINGW32__))

#define JMP_BUF sigjmp_buf

#define SETJMP(x) sigsetjmp(x,~0)

#define LONGJMP(x,y) siglongjmp(x,y)

Line 863 bzero((char *)(q)->b,(w)*sizeof(unsigned int)))

Line 887 bzero((char *)(q)->b,(w)*sizeof(unsigned int)))

#define NEWReal(q) ((q)=(Real)MALLOC_ATOMIC(sizeof(struct oReal)),OID(q)=O_N,NID(q)=N_R)

#define NEWAlg(r) ((r)=(Alg)MALLOC(sizeof(struct oAlg)),OID(r)=O_N,NID(r)=N_A)

#define NEWDAlg(r) ((r)=(DAlg)MALLOC(sizeof(struct oDAlg)),OID(r)=O_N,NID(r)=N_DA)

#define NEWBF(q,l) ((q)=(BF)MALLOC_ATOMIC(TRUESIZE(oBF,(l)-1,long)),OID(q)=O_N,NID(q)=N_B)

#define NEWBF(q) ((q)=(BF)MALLOC(sizeof(struct oBF)),OID(q)=O_N,NID(q)=N_B)

#define NEWPARIBF(q,l) ((q)=(BF)MALLOC_ATOMIC(TRUESIZE(oBF,(l)-1,long)),OID(q)=O_N,NID(q)=N_PARIB)

#define NEWC(r) ((r)=(C)MALLOC(sizeof(struct oC)),OID(r)=O_N,NID(r)=N_C)

#define NEWLM(r) ((r)=(LM)MALLOC(sizeof(struct oLM)),OID(r)=O_N,NID(r)=N_LM)

#define NEWGF2N(r) ((r)=(GF2N)MALLOC(sizeof(struct oGF2N)),OID(r)=O_N,NID(r)=N_GF2N)

Line 963 PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q)

Line 988 PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q)

#define SL(n) ((n)->p)

#define ZALLOC(d) ((Z)MALLOC_ATOMIC(TRUESIZE(oZ,(d)-1,int)))

#if defined(PARI)

#define ToReal(a) (!(a)?(double)0.0:REAL(a)?BDY((Real)a):RATN(a)?RatnToReal((Q)a):BIGFLOAT(a)?mpfrtodbl(BDY((BF)a)):0)

#define ToReal(a) (!(a)?(double)0.0:REAL(a)?BDY((Real)a):RATN(a)?RatnToReal((Q)a):BIGFLOAT(a)?rtodbl(BDY((BF)a)):0)

#else

#define ToReal(a) (!(a)?(double)0.0:REAL(a)?BDY((Real)a):RATN(a)?RatnToReal((Q)a):0.0)

#endif

/* predicates */

#define NUM(p) (OID(p)==O_N)

Line 1239 void rtime_init(void);

Line 1260 void rtime_init(void);

void setmod_gf2n(P);

void mt_sgenrand(unsigned long);

unsigned long mt_genrand(void);

#if defined(VISUAL)

#if defined(VISUAL) || defined(__MINGW32__)

void srandom(unsigned int);

unsigned int random(void);

#endif

Line 1455 double NatToReal(N,int *);

Line 1476 double NatToReal(N,int *);

double RatnToReal(Q);

double pwrreal0(double,int);

double rtodbl(); /* XXX */

double mpfrtodbl(mpfr_t a);

int **almat(int,int);

pointer **almat_pointer(int,int);

int berlecnt(register int,UM);

Line 2649 void mkwcz(int k,int l,Z *t);

Line 2671 void mkwcz(int k,int l,Z *t);

int remzi(Z n,int m);

#if 0 && !defined(VISUAL)

#if 0 && !defined(VISUAL) && !defined(__MINGW32__)

void bzero(const void *,int);

void bcopy(const void *,void *,int);

char *index(char *,int);

Line 2684 int compnbp(VL vl,NBP p1,NBP p2);

Line 2706 int compnbp(VL vl,NBP p1,NBP p2);

#define MPZTOMPQ(z,q) \

(mpq_init(q),mpq_numref(q)[0] = (z)[0],mpz_set_ui(mpq_denref(q),1))

#define MPFRTOBF(g,q) \

(NEWBF(q),BDY(q)[0]=(g)[0],(q))

extern mpz_t ONEMPZ;

extern GZ ONEGZ;

Line 2700 P gptop(P a);

Line 2725 P gptop(P a);

void addgz(GZ n1,GZ n2,GZ *nr);

void subgz(GZ n1,GZ n2,GZ *nr);

void mulgz(GZ n1,GZ n2,GZ *nr);

void mul1gz(GZ n1,int n2,GZ *nr);

void divgz(GZ n1,GZ n2,GZ *nq);

void chsgngz(GZ n,GZ *nr);

void pwrgz(GZ n1,Q n,GZ *nr);

Line 2717 void mkgwc(int k,int l,GZ *t);

Line 2743 void mkgwc(int k,int l,GZ *t);

void gz_ptozp(P p,int sgn,GQ *c,P *pr);

void gz_lgp(P p,GZ *g,GZ *l);

void gz_qltozl(GQ *w,int n,GZ *dvr);

void _mdtodp(DP p,DP *pr);

void add_denomlist(P f);

void algobjtorat(Obj f,Obj *r);

void algtodalg(Alg a,DAlg *r);

void appenduflist(NODE n);

void arf_add(VL,Obj,Obj,Obj *);

void arf_chsgn(Obj,Obj *);

void arf_div(VL,Obj,Obj,Obj *);

void arf_mul(VL,Obj,Obj,Obj *);

void arf_sub(VL,Obj,Obj,Obj *);

void asir_terminate(int);

void check_intr();

void clctalg(P,VL *);

void comm_muld_trunc(VL vl,DP p1,DP p2,DL dl,DP *pr);

void comm_quod(VL vl,DP p1,DP p2,DP *pr);

void create_modorder_spec(int id,LIST shift,struct modorder_spec **s);

void dalgtoalg(DAlg da,Alg *r);

void divsgz(GZ n1,GZ n2,GZ *nq);

void dp_ht(DP p,DP *rp);

void dp_interreduce(LIST f,LIST v,int field,struct order_spec *ord,LIST *rp);

void dp_mbase(NODE,NODE *);

void dp_nf_tab_f(DP p,LIST *tab,DP *rp);

void dp_ptozp(DP,DP *);

void dp_sort(DP p,DP *rp);

void dp_subd(DP,DP,DP *);

void dp_true_nf(NODE,DP,DP *,int,DP *,P *);

void dp_true_nf_marked(NODE b,DP g,DP *ps,DP *hps,DP *rp,P *nmp,P *dnp);

void dp_true_nf_marked_mod(NODE b,DP g,DP *ps,DP *hps,int mod,DP *rp,P *dnp);

void gbcheck_list(NODE f,int n,VECT *gp,LIST *pp);

void gcdsf(VL vl,P *pa,int k,P *r);

void get_algtree(Obj f,VL *r);

void get_vars(Obj,VL *);

void gfstopgfs(GFS a,V v,P *c);

void henmain_incremental(LUM f,LUM *bqlist,ML cqlist, int np, int mod, int start, int bound);

void iftogfs(int n,GFS *c);

void indextogfs(int index,GFS *c);

void init_denomlist();

void inva_chrem(P p0,P p,NODE *pr);

void itogfs(int n,GFS *c);

void lcmn(N n1,N n2,N *nr);

void makevar(char *,P *);

void mdtodp(DP p,DP *pr);

void mfctrsf(VL vl,P f,DCP *dcp);

void mulp_trunc(VL vl,P p1,P p2,VN vn,P *pr);

void nd_det(int mod,MAT f,P *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_postproc(LIST f,LIST v,int m,struct order_spec *ord,int do_check,LIST *rp);

void nd_gr_recompute_trace(LIST f,LIST v,int m,struct order_spec *ord,LIST tlist,LIST *rp);

void nd_gr_trace(LIST f,LIST v,int trace,int homo,int f4,struct order_spec *ord,LIST *rp);

void nd_nf_p(Obj f,LIST g,LIST v,int m,struct order_spec *ord,Obj *rp);

void obj_algtodalg(Obj obj,Obj *r);

void obj_dalgtoalg(Obj obj,Obj *r);

void ox_bcast_102(int root);

void ox_reduce_102(int root,void (*func)());

void print_to_wfep(Obj obj);

void printn(N);

void printz(Z n);

void pthrootgfs(GFS a,GFS *b);

void quop_trunc(VL vl,P p1,P p2,VN vn,P *pr);

void red_by_vect_sf(int m,unsigned int *p,unsigned int *r,unsigned int hc,int len);

void removecont_array(P *c,int n,int full);

void reset_top_weight();

void setfield_dalg(NODE alist);

void setfield_gb(NODE gb,VL vl,struct order_spec *spec);

void setsecuremode(int value);

void sfbfctr_shift(P f,V x,V y,int degbound,GFS *evp,P *sfp,DCP *dcp);

void sfexgcd_by_hensel(BM g,BM h,int dy,BM *ap,BM *bp);

void sfptopsfp(P f,V v,P *gp);

void simpdalg(DAlg da,DAlg *r);

void simple_derivr(VL vl,Obj a,V v,Obj *b);

void substpp(VL vl,P p,V *vvect,P *svect,int nv,P *pr);

void ufctrsf(P p,DCP *dcp);

void vltopl(VL vl,LIST *l);

int arf_comp(VL, Obj, Obj);

int available_mcindex(int ind);

int compare_zero(int n,int *u,int row,int **w);

int create_composite_order_spec(VL vl,LIST order,struct order_spec **specp);

int create_order_spec(VL,Obj,struct order_spec **);

int dalgtoup(DAlg da,P *up,Q *dn);

int dl_equal(int nv,DL dl1,DL dl2);

int dp_redble(DP,DP);

int dpv_hp(DPV p);

int generic_gauss_elim_hensel(MAT mat,MAT *nmmat,Q *dn,int **rindp,int **cindp);

int generic_gauss_elim_hensel_dalg(MAT mat,DP *mb,MAT *nmmat,Q *dn,int **rindp,int **cindp);

int generic_gauss_elim_mod(int **mat0,int row,int col,int md,int *colstat);

int generic_gauss_elim_mod2(int **mat0,int row,int col,int md,int *colstat,int *rowstat);

int get_field_type(P p);

int get_ox_server_id(int index);

int getsecuremode();

int gz_generic_gauss_elim(MAT mat,MAT *nm,GZ *dn,int **rindp,int **cindp);

int gz_generic_gauss_elim_full(MAT mat,MAT *nm,GZ *dn,int **rindp,int **cindp);

int gz_gensolve_check(MAT mat,MAT nm,GZ dn,int *rind,int *cind);

int gz_gensolve_check2(MAT mat,MAT nm,GZ *dn,int *rind,int *cind);

int gz_intmtoratm(MAT mat,GZ md,MAT nm,GZ *dn);

int gz_intmtoratm2(MAT mat,GZ md,MAT nm,GZ *dn,int *stat);

int gz_inttorat(GZ c,GZ m,GZ b,int *sgnp,GZ *nmp,GZ *dnp);

int gz_intvtoratv(GZ *v,int n,GZ md,GZ b,GZ *nm,GZ *dn);

int invdalg(DAlg a,DAlg *c);

int is_eq(Obj a0,Obj a1);

int length(NODE);

int lu_mod(unsigned int **a,int n,unsigned int md,int **rinfo);

int n_bits_gz(GZ a);

int nd_gauss_elim_gz(GZ **mat0,int *sugar,int row,int col,int *colstat);

int poly_is_dependent(P p,V v);

int setsecureflag(char *name,int value);

int sfdegtest(int dy,int bound,UM *d1c,int k,int *in);

int sgnz(Z n);

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>