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