16 #ifndef RUBY_BIG_DECIMAL_H
17 #define RUBY_BIG_DECIMAL_H 1
22 #if defined(__cplusplus)
36 static inline LONG_LONG
47 #if SIZEOF_VALUE <= SIZEOF_INT
49 #elif SIZEOF_VALUE <= SIZEOF_LONG
51 #elif SIZEOF_VALUE <= SIZEOF_LONG_LONG
57 #if 0 || SIZEOF_BDIGITS >= 16
58 # define RMPD_COMPONENT_FIGURES 38
59 # define RMPD_BASE ((BDIGIT)100000000000000000000000000000000000000U)
60 #elif SIZEOF_BDIGITS >= 8
61 # define RMPD_COMPONENT_FIGURES 19
62 # define RMPD_BASE ((BDIGIT)10000000000000000000U)
63 #elif SIZEOF_BDIGITS >= 4
64 # define RMPD_COMPONENT_FIGURES 9
65 # define RMPD_BASE ((BDIGIT)1000000000U)
66 #elif SIZEOF_BDIGITS >= 2
67 # define RMPD_COMPONENT_FIGURES 4
68 # define RMPD_BASE ((BDIGIT)10000U)
70 # define RMPD_COMPONENT_FIGURES 2
71 # define RMPD_BASE ((BDIGIT)100U)
79 #define SZ_INF "Infinity"
80 #define SZ_PINF "+Infinity"
81 #define SZ_NINF "-Infinity"
87 #define VP_EXPORT static
90 #define VP_EXCEPTION_ALL ((unsigned short)0x00FF)
91 #define VP_EXCEPTION_INFINITY ((unsigned short)0x0001)
92 #define VP_EXCEPTION_NaN ((unsigned short)0x0002)
93 #define VP_EXCEPTION_UNDERFLOW ((unsigned short)0x0004)
94 #define VP_EXCEPTION_OVERFLOW ((unsigned short)0x0001)
95 #define VP_EXCEPTION_ZERODIVIDE ((unsigned short)0x0010)
98 #define VP_EXCEPTION_OP ((unsigned short)0x0020)
99 #define VP_EXCEPTION_MEMORY ((unsigned short)0x0040)
101 #define RMPD_EXCEPTION_MODE_DEFAULT 0U
104 #define VP_ROUND_MODE ((unsigned short)0x0100)
105 #define VP_ROUND_UP 1
106 #define VP_ROUND_DOWN 2
107 #define VP_ROUND_HALF_UP 3
108 #define VP_ROUND_HALF_DOWN 4
109 #define VP_ROUND_CEIL 5
110 #define VP_ROUND_FLOOR 6
111 #define VP_ROUND_HALF_EVEN 7
113 #define RMPD_ROUNDING_MODE_DEFAULT VP_ROUND_HALF_UP
115 #define VP_SIGN_NaN 0
116 #define VP_SIGN_POSITIVE_ZERO 1
117 #define VP_SIGN_NEGATIVE_ZERO -1
118 #define VP_SIGN_POSITIVE_FINITE 2
119 #define VP_SIGN_NEGATIVE_FINITE -2
120 #define VP_SIGN_POSITIVE_INFINITE 3
121 #define VP_SIGN_NEGATIVE_INFINITE -3
168 #define VpBaseFig() rmpd_component_figures()
169 #define VpDblFig() rmpd_double_figures()
170 #define VpBaseVal() rmpd_base_value()
227 #define Abs(a) (((a)>= 0)?(a):(-(a)))
228 #define Max(a, b) (((a)>(b))?(a):(b))
229 #define Min(a, b) (((a)>(b))?(b):(a))
231 #define VpMaxPrec(a) ((a)->MaxPrec)
232 #define VpPrec(a) ((a)->Prec)
233 #define VpGetFlag(a) ((a)->flag)
238 #define VpGetSign(a) (((a)->sign>0)?1:(-1))
240 #define VpChangeSign(a,s) {if((s)>0) (a)->sign=(short)Abs((ssize_t)(a)->sign);else (a)->sign=-(short)Abs((ssize_t)(a)->sign);}
242 #define VpSetSign(a,s) {if((s)>0) (a)->sign=(short)VP_SIGN_POSITIVE_FINITE;else (a)->sign=(short)VP_SIGN_NEGATIVE_FINITE;}
245 #define VpSetOne(a) {(a)->Prec=(a)->exponent=(a)->frac[0]=1;(a)->sign=VP_SIGN_POSITIVE_FINITE;}
248 #define VpIsPosZero(a) ((a)->sign==VP_SIGN_POSITIVE_ZERO)
249 #define VpIsNegZero(a) ((a)->sign==VP_SIGN_NEGATIVE_ZERO)
250 #define VpIsZero(a) (VpIsPosZero(a) || VpIsNegZero(a))
251 #define VpSetPosZero(a) ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_POSITIVE_ZERO)
252 #define VpSetNegZero(a) ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_NEGATIVE_ZERO)
253 #define VpSetZero(a,s) ( ((s)>0)?VpSetPosZero(a):VpSetNegZero(a) )
256 #define VpIsNaN(a) ((a)->sign==VP_SIGN_NaN)
257 #define VpSetNaN(a) ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_NaN)
260 #define VpIsPosInf(a) ((a)->sign==VP_SIGN_POSITIVE_INFINITE)
261 #define VpIsNegInf(a) ((a)->sign==VP_SIGN_NEGATIVE_INFINITE)
262 #define VpIsInf(a) (VpIsPosInf(a) || VpIsNegInf(a))
263 #define VpIsDef(a) ( !(VpIsNaN(a)||VpIsInf(a)) )
264 #define VpSetPosInf(a) ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_POSITIVE_INFINITE)
265 #define VpSetNegInf(a) ((a)->frac[0]=0,(a)->Prec=1,(a)->sign=VP_SIGN_NEGATIVE_INFINITE)
266 #define VpSetInf(a,s) ( ((s)>0)?VpSetPosInf(a):VpSetNegInf(a) )
267 #define VpHasVal(a) (a->frac[0])
268 #define VpIsOne(a) ((a->Prec==1)&&(a->frac[0]==1)&&(a->exponent==1))
269 #define VpExponent(a) (a->exponent)
270 #ifdef BIGDECIMAL_DEBUG
271 int VpVarCheck(
Real *
v);
275 #if defined(__cplusplus)
VP_EXPORT size_t VpDivd(Real *c, Real *r, Real *a, Real *b)
VP_EXPORT double VpGetDoublePosInf(void)
VP_EXPORT int VpActiveRound(Real *y, Real *x, unsigned short f, ssize_t il)
VP_EXPORT int VpException(unsigned short f, const char *str, int always)
VP_EXPORT int VpIsRoundMode(unsigned short n)
static long labs(long const x)
VP_EXPORT unsigned short VpGetRoundMode(void)
VP_EXPORT unsigned short VpSetRoundMode(unsigned short n)
VP_EXPORT void VpDtoV(Real *m, double d)
VP_EXPORT void VpFrac(Real *y, Real *x)
VP_EXPORT Real * VpCreateRbObject(size_t mx, const char *str)
VP_EXPORT int VpVtoD(double *d, SIGNED_VALUE *e, Real *m)
VP_EXPORT int VpComp(Real *a, Real *b)
VP_EXPORT int VpCtoV(Real *a, const char *int_chr, size_t ni, const char *frac, size_t nf, const char *exp_chr, size_t ne)
VP_EXPORT Real * VpNewRbClass(size_t mx, const char *str, VALUE klass)
static LONG_LONG llabs(LONG_LONG const x)
static size_t rmpd_component_figures(void)
VP_EXPORT void VpToFString(Real *a, char *psz, size_t fFmt, int fPlus)
VP_EXPORT int VpMidRound(Real *y, unsigned short f, ssize_t nf)
#define RMPD_COMPONENT_FIGURES
VP_EXPORT void VpToString(Real *a, char *psz, size_t fFmt, int fPlus)
VP_EXPORT double VpGetDoubleNaN(void)
VP_EXPORT size_t VpGetPrecLimit(void)
VP_EXPORT void * VpMemAlloc(size_t mb)
VP_EXPORT size_t VpMult(Real *c, Real *a, Real *b)
VP_EXPORT Real * VpOne(void)
VP_EXPORT void VpSzMantissa(Real *a, char *psz)
static VALUE int_chr(int argc, VALUE *argv, VALUE num)
VP_EXPORT void VpFree(Real *pv)
VP_EXPORT Real * VpAlloc(size_t mx, const char *szVal)
VP_EXPORT int VpToSpecialString(Real *a, char *psz, int fPlus)
static BDIGIT rmpd_base_value(void)
VP_EXPORT size_t VpInit(BDIGIT BaseVal)
VP_EXPORT double VpGetDoubleNegInf(void)
VP_EXPORT size_t VpAsgn(Real *c, Real *a, int isw)
static size_t rmpd_double_figures(void)
VP_EXPORT size_t VpAddSub(Real *c, Real *a, Real *b, int operation)
VP_EXPORT int VpLeftRound(Real *y, unsigned short f, ssize_t nf)
VP_EXPORT int VpPower(Real *y, Real *x, SIGNED_VALUE n)
VP_EXPORT size_t VpSetPrecLimit(size_t n)
VP_EXPORT size_t VpNumOfChars(Real *vp, const char *pszFmt)
VP_EXPORT double VpGetDoubleNegZero(void)
VP_EXPORT ssize_t VpExponent10(Real *a)
VP_EXPORT int VpSqrt(Real *y, Real *x)