Initial revision

git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@2244 b2dd03c8-39d4-4d8f-98ff-823fe69b080e

1 files changed, 503 insertions, 0 deletions

diff --git a/ext/bigfloat/bigfloat.h b/ext/bigfloat/bigfloat.h
new file mode 100644
index 0000000000..7a86112c73
--- /dev/null
+++ b/ext/bigfloat/bigfloat.h
@@ -0,0 +1,503 @@
+/*
+ *
+ * Ruby BIGFLOAT(Variable Precision) extension library. 
+ *
+ *  Version 1.1.7(2001/08/27)
+ *  Version 1.1.6(2001/03/28)
+ *
+ * Copyright(C) 1999  by Shigeo Kobayashi(shigeo@tinyforest.gr.jp) 
+ *
+ * You may distribute under the terms of either the GNU General Public 
+ * License or the Artistic License, as specified in the README file 
+ * of this BigFloat distribution. 
+ *
+ */
+
+#ifndef  ____BIG_FLOAT__H____
+#define  ____BIG_FLOAT__H____
+
+#if defined(__cplusplus)
+extern "C" {
+#endif
+
+/*
+ *   #define VP_EXPORT other than static to let VP_ routines 
+ *   be called from outside of this module.
+ */
+#define VP_EXPORT static 
+
+#define U_LONG unsigned long
+#define S_LONG long
+#define U_INT  unsigned int
+#define S_INT  int
+
+/* Exception codes */
+#define VP_EXCEPTION_ALL        ((unsigned short)0xFFFF)
+#define VP_EXCEPTION_INFINITY   ((unsigned short)0x0001)
+#define VP_EXCEPTION_NaN        ((unsigned short)0x0002)
+#define VP_EXCEPTION_UNDERFLOW  ((unsigned short)0x0004)
+#define VP_EXCEPTION_OVERFLOW   ((unsigned short)0x0001) /* 0x0008) */
+#define VP_EXCEPTION_ZERODIVIDE ((unsigned short)0x0001) /* 0x0010) */
+/* Following 2 exceptions cann't controlled by user */
+#define VP_EXCEPTION_OP         ((unsigned short)0x0020)
+#define VP_EXCEPTION_MEMORY     ((unsigned short)0x0040)
+
+#define VP_SIGN_NaN                0 /* NaN                      */
+#define VP_SIGN_POSITIVE_ZERO      1 /* Positive zero            */
+#define VP_SIGN_NEGATIVE_ZERO     -1 /* Negative zero            */
+#define VP_SIGN_POSITIVE_FINITE    2 /* Positive finite number   */
+#define VP_SIGN_NEGATIVE_FINITE   -2 /* Negative finite number   */
+#define VP_SIGN_POSITIVE_INFINITE  3 /* Positive infinite number */
+#define VP_SIGN_NEGATIVE_INFINITE -3 /* Negative infinite number */
+
+/*
+ * VP representation
+ *  r = 0.xxxxxxxxx *BASE**exponent
+ */
+typedef struct {
+	VALUE  obj;     /* Back pointer(VALUE) for Ruby object.     */
+	U_LONG MaxPrec;	/* Maximum precision size                   */
+					/* This is the actual size of pfrac[]       */
+					/*(frac[0] to frac[MaxPrec] are available). */
+	U_LONG Prec;	/* Current precision size.                  */
+					/* This indicates how much the.             */
+					/* the array frac[] is actually used.       */
+	S_INT  exponent;/* Exponent part.                           */
+	short  sign;	/* Attributes of the value.                 */
+					/*
+					 *		==0 : NaN
+					 *		  1 : Positive zero
+					 *		 -1 : Negative zero
+					 *		  2 : Positive number
+					 *		 -2 : Negative number
+					 *		  3 : Positive infinite number
+					 *		 -3 : Negative infinite number
+					 */
+	short  flag;   /* Not used in vp_routines,space for user.  */
+	U_LONG frac[1];	/* Pointer to array of fraction part.       */
+} Real;
+
+/*  
+ *  ------------------
+ *   EXPORTables.
+ *  ------------------
+ */
+
+VP_EXPORT  Real *
+#ifdef HAVE_STDARG_PROTOTYPES
+VpNewRbClass(U_LONG mx,char *str,VALUE klass);
+#else
+VpNewRbClass();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT  Real *
+#ifdef HAVE_STDARG_PROTOTYPES
+VpCreateRbObject(U_LONG mx,char *str);
+#else
+VpCreateRbObject();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT S_LONG VpBaseFig();
+VP_EXPORT S_LONG VpDblFig();
+VP_EXPORT U_LONG VpBaseVal();
+
+/* Zero,Inf,NaN (isinf(),isnan() used to check) */
+VP_EXPORT double VpGetDoubleNaN();
+VP_EXPORT double VpGetDoublePosInf();
+VP_EXPORT double VpGetDoubleNegInf();
+VP_EXPORT double VpGetDoubleNegZero();
+
+/* These 2 functions added at v1.1.7 */
+VP_EXPORT U_LONG VpGetPrecLimit();
+#ifdef HAVE_STDARG_PROTOTYPES
+VP_EXPORT U_LONG VpSetPrecLimit(U_LONG n);
+#else
+VP_EXPORT U_LONG VpSetPrecLimit();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpException(unsigned short f,char *str,int always);
+#else
+VpException();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpIsNegDoubleZero(double v);
+#else
+VpIsNegDoubleZero();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT U_LONG
+#ifdef HAVE_STDARG_PROTOTYPES
+VpNumOfChars(Real *vp);
+#else
+VpNumOfChars();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT U_LONG
+#ifdef HAVE_STDARG_PROTOTYPES
+VpInit(U_LONG BaseVal);
+#else
+VpInit();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpFree(Real *pv);
+#else
+VpFree();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT Real *
+#ifdef HAVE_STDARG_PROTOTYPES
+VpAlloc(U_LONG mx, char *szVal);
+#else
+VpAlloc();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpAsgn(Real *c,Real *a,int isw);
+#else
+VpAsgn();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpAddSub(Real *c,Real *a,Real *b,int operation);
+#else
+VpAddSub();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpMult(Real *c,Real *a,Real *b);
+#else
+VpMult();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpDivd(Real *c,Real *r,Real *a,Real *b);
+#else
+VpDivd(c, r, a, b);
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpComp(Real *a,Real *b);
+#else
+VpComp();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT S_LONG
+#ifdef HAVE_STDARG_PROTOTYPES
+VpExponent10(Real *a);
+#else
+VpExponent10();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpSzMantissa(Real *a,char *psz);
+#else
+VpSzMantissa();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpToString(Real *a,char *psz,int fFmt);
+#else
+VpToString();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpCtoV(Real *a,char *int_chr,U_LONG ni,char *frac,U_LONG nf,char *exp_chr,U_LONG ne);
+#else
+VpCtoV();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpVtoD(double *d,U_LONG *e,Real *m);
+#else
+VpVtoD();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpDtoV(Real *m,double d);
+#else
+VpDtoV();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpItoV(Real *m,S_INT ival);
+#else
+VpItoV();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpSqrt(Real *y,Real *x);
+#else
+VpSqrt();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpRound(Real *y,Real *x,int sw,int f,int il);
+#else
+VpRound();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpFrac(Real *y,Real *x);
+#else
+VpFrac();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpPower(Real *y,Real *x,S_INT n);
+#else
+VpPower();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpPai(Real *y);
+#else
+VpPai();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpExp1(Real *y);
+#else
+VpExp1();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpExp(Real *y,Real *x);
+#else
+VpExp();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpSinCos(Real *psin,Real *pcos,Real *x);
+#else
+VpSinCos();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+VP_EXPORT int
+#ifdef HAVE_STDARG_PROTOTYPES
+VPrint(FILE *fp,char *cntl_chr,Real *a);
+#else
+VPrint();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+/*  
+ *  ------------------
+ *   static routines.
+ *  ------------------
+ */
+
+static int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpIsDefOP(Real *c,Real *a,Real *b,int sw);
+#else
+VpIsDefOP();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static int
+#ifdef HAVE_STDARG_PROTOTYPES
+AddExponent(Real *a,S_INT n);
+#else
+AddExponent();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static void *
+#ifdef HAVE_STDARG_PROTOTYPES
+VpMemAlloc(U_LONG mb);
+#else
+VpMemAlloc();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static unsigned short VpGetException();
+
+static void 
+#ifdef HAVE_STDARG_PROTOTYPES
+VpSetException(unsigned short f);
+#else
+VpSetException();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static int
+#ifdef HAVE_STDARG_PROTOTYPES
+MemCmp(	unsigned char *a,unsigned char *b,int n);
+#else
+MemCmp(a,b,n);
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpAddAbs(Real *a,Real *b,Real *c);
+#else
+VpAddAbs();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpSubAbs(Real *a,Real *b,Real *c);
+#else
+VpSubAbs();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static U_LONG
+#ifdef HAVE_STDARG_PROTOTYPES
+VpSetPTR(Real *a,Real *b,Real *c,U_LONG *a_pos,U_LONG *b_pos,U_LONG *c_pos,U_LONG *av,U_LONG *bv);
+#else
+VpSetPTR();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpNmlz(Real *a);
+#else
+VpNmlz();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static void
+#ifdef HAVE_STDARG_PROTOTYPES
+VpFormatSt(char *psz,S_INT fFmt);
+#else
+VpFormatSt();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static int
+#ifdef HAVE_STDARG_PROTOTYPES
+VpRdup(Real *m);
+#else
+VpRdup();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static U_LONG
+#ifdef HAVE_STDARG_PROTOTYPES
+SkipWhiteChar(char *szVal);
+#else
+SkipWhiteChar();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static U_LONG 
+#ifdef HAVE_STDARG_PROTOTYPES
+GetAddSubPrec(Real *a,Real *b);
+#else
+GetAddSubPrec();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static S_INT
+#ifdef HAVE_STDARG_PROTOTYPES
+GetPositiveInt(VALUE v);
+#else
+GetPositiveInt();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static Real *
+#ifdef HAVE_STDARG_PROTOTYPES
+GetVpValue(VALUE v,int must);
+#else
+GetVpValue();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static VALUE
+#ifdef HAVE_STDARG_PROTOTYPES
+ToValue(Real *p);
+#else
+ToValue();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static S_INT 
+#ifdef HAVE_STDARG_PROTOTYPES
+BigFloatCmp(VALUE self,VALUE r);
+#else
+BigFloatCmp();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static VALUE
+#ifdef HAVE_STDARG_PROTOTYPES
+BigFloat_divide(Real **c,Real **res,Real **div,VALUE self,VALUE r);
+#else
+BigFloat_divide();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+static VALUE
+#ifdef HAVE_STDARG_PROTOTYPES
+BigFloat_DoDivmod(VALUE self,VALUE r,Real **div,Real **mod);
+#else
+BigFloat_DoDivmod();
+#endif /* HAVE_STDARG_PROTOTYPES */
+
+/*  
+ *  ------------------
+ *  MACRO definitions.
+ *  ------------------
+ */
+#define Abs(a)     (((a)>= 0)?(a):(-(a)))
+#define Max(a, b)  (((a)>(b))?(a):(b))
+#define Min(a, b)  (((a)>(b))?(b):(a))
+
+#define IsWhiteChar(ch) (((ch==' ')||(ch=='\n')||(ch=='\t')||(ch=='\b'))?1:0)
+
+#define VpMaxPrec(a)   ((a)->MaxPrec)
+#define VpPrec(a)      ((a)->Prec)
+#define VpGetFlag(a)   ((a)->flag)
+
+/* Sign */
+
+/* VpGetSign(a) returns 1,-1 if a>0,a<0 respectively */
+#define VpGetSign(a) (((a)->sign>0)?1:(-1))
+/* Change sign of a to a>0,a<0 if s = 1,-1 respectively */
+#define VpChangeSign(a,s) {if((s)>0) (a)->sign=(short)Abs((S_LONG)(a)->sign);else (a)->sign=-(short)Abs((S_LONG)(a)->sign);}
+/* Sets sign of a to a>0,a<0 if s = 1,-1 respectively */
+#define VpSetSign(a,s)    {if((s)>0) (a)->sign=(short)VP_SIGN_POSITIVE_FINITE;else (a)->sign=(short)VP_SIGN_NEGATIVE_FINITE;}
+
+/* 1 */
+#define VpSetOne(a)       {(a)->frac[0]=(a)->exponent=(a)->Prec=1;(a)->sign=VP_SIGN_POSITIVE_FINITE;}
+
+/* ZEROs */
+#define VpIsPosZero(a)  ((a)->sign==VP_SIGN_POSITIVE_ZERO)
+#define VpIsNegZero(a)  ((a)->sign==VP_SIGN_NEGATIVE_ZERO)
+#define VpIsZero(a)     (VpIsPosZero(a) || VpIsNegZero(a))
+#define VpSetPosZero(a) ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_POSITIVE_ZERO)
+#define VpSetNegZero(a) ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_NEGATIVE_ZERO)
+#define VpSetZero(a,s)  ( ((s)>0)?VpSetPosZero(a):VpSetNegZero(a) )
+
+/* NaN */
+#define VpIsNaN(a)      ((a)->sign==VP_SIGN_NaN)
+#define VpSetNaN(a)     ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_NaN)
+
+/* Infinity */
+#define VpIsPosInf(a)   ((a)->sign==VP_SIGN_POSITIVE_INFINITE)
+#define VpIsNegInf(a)   ((a)->sign==VP_SIGN_NEGATIVE_INFINITE)
+#define VpIsInf(a)      (VpIsPosInf(a) || VpIsNegInf(a))
+#define VpIsDef(a)      ( !(VpIsNaN(a)||VpIsInf(a)) )
+#define VpSetPosInf(a)  ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_POSITIVE_INFINITE)
+#define VpSetNegInf(a)  ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_NEGATIVE_INFINITE)
+#define VpSetInf(a,s)   ( ((s)>0)?VpSetPosInf(a):VpSetNegInf(a) )
+
+#ifdef _DEBUG
+int VpVarCheck(Real * v);
+#endif /* _DEBUG */
+
+#if defined(__cplusplus)
+}  /* extern "C" { */
+#endif
+#endif //____BIG_FLOAT__H____