version 1.51, 2004/08/18 00:17:02 |
version 1.55, 2004/12/18 16:50:10 |
|
|
* 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/asir2000/include/ca.h,v 1.50 2004/07/13 07:59:53 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/include/ca.h,v 1.54 2004/12/10 07:36:34 noro Exp $ |
*/ |
*/ |
#include <stdio.h> |
#include <stdio.h> |
#include <stdlib.h> |
#include <stdlib.h> |
Line 125 typedef void * pointer; |
|
Line 125 typedef void * pointer; |
|
#define O_DPV 22 |
#define O_DPV 22 |
#define O_QUOTEARG 23 |
#define O_QUOTEARG 23 |
#define O_VOID -1 |
#define O_VOID -1 |
|
/* IMAT */ |
|
#define O_IMAT 24 |
|
/* IMAT */ |
|
|
#define N_Q 0 |
#define N_Q 0 |
#define N_R 1 |
#define N_R 1 |
Line 137 typedef void * pointer; |
|
Line 140 typedef void * pointer; |
|
#define N_GFPN 8 |
#define N_GFPN 8 |
#define N_GFS 9 |
#define N_GFS 9 |
#define N_GFSN 10 |
#define N_GFSN 10 |
|
#define N_DA 11 |
|
|
#define ORD_REVGRADLEX 0 |
#define ORD_REVGRADLEX 0 |
#define ORD_GRADLEX 1 |
#define ORD_GRADLEX 1 |
Line 164 typedef struct oN { |
|
Line 168 typedef struct oN { |
|
unsigned int b[1]; |
unsigned int b[1]; |
} *N; |
} *N; |
|
|
|
typedef struct oZ { |
|
int p; |
|
unsigned int b[1]; |
|
} *Z; |
|
|
typedef struct oUP2 { |
typedef struct oUP2 { |
int w; |
int w; |
unsigned int b[1]; |
unsigned int b[1]; |
Line 246 typedef struct oGFPN { |
|
Line 255 typedef struct oGFPN { |
|
struct oUP *body; |
struct oUP *body; |
} *GFPN; |
} *GFPN; |
|
|
|
typedef struct oDAlg { |
|
short id; |
|
char nid; |
|
char pad; |
|
struct oDP *nm; |
|
struct oQ *dn; |
|
} *DAlg; |
|
|
|
|
typedef struct oNum { |
typedef struct oNum { |
short id; |
short id; |
char nid; |
char nid; |
Line 308 typedef struct oGF2MAT { |
|
Line 326 typedef struct oGF2MAT { |
|
unsigned int **body; |
unsigned int **body; |
} *GF2MAT, *GFMMAT; |
} *GF2MAT, *GFMMAT; |
|
|
|
/* IMAT */ |
|
#define IMATCH 64 |
|
|
|
typedef struct oIENT { |
|
int cr; |
|
int row, col; |
|
pointer *body; |
|
} IENT; |
|
|
|
typedef struct oIMATC { |
|
pointer *fore; |
|
pointer *next; |
|
IENT ient[IMATCH]; |
|
} *IMATC; |
|
|
|
typedef struct oIMAT { |
|
short id; |
|
int row, col, clen; |
|
pointer *root; |
|
pointer *toor; |
|
} *IMAT; |
|
/* IMAT */ |
typedef struct oLIST { |
typedef struct oLIST { |
short id; |
short id; |
short pad; |
short pad; |
Line 583 struct modorder_spec { |
|
Line 623 struct modorder_spec { |
|
int *degree_shift; |
int *degree_shift; |
}; |
}; |
|
|
|
typedef struct oNumberField { |
|
int n; |
|
int dim; |
|
VL vl; |
|
P *defpoly; |
|
DP *mb; |
|
DP *ps; |
|
struct oDAlg *one; |
|
NODE ind; |
|
struct order_spec *spec; |
|
} *NumberField; |
|
|
/* structure for cputime */ |
/* structure for cputime */ |
|
|
struct oEGT { |
struct oEGT { |
Line 761 bzero((char *)(q)->b,(w)*sizeof(unsigned int))) |
|
Line 813 bzero((char *)(q)->b,(w)*sizeof(unsigned int))) |
|
#define NEWMAT(l) ((l)=(MAT)MALLOC(sizeof(struct oMAT)),OID(l)=O_MAT) |
#define NEWMAT(l) ((l)=(MAT)MALLOC(sizeof(struct oMAT)),OID(l)=O_MAT) |
#define NEWGF2MAT(l) ((l)=(GF2MAT)MALLOC(sizeof(struct oGF2MAT)),OID(l)=O_GF2MAT) |
#define NEWGF2MAT(l) ((l)=(GF2MAT)MALLOC(sizeof(struct oGF2MAT)),OID(l)=O_GF2MAT) |
#define NEWGFMMAT(l) ((l)=(GFMMAT)MALLOC(sizeof(struct oGF2MAT)),OID(l)=O_GFMMAT) |
#define NEWGFMMAT(l) ((l)=(GFMMAT)MALLOC(sizeof(struct oGF2MAT)),OID(l)=O_GFMMAT) |
|
/* IMAT */ |
|
#define NEWIMAT(l) ((l)=(IMAT)MALLOC(sizeof(struct oIMAT)),OID(l)=O_IMAT,l->clen=0,l->root=0,l->toor=0) |
|
#define NEWIENT(l) ((l)=(IMATC)MALLOC(sizeof(struct oIMATC)),l->fore=0,l->next=0) |
|
/* IMAT */ |
#define NEWReal(q) ((q)=(Real)MALLOC_ATOMIC(sizeof(struct oReal)),OID(q)=O_N,NID(q)=N_R) |
#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 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,l) ((q)=(BF)MALLOC_ATOMIC(TRUESIZE(oBF,(l)-1,long)),OID(q)=O_N,NID(q)=N_B) |
#define NEWC(r) ((r)=(C)MALLOC(sizeof(struct oC)),OID(r)=O_N,NID(r)=N_C) |
#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 NEWLM(r) ((r)=(LM)MALLOC(sizeof(struct oLM)),OID(r)=O_N,NID(r)=N_LM) |
Line 855 PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q) |
|
Line 912 PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q) |
|
#define MKReal(a,b) (!(a)?((b)=0):(NEWReal(b),BDY(b)=(a),(b))) |
#define MKReal(a,b) (!(a)?((b)=0):(NEWReal(b),BDY(b)=(a),(b))) |
#define MKAlg(b,r) \ |
#define MKAlg(b,r) \ |
(!(b)?((r)=0):NUM(b)?((r)=(Alg)(b)):(NEWAlg(r),BDY(r)=(Obj)(b),(r))) |
(!(b)?((r)=0):NUM(b)?((r)=(Alg)(b)):(NEWAlg(r),BDY(r)=(Obj)(b),(r))) |
|
#define MKDAlg(dp,den,r) (NEWDAlg(r),(r)->nm = (dp),(r)->dn=(den)) |
|
|
|
#define IMM_MAX 1073741823 |
|
#define IMM_MIN -1073741823 |
|
|
|
#define SL(n) ((n)->p) |
|
#define ZALLOC(d) ((Z)MALLOC_ATOMIC(TRUESIZE(oZ,(d)-1,int))) |
|
|
#if defined(PARI) |
#if defined(PARI) |
#define ToReal(a) (!(a)?(double)0.0:REAL(a)?BDY((Real)a):RATN(a)?RatnToReal((Q)a):BIGFLOAT(a)?rtodbl(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 |
#else |
Line 1439 void chsgngfsn(GFSN,GFSN *); |
|
Line 1503 void chsgngfsn(GFSN,GFSN *); |
|
void pwrgfsn(GFSN,Q, GFSN *); |
void pwrgfsn(GFSN,Q, GFSN *); |
int cmpgfsn(GFSN,GFSN); |
int cmpgfsn(GFSN,GFSN); |
|
|
|
void adddalg(DAlg,DAlg,DAlg *); |
|
void subdalg(DAlg,DAlg,DAlg *); |
|
void muldalg(DAlg,DAlg,DAlg *); |
|
void divdalg(DAlg,DAlg,DAlg *); |
|
void chsgndalg(DAlg,DAlg *); |
|
void pwrdalg(DAlg,Q, DAlg *); |
|
int cmpdalg(DAlg,DAlg); |
|
|
void addalg(Num,Num,Num *); |
void addalg(Num,Num,Num *); |
void addbf(Num,Num,Num *); |
void addbf(Num,Num,Num *); |
void addcplx(Num,Num,Num *); |
void addcplx(Num,Num,Num *); |
Line 1760 void savep(FILE *,P); |
|
Line 1832 void savep(FILE *,P); |
|
void savegf2n(FILE *,GF2N); |
void savegf2n(FILE *,GF2N); |
void savegfpn(FILE *,GFPN); |
void savegfpn(FILE *,GFPN); |
void savegfs(FILE *,GFS); |
void savegfs(FILE *,GFS); |
|
void savedalg(FILE *,DAlg); |
void savelm(FILE *,LM); |
void savelm(FILE *,LM); |
void savemi(FILE *,MQ); |
void savemi(FILE *,MQ); |
void savecplx(FILE *,C); |
void savecplx(FILE *,C); |
Line 1785 void loadgf2n(FILE *,GF2N *); |
|
Line 1858 void loadgf2n(FILE *,GF2N *); |
|
void loadgfpn(FILE *,GFPN *); |
void loadgfpn(FILE *,GFPN *); |
void loadgfs(FILE *,GFS *); |
void loadgfs(FILE *,GFS *); |
void loadgfsn(FILE *,GFSN *); |
void loadgfsn(FILE *,GFSN *); |
|
void loaddalg(FILE *,DAlg *); |
void loadlm(FILE *,LM *); |
void loadlm(FILE *,LM *); |
void loadmi(FILE *,MQ *); |
void loadmi(FILE *,MQ *); |
void loadcplx(FILE *,C *); |
void loadcplx(FILE *,C *); |
Line 2486 void afctrmain(VL vl,P p0,P p,int init,DCP *dcp); |
|
Line 2560 void afctrmain(VL vl,P p0,P p,int init,DCP *dcp); |
|
int divtmp(VL vl,int mod,P p1,P p2,P *q); |
int divtmp(VL vl,int mod,P p1,P p2,P *q); |
int divtdcmp(VL vl,int mod,P p1,P p2,P *q); |
int divtdcmp(VL vl,int mod,P p1,P p2,P *q); |
void GC_gcollect(); |
void GC_gcollect(); |
|
|
|
/* IMAT */ |
|
void Pnewimat(NODE, IMAT *); |
|
void PChsgnI(NODE, IMAT *); |
|
void Pm2Im(NODE, IMAT *); |
|
void PIm2m(NODE, MAT *); |
|
|
|
void AddMatI(VL, IMAT, IMAT, IMAT *); |
|
void MulMatI(VL, IMAT, IMAT, IMAT *); |
|
void MulMatG(VL, Obj, Obj, Obj *); |
|
void MulrMatI(VL, Obj, Obj, Obj *); |
|
void MulMatS(VL, IMAT, IMAT, IMAT *); |
|
void PutIent(IMAT, int, int, Obj); |
|
void GetIent(IMAT, int, int, Obj); |
|
void GetIbody(IMAT, int, int, Obj *); |
|
void ChsgnI(IMAT, IMAT *c); |
|
void AppendIent(IMAT, int, int, Obj); |
|
void MEnt(int, int, int, Obj, IENT *); |
|
void GetForeIent(IMATC *, IENT *, int *); |
|
void GetNextIent(IMATC *, IENT *, int *); |
|
void SubMatI(VL, IMAT, IMAT, IMAT *); |
|
/* IMAT */ |
|
|
|
Z stoz(int c); |
|
Z utoz(unsigned int c); |
|
Z qtoz(Q n); |
|
Q ztoq(Z n); |
|
Z chsgnz(Z n); |
|
Z simpz(Z n); |
|
Z dupz(Z n); |
|
Z absz(Z n); |
|
Z addz(Z n1,Z n2); |
|
Z subz(Z n1,Z n2); |
|
Z mulz(Z n1,Z n2); |
|
Z divsz(Z n1,Z n2); |
|
Z divz(Z n1,Z n2,Z *rem); |
|
Z gcdz(Z n1,Z n2); |
|
Z gcdz_cofactor(Z n1,Z n2,Z *c1,Z *c2); |
|
Z estimate_array_gcdz(Z *a,int n); |
|
Z array_gcdz(Z *a,int n); |
|
void mkwcz(int k,int l,Z *t); |
|
int remzi(Z n,int m); |
|
|
|
|
#if 0 && !defined(VISUAL) |
#if 0 && !defined(VISUAL) |
void bzero(const void *,int); |
void bzero(const void *,int); |