1*5796c8dcSSimon Schubert /* Decimal number arithmetic module header for the decNumber C Library. 2*5796c8dcSSimon Schubert Copyright (C) 2005, 2007, 2009 Free Software Foundation, Inc. 3*5796c8dcSSimon Schubert Contributed by IBM Corporation. Author Mike Cowlishaw. 4*5796c8dcSSimon Schubert 5*5796c8dcSSimon Schubert This file is part of GCC. 6*5796c8dcSSimon Schubert 7*5796c8dcSSimon Schubert GCC is free software; you can redistribute it and/or modify it under 8*5796c8dcSSimon Schubert the terms of the GNU General Public License as published by the Free 9*5796c8dcSSimon Schubert Software Foundation; either version 3, or (at your option) any later 10*5796c8dcSSimon Schubert version. 11*5796c8dcSSimon Schubert 12*5796c8dcSSimon Schubert GCC is distributed in the hope that it will be useful, but WITHOUT ANY 13*5796c8dcSSimon Schubert WARRANTY; without even the implied warranty of MERCHANTABILITY or 14*5796c8dcSSimon Schubert FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 15*5796c8dcSSimon Schubert for more details. 16*5796c8dcSSimon Schubert 17*5796c8dcSSimon Schubert Under Section 7 of GPL version 3, you are granted additional 18*5796c8dcSSimon Schubert permissions described in the GCC Runtime Library Exception, version 19*5796c8dcSSimon Schubert 3.1, as published by the Free Software Foundation. 20*5796c8dcSSimon Schubert 21*5796c8dcSSimon Schubert You should have received a copy of the GNU General Public License and 22*5796c8dcSSimon Schubert a copy of the GCC Runtime Library Exception along with this program; 23*5796c8dcSSimon Schubert see the files COPYING3 and COPYING.RUNTIME respectively. If not, see 24*5796c8dcSSimon Schubert <http://www.gnu.org/licenses/>. */ 25*5796c8dcSSimon Schubert 26*5796c8dcSSimon Schubert /* ------------------------------------------------------------------ */ 27*5796c8dcSSimon Schubert /* Decimal Number arithmetic module header */ 28*5796c8dcSSimon Schubert /* ------------------------------------------------------------------ */ 29*5796c8dcSSimon Schubert 30*5796c8dcSSimon Schubert #if !defined(DECNUMBER) 31*5796c8dcSSimon Schubert #define DECNUMBER 32*5796c8dcSSimon Schubert #define DECNAME "decNumber" /* Short name */ 33*5796c8dcSSimon Schubert #define DECFULLNAME "Decimal Number Module" /* Verbose name */ 34*5796c8dcSSimon Schubert #define DECAUTHOR "Mike Cowlishaw" /* Who to blame */ 35*5796c8dcSSimon Schubert 36*5796c8dcSSimon Schubert #if !defined(DECCONTEXT) 37*5796c8dcSSimon Schubert #include "decContext.h" 38*5796c8dcSSimon Schubert #endif 39*5796c8dcSSimon Schubert 40*5796c8dcSSimon Schubert /* Bit settings for decNumber.bits */ 41*5796c8dcSSimon Schubert #define DECNEG 0x80 /* Sign; 1=negative, 0=positive or zero */ 42*5796c8dcSSimon Schubert #define DECINF 0x40 /* 1=Infinity */ 43*5796c8dcSSimon Schubert #define DECNAN 0x20 /* 1=NaN */ 44*5796c8dcSSimon Schubert #define DECSNAN 0x10 /* 1=sNaN */ 45*5796c8dcSSimon Schubert /* The remaining bits are reserved; they must be 0 */ 46*5796c8dcSSimon Schubert #define DECSPECIAL (DECINF|DECNAN|DECSNAN) /* any special value */ 47*5796c8dcSSimon Schubert 48*5796c8dcSSimon Schubert /* Define the decNumber data structure. The size and shape of the */ 49*5796c8dcSSimon Schubert /* units array in the structure is determined by the following */ 50*5796c8dcSSimon Schubert /* constant. This must not be changed without recompiling the */ 51*5796c8dcSSimon Schubert /* decNumber library modules. */ 52*5796c8dcSSimon Schubert 53*5796c8dcSSimon Schubert #define DECDPUN 3 /* DECimal Digits Per UNit [must be >0 */ 54*5796c8dcSSimon Schubert /* and <10; 3 or powers of 2 are best]. */ 55*5796c8dcSSimon Schubert 56*5796c8dcSSimon Schubert /* DECNUMDIGITS is the default number of digits that can be held in */ 57*5796c8dcSSimon Schubert /* the structure. If undefined, 1 is assumed and it is assumed */ 58*5796c8dcSSimon Schubert /* that the structure will be immediately followed by extra space, */ 59*5796c8dcSSimon Schubert /* as required. DECNUMDIGITS is always >0. */ 60*5796c8dcSSimon Schubert #if !defined(DECNUMDIGITS) 61*5796c8dcSSimon Schubert #define DECNUMDIGITS 1 62*5796c8dcSSimon Schubert #endif 63*5796c8dcSSimon Schubert 64*5796c8dcSSimon Schubert /* The size (integer data type) of each unit is determined by the */ 65*5796c8dcSSimon Schubert /* number of digits it will hold. */ 66*5796c8dcSSimon Schubert #if DECDPUN<=2 67*5796c8dcSSimon Schubert #define decNumberUnit uint8_t 68*5796c8dcSSimon Schubert #elif DECDPUN<=4 69*5796c8dcSSimon Schubert #define decNumberUnit uint16_t 70*5796c8dcSSimon Schubert #else 71*5796c8dcSSimon Schubert #define decNumberUnit uint32_t 72*5796c8dcSSimon Schubert #endif 73*5796c8dcSSimon Schubert /* The number of units needed is ceil(DECNUMDIGITS/DECDPUN) */ 74*5796c8dcSSimon Schubert #define DECNUMUNITS ((DECNUMDIGITS+DECDPUN-1)/DECDPUN) 75*5796c8dcSSimon Schubert 76*5796c8dcSSimon Schubert /* The data structure... */ 77*5796c8dcSSimon Schubert typedef struct { 78*5796c8dcSSimon Schubert int32_t digits; /* Count of digits in the coefficient; >0 */ 79*5796c8dcSSimon Schubert int32_t exponent; /* Unadjusted exponent, unbiased, in */ 80*5796c8dcSSimon Schubert /* range: -1999999997 through 999999999 */ 81*5796c8dcSSimon Schubert uint8_t bits; /* Indicator bits (see above) */ 82*5796c8dcSSimon Schubert /* Coefficient, from least significant unit */ 83*5796c8dcSSimon Schubert decNumberUnit lsu[DECNUMUNITS]; 84*5796c8dcSSimon Schubert } decNumber; 85*5796c8dcSSimon Schubert 86*5796c8dcSSimon Schubert /* Notes: */ 87*5796c8dcSSimon Schubert /* 1. If digits is > DECDPUN then there will one or more */ 88*5796c8dcSSimon Schubert /* decNumberUnits immediately following the first element of lsu.*/ 89*5796c8dcSSimon Schubert /* These contain the remaining (more significant) digits of the */ 90*5796c8dcSSimon Schubert /* number, and may be in the lsu array, or may be guaranteed by */ 91*5796c8dcSSimon Schubert /* some other mechanism (such as being contained in another */ 92*5796c8dcSSimon Schubert /* structure, or being overlaid on dynamically allocated */ 93*5796c8dcSSimon Schubert /* storage). */ 94*5796c8dcSSimon Schubert /* */ 95*5796c8dcSSimon Schubert /* Each integer of the coefficient (except potentially the last) */ 96*5796c8dcSSimon Schubert /* contains DECDPUN digits (e.g., a value in the range 0 through */ 97*5796c8dcSSimon Schubert /* 99999999 if DECDPUN is 8, or 0 through 999 if DECDPUN is 3). */ 98*5796c8dcSSimon Schubert /* */ 99*5796c8dcSSimon Schubert /* 2. A decNumber converted to a string may need up to digits+14 */ 100*5796c8dcSSimon Schubert /* characters. The worst cases (non-exponential and exponential */ 101*5796c8dcSSimon Schubert /* formats) are -0.00000{9...}# and -9.{9...}E+999999999# */ 102*5796c8dcSSimon Schubert /* (where # is '\0') */ 103*5796c8dcSSimon Schubert 104*5796c8dcSSimon Schubert 105*5796c8dcSSimon Schubert /* ---------------------------------------------------------------- */ 106*5796c8dcSSimon Schubert /* decNumber public functions and macros */ 107*5796c8dcSSimon Schubert /* ---------------------------------------------------------------- */ 108*5796c8dcSSimon Schubert 109*5796c8dcSSimon Schubert #include "decNumberSymbols.h" 110*5796c8dcSSimon Schubert 111*5796c8dcSSimon Schubert #ifdef __cplusplus 112*5796c8dcSSimon Schubert extern "C" { 113*5796c8dcSSimon Schubert #endif 114*5796c8dcSSimon Schubert 115*5796c8dcSSimon Schubert /* Conversions */ 116*5796c8dcSSimon Schubert decNumber * decNumberFromInt32(decNumber *, int32_t); 117*5796c8dcSSimon Schubert decNumber * decNumberFromUInt32(decNumber *, uint32_t); 118*5796c8dcSSimon Schubert decNumber * decNumberFromString(decNumber *, const char *, decContext *); 119*5796c8dcSSimon Schubert char * decNumberToString(const decNumber *, char *); 120*5796c8dcSSimon Schubert char * decNumberToEngString(const decNumber *, char *); 121*5796c8dcSSimon Schubert uint32_t decNumberToUInt32(const decNumber *, decContext *); 122*5796c8dcSSimon Schubert int32_t decNumberToInt32(const decNumber *, decContext *); 123*5796c8dcSSimon Schubert uint8_t * decNumberGetBCD(const decNumber *, uint8_t *); 124*5796c8dcSSimon Schubert decNumber * decNumberSetBCD(decNumber *, const uint8_t *, uint32_t); 125*5796c8dcSSimon Schubert 126*5796c8dcSSimon Schubert /* Operators and elementary functions */ 127*5796c8dcSSimon Schubert decNumber * decNumberAbs(decNumber *, const decNumber *, decContext *); 128*5796c8dcSSimon Schubert decNumber * decNumberAdd(decNumber *, const decNumber *, const decNumber *, decContext *); 129*5796c8dcSSimon Schubert decNumber * decNumberAnd(decNumber *, const decNumber *, const decNumber *, decContext *); 130*5796c8dcSSimon Schubert decNumber * decNumberCompare(decNumber *, const decNumber *, const decNumber *, decContext *); 131*5796c8dcSSimon Schubert decNumber * decNumberCompareSignal(decNumber *, const decNumber *, const decNumber *, decContext *); 132*5796c8dcSSimon Schubert decNumber * decNumberCompareTotal(decNumber *, const decNumber *, const decNumber *, decContext *); 133*5796c8dcSSimon Schubert decNumber * decNumberCompareTotalMag(decNumber *, const decNumber *, const decNumber *, decContext *); 134*5796c8dcSSimon Schubert decNumber * decNumberDivide(decNumber *, const decNumber *, const decNumber *, decContext *); 135*5796c8dcSSimon Schubert decNumber * decNumberDivideInteger(decNumber *, const decNumber *, const decNumber *, decContext *); 136*5796c8dcSSimon Schubert decNumber * decNumberExp(decNumber *, const decNumber *, decContext *); 137*5796c8dcSSimon Schubert decNumber * decNumberFMA(decNumber *, const decNumber *, const decNumber *, const decNumber *, decContext *); 138*5796c8dcSSimon Schubert decNumber * decNumberInvert(decNumber *, const decNumber *, decContext *); 139*5796c8dcSSimon Schubert decNumber * decNumberLn(decNumber *, const decNumber *, decContext *); 140*5796c8dcSSimon Schubert decNumber * decNumberLogB(decNumber *, const decNumber *, decContext *); 141*5796c8dcSSimon Schubert decNumber * decNumberLog10(decNumber *, const decNumber *, decContext *); 142*5796c8dcSSimon Schubert decNumber * decNumberMax(decNumber *, const decNumber *, const decNumber *, decContext *); 143*5796c8dcSSimon Schubert decNumber * decNumberMaxMag(decNumber *, const decNumber *, const decNumber *, decContext *); 144*5796c8dcSSimon Schubert decNumber * decNumberMin(decNumber *, const decNumber *, const decNumber *, decContext *); 145*5796c8dcSSimon Schubert decNumber * decNumberMinMag(decNumber *, const decNumber *, const decNumber *, decContext *); 146*5796c8dcSSimon Schubert decNumber * decNumberMinus(decNumber *, const decNumber *, decContext *); 147*5796c8dcSSimon Schubert decNumber * decNumberMultiply(decNumber *, const decNumber *, const decNumber *, decContext *); 148*5796c8dcSSimon Schubert decNumber * decNumberNormalize(decNumber *, const decNumber *, decContext *); 149*5796c8dcSSimon Schubert decNumber * decNumberOr(decNumber *, const decNumber *, const decNumber *, decContext *); 150*5796c8dcSSimon Schubert decNumber * decNumberPlus(decNumber *, const decNumber *, decContext *); 151*5796c8dcSSimon Schubert decNumber * decNumberPower(decNumber *, const decNumber *, const decNumber *, decContext *); 152*5796c8dcSSimon Schubert decNumber * decNumberQuantize(decNumber *, const decNumber *, const decNumber *, decContext *); 153*5796c8dcSSimon Schubert decNumber * decNumberReduce(decNumber *, const decNumber *, decContext *); 154*5796c8dcSSimon Schubert decNumber * decNumberRemainder(decNumber *, const decNumber *, const decNumber *, decContext *); 155*5796c8dcSSimon Schubert decNumber * decNumberRemainderNear(decNumber *, const decNumber *, const decNumber *, decContext *); 156*5796c8dcSSimon Schubert decNumber * decNumberRescale(decNumber *, const decNumber *, const decNumber *, decContext *); 157*5796c8dcSSimon Schubert decNumber * decNumberRotate(decNumber *, const decNumber *, const decNumber *, decContext *); 158*5796c8dcSSimon Schubert decNumber * decNumberSameQuantum(decNumber *, const decNumber *, const decNumber *); 159*5796c8dcSSimon Schubert decNumber * decNumberScaleB(decNumber *, const decNumber *, const decNumber *, decContext *); 160*5796c8dcSSimon Schubert decNumber * decNumberShift(decNumber *, const decNumber *, const decNumber *, decContext *); 161*5796c8dcSSimon Schubert decNumber * decNumberSquareRoot(decNumber *, const decNumber *, decContext *); 162*5796c8dcSSimon Schubert decNumber * decNumberSubtract(decNumber *, const decNumber *, const decNumber *, decContext *); 163*5796c8dcSSimon Schubert decNumber * decNumberToIntegralExact(decNumber *, const decNumber *, decContext *); 164*5796c8dcSSimon Schubert decNumber * decNumberToIntegralValue(decNumber *, const decNumber *, decContext *); 165*5796c8dcSSimon Schubert decNumber * decNumberXor(decNumber *, const decNumber *, const decNumber *, decContext *); 166*5796c8dcSSimon Schubert 167*5796c8dcSSimon Schubert /* Utilities */ 168*5796c8dcSSimon Schubert enum decClass decNumberClass(const decNumber *, decContext *); 169*5796c8dcSSimon Schubert const char * decNumberClassToString(enum decClass); 170*5796c8dcSSimon Schubert decNumber * decNumberCopy(decNumber *, const decNumber *); 171*5796c8dcSSimon Schubert decNumber * decNumberCopyAbs(decNumber *, const decNumber *); 172*5796c8dcSSimon Schubert decNumber * decNumberCopyNegate(decNumber *, const decNumber *); 173*5796c8dcSSimon Schubert decNumber * decNumberCopySign(decNumber *, const decNumber *, const decNumber *); 174*5796c8dcSSimon Schubert decNumber * decNumberNextMinus(decNumber *, const decNumber *, decContext *); 175*5796c8dcSSimon Schubert decNumber * decNumberNextPlus(decNumber *, const decNumber *, decContext *); 176*5796c8dcSSimon Schubert decNumber * decNumberNextToward(decNumber *, const decNumber *, const decNumber *, decContext *); 177*5796c8dcSSimon Schubert decNumber * decNumberTrim(decNumber *); 178*5796c8dcSSimon Schubert const char * decNumberVersion(void); 179*5796c8dcSSimon Schubert decNumber * decNumberZero(decNumber *); 180*5796c8dcSSimon Schubert 181*5796c8dcSSimon Schubert /* Functions for testing decNumbers (normality depends on context) */ 182*5796c8dcSSimon Schubert int32_t decNumberIsNormal(const decNumber *, decContext *); 183*5796c8dcSSimon Schubert int32_t decNumberIsSubnormal(const decNumber *, decContext *); 184*5796c8dcSSimon Schubert 185*5796c8dcSSimon Schubert /* Macros for testing decNumber *dn */ 186*5796c8dcSSimon Schubert #define decNumberIsCanonical(dn) (1) /* All decNumbers are saintly */ 187*5796c8dcSSimon Schubert #define decNumberIsFinite(dn) (((dn)->bits&DECSPECIAL)==0) 188*5796c8dcSSimon Schubert #define decNumberIsInfinite(dn) (((dn)->bits&DECINF)!=0) 189*5796c8dcSSimon Schubert #define decNumberIsNaN(dn) (((dn)->bits&(DECNAN|DECSNAN))!=0) 190*5796c8dcSSimon Schubert #define decNumberIsNegative(dn) (((dn)->bits&DECNEG)!=0) 191*5796c8dcSSimon Schubert #define decNumberIsQNaN(dn) (((dn)->bits&(DECNAN))!=0) 192*5796c8dcSSimon Schubert #define decNumberIsSNaN(dn) (((dn)->bits&(DECSNAN))!=0) 193*5796c8dcSSimon Schubert #define decNumberIsSpecial(dn) (((dn)->bits&DECSPECIAL)!=0) 194*5796c8dcSSimon Schubert #define decNumberIsZero(dn) (*(dn)->lsu==0 \ 195*5796c8dcSSimon Schubert && (dn)->digits==1 \ 196*5796c8dcSSimon Schubert && (((dn)->bits&DECSPECIAL)==0)) 197*5796c8dcSSimon Schubert #define decNumberRadix(dn) (10) 198*5796c8dcSSimon Schubert 199*5796c8dcSSimon Schubert #ifdef __cplusplus 200*5796c8dcSSimon Schubert } 201*5796c8dcSSimon Schubert #endif 202*5796c8dcSSimon Schubert 203*5796c8dcSSimon Schubert #endif 204