504 lines
12 KiB
C
504 lines
12 KiB
C
/*
|
|
*
|
|
* 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____
|