CoCalc -- mpdecimal.h

GitHub Repository: allendowney/cpython
Path: blob/main/Modules/_decimal/libmpdec/mpdecimal.h
¹² views
1
/*
2
 * Copyright (c) 2008-2020 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
 *
11
 * 2. Redistributions in binary form must reproduce the above copyright
12
 *    notice, this list of conditions and the following disclaimer in the
13
 *    documentation and/or other materials provided with the distribution.
14
 *
15
 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS" AND
16
 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17
 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18
 * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
19
 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
20
 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
21
 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
22
 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
23
 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
24
 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
25
 * SUCH DAMAGE.
26
 */
27

28

29
#ifndef LIBMPDEC_MPDECIMAL_H_
30
#define LIBMPDEC_MPDECIMAL_H_
31

32

33
#ifndef _MSC_VER
34
  #include "pyconfig.h"
35
#endif
36

37
#ifdef __cplusplus
38
  #include <cinttypes>
39
  #include <climits>
40
  #include <cstdint>
41
  #include <cstdio>
42
  #include <cstdlib>
43
  #define MPD_UINT8_C(x) (static_cast<uint8_t>(x))
44
extern "C" {
45
#else
46
  #include <inttypes.h>
47
  #include <limits.h>
48
  #include <stdint.h>
49
  #include <stdio.h>
50
  #include <stdlib.h>
51
  #define MPD_UINT8_C(x) ((uint8_t)x)
52
#endif
53

54

55
#if (defined(__linux__) || defined(__FreeBSD__) || defined(__APPLE__)) && \
56
    defined(__GNUC__) && __GNUC__ >= 4 && !defined(__INTEL_COMPILER)
57
  #define MPD_PRAGMA(x) _Pragma(x)
58
  #define MPD_HIDE_SYMBOLS_START "GCC visibility push(hidden)"
59
  #define MPD_HIDE_SYMBOLS_END "GCC visibility pop"
60
#else
61
  #define MPD_PRAGMA(x)
62
  #define MPD_HIDE_SYMBOLS_START
63
  #define MPD_HIDE_SYMBOLS_END
64
#endif
65

66
#if defined(_MSC_VER)
67
  #define EXTINLINE extern inline
68
#else
69
  #define EXTINLINE
70
#endif
71

72

73
/* This header file is internal for the purpose of building _decimal.so.
74
 * All symbols should have local scope in the DSO. */
75
MPD_PRAGMA(MPD_HIDE_SYMBOLS_START)
76

77

78
/******************************************************************************/
79
/*                                  Version                                   */
80
/******************************************************************************/
81

82
#define MPD_MAJOR_VERSION 2
83
#define MPD_MINOR_VERSION 5
84
#define MPD_MICRO_VERSION 1
85

86
#define MPD_VERSION "2.5.1"
87

88
#define MPD_VERSION_HEX ((MPD_MAJOR_VERSION << 24) | \
89
                         (MPD_MINOR_VERSION << 16) | \
90
                         (MPD_MICRO_VERSION <<  8))
91

92
const char *mpd_version(void);
93

94

95
/******************************************************************************/
96
/*                              Configuration                                 */
97
/******************************************************************************/
98

99
#if defined(UNIVERSAL)
100
  #if defined(CONFIG_64) || defined(CONFIG_32)
101
    #error "cannot use CONFIG_64 or CONFIG_32 with UNIVERSAL."
102
  #endif
103
  #if defined(__ppc__)
104
    #define CONFIG_32
105
    #define ANSI
106
  #elif defined(__ppc64__)
107
    #define CONFIG_64
108
    #define ANSI
109
  #elif defined(__i386__)
110
    #define CONFIG_32
111
    #define ANSI
112
  #elif defined(__x86_64__)
113
    #define CONFIG_64
114
    #define ASM
115
  #elif defined(__arm64__)
116
    #define CONFIG_64
117
    #define ANSI
118
  #else
119
    #error "unknown architecture for universal build."
120
  #endif
121
#endif
122

123

124
/* BEGIN CONFIG_64 */
125
#if defined(CONFIG_64)
126
/* types for modular and base arithmetic */
127
#define MPD_UINT_MAX UINT64_MAX
128
#define MPD_BITS_PER_UINT 64
129
typedef uint64_t mpd_uint_t;  /* unsigned mod type */
130

131
#define MPD_SIZE_MAX SIZE_MAX
132
typedef size_t mpd_size_t; /* unsigned size type */
133

134
/* type for exp, digits, len, prec */
135
#define MPD_SSIZE_MAX INT64_MAX
136
#define MPD_SSIZE_MIN INT64_MIN
137
typedef int64_t mpd_ssize_t;
138
#define _mpd_strtossize strtoll
139

140
/* decimal arithmetic */
141
#define MPD_RADIX 10000000000000000000ULL  /* 10**19 */
142
#define MPD_RDIGITS 19
143
#define MPD_MAX_POW10 19
144
#define MPD_EXPDIGITS 19  /* MPD_EXPDIGITS <= MPD_RDIGITS+1 */
145

146
#define MPD_MAXTRANSFORM_2N 4294967296ULL      /* 2**32 */
147
#define MPD_MAX_PREC 999999999999999999LL
148
#define MPD_MAX_PREC_LOG2 64
149
#define MPD_ELIMIT  1000000000000000000LL
150
#define MPD_MAX_EMAX   999999999999999999LL    /* ELIMIT-1 */
151
#define MPD_MIN_EMIN  (-999999999999999999LL)  /* -EMAX */
152
#define MPD_MIN_ETINY (MPD_MIN_EMIN-(MPD_MAX_PREC-1))
153
#define MPD_EXP_INF 2000000000000000001LL
154
#define MPD_EXP_CLAMP (-4000000000000000001LL)
155
#define MPD_MAXIMPORT 105263157894736842L /* ceil((2*MPD_MAX_PREC)/MPD_RDIGITS) */
156
#define MPD_IEEE_CONTEXT_MAX_BITS 512     /* 16*(log2(MPD_MAX_EMAX / 3)-3) */
157

158
/* conversion specifiers */
159
#define PRI_mpd_uint_t PRIu64
160
#define PRI_mpd_ssize_t PRIi64
161
/* END CONFIG_64 */
162

163

164
/* BEGIN CONFIG_32 */
165
#elif defined(CONFIG_32)
166
/* types for modular and base arithmetic */
167
#define MPD_UINT_MAX UINT32_MAX
168
#define MPD_BITS_PER_UINT 32
169
typedef uint32_t mpd_uint_t;  /* unsigned mod type */
170

171
#ifndef LEGACY_COMPILER
172
#define MPD_UUINT_MAX UINT64_MAX
173
typedef uint64_t mpd_uuint_t; /* double width unsigned mod type */
174
#endif
175

176
#define MPD_SIZE_MAX SIZE_MAX
177
typedef size_t mpd_size_t; /* unsigned size type */
178

179
/* type for dec->len, dec->exp, ctx->prec */
180
#define MPD_SSIZE_MAX INT32_MAX
181
#define MPD_SSIZE_MIN INT32_MIN
182
typedef int32_t mpd_ssize_t;
183
#define _mpd_strtossize strtol
184

185
/* decimal arithmetic */
186
#define MPD_RADIX 1000000000UL  /* 10**9 */
187
#define MPD_RDIGITS 9
188
#define MPD_MAX_POW10 9
189
#define MPD_EXPDIGITS 10 /* MPD_EXPDIGITS <= MPD_RDIGITS+1 */
190

191
#define MPD_MAXTRANSFORM_2N 33554432UL /* 2**25 */
192
#define MPD_MAX_PREC 425000000L
193
#define MPD_MAX_PREC_LOG2 32
194
#define MPD_ELIMIT 425000001L
195
#define MPD_MAX_EMAX 425000000L        /* ELIMIT-1 */
196
#define MPD_MIN_EMIN (-425000000L)     /* -EMAX */
197
#define MPD_MIN_ETINY (MPD_MIN_EMIN-(MPD_MAX_PREC-1))
198
#define MPD_EXP_INF 1000000001L       /* allows for emax=999999999 in the tests */
199
#define MPD_EXP_CLAMP (-2000000001L)  /* allows for emin=-999999999 in the tests */
200
#define MPD_MAXIMPORT 94444445L       /* ceil((2*MPD_MAX_PREC)/MPD_RDIGITS) */
201
#define MPD_IEEE_CONTEXT_MAX_BITS 256 /* 16*(log2(MPD_MAX_EMAX / 3)-3) */
202

203
/* conversion specifiers */
204
#define PRI_mpd_uint_t PRIu32
205
#define PRI_mpd_ssize_t PRIi32
206
/* END CONFIG_32 */
207

208
#else
209
  #error "define CONFIG_64 or CONFIG_32"
210
#endif
211
/* END CONFIG */
212

213

214
#if MPD_SIZE_MAX != MPD_UINT_MAX
215
  #error "unsupported platform: need mpd_size_t == mpd_uint_t"
216
#endif
217

218

219
/******************************************************************************/
220
/*                                Context                                     */
221
/******************************************************************************/
222

223
enum {
224
    MPD_ROUND_UP,          /* round away from 0               */
225
    MPD_ROUND_DOWN,        /* round toward 0 (truncate)       */
226
    MPD_ROUND_CEILING,     /* round toward +infinity          */
227
    MPD_ROUND_FLOOR,       /* round toward -infinity          */
228
    MPD_ROUND_HALF_UP,     /* 0.5 is rounded up               */
229
    MPD_ROUND_HALF_DOWN,   /* 0.5 is rounded down             */
230
    MPD_ROUND_HALF_EVEN,   /* 0.5 is rounded to even          */
231
    MPD_ROUND_05UP,        /* round zero or five away from 0  */
232
    MPD_ROUND_TRUNC,       /* truncate, but set infinity      */
233
    MPD_ROUND_GUARD
234
};
235

236
enum { MPD_CLAMP_DEFAULT, MPD_CLAMP_IEEE_754, MPD_CLAMP_GUARD };
237

238
extern const char * const mpd_round_string[MPD_ROUND_GUARD];
239
extern const char * const mpd_clamp_string[MPD_CLAMP_GUARD];
240

241

242
typedef struct mpd_context_t {
243
    mpd_ssize_t prec;   /* precision */
244
    mpd_ssize_t emax;   /* max positive exp */
245
    mpd_ssize_t emin;   /* min negative exp */
246
    uint32_t traps;     /* status events that should be trapped */
247
    uint32_t status;    /* status flags */
248
    uint32_t newtrap;   /* set by mpd_addstatus_raise() */
249
    int      round;     /* rounding mode */
250
    int      clamp;     /* clamp mode */
251
    int      allcr;     /* all functions correctly rounded */
252
} mpd_context_t;
253

254

255
/* Status flags */
256
#define MPD_Clamped             0x00000001U
257
#define MPD_Conversion_syntax   0x00000002U
258
#define MPD_Division_by_zero    0x00000004U
259
#define MPD_Division_impossible 0x00000008U
260
#define MPD_Division_undefined  0x00000010U
261
#define MPD_Fpu_error           0x00000020U
262
#define MPD_Inexact             0x00000040U
263
#define MPD_Invalid_context     0x00000080U
264
#define MPD_Invalid_operation   0x00000100U
265
#define MPD_Malloc_error        0x00000200U
266
#define MPD_Not_implemented     0x00000400U
267
#define MPD_Overflow            0x00000800U
268
#define MPD_Rounded             0x00001000U
269
#define MPD_Subnormal           0x00002000U
270
#define MPD_Underflow           0x00004000U
271
#define MPD_Max_status         (0x00008000U-1U)
272

273
/* Conditions that result in an IEEE 754 exception */
274
#define MPD_IEEE_Invalid_operation (MPD_Conversion_syntax |   \
275
                                    MPD_Division_impossible | \
276
                                    MPD_Division_undefined |  \
277
                                    MPD_Fpu_error |           \
278
                                    MPD_Invalid_context |     \
279
                                    MPD_Invalid_operation |   \
280
                                    MPD_Malloc_error)         \
281

282
/* Errors that require the result of an operation to be set to NaN */
283
#define MPD_Errors (MPD_IEEE_Invalid_operation | \
284
                    MPD_Division_by_zero)
285

286
/* Default traps */
287
#define MPD_Traps (MPD_IEEE_Invalid_operation | \
288
                   MPD_Division_by_zero |       \
289
                   MPD_Overflow |               \
290
                   MPD_Underflow)
291

292
/* Official name */
293
#define MPD_Insufficient_storage MPD_Malloc_error
294

295
/* IEEE 754 interchange format contexts */
296
#define MPD_DECIMAL32 32
297
#define MPD_DECIMAL64 64
298
#define MPD_DECIMAL128 128
299

300

301
#define MPD_MINALLOC_MIN 2
302
#define MPD_MINALLOC_MAX 64
303
extern mpd_ssize_t MPD_MINALLOC;
304
extern void (* mpd_traphandler)(mpd_context_t *);
305
void mpd_dflt_traphandler(mpd_context_t *);
306

307
void mpd_setminalloc(mpd_ssize_t n);
308
void mpd_init(mpd_context_t *ctx, mpd_ssize_t prec);
309

310
void mpd_maxcontext(mpd_context_t *ctx);
311
void mpd_defaultcontext(mpd_context_t *ctx);
312
void mpd_basiccontext(mpd_context_t *ctx);
313
int mpd_ieee_context(mpd_context_t *ctx, int bits);
314

315
mpd_ssize_t mpd_getprec(const mpd_context_t *ctx);
316
mpd_ssize_t mpd_getemax(const mpd_context_t *ctx);
317
mpd_ssize_t mpd_getemin(const mpd_context_t *ctx);
318
int mpd_getround(const mpd_context_t *ctx);
319
uint32_t mpd_gettraps(const mpd_context_t *ctx);
320
uint32_t mpd_getstatus(const mpd_context_t *ctx);
321
int mpd_getclamp(const mpd_context_t *ctx);
322
int mpd_getcr(const mpd_context_t *ctx);
323

324
int mpd_qsetprec(mpd_context_t *ctx, mpd_ssize_t prec);
325
int mpd_qsetemax(mpd_context_t *ctx, mpd_ssize_t emax);
326
int mpd_qsetemin(mpd_context_t *ctx, mpd_ssize_t emin);
327
int mpd_qsetround(mpd_context_t *ctx, int newround);
328
int mpd_qsettraps(mpd_context_t *ctx, uint32_t flags);
329
int mpd_qsetstatus(mpd_context_t *ctx, uint32_t flags);
330
int mpd_qsetclamp(mpd_context_t *ctx, int c);
331
int mpd_qsetcr(mpd_context_t *ctx, int c);
332
void mpd_addstatus_raise(mpd_context_t *ctx, uint32_t flags);
333

334

335
/******************************************************************************/
336
/*                           Decimal Arithmetic                               */
337
/******************************************************************************/
338

339
/* mpd_t flags */
340
#define MPD_POS                 MPD_UINT8_C(0)
341
#define MPD_NEG                 MPD_UINT8_C(1)
342
#define MPD_INF                 MPD_UINT8_C(2)
343
#define MPD_NAN                 MPD_UINT8_C(4)
344
#define MPD_SNAN                MPD_UINT8_C(8)
345
#define MPD_SPECIAL (MPD_INF|MPD_NAN|MPD_SNAN)
346
#define MPD_STATIC              MPD_UINT8_C(16)
347
#define MPD_STATIC_DATA         MPD_UINT8_C(32)
348
#define MPD_SHARED_DATA         MPD_UINT8_C(64)
349
#define MPD_CONST_DATA          MPD_UINT8_C(128)
350
#define MPD_DATAFLAGS (MPD_STATIC_DATA|MPD_SHARED_DATA|MPD_CONST_DATA)
351

352
/* mpd_t */
353
typedef struct mpd_t {
354
    uint8_t flags;
355
    mpd_ssize_t exp;
356
    mpd_ssize_t digits;
357
    mpd_ssize_t len;
358
    mpd_ssize_t alloc;
359
    mpd_uint_t *data;
360
} mpd_t;
361

362

363
/******************************************************************************/
364
/*                                    Triple                                  */
365
/******************************************************************************/
366

367
/* status cases for getting a triple */
368
enum mpd_triple_class {
369
  MPD_TRIPLE_NORMAL,
370
  MPD_TRIPLE_INF,
371
  MPD_TRIPLE_QNAN,
372
  MPD_TRIPLE_SNAN,
373
  MPD_TRIPLE_ERROR,
374
};
375

376
typedef struct {
377
  enum mpd_triple_class tag;
378
  uint8_t sign;
379
  uint64_t hi;
380
  uint64_t lo;
381
  int64_t exp;
382
} mpd_uint128_triple_t;
383

384
int mpd_from_uint128_triple(mpd_t *result, const mpd_uint128_triple_t *triple, uint32_t *status);
385
mpd_uint128_triple_t mpd_as_uint128_triple(const mpd_t *a);
386

387

388
/******************************************************************************/
389
/*                       Quiet, thread-safe functions                         */
390
/******************************************************************************/
391

392
/* format specification */
393
typedef struct mpd_spec_t {
394
    mpd_ssize_t min_width; /* minimum field width */
395
    mpd_ssize_t prec;      /* fraction digits or significant digits */
396
    char type;             /* conversion specifier */
397
    char align;            /* alignment */
398
    char sign;             /* sign printing/alignment */
399
    char fill[5];          /* fill character */
400
    const char *dot;       /* decimal point */
401
    const char *sep;       /* thousands separator */
402
    const char *grouping;  /* grouping of digits */
403
} mpd_spec_t;
404

405
/* output to a string */
406
char *mpd_to_sci(const mpd_t *dec, int fmt);
407
char *mpd_to_eng(const mpd_t *dec, int fmt);
408
mpd_ssize_t mpd_to_sci_size(char **res, const mpd_t *dec, int fmt);
409
mpd_ssize_t mpd_to_eng_size(char **res, const mpd_t *dec, int fmt);
410
int mpd_validate_lconv(mpd_spec_t *spec);
411
int mpd_parse_fmt_str(mpd_spec_t *spec, const char *fmt, int caps);
412
char *mpd_qformat_spec(const mpd_t *dec, const mpd_spec_t *spec, const mpd_context_t *ctx, uint32_t *status);
413
char *mpd_qformat(const mpd_t *dec, const char *fmt, const mpd_context_t *ctx, uint32_t *status);
414

415
#define MPD_NUM_FLAGS 15
416
#define MPD_MAX_FLAG_STRING 208
417
#define MPD_MAX_FLAG_LIST (MPD_MAX_FLAG_STRING+18)
418
#define MPD_MAX_SIGNAL_LIST 121
419
int mpd_snprint_flags(char *dest, int nmemb, uint32_t flags);
420
int mpd_lsnprint_flags(char *dest, int nmemb, uint32_t flags, const char *flag_string[]);
421
int mpd_lsnprint_signals(char *dest, int nmemb, uint32_t flags, const char *signal_string[]);
422

423
/* output to a file */
424
void mpd_fprint(FILE *file, const mpd_t *dec);
425
void mpd_print(const mpd_t *dec);
426

427
/* assignment from a string */
428
void mpd_qset_string(mpd_t *dec, const char *s, const mpd_context_t *ctx, uint32_t *status);
429
void mpd_qset_string_exact(mpd_t *dec, const char *s, uint32_t *status);
430

431
/* set to NaN with error flags */
432
void mpd_seterror(mpd_t *result, uint32_t flags, uint32_t *status);
433
/* set a special with sign and type */
434
void mpd_setspecial(mpd_t *result, uint8_t sign, uint8_t type);
435
/* set coefficient to zero or all nines */
436
void mpd_zerocoeff(mpd_t *result);
437
void mpd_qmaxcoeff(mpd_t *result, const mpd_context_t *ctx, uint32_t *status);
438

439
/* quietly assign a C integer type to an mpd_t */
440
void mpd_qset_ssize(mpd_t *result, mpd_ssize_t a, const mpd_context_t *ctx, uint32_t *status);
441
void mpd_qset_i32(mpd_t *result, int32_t a, const mpd_context_t *ctx, uint32_t *status);
442
void mpd_qset_uint(mpd_t *result, mpd_uint_t a, const mpd_context_t *ctx, uint32_t *status);
443
void mpd_qset_u32(mpd_t *result, uint32_t a, const mpd_context_t *ctx, uint32_t *status);
444
#ifndef LEGACY_COMPILER
445
void mpd_qset_i64(mpd_t *result, int64_t a, const mpd_context_t *ctx, uint32_t *status);
446
void mpd_qset_u64(mpd_t *result, uint64_t a, const mpd_context_t *ctx, uint32_t *status);
447
void mpd_qset_i64_exact(mpd_t *result, int64_t a, uint32_t *status);
448
void mpd_qset_u64_exact(mpd_t *result, uint64_t a, uint32_t *status);
449
#endif
450

451
/* quietly assign a C integer type to an mpd_t with a static coefficient */
452
void mpd_qsset_ssize(mpd_t *result, mpd_ssize_t a, const mpd_context_t *ctx, uint32_t *status);
453
void mpd_qsset_i32(mpd_t *result, int32_t a, const mpd_context_t *ctx, uint32_t *status);
454
void mpd_qsset_uint(mpd_t *result, mpd_uint_t a, const mpd_context_t *ctx, uint32_t *status);
455
void mpd_qsset_u32(mpd_t *result, uint32_t a, const mpd_context_t *ctx, uint32_t *status);
456

457
/* quietly get a C integer type from an mpd_t */
458
mpd_ssize_t mpd_qget_ssize(const mpd_t *dec, uint32_t *status);
459
mpd_uint_t mpd_qget_uint(const mpd_t *dec, uint32_t *status);
460
mpd_uint_t mpd_qabs_uint(const mpd_t *dec, uint32_t *status);
461

462
int32_t mpd_qget_i32(const mpd_t *dec, uint32_t *status);
463
uint32_t mpd_qget_u32(const mpd_t *dec, uint32_t *status);
464
#ifndef LEGACY_COMPILER
465
int64_t mpd_qget_i64(const mpd_t *dec, uint32_t *status);
466
uint64_t mpd_qget_u64(const mpd_t *dec, uint32_t *status);
467
#endif
468

469
/* quiet functions */
470
int mpd_qcheck_nan(mpd_t *nanresult, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
471
int mpd_qcheck_nans(mpd_t *nanresult, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
472
void mpd_qfinalize(mpd_t *result, const mpd_context_t *ctx, uint32_t *status);
473

474
const char *mpd_class(const mpd_t *a, const mpd_context_t *ctx);
475

476
int mpd_qcopy(mpd_t *result, const mpd_t *a, uint32_t *status);
477
int mpd_qcopy_cxx(mpd_t *result, const mpd_t *a);
478
mpd_t *mpd_qncopy(const mpd_t *a);
479
int mpd_qcopy_abs(mpd_t *result, const mpd_t *a, uint32_t *status);
480
int mpd_qcopy_negate(mpd_t *result, const mpd_t *a, uint32_t *status);
481
int mpd_qcopy_sign(mpd_t *result, const mpd_t *a, const mpd_t *b, uint32_t *status);
482

483
void mpd_qand(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
484
void mpd_qinvert(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
485
void mpd_qlogb(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
486
void mpd_qor(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
487
void mpd_qscaleb(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
488
void mpd_qxor(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
489
int mpd_same_quantum(const mpd_t *a, const mpd_t *b);
490

491
void mpd_qrotate(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
492
int mpd_qshiftl(mpd_t *result, const mpd_t *a, mpd_ssize_t n, uint32_t *status);
493
mpd_uint_t mpd_qshiftr(mpd_t *result, const mpd_t *a, mpd_ssize_t n, uint32_t *status);
494
mpd_uint_t mpd_qshiftr_inplace(mpd_t *result, mpd_ssize_t n);
495
void mpd_qshift(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
496
void mpd_qshiftn(mpd_t *result, const mpd_t *a, mpd_ssize_t n, const mpd_context_t *ctx, uint32_t *status);
497

498
int mpd_qcmp(const mpd_t *a, const mpd_t *b, uint32_t *status);
499
int mpd_qcompare(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
500
int mpd_qcompare_signal(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
501
int mpd_cmp_total(const mpd_t *a, const mpd_t *b);
502
int mpd_cmp_total_mag(const mpd_t *a, const mpd_t *b);
503
int mpd_compare_total(mpd_t *result, const mpd_t *a, const mpd_t *b);
504
int mpd_compare_total_mag(mpd_t *result, const mpd_t *a, const mpd_t *b);
505

506
void mpd_qround_to_intx(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
507
void mpd_qround_to_int(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
508
void mpd_qtrunc(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
509
void mpd_qfloor(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
510
void mpd_qceil(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
511

512
void mpd_qabs(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
513
void mpd_qmax(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
514
void mpd_qmax_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
515
void mpd_qmin(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
516
void mpd_qmin_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
517
void mpd_qminus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
518
void mpd_qplus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
519
void mpd_qnext_minus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
520
void mpd_qnext_plus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
521
void mpd_qnext_toward(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
522
void mpd_qquantize(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
523
void mpd_qrescale(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, const mpd_context_t *ctx, uint32_t *status);
524
void mpd_qrescale_fmt(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, const mpd_context_t *ctx, uint32_t *status);
525
void mpd_qreduce(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
526
void mpd_qadd(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
527
void mpd_qadd_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
528
void mpd_qadd_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
529
void mpd_qadd_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
530
void mpd_qadd_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
531
void mpd_qsub(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
532
void mpd_qsub_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
533
void mpd_qsub_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
534
void mpd_qsub_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
535
void mpd_qsub_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
536
void mpd_qmul(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
537
void mpd_qmul_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
538
void mpd_qmul_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
539
void mpd_qmul_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
540
void mpd_qmul_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
541
void 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);
542
void mpd_qdiv(mpd_t *q, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
543
void mpd_qdiv_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
544
void mpd_qdiv_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
545
void mpd_qdiv_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
546
void mpd_qdiv_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
547
void mpd_qdivint(mpd_t *q, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
548
void mpd_qrem(mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
549
void mpd_qrem_near(mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
550
void mpd_qdivmod(mpd_t *q, mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
551
void mpd_qpow(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_context_t *ctx, uint32_t *status);
552
void 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);
553
void mpd_qexp(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
554
void mpd_qln10(mpd_t *result, mpd_ssize_t prec, uint32_t *status);
555
void mpd_qln(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
556
void mpd_qlog10(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
557
void mpd_qsqrt(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
558
void mpd_qinvroot(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
559

560
#ifndef LEGACY_COMPILER
561
void mpd_qadd_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
562
void mpd_qadd_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
563
void mpd_qsub_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
564
void mpd_qsub_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
565
void mpd_qmul_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
566
void mpd_qmul_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
567
void mpd_qdiv_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
568
void mpd_qdiv_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
569
#endif
570

571

572
size_t mpd_sizeinbase(const mpd_t *a, uint32_t base);
573
void mpd_qimport_u16(mpd_t *result, const uint16_t *srcdata, size_t srclen,
574
                     uint8_t srcsign, uint32_t srcbase,
575
                     const mpd_context_t *ctx, uint32_t *status);
576
void mpd_qimport_u32(mpd_t *result, const uint32_t *srcdata, size_t srclen,
577
                     uint8_t srcsign, uint32_t srcbase,
578
                     const mpd_context_t *ctx, uint32_t *status);
579
size_t mpd_qexport_u16(uint16_t **rdata, size_t rlen, uint32_t base,
580
                       const mpd_t *src, uint32_t *status);
581
size_t mpd_qexport_u32(uint32_t **rdata, size_t rlen, uint32_t base,
582
                       const mpd_t *src, uint32_t *status);
583

584

585
/******************************************************************************/
586
/*                           Signalling functions                             */
587
/******************************************************************************/
588

589
char *mpd_format(const mpd_t *dec, const char *fmt, mpd_context_t *ctx);
590
void mpd_import_u16(mpd_t *result, const uint16_t *srcdata, size_t srclen, uint8_t srcsign, uint32_t base, mpd_context_t *ctx);
591
void mpd_import_u32(mpd_t *result, const uint32_t *srcdata, size_t srclen, uint8_t srcsign, uint32_t base, mpd_context_t *ctx);
592
size_t mpd_export_u16(uint16_t **rdata, size_t rlen, uint32_t base, const mpd_t *src, mpd_context_t *ctx);
593
size_t mpd_export_u32(uint32_t **rdata, size_t rlen, uint32_t base, const mpd_t *src, mpd_context_t *ctx);
594
void mpd_finalize(mpd_t *result, mpd_context_t *ctx);
595
int mpd_check_nan(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
596
int mpd_check_nans(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
597
void mpd_set_string(mpd_t *result, const char *s, mpd_context_t *ctx);
598
void mpd_maxcoeff(mpd_t *result, mpd_context_t *ctx);
599
void mpd_sset_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx);
600
void mpd_sset_i32(mpd_t *result, int32_t a, mpd_context_t *ctx);
601
void mpd_sset_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx);
602
void mpd_sset_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx);
603
void mpd_set_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx);
604
void mpd_set_i32(mpd_t *result, int32_t a, mpd_context_t *ctx);
605
void mpd_set_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx);
606
void mpd_set_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx);
607
#ifndef LEGACY_COMPILER
608
void mpd_set_i64(mpd_t *result, int64_t a, mpd_context_t *ctx);
609
void mpd_set_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx);
610
#endif
611
mpd_ssize_t mpd_get_ssize(const mpd_t *a, mpd_context_t *ctx);
612
mpd_uint_t mpd_get_uint(const mpd_t *a, mpd_context_t *ctx);
613
mpd_uint_t mpd_abs_uint(const mpd_t *a, mpd_context_t *ctx);
614
int32_t mpd_get_i32(const mpd_t *a, mpd_context_t *ctx);
615
uint32_t mpd_get_u32(const mpd_t *a, mpd_context_t *ctx);
616
#ifndef LEGACY_COMPILER
617
int64_t mpd_get_i64(const mpd_t *a, mpd_context_t *ctx);
618
uint64_t mpd_get_u64(const mpd_t *a, mpd_context_t *ctx);
619
#endif
620
void mpd_and(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
621
void mpd_copy(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
622
void mpd_canonical(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
623
void mpd_copy_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
624
void mpd_copy_negate(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
625
void mpd_copy_sign(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
626
void mpd_invert(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
627
void mpd_logb(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
628
void mpd_or(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
629
void mpd_rotate(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
630
void mpd_scaleb(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
631
void mpd_shiftl(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx);
632
mpd_uint_t mpd_shiftr(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx);
633
void mpd_shiftn(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx);
634
void mpd_shift(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
635
void mpd_xor(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
636
void mpd_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
637
int mpd_cmp(const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
638
int mpd_compare(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
639
int mpd_compare_signal(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
640
void mpd_add(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
641
void mpd_add_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
642
void mpd_add_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
643
void mpd_add_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
644
void mpd_add_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
645
void mpd_sub(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
646
void mpd_sub_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
647
void mpd_sub_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
648
void mpd_sub_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
649
void mpd_sub_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
650
void mpd_div(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
651
void mpd_div_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
652
void mpd_div_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
653
void mpd_div_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
654
void mpd_div_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
655
void mpd_divmod(mpd_t *q, mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
656
void mpd_divint(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
657
void mpd_exp(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
658
void mpd_fma(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_t *c, mpd_context_t *ctx);
659
void mpd_ln(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
660
void mpd_log10(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
661
void mpd_max(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
662
void mpd_max_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
663
void mpd_min(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
664
void mpd_min_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
665
void mpd_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
666
void mpd_mul(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
667
void mpd_mul_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
668
void mpd_mul_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
669
void mpd_mul_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
670
void mpd_mul_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
671
void mpd_next_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
672
void mpd_next_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
673
void mpd_next_toward(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
674
void mpd_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
675
void mpd_pow(mpd_t *result, const mpd_t *base, const mpd_t *exp, mpd_context_t *ctx);
676
void mpd_powmod(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_t *mod, mpd_context_t *ctx);
677
void mpd_quantize(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
678
void mpd_rescale(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, mpd_context_t *ctx);
679
void mpd_reduce(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
680
void mpd_rem(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
681
void mpd_rem_near(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
682
void mpd_round_to_intx(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
683
void mpd_round_to_int(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
684
void mpd_trunc(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
685
void mpd_floor(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
686
void mpd_ceil(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
687
void mpd_sqrt(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
688
void mpd_invroot(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
689

690
#ifndef LEGACY_COMPILER
691
void mpd_add_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
692
void mpd_add_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
693
void mpd_sub_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
694
void mpd_sub_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
695
void mpd_div_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
696
void mpd_div_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
697
void mpd_mul_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
698
void mpd_mul_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
699
#endif
700

701

702
/******************************************************************************/
703
/*                          Configuration specific                            */
704
/******************************************************************************/
705

706
#ifdef CONFIG_64
707
void mpd_qsset_i64(mpd_t *result, int64_t a, const mpd_context_t *ctx, uint32_t *status);
708
void mpd_qsset_u64(mpd_t *result, uint64_t a, const mpd_context_t *ctx, uint32_t *status);
709
void mpd_sset_i64(mpd_t *result, int64_t a, mpd_context_t *ctx);
710
void mpd_sset_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx);
711
#endif
712

713

714
/******************************************************************************/
715
/*                       Get attributes of a decimal                          */
716
/******************************************************************************/
717

718
EXTINLINE mpd_ssize_t mpd_adjexp(const mpd_t *dec);
719
EXTINLINE mpd_ssize_t mpd_etiny(const mpd_context_t *ctx);
720
EXTINLINE mpd_ssize_t mpd_etop(const mpd_context_t *ctx);
721
EXTINLINE mpd_uint_t mpd_msword(const mpd_t *dec);
722
EXTINLINE int mpd_word_digits(mpd_uint_t word);
723
/* most significant digit of a word */
724
EXTINLINE mpd_uint_t mpd_msd(mpd_uint_t word);
725
/* least significant digit of a word */
726
EXTINLINE mpd_uint_t mpd_lsd(mpd_uint_t word);
727
/* coefficient size needed to store 'digits' */
728
EXTINLINE mpd_ssize_t mpd_digits_to_size(mpd_ssize_t digits);
729
/* number of digits in the exponent, undefined for MPD_SSIZE_MIN */
730
EXTINLINE int mpd_exp_digits(mpd_ssize_t exp);
731
EXTINLINE int mpd_iscanonical(const mpd_t *dec);
732
EXTINLINE int mpd_isfinite(const mpd_t *dec);
733
EXTINLINE int mpd_isinfinite(const mpd_t *dec);
734
EXTINLINE int mpd_isinteger(const mpd_t *dec);
735
EXTINLINE int mpd_isnan(const mpd_t *dec);
736
EXTINLINE int mpd_isnegative(const mpd_t *dec);
737
EXTINLINE int mpd_ispositive(const mpd_t *dec);
738
EXTINLINE int mpd_isqnan(const mpd_t *dec);
739
EXTINLINE int mpd_issigned(const mpd_t *dec);
740
EXTINLINE int mpd_issnan(const mpd_t *dec);
741
EXTINLINE int mpd_isspecial(const mpd_t *dec);
742
EXTINLINE int mpd_iszero(const mpd_t *dec);
743
/* undefined for special numbers */
744
EXTINLINE int mpd_iszerocoeff(const mpd_t *dec);
745
EXTINLINE int mpd_isnormal(const mpd_t *dec, const mpd_context_t *ctx);
746
EXTINLINE int mpd_issubnormal(const mpd_t *dec, const mpd_context_t *ctx);
747
/* odd word */
748
EXTINLINE int mpd_isoddword(mpd_uint_t word);
749
/* odd coefficient */
750
EXTINLINE int mpd_isoddcoeff(const mpd_t *dec);
751
/* odd decimal, only defined for integers */
752
int mpd_isodd(const mpd_t *dec);
753
/* even decimal, only defined for integers */
754
int mpd_iseven(const mpd_t *dec);
755
/* 0 if dec is positive, 1 if dec is negative */
756
EXTINLINE uint8_t mpd_sign(const mpd_t *dec);
757
/* 1 if dec is positive, -1 if dec is negative */
758
EXTINLINE int mpd_arith_sign(const mpd_t *dec);
759
EXTINLINE long mpd_radix(void);
760
EXTINLINE int mpd_isdynamic(const mpd_t *dec);
761
EXTINLINE int mpd_isstatic(const mpd_t *dec);
762
EXTINLINE int mpd_isdynamic_data(const mpd_t *dec);
763
EXTINLINE int mpd_isstatic_data(const mpd_t *dec);
764
EXTINLINE int mpd_isshared_data(const mpd_t *dec);
765
EXTINLINE int mpd_isconst_data(const mpd_t *dec);
766
EXTINLINE mpd_ssize_t mpd_trail_zeros(const mpd_t *dec);
767

768

769
/******************************************************************************/
770
/*                       Set attributes of a decimal                          */
771
/******************************************************************************/
772

773
/* set number of decimal digits in the coefficient */
774
EXTINLINE void mpd_setdigits(mpd_t *result);
775
EXTINLINE void mpd_set_sign(mpd_t *result, uint8_t sign);
776
/* copy sign from another decimal */
777
EXTINLINE void mpd_signcpy(mpd_t *result, const mpd_t *a);
778
EXTINLINE void mpd_set_infinity(mpd_t *result);
779
EXTINLINE void mpd_set_qnan(mpd_t *result);
780
EXTINLINE void mpd_set_snan(mpd_t *result);
781
EXTINLINE void mpd_set_negative(mpd_t *result);
782
EXTINLINE void mpd_set_positive(mpd_t *result);
783
EXTINLINE void mpd_set_dynamic(mpd_t *result);
784
EXTINLINE void mpd_set_static(mpd_t *result);
785
EXTINLINE void mpd_set_dynamic_data(mpd_t *result);
786
EXTINLINE void mpd_set_static_data(mpd_t *result);
787
EXTINLINE void mpd_set_shared_data(mpd_t *result);
788
EXTINLINE void mpd_set_const_data(mpd_t *result);
789
EXTINLINE void mpd_clear_flags(mpd_t *result);
790
EXTINLINE void mpd_set_flags(mpd_t *result, uint8_t flags);
791
EXTINLINE void mpd_copy_flags(mpd_t *result, const mpd_t *a);
792

793

794
/******************************************************************************/
795
/*                              Error Macros                                  */
796
/******************************************************************************/
797

798
#define mpd_err_fatal(...) \
799
    do {fprintf(stderr, "%s:%d: error: ", __FILE__, __LINE__); \
800
        fprintf(stderr, __VA_ARGS__);  fputc('\n', stderr);    \
801
        abort();                                               \
802
    } while (0)
803
#define mpd_err_warn(...) \
804
    do {fprintf(stderr, "%s:%d: warning: ", __FILE__, __LINE__); \
805
        fprintf(stderr, __VA_ARGS__); fputc('\n', stderr);       \
806
    } while (0)
807

808

809
/******************************************************************************/
810
/*                            Memory handling                                 */
811
/******************************************************************************/
812

813
extern void *(* mpd_mallocfunc)(size_t size);
814
extern void *(* mpd_callocfunc)(size_t nmemb, size_t size);
815
extern void *(* mpd_reallocfunc)(void *ptr, size_t size);
816
extern void (* mpd_free)(void *ptr);
817

818
void *mpd_callocfunc_em(size_t nmemb, size_t size);
819

820
void *mpd_alloc(mpd_size_t nmemb, mpd_size_t size);
821
void *mpd_calloc(mpd_size_t nmemb, mpd_size_t size);
822
void *mpd_realloc(void *ptr, mpd_size_t nmemb, mpd_size_t size, uint8_t *err);
823
void *mpd_sh_alloc(mpd_size_t struct_size, mpd_size_t nmemb, mpd_size_t size);
824

825
mpd_t *mpd_qnew(void);
826
mpd_t *mpd_new(mpd_context_t *ctx);
827
mpd_t *mpd_qnew_size(mpd_ssize_t nwords);
828
EXTINLINE void mpd_del(mpd_t *dec);
829

830
EXTINLINE void mpd_uint_zero(mpd_uint_t *dest, mpd_size_t len);
831
EXTINLINE int mpd_qresize(mpd_t *result, mpd_ssize_t nwords, uint32_t *status);
832
EXTINLINE int mpd_qresize_zero(mpd_t *result, mpd_ssize_t nwords, uint32_t *status);
833
EXTINLINE void mpd_minalloc(mpd_t *result);
834

835
int mpd_resize(mpd_t *result, mpd_ssize_t nwords, mpd_context_t *ctx);
836
int mpd_resize_zero(mpd_t *result, mpd_ssize_t nwords, mpd_context_t *ctx);
837

838

839
MPD_PRAGMA(MPD_HIDE_SYMBOLS_END) /* restore previous scope rules */
840

841

842
#ifdef __cplusplus
843
} /* END extern "C" */
844
#endif
845

846

847
#endif /* LIBMPDEC_MPDECIMAL_H_ */
848

849
Product

Resources

Company
