[BACK]Return to ca.h CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / include

Diff for /OpenXM_contrib2/asir2000/include/ca.h between version 1.13 and 1.26

version 1.13, 2001/04/20 02:34:22 version 1.26, 2001/09/17 07:16:59
Line 45 
Line 45 
  * 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.12 2001/03/13 01:10:26 noro Exp $   * $OpenXM: OpenXM_contrib2/asir2000/include/ca.h,v 1.25 2001/09/13 03:04:28 noro Exp $
 */  */
 #include <stdio.h>  #include <stdio.h>
   
Line 121  typedef void * pointer;
Line 121  typedef void * pointer;
 #define N_GF2N 7  #define N_GF2N 7
 #define N_GFPN 8  #define N_GFPN 8
 #define N_GFS 9  #define N_GFS 9
   #define N_GFSN 10
   
 #define ORD_REVGRADLEX 0  #define ORD_REVGRADLEX 0
 #define ORD_GRADLEX 1  #define ORD_GRADLEX 1
Line 246  typedef struct oGFS {
Line 247  typedef struct oGFS {
         int cont;          int cont;
 } *GFS;  } *GFS;
   
   typedef struct oGFSN {
           short id;
           char nid;
           char pad;
           struct oUM *body;
   } *GFSN;
   
 typedef struct oP {  typedef struct oP {
         short id;          short id;
Line 360  typedef struct oDL {
Line 367  typedef struct oDL {
         int d[1];          int d[1];
 } *DL;  } *DL;
   
   /*
    * compressed DP
    */
   
   typedef struct oCDP {
           int len;
           int psindex;
           unsigned short *body;
   } *CDP;
   
   typedef struct oCM {
           int index;
           int c;
   } *CM;
   
   /* bucket list for DL */
   
   typedef struct oDLBUCKET {
           int td;
           struct oNODE *body;
           struct oDLBUCKET *next;
   } *DLBUCKET;
   
 typedef struct oVL {  typedef struct oVL {
         V v;          V v;
         struct oVL *next;          struct oVL *next;
Line 370  typedef struct oNODE {
Line 400  typedef struct oNODE {
         struct oNODE *next;          struct oNODE *next;
 } *NODE;  } *NODE;
   
   /* univariate poly over small finite field; dense */
 typedef struct oUM {  typedef struct oUM {
         int d;          int d;
         int c[1];          int c[1];
 } *UM;  } *UM;
   
   /* univariate poly with padic coeff */
 typedef struct oLUM {  typedef struct oLUM {
         int d;          int d;
         int *c[1];          int *c[1];
 } *LUM;  } *LUM;
   
   /* bivariate poly over small finite field; dense */
   
   typedef struct oBM {
           int d;
           UM c[1];
   } *BM;
   
 typedef struct oML {  typedef struct oML {
         int n;          int n;
         int mod;          int mod;
Line 442  struct oEGT {
Line 481  struct oEGT {
 #define FF_GF2N 2  #define FF_GF2N 2
 #define FF_GFPN 3  #define FF_GFPN 3
 #define FF_GFS 4  #define FF_GFS 4
   #define FF_GFSN 5
   
 /* include interval.h */  /* include interval.h */
 #include "interval.h"  #include "interval.h"
Line 529  bzero((char *)(p),(int)(((n)+1)*sizeof(type))))
Line 569  bzero((char *)(p),(int)(((n)+1)*sizeof(type))))
         (p) = ___q___;\          (p) = ___q___;\
 }  }
   
   #define W_BMALLOC(dx,dy,p)\
   {\
           BM ___q___;\
           int ___i___;\
           UM *___c___;\
           (___q___) = (BM)ALLOCA(TRUESIZE(oBM,(dy),UM));\
           DEG(___q___) = dy;\
           ___c___ = (UM *)COEF(___q___);\
           for ( ___i___ = 0; ___i___ <= dy; ___i___++ ) {\
                   ___c___[___i___] = W_UMALLOC(dx);\
                   clearum(___c___[___i___],dx);\
           }\
           (p) = ___q___;\
   }
   
 #define NEWUP2(q,w)\  #define NEWUP2(q,w)\
 ((q)=(UP2)MALLOC_ATOMIC(TRUESIZE(oUP2,(w)-1,unsigned int)),\  ((q)=(UP2)MALLOC_ATOMIC(TRUESIZE(oUP2,(w)-1,unsigned int)),\
 bzero((char *)(q)->b,(w)*sizeof(unsigned int)))  bzero((char *)(q)->b,(w)*sizeof(unsigned int)))
Line 542  bzero((char *)(q)->b,(w)*sizeof(unsigned int)))
Line 597  bzero((char *)(q)->b,(w)*sizeof(unsigned int)))
 #define NEWQ(q) ((q)=(Q)MALLOC(sizeof(struct oQ)),OID(q)=O_N,NID(q)=N_Q)  #define NEWQ(q) ((q)=(Q)MALLOC(sizeof(struct oQ)),OID(q)=O_N,NID(q)=N_Q)
 #define NEWMQ(q) ((q)=(MQ)MALLOC_ATOMIC(sizeof(struct oMQ)),OID(q)=O_N,NID(q)=N_M)  #define NEWMQ(q) ((q)=(MQ)MALLOC_ATOMIC(sizeof(struct oMQ)),OID(q)=O_N,NID(q)=N_M)
 #define NEWGFS(q) ((q)=(GFS)MALLOC_ATOMIC(sizeof(struct oGFS)),OID(q)=O_N,NID(q)=N_GFS)  #define NEWGFS(q) ((q)=(GFS)MALLOC_ATOMIC(sizeof(struct oGFS)),OID(q)=O_N,NID(q)=N_GFS)
   #define NEWGFSN(q) ((q)=(GFSN)MALLOC(sizeof(struct oGFSN)),OID(q)=O_N,NID(q)=N_GFSN)
 #define NEWP(p) ((p)=(P)MALLOC(sizeof(struct oP)),OID(p)=O_P)  #define NEWP(p) ((p)=(P)MALLOC(sizeof(struct oP)),OID(p)=O_P)
 #define NEWR(r) ((r)=(R)MALLOC(sizeof(struct oR)),OID(r)=O_R,(r)->reduced=0)  #define NEWR(r) ((r)=(R)MALLOC(sizeof(struct oR)),OID(r)=O_R,(r)->reduced=0)
 #define NEWLIST(l) ((l)=(LIST)MALLOC(sizeof(struct oLIST)),OID(l)=O_LIST)  #define NEWLIST(l) ((l)=(LIST)MALLOC(sizeof(struct oLIST)),OID(l)=O_LIST)
Line 560  bzero((char *)(q)->b,(w)*sizeof(unsigned int)))
Line 616  bzero((char *)(q)->b,(w)*sizeof(unsigned int)))
 #define NEWV(v) ((v)=(V)MALLOC(sizeof(struct oV)))  #define NEWV(v) ((v)=(V)MALLOC(sizeof(struct oV)))
 #define NEWVL(vl) ((vl)=(VL)MALLOC(sizeof(struct oVL)))  #define NEWVL(vl) ((vl)=(VL)MALLOC(sizeof(struct oVL)))
 #define NEWMP(m) ((m)=(MP)MALLOC(sizeof(struct oMP)))  #define NEWMP(m) ((m)=(MP)MALLOC(sizeof(struct oMP)))
   #define NEWDLBUCKET(a) ((a)=(DLBUCKET)MALLOC(sizeof(struct oDLBUCKET)))
   
 #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)
Line 615  if(!(r)){NEWNODE(r);(c)=(r);}else{NEWNODE(NEXT(c));(c)
Line 672  if(!(r)){NEWNODE(r);(c)=(r);}else{NEWNODE(NEXT(c));(c)
 if(!(r)){NEWMP(r);(c)=(r);}else{NEWMP(NEXT(c));(c)=NEXT(c);}  if(!(r)){NEWMP(r);(c)=(r);}else{NEWMP(NEXT(c));(c)=NEXT(c);}
 #define NEXTMP2(r,c,s) \  #define NEXTMP2(r,c,s) \
 if(!(r)){(c)=(r)=(s);}else{NEXT(c)=(s);(c)=(s);}  if(!(r)){(c)=(r)=(s);}else{NEXT(c)=(s);(c)=(s);}
   #define NEXTDLBUCKET(r,c) \
   if(!(r)){NEWDLBUCKET(r);(c)=(r);}else{NEWDLBUCKET(NEXT(c));(c)=NEXT(c);}
   
 /* convertors */  /* convertors */
 #define NTOQ(n,s,q) \  #define NTOQ(n,s,q) \
Line 631  PL(NM(q))=1,BD(NM(q))[0]=ABS(n),DN(q)=0,(q)))
Line 690  PL(NM(q))=1,BD(NM(q))[0]=ABS(n),DN(q)=0,(q)))
 ((a)?(NEWMQ(b),CONT(b)=(unsigned int)(a),(b)):((b)=0))  ((a)?(NEWMQ(b),CONT(b)=(unsigned int)(a),(b)):((b)=0))
 #define MKGFS(a,b) \  #define MKGFS(a,b) \
 ((NEWGFS(b),CONT(b)=(a),(b)))  ((NEWGFS(b),CONT(b)=(a),(b)))
   #define MKGFSN(a,b) \
   ((DEG(a)>=0)?(NEWGFSN(b),BDY(b)=(a),(b)):((b)=0))
 #define STOMQ(a,b) \  #define STOMQ(a,b) \
 ((a)?(NEWMQ(b),CONT(b)=(a),(b)):((b)=0))  ((a)?(NEWMQ(b),CONT(b)=(a),(b)):((b)=0))
 #define UTON(u,n) \  #define UTON(u,n) \
Line 658  PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q)
Line 719  PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q)
 /* predicates */  /* predicates */
 #define NUM(p) (OID(p)==O_N)  #define NUM(p) (OID(p)==O_N)
 #define RAT(p) (OID(p)==O_R)  #define RAT(p) (OID(p)==O_R)
 #define INT(q) (!DN((Q)q))  
 #define RATN(a) (NID(a)==N_Q)  #define RATN(a) (NID(a)==N_Q)
   #define INT(q) (!(q)||(NUM(q)&&RATN((Num)q)&&!DN((Q)q)))
 #define REAL(a) (NID(a)==N_R)  #define REAL(a) (NID(a)==N_R)
 #define BIGFLOAT(a) (NID(a)==N_B)  #define BIGFLOAT(a) (NID(a)==N_B)
 #define SFF(a) (NID(a)==N_M)  #define SFF(a) (NID(a)==N_M)
Line 677  PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q)
Line 738  PL(NM(q))=1,BD(NM(q))[0]=(unsigned int)(n),DN(q)=0,(q)
 #define ITOS(p) (((unsigned int)(p))&0x7fffffff)  #define ITOS(p) (((unsigned int)(p))&0x7fffffff)
 #define STOI(i) ((P)((unsigned int)(i)|0x80000000))  #define STOI(i) ((P)((unsigned int)(i)|0x80000000))
   
   /* immediate GFS representation */
   
   #define IFTOF(p) ((int)(((unsigned int)(p))&0x7fffffff))
   #define FTOIF(i) ((int)(((unsigned int)(i)|0x80000000)))
   
 struct cdl {  struct cdl {
         P c;          P c;
         DL d;          DL d;
Line 734  extern UP2 ONEUP2;
Line 800  extern UP2 ONEUP2;
 extern FILE *asir_out;  extern FILE *asir_out;
 #if defined(__GNUC__)  #if defined(__GNUC__)
 extern const int sprime[];  extern const int sprime[];
 extern const int lprime[];  
 #else  #else
 extern int sprime[];  extern int sprime[];
 extern int lprime[];  
 #endif  #endif
   
   extern int lprime_size;
   extern int *lprime;
   
 extern void (*addnumt[])();  extern void (*addnumt[])();
 extern void (*subnumt[])();  extern void (*subnumt[])();
 extern void (*mulnumt[])();  extern void (*mulnumt[])();
Line 1089  int int_bits(int);
Line 1156  int int_bits(int);
   
   
 LUM LUMALLOC(int, int);  LUM LUMALLOC(int, int);
   BM BMALLOC(int, int);
 Obj ToAlg(Num);  Obj ToAlg(Num);
 UM *berlemain(register int, UM, UM *);  UM *berlemain(register int, UM, UM *);
 void *Risa_GC_malloc(size_t);  void *Risa_GC_malloc(size_t);
Line 1185  void chsgngfs(GFS,GFS *);
Line 1253  void chsgngfs(GFS,GFS *);
 void pwrgfs(GFS,Q, GFS *);  void pwrgfs(GFS,Q, GFS *);
 int cmpgfs(GFS,GFS);  int cmpgfs(GFS,GFS);
   
   void addgfsn(GFSN,GFSN,GFSN *);
   void subgfsn(GFSN,GFSN,GFSN *);
   void mulgfsn(GFSN,GFSN,GFSN *);
   void divgfsn(GFSN,GFSN,GFSN *);
   void chsgngfsn(GFSN,GFSN *);
   void pwrgfsn(GFSN,Q, GFSN *);
   int cmpgfsn(GFSN,GFSN);
   
 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 1209  void addlm(LM,LM,LM *);
Line 1285  void addlm(LM,LM,LM *);
 void addstr(VL,STRING,STRING,STRING *);  void addstr(VL,STRING,STRING,STRING *);
 void addum(int,UM,UM,UM);  void addum(int,UM,UM,UM);
 void addvect(VL,VECT,VECT,VECT *);  void addvect(VL,VECT,VECT,VECT *);
   void addquote(VL,QUOTE,QUOTE,QUOTE *);
 void adjc(VL,P,P,P,Q,P *,P *);  void adjc(VL,P,P,P,Q,P *,P *);
 void afctr(VL,P,P,DCP *);  void afctr(VL,P,P,DCP *);
 void afctrmain(VL,P,P,int,DCP *);  void afctrmain(VL,P,P,int,DCP *);
Line 1232  void chsgnreal(Num,Num *);
Line 1309  void chsgnreal(Num,Num *);
 void chsgnmi(MQ,MQ *);  void chsgnmi(MQ,MQ *);
 void chsgnlm(LM,LM *);  void chsgnlm(LM,LM *);
 void chsgnvect(VECT,VECT *);  void chsgnvect(VECT,VECT *);
   void chsgnquote(QUOTE,QUOTE *);
 void clctv(VL,P,VL *);  void clctv(VL,P,VL *);
 void clctvr(VL,Obj,VL *);  void clctvr(VL,Obj,VL *);
 void cm2p(Q,Q,P,P *);  void cm2p(Q,Q,P,P *);
Line 1273  void divsrdcp(VL,P,P,P *,P *);
Line 1351  void divsrdcp(VL,P,P,P *,P *);
 void divsrmp(VL,int,P,P,P *,P *);  void divsrmp(VL,int,P,P,P *,P *);
 void divsrp(VL,P,P,P *,P *);  void divsrp(VL,P,P,P *,P *);
 void divvect(VL,Obj,Obj,Obj *);  void divvect(VL,Obj,Obj,Obj *);
   void divquote(VL,QUOTE,QUOTE,QUOTE *);
 void dtest(P,ML,int,DCP *);  void dtest(P,ML,int,DCP *);
 void dtestroot(int,int,P,LUM,struct oDUM *,DCP *);  void dtestroot(int,int,P,LUM,struct oDUM *,DCP *);
 void dtestroot1(int,int,P,LUM,P *);  void dtestroot1(int,int,P,LUM,P *);
Line 1368  void mulsgn(VN,VN,int,VN);
Line 1447  void mulsgn(VN,VN,int,VN);
 void mulsum(register int,UM,register int,UM);  void mulsum(register int,UM,register int,UM);
 void mulum(register int,UM,UM,UM);  void mulum(register int,UM,UM,UM);
 void mulvect(VL,Obj,Obj,Obj *);  void mulvect(VL,Obj,Obj,Obj *);
   void mulquote(VL,QUOTE,QUOTE,QUOTE *);
 void mulvectmat(VL,VECT,MAT,VECT *);  void mulvectmat(VL,VECT,MAT,VECT *);
 void next(VN);  void next(VN);
 void nezgcdnp_sqfr_primitive(VL,P,P *,int,P *);  void nezgcdnp_sqfr_primitive(VL,P,P *,int,P *);
Line 1413  void pwrnum(VL,Num,Num,Num *);
Line 1493  void pwrnum(VL,Num,Num,Num *);
 void pwrp(VL,P,Q,P *);  void pwrp(VL,P,Q,P *);
 void pwrq(Q,Q,Q *);  void pwrq(Q,Q,Q *);
 void pwrr(VL,Obj,Obj,Obj *);  void pwrr(VL,Obj,Obj,Obj *);
   void pwrquote(VL,QUOTE,QUOTE,QUOTE *);
 void pwrreal(Num,Num,Real *);  void pwrreal(Num,Num,Real *);
 void pwrmi(MQ,Q,MQ *);  void pwrmi(MQ,Q,MQ *);
 void pwrlm(LM,Q,LM *);  void pwrlm(LM,Q,LM *);
Line 1471  void substp(VL,P,V,P,P *);
Line 1552  void substp(VL,P,V,P,P *);
 void substvp(VL,P,VN,P *);  void substvp(VL,P,VN,P *);
 void subum(int,UM,UM,UM);  void subum(int,UM,UM,UM);
 void subvect(VL,VECT,VECT,VECT *);  void subvect(VL,VECT,VECT,VECT *);
   void subquote(VL,QUOTE,QUOTE,QUOTE *);
 void toreim(Num,Num *,Num *);  void toreim(Num,Num *,Num *);
 void ucsump(P,Q *);  void ucsump(P,Q *);
 void udivpwm(Q,P,P,P *,P *);  void udivpwm(Q,P,P,P *,P *);
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.26
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>