1/* 2 * Copyright (c) 2008-2024 Stefan Krah. All rights reserved. 3 * 4 * Redistribution and use in source and binary forms, with or without 5 * modification, are permitted provided that the following conditions 6 * are met: 7 * 8 * 1. Redistributions of source code must retain the above copyright 9 * notice, this list of conditions and the following disclaimer. 10 * 2. Redistributions in binary form must reproduce the above copyright 11 * notice, this list of conditions and the following disclaimer in the 12 * documentation and/or other materials provided with the distribution. 13 * 14 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND 15 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE 16 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE 17 * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE 18 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL 19 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS 20 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) 21 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 22 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY 23 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF 24 * SUCH DAMAGE. 25 */ 26 27 28#ifndef LIBMPDEC_MPDECIMAL_H_ 29#define LIBMPDEC_MPDECIMAL_H_ 30 31 32#ifdef __cplusplus 33 #include <cinttypes> 34 #include <climits> 35 #include <cstdint> 36 #include <cstdio> 37 #include <cstdlib> 38 #define MPD_UINT8_C(x) (static_cast<uint8_t>(x)) 39extern "C" { 40#else 41 #include <inttypes.h> 42 #include <limits.h> 43 #include <stdint.h> 44 #include <stdio.h> 45 #include <stdlib.h> 46 #define MPD_UINT8_C(x) ((uint8_t)x) 47#endif 48 49 50#if (defined(__linux__) || defined(__FreeBSD__) || defined(__APPLE__)) && \ 51 defined(__GNUC__) && __GNUC__ >= 4 && !defined(__INTEL_COMPILER) 52 #define MPD_PRAGMA(x) _Pragma(x) 53 #define MPD_HIDE_SYMBOLS_START "GCC visibility push(hidden)" 54 #define MPD_HIDE_SYMBOLS_END "GCC visibility pop" 55#else 56 #define MPD_PRAGMA(x) 57 #define MPD_HIDE_SYMBOLS_START 58 #define MPD_HIDE_SYMBOLS_END 59#endif 60 61 62/******************************************************************************/ 63/* Version */ 64/******************************************************************************/ 65 66#define MPD_MAJOR_VERSION 4 67#define MPD_MINOR_VERSION 0 68#define MPD_MICRO_VERSION 0 69 70#define MPD_VERSION "4.0.0" 71 72#define MPD_VERSION_HEX ((MPD_MAJOR_VERSION << 24) | \ 73 (MPD_MINOR_VERSION << 16) | \ 74 (MPD_MICRO_VERSION << 8)) 75 76const char *mpd_version(void); 77 78 79/******************************************************************************/ 80/* Configuration */ 81/******************************************************************************/ 82@MPD_HEADER_CONFIG@ 83 84 85/* BEGIN MPD_CONFIG_64 */ 86#if defined(MPD_CONFIG_64) 87/* types for modular and base arithmetic */ 88#define MPD_UINT_MAX UINT64_MAX 89#define MPD_BITS_PER_UINT 64 90typedef uint64_t mpd_uint_t; /* unsigned mod type */ 91 92#define MPD_SIZE_MAX SIZE_MAX 93typedef size_t mpd_size_t; /* unsigned size type */ 94 95/* type for exp, digits, len, prec */ 96#define MPD_SSIZE_MAX INT64_MAX 97#define MPD_SSIZE_MIN INT64_MIN 98typedef int64_t mpd_ssize_t; 99#define _mpd_strtossize strtoll 100 101/* decimal arithmetic */ 102#define MPD_RADIX 10000000000000000000ULL /* 10**19 */ 103#define MPD_RDIGITS 19 104#define MPD_MAX_POW10 19 105#define MPD_EXPDIGITS 19 /* MPD_EXPDIGITS <= MPD_RDIGITS+1 */ 106 107#define MPD_MAXTRANSFORM_2N 4294967296ULL /* 2**32 */ 108#define MPD_MAX_PREC 999999999999999999LL 109#define MPD_MAX_PREC_LOG2 64 110#define MPD_ELIMIT 1000000000000000000LL 111#define MPD_MAX_EMAX 999999999999999999LL /* ELIMIT-1 */ 112#define MPD_MIN_EMIN (-999999999999999999LL) /* -EMAX */ 113#define MPD_MIN_ETINY (MPD_MIN_EMIN-(MPD_MAX_PREC-1)) 114#define MPD_EXP_INF 2000000000000000001LL 115#define MPD_EXP_CLAMP (-4000000000000000001LL) 116#define MPD_MAXIMPORT 105263157894736842L /* ceil((2*MPD_MAX_PREC)/MPD_RDIGITS) */ 117#define MPD_IEEE_CONTEXT_MAX_BITS 512 /* 16*(log2(MPD_MAX_EMAX / 3)-3) */ 118 119/* conversion specifiers */ 120#define PRI_mpd_uint_t PRIu64 121#define PRI_mpd_ssize_t PRIi64 122/* END MPD_CONFIG_64 */ 123 124 125/* BEGIN MPD_CONFIG_32 */ 126#elif defined(MPD_CONFIG_32) 127/* types for modular and base arithmetic */ 128#define MPD_UINT_MAX UINT32_MAX 129#define MPD_BITS_PER_UINT 32 130typedef uint32_t mpd_uint_t; /* unsigned mod type */ 131 132#ifndef MPD_LEGACY_COMPILER 133#define MPD_UUINT_MAX UINT64_MAX 134typedef uint64_t mpd_uuint_t; /* double width unsigned mod type */ 135#endif 136 137#define MPD_SIZE_MAX SIZE_MAX 138typedef size_t mpd_size_t; /* unsigned size type */ 139 140/* type for dec->len, dec->exp, ctx->prec */ 141#define MPD_SSIZE_MAX INT32_MAX 142#define MPD_SSIZE_MIN INT32_MIN 143typedef int32_t mpd_ssize_t; 144#define _mpd_strtossize strtol 145 146/* decimal arithmetic */ 147#define MPD_RADIX 1000000000UL /* 10**9 */ 148#define MPD_RDIGITS 9 149#define MPD_MAX_POW10 9 150#define MPD_EXPDIGITS 10 /* MPD_EXPDIGITS <= MPD_RDIGITS+1 */ 151 152#define MPD_MAXTRANSFORM_2N 33554432UL /* 2**25 */ 153#define MPD_MAX_PREC 425000000L 154#define MPD_MAX_PREC_LOG2 32 155#define MPD_ELIMIT 425000001L 156#define MPD_MAX_EMAX 425000000L /* ELIMIT-1 */ 157#define MPD_MIN_EMIN (-425000000L) /* -EMAX */ 158#define MPD_MIN_ETINY (MPD_MIN_EMIN-(MPD_MAX_PREC-1)) 159#define MPD_EXP_INF 1000000001L /* allows for emax=999999999 in the tests */ 160#define MPD_EXP_CLAMP (-2000000001L) /* allows for emin=-999999999 in the tests */ 161#define MPD_MAXIMPORT 94444445L /* ceil((2*MPD_MAX_PREC)/MPD_RDIGITS) */ 162#define MPD_IEEE_CONTEXT_MAX_BITS 256 /* 16*(log2(MPD_MAX_EMAX / 3)-3) */ 163 164/* conversion specifiers */ 165#define PRI_mpd_uint_t PRIu32 166#define PRI_mpd_ssize_t PRIi32 167/* END MPD_CONFIG_32 */ 168 169#else 170 #error "define MPD_CONFIG_64 or MPD_CONFIG_32" 171#endif 172/* END CONFIG */ 173 174 175#if MPD_SIZE_MAX != MPD_UINT_MAX 176 #error "unsupported platform: need mpd_size_t == mpd_uint_t" 177#endif 178 179 180/******************************************************************************/ 181/* Context */ 182/******************************************************************************/ 183 184enum { 185 MPD_ROUND_UP, /* round away from 0 */ 186 MPD_ROUND_DOWN, /* round toward 0 (truncate) */ 187 MPD_ROUND_CEILING, /* round toward +infinity */ 188 MPD_ROUND_FLOOR, /* round toward -infinity */ 189 MPD_ROUND_HALF_UP, /* 0.5 is rounded up */ 190 MPD_ROUND_HALF_DOWN, /* 0.5 is rounded down */ 191 MPD_ROUND_HALF_EVEN, /* 0.5 is rounded to even */ 192 MPD_ROUND_05UP, /* round zero or five away from 0 */ 193 MPD_ROUND_TRUNC, /* truncate, but set infinity */ 194 MPD_ROUND_GUARD 195}; 196 197enum { MPD_CLAMP_DEFAULT, MPD_CLAMP_IEEE_754, MPD_CLAMP_GUARD }; 198 199extern const char * const mpd_round_string[MPD_ROUND_GUARD]; 200extern const char * const mpd_clamp_string[MPD_CLAMP_GUARD]; 201 202 203typedef struct mpd_context_t { 204 mpd_ssize_t prec; /* precision */ 205 mpd_ssize_t emax; /* max positive exp */ 206 mpd_ssize_t emin; /* min negative exp */ 207 uint32_t traps; /* status events that should be trapped */ 208 uint32_t status; /* status flags */ 209 uint32_t newtrap; /* set by mpd_addstatus_raise() */ 210 int round; /* rounding mode */ 211 int clamp; /* clamp mode */ 212 int allcr; /* all functions correctly rounded */ 213} mpd_context_t; 214 215 216/* Status flags */ 217#define MPD_Clamped 0x00000001U 218#define MPD_Conversion_syntax 0x00000002U 219#define MPD_Division_by_zero 0x00000004U 220#define MPD_Division_impossible 0x00000008U 221#define MPD_Division_undefined 0x00000010U 222#define MPD_Fpu_error 0x00000020U 223#define MPD_Inexact 0x00000040U 224#define MPD_Invalid_context 0x00000080U 225#define MPD_Invalid_operation 0x00000100U 226#define MPD_Malloc_error 0x00000200U 227#define MPD_Not_implemented 0x00000400U 228#define MPD_Overflow 0x00000800U 229#define MPD_Rounded 0x00001000U 230#define MPD_Subnormal 0x00002000U 231#define MPD_Underflow 0x00004000U 232#define MPD_Max_status (0x00008000U-1U) 233 234/* Conditions that result in an IEEE 754 exception */ 235#define MPD_IEEE_Invalid_operation (MPD_Conversion_syntax | \ 236 MPD_Division_impossible | \ 237 MPD_Division_undefined | \ 238 MPD_Fpu_error | \ 239 MPD_Invalid_context | \ 240 MPD_Invalid_operation | \ 241 MPD_Malloc_error) \ 242 243/* Errors that require the result of an operation to be set to NaN */ 244#define MPD_Errors (MPD_IEEE_Invalid_operation | \ 245 MPD_Division_by_zero) 246 247/* Default traps */ 248#define MPD_Traps (MPD_IEEE_Invalid_operation | \ 249 MPD_Division_by_zero | \ 250 MPD_Overflow | \ 251 MPD_Underflow) 252 253/* Official name */ 254#define MPD_Insufficient_storage MPD_Malloc_error 255 256/* IEEE 754 interchange format contexts */ 257#define MPD_DECIMAL32 32 258#define MPD_DECIMAL64 64 259#define MPD_DECIMAL128 128 260 261 262#define MPD_MINALLOC_MIN 2 263#define MPD_MINALLOC_MAX 64 264extern mpd_ssize_t MPD_MINALLOC; 265extern void (* mpd_traphandler)(mpd_context_t *); 266void mpd_dflt_traphandler(mpd_context_t *); 267 268void mpd_setminalloc(mpd_ssize_t n); 269void mpd_init(mpd_context_t *ctx, mpd_ssize_t prec); 270 271void mpd_maxcontext(mpd_context_t *ctx); 272void mpd_defaultcontext(mpd_context_t *ctx); 273void mpd_basiccontext(mpd_context_t *ctx); 274int mpd_ieee_context(mpd_context_t *ctx, int bits); 275 276mpd_ssize_t mpd_getprec(const mpd_context_t *ctx); 277mpd_ssize_t mpd_getemax(const mpd_context_t *ctx); 278mpd_ssize_t mpd_getemin(const mpd_context_t *ctx); 279int mpd_getround(const mpd_context_t *ctx); 280uint32_t mpd_gettraps(const mpd_context_t *ctx); 281uint32_t mpd_getstatus(const mpd_context_t *ctx); 282int mpd_getclamp(const mpd_context_t *ctx); 283int mpd_getcr(const mpd_context_t *ctx); 284 285int mpd_qsetprec(mpd_context_t *ctx, mpd_ssize_t prec); 286int mpd_qsetemax(mpd_context_t *ctx, mpd_ssize_t emax); 287int mpd_qsetemin(mpd_context_t *ctx, mpd_ssize_t emin); 288int mpd_qsetround(mpd_context_t *ctx, int newround); 289int mpd_qsettraps(mpd_context_t *ctx, uint32_t flags); 290int mpd_qsetstatus(mpd_context_t *ctx, uint32_t flags); 291int mpd_qsetclamp(mpd_context_t *ctx, int c); 292int mpd_qsetcr(mpd_context_t *ctx, int c); 293void mpd_addstatus_raise(mpd_context_t *ctx, uint32_t flags); 294 295 296/******************************************************************************/ 297/* Decimal Arithmetic */ 298/******************************************************************************/ 299 300/* mpd_t flags */ 301#define MPD_POS MPD_UINT8_C(0) 302#define MPD_NEG MPD_UINT8_C(1) 303#define MPD_INF MPD_UINT8_C(2) 304#define MPD_NAN MPD_UINT8_C(4) 305#define MPD_SNAN MPD_UINT8_C(8) 306#define MPD_SPECIAL (MPD_INF|MPD_NAN|MPD_SNAN) 307#define MPD_STATIC MPD_UINT8_C(16) 308#define MPD_STATIC_DATA MPD_UINT8_C(32) 309#define MPD_SHARED_DATA MPD_UINT8_C(64) 310#define MPD_CONST_DATA MPD_UINT8_C(128) 311#define MPD_DATAFLAGS (MPD_STATIC_DATA|MPD_SHARED_DATA|MPD_CONST_DATA) 312 313/* mpd_t */ 314typedef struct mpd_t { 315 uint8_t flags; 316 mpd_ssize_t exp; 317 mpd_ssize_t digits; 318 mpd_ssize_t len; 319 mpd_ssize_t alloc; 320 mpd_uint_t *data; 321} mpd_t; 322 323 324/******************************************************************************/ 325/* Triple */ 326/******************************************************************************/ 327 328/* status cases for getting a triple */ 329enum mpd_triple_class { 330 MPD_TRIPLE_NORMAL, 331 MPD_TRIPLE_INF, 332 MPD_TRIPLE_QNAN, 333 MPD_TRIPLE_SNAN, 334 MPD_TRIPLE_ERROR, 335}; 336 337typedef struct { 338 enum mpd_triple_class tag; 339 uint8_t sign; 340 uint64_t hi; 341 uint64_t lo; 342 int64_t exp; 343} mpd_uint128_triple_t; 344 345int mpd_from_uint128_triple(mpd_t *result, const mpd_uint128_triple_t *triple, uint32_t *status); 346mpd_uint128_triple_t mpd_as_uint128_triple(const mpd_t *a); 347 348 349/******************************************************************************/ 350/* Quiet, thread-safe functions */ 351/******************************************************************************/ 352 353/* format specification */ 354typedef struct mpd_spec_t { 355 mpd_ssize_t min_width; /* minimum field width */ 356 mpd_ssize_t prec; /* fraction digits or significant digits */ 357 char type; /* conversion specifier */ 358 char align; /* alignment */ 359 char sign; /* sign printing/alignment */ 360 char sign_coerce; /* coerce to positive zero */ 361 char fill[5]; /* fill character */ 362 const char *dot; /* decimal point */ 363 const char *sep; /* thousands separator */ 364 const char *grouping; /* grouping of digits */ 365} mpd_spec_t; 366 367/* output to a string */ 368char *mpd_to_sci(const mpd_t *dec, int fmt); 369char *mpd_to_eng(const mpd_t *dec, int fmt); 370mpd_ssize_t mpd_to_sci_size(char **res, const mpd_t *dec, int fmt); 371mpd_ssize_t mpd_to_eng_size(char **res, const mpd_t *dec, int fmt); 372int mpd_validate_lconv(mpd_spec_t *spec); 373int mpd_parse_fmt_str(mpd_spec_t *spec, const char *fmt, int caps); 374char *mpd_qformat_spec(const mpd_t *dec, const mpd_spec_t *spec, const mpd_context_t *ctx, uint32_t *status); 375char *mpd_qformat(const mpd_t *dec, const char *fmt, const mpd_context_t *ctx, uint32_t *status); 376 377#define MPD_NUM_FLAGS 15 378#define MPD_MAX_FLAG_STRING 208 379#define MPD_MAX_FLAG_LIST (MPD_MAX_FLAG_STRING+18) 380#define MPD_MAX_SIGNAL_LIST 121 381int mpd_snprint_flags(char *dest, int nmemb, uint32_t flags); 382int mpd_lsnprint_flags(char *dest, int nmemb, uint32_t flags, const char *flag_string[]); 383int mpd_lsnprint_signals(char *dest, int nmemb, uint32_t flags, const char *signal_string[]); 384 385/* output to a file */ 386void mpd_fprint(FILE *file, const mpd_t *dec); 387void mpd_print(const mpd_t *dec); 388 389/* assignment from a string */ 390void mpd_qset_string(mpd_t *dec, const char *s, const mpd_context_t *ctx, uint32_t *status); 391void mpd_qset_string_exact(mpd_t *dec, const char *s, uint32_t *status); 392 393/* set to NaN with error flags */ 394void mpd_seterror(mpd_t *result, uint32_t flags, uint32_t *status); 395/* set a special with sign and type */ 396void mpd_setspecial(mpd_t *result, uint8_t sign, uint8_t type); 397/* set coefficient to zero or all nines */ 398void mpd_zerocoeff(mpd_t *result); 399void mpd_qmaxcoeff(mpd_t *result, const mpd_context_t *ctx, uint32_t *status); 400 401/* quietly assign a C integer type to an mpd_t */ 402void mpd_qset_ssize(mpd_t *result, mpd_ssize_t a, const mpd_context_t *ctx, uint32_t *status); 403void mpd_qset_i32(mpd_t *result, int32_t a, const mpd_context_t *ctx, uint32_t *status); 404void mpd_qset_uint(mpd_t *result, mpd_uint_t a, const mpd_context_t *ctx, uint32_t *status); 405void mpd_qset_u32(mpd_t *result, uint32_t a, const mpd_context_t *ctx, uint32_t *status); 406#ifndef MPD_LEGACY_COMPILER 407void mpd_qset_i64(mpd_t *result, int64_t a, const mpd_context_t *ctx, uint32_t *status); 408void mpd_qset_u64(mpd_t *result, uint64_t a, const mpd_context_t *ctx, uint32_t *status); 409void mpd_qset_i64_exact(mpd_t *result, int64_t a, uint32_t *status); 410void mpd_qset_u64_exact(mpd_t *result, uint64_t a, uint32_t *status); 411#endif 412 413/* quietly assign a C integer type to an mpd_t with a static coefficient */ 414void mpd_qsset_ssize(mpd_t *result, mpd_ssize_t a, const mpd_context_t *ctx, uint32_t *status); 415void mpd_qsset_i32(mpd_t *result, int32_t a, const mpd_context_t *ctx, uint32_t *status); 416void mpd_qsset_uint(mpd_t *result, mpd_uint_t a, const mpd_context_t *ctx, uint32_t *status); 417void mpd_qsset_u32(mpd_t *result, uint32_t a, const mpd_context_t *ctx, uint32_t *status); 418 419/* quietly get a C integer type from an mpd_t */ 420mpd_ssize_t mpd_qget_ssize(const mpd_t *dec, uint32_t *status); 421mpd_uint_t mpd_qget_uint(const mpd_t *dec, uint32_t *status); 422mpd_uint_t mpd_qabs_uint(const mpd_t *dec, uint32_t *status); 423 424int32_t mpd_qget_i32(const mpd_t *dec, uint32_t *status); 425uint32_t mpd_qget_u32(const mpd_t *dec, uint32_t *status); 426#ifndef MPD_LEGACY_COMPILER 427int64_t mpd_qget_i64(const mpd_t *dec, uint32_t *status); 428uint64_t mpd_qget_u64(const mpd_t *dec, uint32_t *status); 429#endif 430 431/* quiet functions */ 432int mpd_qcheck_nan(mpd_t *nanresult, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 433int mpd_qcheck_nans(mpd_t *nanresult, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 434void mpd_qfinalize(mpd_t *result, const mpd_context_t *ctx, uint32_t *status); 435 436const char *mpd_class(const mpd_t *a, const mpd_context_t *ctx); 437 438int mpd_qcopy(mpd_t *result, const mpd_t *a, uint32_t *status); 439int mpd_qcopy_cxx(mpd_t *result, const mpd_t *a); 440mpd_t *mpd_qncopy(const mpd_t *a); 441int mpd_qcopy_abs(mpd_t *result, const mpd_t *a, uint32_t *status); 442int mpd_qcopy_negate(mpd_t *result, const mpd_t *a, uint32_t *status); 443int mpd_qcopy_sign(mpd_t *result, const mpd_t *a, const mpd_t *b, uint32_t *status); 444 445void mpd_qand(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 446void mpd_qinvert(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 447void mpd_qlogb(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 448void mpd_qor(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 449void mpd_qscaleb(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 450void mpd_qxor(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 451int mpd_same_quantum(const mpd_t *a, const mpd_t *b); 452 453void mpd_qrotate(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 454int mpd_qshiftl(mpd_t *result, const mpd_t *a, mpd_ssize_t n, uint32_t *status); 455mpd_uint_t mpd_qshiftr(mpd_t *result, const mpd_t *a, mpd_ssize_t n, uint32_t *status); 456mpd_uint_t mpd_qshiftr_inplace(mpd_t *result, mpd_ssize_t n); 457void mpd_qshift(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 458void mpd_qshiftn(mpd_t *result, const mpd_t *a, mpd_ssize_t n, const mpd_context_t *ctx, uint32_t *status); 459 460int mpd_qcmp(const mpd_t *a, const mpd_t *b, uint32_t *status); 461int mpd_qcompare(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 462int mpd_qcompare_signal(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 463int mpd_cmp_total(const mpd_t *a, const mpd_t *b); 464int mpd_cmp_total_mag(const mpd_t *a, const mpd_t *b); 465int mpd_compare_total(mpd_t *result, const mpd_t *a, const mpd_t *b); 466int mpd_compare_total_mag(mpd_t *result, const mpd_t *a, const mpd_t *b); 467 468void mpd_qround_to_intx(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 469void mpd_qround_to_int(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 470void mpd_qtrunc(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 471void mpd_qfloor(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 472void mpd_qceil(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 473 474void mpd_qabs(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 475void mpd_qmax(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 476void mpd_qmax_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 477void mpd_qmin(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 478void mpd_qmin_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 479void mpd_qminus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 480void mpd_qplus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 481void mpd_qnext_minus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 482void mpd_qnext_plus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 483void mpd_qnext_toward(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 484void mpd_qquantize(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 485void mpd_qrescale(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, const mpd_context_t *ctx, uint32_t *status); 486void mpd_qrescale_fmt(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, const mpd_context_t *ctx, uint32_t *status); 487void mpd_qreduce(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 488void mpd_qadd(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 489void mpd_qadd_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status); 490void mpd_qadd_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status); 491void mpd_qadd_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status); 492void mpd_qadd_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status); 493void mpd_qsub(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 494void mpd_qsub_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status); 495void mpd_qsub_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status); 496void mpd_qsub_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status); 497void mpd_qsub_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status); 498void mpd_qmul(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 499void mpd_qmul_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status); 500void mpd_qmul_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status); 501void mpd_qmul_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status); 502void mpd_qmul_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status); 503void mpd_qfma(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_t *c, const mpd_context_t *ctx, uint32_t *status); 504void mpd_qdiv(mpd_t *q, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 505void mpd_qdiv_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status); 506void mpd_qdiv_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status); 507void mpd_qdiv_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status); 508void mpd_qdiv_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status); 509void mpd_qdivint(mpd_t *q, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 510void mpd_qrem(mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 511void mpd_qrem_near(mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 512void mpd_qdivmod(mpd_t *q, mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status); 513void mpd_qpow(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_context_t *ctx, uint32_t *status); 514void mpd_qpowmod(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_t *mod, const mpd_context_t *ctx, uint32_t *status); 515void mpd_qexp(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 516void mpd_qln10(mpd_t *result, mpd_ssize_t prec, uint32_t *status); 517void mpd_qln(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 518void mpd_qlog10(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 519void mpd_qsqrt(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 520void mpd_qinvroot(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status); 521#ifndef MPD_LEGACY_COMPILER 522void mpd_qadd_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status); 523void mpd_qadd_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status); 524void mpd_qsub_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status); 525void mpd_qsub_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status); 526void mpd_qmul_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status); 527void mpd_qmul_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status); 528void mpd_qdiv_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status); 529void mpd_qdiv_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status); 530#endif 531 532 533size_t mpd_sizeinbase(const mpd_t *a, uint32_t base); 534void mpd_qimport_u16(mpd_t *result, const uint16_t *srcdata, size_t srclen, 535 uint8_t srcsign, uint32_t srcbase, 536 const mpd_context_t *ctx, uint32_t *status); 537void mpd_qimport_u32(mpd_t *result, const uint32_t *srcdata, size_t srclen, 538 uint8_t srcsign, uint32_t srcbase, 539 const mpd_context_t *ctx, uint32_t *status); 540size_t mpd_qexport_u16(uint16_t **rdata, size_t rlen, uint32_t base, 541 const mpd_t *src, uint32_t *status); 542size_t mpd_qexport_u32(uint32_t **rdata, size_t rlen, uint32_t base, 543 const mpd_t *src, uint32_t *status); 544 545 546/******************************************************************************/ 547/* Signalling functions */ 548/******************************************************************************/ 549 550char *mpd_format(const mpd_t *dec, const char *fmt, mpd_context_t *ctx); 551void mpd_import_u16(mpd_t *result, const uint16_t *srcdata, size_t srclen, uint8_t srcsign, uint32_t base, mpd_context_t *ctx); 552void mpd_import_u32(mpd_t *result, const uint32_t *srcdata, size_t srclen, uint8_t srcsign, uint32_t base, mpd_context_t *ctx); 553size_t mpd_export_u16(uint16_t **rdata, size_t rlen, uint32_t base, const mpd_t *src, mpd_context_t *ctx); 554size_t mpd_export_u32(uint32_t **rdata, size_t rlen, uint32_t base, const mpd_t *src, mpd_context_t *ctx); 555void mpd_finalize(mpd_t *result, mpd_context_t *ctx); 556int mpd_check_nan(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 557int mpd_check_nans(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 558void mpd_set_string(mpd_t *result, const char *s, mpd_context_t *ctx); 559void mpd_maxcoeff(mpd_t *result, mpd_context_t *ctx); 560void mpd_sset_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx); 561void mpd_sset_i32(mpd_t *result, int32_t a, mpd_context_t *ctx); 562void mpd_sset_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx); 563void mpd_sset_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx); 564void mpd_set_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx); 565void mpd_set_i32(mpd_t *result, int32_t a, mpd_context_t *ctx); 566void mpd_set_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx); 567void mpd_set_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx); 568#ifndef MPD_LEGACY_COMPILER 569void mpd_set_i64(mpd_t *result, int64_t a, mpd_context_t *ctx); 570void mpd_set_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx); 571#endif 572mpd_ssize_t mpd_get_ssize(const mpd_t *a, mpd_context_t *ctx); 573mpd_uint_t mpd_get_uint(const mpd_t *a, mpd_context_t *ctx); 574mpd_uint_t mpd_abs_uint(const mpd_t *a, mpd_context_t *ctx); 575int32_t mpd_get_i32(const mpd_t *a, mpd_context_t *ctx); 576uint32_t mpd_get_u32(const mpd_t *a, mpd_context_t *ctx); 577#ifndef MPD_LEGACY_COMPILER 578int64_t mpd_get_i64(const mpd_t *a, mpd_context_t *ctx); 579uint64_t mpd_get_u64(const mpd_t *a, mpd_context_t *ctx); 580#endif 581void mpd_and(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 582void mpd_copy(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 583void mpd_canonical(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 584void mpd_copy_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 585void mpd_copy_negate(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 586void mpd_copy_sign(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 587void mpd_invert(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 588void mpd_logb(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 589void mpd_or(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 590void mpd_rotate(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 591void mpd_scaleb(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 592void mpd_shiftl(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx); 593mpd_uint_t mpd_shiftr(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx); 594void mpd_shiftn(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx); 595void mpd_shift(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 596void mpd_xor(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 597void mpd_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 598int mpd_cmp(const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 599int mpd_compare(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 600int mpd_compare_signal(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 601void mpd_add(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 602void mpd_add_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx); 603void mpd_add_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx); 604void mpd_add_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx); 605void mpd_add_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx); 606void mpd_sub(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 607void mpd_sub_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx); 608void mpd_sub_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx); 609void mpd_sub_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx); 610void mpd_sub_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx); 611void mpd_div(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 612void mpd_div_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx); 613void mpd_div_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx); 614void mpd_div_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx); 615void mpd_div_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx); 616void mpd_divmod(mpd_t *q, mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 617void mpd_divint(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 618void mpd_exp(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 619void mpd_fma(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_t *c, mpd_context_t *ctx); 620void mpd_ln(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 621void mpd_log10(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 622void mpd_max(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 623void mpd_max_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 624void mpd_min(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 625void mpd_min_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 626void mpd_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 627void mpd_mul(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 628void mpd_mul_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx); 629void mpd_mul_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx); 630void mpd_mul_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx); 631void mpd_mul_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx); 632void mpd_next_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 633void mpd_next_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 634void mpd_next_toward(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 635void mpd_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 636void mpd_pow(mpd_t *result, const mpd_t *base, const mpd_t *exp, mpd_context_t *ctx); 637void mpd_powmod(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_t *mod, mpd_context_t *ctx); 638void mpd_quantize(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 639void mpd_rescale(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, mpd_context_t *ctx); 640void mpd_reduce(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 641void mpd_rem(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 642void mpd_rem_near(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx); 643void mpd_round_to_intx(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 644void mpd_round_to_int(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 645void mpd_trunc(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 646void mpd_floor(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 647void mpd_ceil(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 648void mpd_sqrt(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 649void mpd_invroot(mpd_t *result, const mpd_t *a, mpd_context_t *ctx); 650#ifndef MPD_LEGACY_COMPILER 651void mpd_add_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx); 652void mpd_add_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx); 653void mpd_sub_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx); 654void mpd_sub_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx); 655void mpd_div_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx); 656void mpd_div_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx); 657void mpd_mul_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx); 658void mpd_mul_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx); 659#endif 660 661 662/******************************************************************************/ 663/* Configuration specific */ 664/******************************************************************************/ 665 666#ifdef MPD_CONFIG_64 667void mpd_qsset_i64(mpd_t *result, int64_t a, const mpd_context_t *ctx, uint32_t *status); 668void mpd_qsset_u64(mpd_t *result, uint64_t a, const mpd_context_t *ctx, uint32_t *status); 669void mpd_sset_i64(mpd_t *result, int64_t a, mpd_context_t *ctx); 670void mpd_sset_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx); 671#endif 672 673 674/******************************************************************************/ 675/* Get attributes of a decimal */ 676/******************************************************************************/ 677 678mpd_ssize_t mpd_adjexp(const mpd_t *dec); 679mpd_ssize_t mpd_etiny(const mpd_context_t *ctx); 680mpd_ssize_t mpd_etop(const mpd_context_t *ctx); 681mpd_uint_t mpd_msword(const mpd_t *dec); 682int mpd_word_digits(mpd_uint_t word); 683/* most significant digit of a word */ 684mpd_uint_t mpd_msd(mpd_uint_t word); 685/* least significant digit of a word */ 686mpd_uint_t mpd_lsd(mpd_uint_t word); 687/* coefficient size needed to store 'digits' */ 688mpd_ssize_t mpd_digits_to_size(mpd_ssize_t digits); 689/* number of digits in the exponent, undefined for MPD_SSIZE_MIN */ 690int mpd_exp_digits(mpd_ssize_t exp); 691int mpd_iscanonical(const mpd_t *dec); 692int mpd_isfinite(const mpd_t *dec); 693int mpd_isinfinite(const mpd_t *dec); 694int mpd_isinteger(const mpd_t *dec); 695int mpd_isnan(const mpd_t *dec); 696int mpd_isnegative(const mpd_t *dec); 697int mpd_ispositive(const mpd_t *dec); 698int mpd_isqnan(const mpd_t *dec); 699int mpd_issigned(const mpd_t *dec); 700int mpd_issnan(const mpd_t *dec); 701int mpd_isspecial(const mpd_t *dec); 702int mpd_iszero(const mpd_t *dec); 703/* undefined for special numbers */ 704int mpd_iszerocoeff(const mpd_t *dec); 705int mpd_isnormal(const mpd_t *dec, const mpd_context_t *ctx); 706int mpd_issubnormal(const mpd_t *dec, const mpd_context_t *ctx); 707/* odd word */ 708int mpd_isoddword(mpd_uint_t word); 709/* odd coefficient */ 710int mpd_isoddcoeff(const mpd_t *dec); 711/* odd decimal, only defined for integers */ 712int mpd_isodd(const mpd_t *dec); 713/* even decimal, only defined for integers */ 714int mpd_iseven(const mpd_t *dec); 715/* 0 if dec is positive, 1 if dec is negative */ 716uint8_t mpd_sign(const mpd_t *dec); 717/* 1 if dec is positive, -1 if dec is negative */ 718int mpd_arith_sign(const mpd_t *dec); 719long mpd_radix(void); 720int mpd_isdynamic(const mpd_t *dec); 721int mpd_isstatic(const mpd_t *dec); 722int mpd_isdynamic_data(const mpd_t *dec); 723int mpd_isstatic_data(const mpd_t *dec); 724int mpd_isshared_data(const mpd_t *dec); 725int mpd_isconst_data(const mpd_t *dec); 726mpd_ssize_t mpd_trail_zeros(const mpd_t *dec); 727 728 729/******************************************************************************/ 730/* Set attributes of a decimal */ 731/******************************************************************************/ 732 733/* set number of decimal digits in the coefficient */ 734void mpd_setdigits(mpd_t *result); 735void mpd_set_sign(mpd_t *result, uint8_t sign); 736/* copy sign from another decimal */ 737void mpd_signcpy(mpd_t *result, const mpd_t *a); 738void mpd_set_infinity(mpd_t *result); 739void mpd_set_qnan(mpd_t *result); 740void mpd_set_snan(mpd_t *result); 741void mpd_set_negative(mpd_t *result); 742void mpd_set_positive(mpd_t *result); 743void mpd_set_dynamic(mpd_t *result); 744void mpd_set_static(mpd_t *result); 745void mpd_set_dynamic_data(mpd_t *result); 746void mpd_set_static_data(mpd_t *result); 747void mpd_set_shared_data(mpd_t *result); 748void mpd_set_const_data(mpd_t *result); 749void mpd_clear_flags(mpd_t *result); 750void mpd_set_flags(mpd_t *result, uint8_t flags); 751void mpd_copy_flags(mpd_t *result, const mpd_t *a); 752 753 754/******************************************************************************/ 755/* Error Macros */ 756/******************************************************************************/ 757 758#define mpd_err_fatal(...) \ 759 do {fprintf(stderr, "%s:%d: error: ", __FILE__, __LINE__); \ 760 fprintf(stderr, __VA_ARGS__); fputc('\n', stderr); \ 761 abort(); \ 762 } while (0) 763#define mpd_err_warn(...) \ 764 do {fprintf(stderr, "%s:%d: warning: ", __FILE__, __LINE__); \ 765 fprintf(stderr, __VA_ARGS__); fputc('\n', stderr); \ 766 } while (0) 767 768 769/******************************************************************************/ 770/* Memory handling */ 771/******************************************************************************/ 772 773extern void *(* mpd_mallocfunc)(size_t size); 774extern void *(* mpd_callocfunc)(size_t nmemb, size_t size); 775extern void *(* mpd_reallocfunc)(void *ptr, size_t size); 776extern void (* mpd_free)(void *ptr); 777 778void *mpd_callocfunc_em(size_t nmemb, size_t size); 779 780void *mpd_alloc(mpd_size_t nmemb, mpd_size_t size); 781void *mpd_calloc(mpd_size_t nmemb, mpd_size_t size); 782void *mpd_realloc(void *ptr, mpd_size_t nmemb, mpd_size_t size, uint8_t *err); 783void *mpd_sh_alloc(mpd_size_t struct_size, mpd_size_t nmemb, mpd_size_t size); 784 785mpd_t *mpd_qnew(void); 786mpd_t *mpd_new(mpd_context_t *ctx); 787mpd_t *mpd_qnew_size(mpd_ssize_t nwords); 788void mpd_del(mpd_t *dec); 789 790void mpd_uint_zero(mpd_uint_t *dest, mpd_size_t len); 791int mpd_qresize(mpd_t *result, mpd_ssize_t nwords, uint32_t *status); 792int mpd_qresize_zero(mpd_t *result, mpd_ssize_t nwords, uint32_t *status); 793void mpd_minalloc(mpd_t *result); 794 795int mpd_resize(mpd_t *result, mpd_ssize_t nwords, mpd_context_t *ctx); 796int mpd_resize_zero(mpd_t *result, mpd_ssize_t nwords, mpd_context_t *ctx); 797 798 799#ifdef __cplusplus 800} /* END extern "C" */ 801#endif 802 803 804#endif /* LIBMPDEC_MPDECIMAL_H_ */ 805