/*
 * Copyright (c) 2008-2016 Stefan Krah. All rights reserved.
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions
 * are met:
 *
 * 1. Redistributions of source code must retain the above copyright
 *    notice, this list of conditions and the following disclaimer.
 *
 * 2. Redistributions in binary form must reproduce the above copyright
 *    notice, this list of conditions and the following disclaimer in the
 *    documentation and/or other materials provided with the distribution.
 *
 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS "AS IS" AND
 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
 * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
 * SUCH DAMAGE.
 */


#ifndef MPDECIMAL_H
#define MPDECIMAL_H


#ifdef __cplusplus
extern "C" {
  #ifndef __STDC_LIMIT_MACROS
    #define __STDC_LIMIT_MACROS
    #define MPD_CLEAR_STDC_LIMIT_MACROS
  #endif
#endif


#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <limits.h>
#include <assert.h>
#include <stdint.h>
#include <inttypes.h>


#ifndef __GNUC_STDC_INLINE__
  #define __GNUC_STDC_INLINE__ 1
#endif
#if defined(__GNUC__) && !defined(__INTEL_COMPILER)
  #define UNUSED __attribute__((unused))
#else
  #define UNUSED
#endif
#if (defined(__linux__) || defined(__FreeBSD__) || defined(__APPLE__)) && \
    defined(__GNUC__) && __GNUC__ >= 4 && !defined(__INTEL_COMPILER)
  #define MPD_PRAGMA(x) _Pragma(x)
  #define MPD_HIDE_SYMBOLS_START "GCC visibility push(hidden)"
  #define MPD_HIDE_SYMBOLS_END "GCC visibility pop"
#else
  #define MPD_PRAGMA(x)
  #define MPD_HIDE_SYMBOLS_START
  #define MPD_HIDE_SYMBOLS_END
#endif


#if !defined(LEGACY_COMPILER)
  #if !defined(UINT64_MAX)
    /* The following #error is just a warning. If the compiler indeed does
     * not have uint64_t, it is perfectly safe to comment out the #error. */
    #error "Warning: Compiler without uint64_t. Comment out this line."
    #define LEGACY_COMPILER
  #endif
#endif

#if _CRTBLD
#define __MANGLE_ME(x) __mingw_dfp_##x
#define _mpd_baseadd __MANGLE_ME(_mpd_baseadd)
#define _mpd_baseaddto __MANGLE_ME(_mpd_baseaddto)
#define _mpd_basedivmod __MANGLE_ME(_mpd_basedivmod)
#define _mpd_baseincr __MANGLE_ME(_mpd_baseincr)
#define _mpd_basemul __MANGLE_ME(_mpd_basemul)
#define _mpd_baseshiftl __MANGLE_ME(_mpd_baseshiftl)
#define _mpd_baseshiftr __MANGLE_ME(_mpd_baseshiftr)
#define _mpd_basesub __MANGLE_ME(_mpd_basesub)
#define _mpd_basesubfrom __MANGLE_ME(_mpd_basesubfrom)
#define _mpd_getkernel __MANGLE_ME(_mpd_getkernel)
#define _mpd_init_fnt_params __MANGLE_ME(_mpd_init_fnt_params)
#define _mpd_init_w3table __MANGLE_ME(_mpd_init_w3table)
#define _mpd_shortadd __MANGLE_ME(_mpd_shortadd)
#define _mpd_shortadd_b __MANGLE_ME(_mpd_shortadd_b)
#define _mpd_shortdiv __MANGLE_ME(_mpd_shortdiv)
#define _mpd_shortdiv_b __MANGLE_ME(_mpd_shortdiv_b)
#define _mpd_shortmul __MANGLE_ME(_mpd_shortmul)
#define _mpd_shortmul_b __MANGLE_ME(_mpd_shortmul_b)
#define _mpd_shortmul_c __MANGLE_ME(_mpd_shortmul_c)
#define crt3 __MANGLE_ME(crt3)
#define fnt_autoconvolute __MANGLE_ME(fnt_autoconvolute)
#define fnt_convolute __MANGLE_ME(fnt_convolute)
#define fnt_dif2 __MANGLE_ME(fnt_dif2)
#define four_step_fnt __MANGLE_ME(four_step_fnt)
#define inv_four_step_fnt __MANGLE_ME(inv_four_step_fnt)
#define inv_six_step_fnt __MANGLE_ME(inv_six_step_fnt)
#define mpd_abs __MANGLE_ME(mpd_abs)
#define mpd_abs_uint __MANGLE_ME(mpd_abs_uint)
#define mpd_add __MANGLE_ME(mpd_add)
#define mpd_add_i32 __MANGLE_ME(mpd_add_i32)
#define mpd_add_i64 __MANGLE_ME(mpd_add_i64)
#define mpd_add_ssize __MANGLE_ME(mpd_add_ssize)
#define mpd_add_u32 __MANGLE_ME(mpd_add_u32)
#define mpd_add_u64 __MANGLE_ME(mpd_add_u64)
#define mpd_add_uint __MANGLE_ME(mpd_add_uint)
#define mpd_addstatus_raise __MANGLE_ME(mpd_addstatus_raise)
#define mpd_adjexp __MANGLE_ME(mpd_adjexp)
#define mpd_alloc __MANGLE_ME(mpd_alloc)
#define mpd_and __MANGLE_ME(mpd_and)
#define mpd_arith_sign __MANGLE_ME(mpd_arith_sign)
#define mpd_basiccontext __MANGLE_ME(mpd_basiccontext)
#define mpd_calloc __MANGLE_ME(mpd_calloc)
#define mpd_callocfunc_em __MANGLE_ME(mpd_callocfunc_em)
#define mpd_canonical __MANGLE_ME(mpd_canonical)
#define mpd_ceil __MANGLE_ME(mpd_ceil)
#define mpd_check_nan __MANGLE_ME(mpd_check_nan)
#define mpd_check_nans __MANGLE_ME(mpd_check_nans)
#define mpd_class __MANGLE_ME(mpd_class)
#define mpd_clear_flags __MANGLE_ME(mpd_clear_flags)
#define mpd_cmp __MANGLE_ME(mpd_cmp)
#define mpd_cmp_total __MANGLE_ME(mpd_cmp_total)
#define mpd_cmp_total_mag __MANGLE_ME(mpd_cmp_total_mag)
#define mpd_compare __MANGLE_ME(mpd_compare)
#define mpd_compare_signal __MANGLE_ME(mpd_compare_signal)
#define mpd_compare_total __MANGLE_ME(mpd_compare_total)
#define mpd_compare_total_mag __MANGLE_ME(mpd_compare_total_mag)
#define mpd_copy __MANGLE_ME(mpd_copy)
#define mpd_copy_abs __MANGLE_ME(mpd_copy_abs)
#define mpd_copy_flags __MANGLE_ME(mpd_copy_flags)
#define mpd_copy_negate __MANGLE_ME(mpd_copy_negate)
#define mpd_copy_sign __MANGLE_ME(mpd_copy_sign)
#define mpd_defaultcontext __MANGLE_ME(mpd_defaultcontext)
#define mpd_del __MANGLE_ME(mpd_del)
#define mpd_dflt_traphandler __MANGLE_ME(mpd_dflt_traphandler)
#define mpd_digits_to_size __MANGLE_ME(mpd_digits_to_size)
#define mpd_div __MANGLE_ME(mpd_div)
#define mpd_div_i32 __MANGLE_ME(mpd_div_i32)
#define mpd_div_i64 __MANGLE_ME(mpd_div_i64)
#define mpd_div_ssize __MANGLE_ME(mpd_div_ssize)
#define mpd_div_u32 __MANGLE_ME(mpd_div_u32)
#define mpd_div_u64 __MANGLE_ME(mpd_div_u64)
#define mpd_div_uint __MANGLE_ME(mpd_div_uint)
#define mpd_divint __MANGLE_ME(mpd_divint)
#define mpd_divmod __MANGLE_ME(mpd_divmod)
#define mpd_etiny __MANGLE_ME(mpd_etiny)
#define mpd_etop __MANGLE_ME(mpd_etop)
#define mpd_exp __MANGLE_ME(mpd_exp)
#define mpd_exp_digits __MANGLE_ME(mpd_exp_digits)
#define mpd_export_u16 __MANGLE_ME(mpd_export_u16)
#define mpd_export_u32 __MANGLE_ME(mpd_export_u32)
#define mpd_finalize __MANGLE_ME(mpd_finalize)
#define mpd_floor __MANGLE_ME(mpd_floor)
#define mpd_fma __MANGLE_ME(mpd_fma)
#define mpd_format __MANGLE_ME(mpd_format)
#define mpd_fprint __MANGLE_ME(mpd_fprint)
#define mpd_get_i32 __MANGLE_ME(mpd_get_i32)
#define mpd_get_i64 __MANGLE_ME(mpd_get_i64)
#define mpd_get_ssize __MANGLE_ME(mpd_get_ssize)
#define mpd_get_u32 __MANGLE_ME(mpd_get_u32)
#define mpd_get_u64 __MANGLE_ME(mpd_get_u64)
#define mpd_get_uint __MANGLE_ME(mpd_get_uint)
#define mpd_getclamp __MANGLE_ME(mpd_getclamp)
#define mpd_getcr __MANGLE_ME(mpd_getcr)
#define mpd_getemax __MANGLE_ME(mpd_getemax)
#define mpd_getemin __MANGLE_ME(mpd_getemin)
#define mpd_getprec __MANGLE_ME(mpd_getprec)
#define mpd_getround __MANGLE_ME(mpd_getround)
#define mpd_getstatus __MANGLE_ME(mpd_getstatus)
#define mpd_gettraps __MANGLE_ME(mpd_gettraps)
#define mpd_ieee_context __MANGLE_ME(mpd_ieee_context)
#define mpd_import_u16 __MANGLE_ME(mpd_import_u16)
#define mpd_import_u32 __MANGLE_ME(mpd_import_u32)
#define mpd_init __MANGLE_ME(mpd_init)
#define mpd_invert __MANGLE_ME(mpd_invert)
#define mpd_invroot __MANGLE_ME(mpd_invroot)
#define mpd_iscanonical __MANGLE_ME(mpd_iscanonical)
#define mpd_isconst_data __MANGLE_ME(mpd_isconst_data)
#define mpd_isdynamic __MANGLE_ME(mpd_isdynamic)
#define mpd_isdynamic_data __MANGLE_ME(mpd_isdynamic_data)
#define mpd_iseven __MANGLE_ME(mpd_iseven)
#define mpd_isfinite __MANGLE_ME(mpd_isfinite)
#define mpd_isinfinite __MANGLE_ME(mpd_isinfinite)
#define mpd_isinteger __MANGLE_ME(mpd_isinteger)
#define mpd_isnan __MANGLE_ME(mpd_isnan)
#define mpd_isnegative __MANGLE_ME(mpd_isnegative)
#define mpd_isnormal __MANGLE_ME(mpd_isnormal)
#define mpd_isodd __MANGLE_ME(mpd_isodd)
#define mpd_isoddcoeff __MANGLE_ME(mpd_isoddcoeff)
#define mpd_isoddword __MANGLE_ME(mpd_isoddword)
#define mpd_ispositive __MANGLE_ME(mpd_ispositive)
#define mpd_isqnan __MANGLE_ME(mpd_isqnan)
#define mpd_isshared_data __MANGLE_ME(mpd_isshared_data)
#define mpd_issigned __MANGLE_ME(mpd_issigned)
#define mpd_issnan __MANGLE_ME(mpd_issnan)
#define mpd_isspecial __MANGLE_ME(mpd_isspecial)
#define mpd_isstatic __MANGLE_ME(mpd_isstatic)
#define mpd_isstatic_data __MANGLE_ME(mpd_isstatic_data)
#define mpd_issubnormal __MANGLE_ME(mpd_issubnormal)
#define mpd_iszero __MANGLE_ME(mpd_iszero)
#define mpd_iszerocoeff __MANGLE_ME(mpd_iszerocoeff)
#define mpd_ln __MANGLE_ME(mpd_ln)
#define mpd_log10 __MANGLE_ME(mpd_log10)
#define mpd_logb __MANGLE_ME(mpd_logb)
#define mpd_lsd __MANGLE_ME(mpd_lsd)
#define mpd_lsnprint_flags __MANGLE_ME(mpd_lsnprint_flags)
#define mpd_lsnprint_signals __MANGLE_ME(mpd_lsnprint_signals)
#define mpd_max __MANGLE_ME(mpd_max)
#define mpd_max_mag __MANGLE_ME(mpd_max_mag)
#define mpd_maxcoeff __MANGLE_ME(mpd_maxcoeff)
#define mpd_maxcontext __MANGLE_ME(mpd_maxcontext)
#define mpd_min __MANGLE_ME(mpd_min)
#define mpd_min_mag __MANGLE_ME(mpd_min_mag)
#define mpd_minalloc __MANGLE_ME(mpd_minalloc)
#define mpd_minus __MANGLE_ME(mpd_minus)
#define mpd_msd __MANGLE_ME(mpd_msd)
#define mpd_msword __MANGLE_ME(mpd_msword)
#define mpd_mul __MANGLE_ME(mpd_mul)
#define mpd_mul_i32 __MANGLE_ME(mpd_mul_i32)
#define mpd_mul_i64 __MANGLE_ME(mpd_mul_i64)
#define mpd_mul_ssize __MANGLE_ME(mpd_mul_ssize)
#define mpd_mul_u32 __MANGLE_ME(mpd_mul_u32)
#define mpd_mul_u64 __MANGLE_ME(mpd_mul_u64)
#define mpd_mul_uint __MANGLE_ME(mpd_mul_uint)
#define mpd_new __MANGLE_ME(mpd_new)
#define mpd_next_minus __MANGLE_ME(mpd_next_minus)
#define mpd_next_plus __MANGLE_ME(mpd_next_plus)
#define mpd_next_toward __MANGLE_ME(mpd_next_toward)
#define mpd_or __MANGLE_ME(mpd_or)
#define mpd_parse_fmt_str __MANGLE_ME(mpd_parse_fmt_str)
#define mpd_plus __MANGLE_ME(mpd_plus)
#define mpd_pow __MANGLE_ME(mpd_pow)
#define mpd_powmod __MANGLE_ME(mpd_powmod)
#define mpd_print __MANGLE_ME(mpd_print)
#define mpd_qabs __MANGLE_ME(mpd_qabs)
#define mpd_qabs_uint __MANGLE_ME(mpd_qabs_uint)
#define mpd_qadd __MANGLE_ME(mpd_qadd)
#define mpd_qadd_i32 __MANGLE_ME(mpd_qadd_i32)
#define mpd_qadd_i64 __MANGLE_ME(mpd_qadd_i64)
#define mpd_qadd_ssize __MANGLE_ME(mpd_qadd_ssize)
#define mpd_qadd_u32 __MANGLE_ME(mpd_qadd_u32)
#define mpd_qadd_u64 __MANGLE_ME(mpd_qadd_u64)
#define mpd_qadd_uint __MANGLE_ME(mpd_qadd_uint)
#define mpd_qand __MANGLE_ME(mpd_qand)
#define mpd_qceil __MANGLE_ME(mpd_qceil)
#define mpd_qcheck_nan __MANGLE_ME(mpd_qcheck_nan)
#define mpd_qcheck_nans __MANGLE_ME(mpd_qcheck_nans)
#define mpd_qcmp __MANGLE_ME(mpd_qcmp)
#define mpd_qcompare __MANGLE_ME(mpd_qcompare)
#define mpd_qcompare_signal __MANGLE_ME(mpd_qcompare_signal)
#define mpd_qcopy __MANGLE_ME(mpd_qcopy)
#define mpd_qcopy_abs __MANGLE_ME(mpd_qcopy_abs)
#define mpd_qcopy_negate __MANGLE_ME(mpd_qcopy_negate)
#define mpd_qcopy_sign __MANGLE_ME(mpd_qcopy_sign)
#define mpd_qdiv __MANGLE_ME(mpd_qdiv)
#define mpd_qdiv_i32 __MANGLE_ME(mpd_qdiv_i32)
#define mpd_qdiv_i64 __MANGLE_ME(mpd_qdiv_i64)
#define mpd_qdiv_ssize __MANGLE_ME(mpd_qdiv_ssize)
#define mpd_qdiv_u32 __MANGLE_ME(mpd_qdiv_u32)
#define mpd_qdiv_u64 __MANGLE_ME(mpd_qdiv_u64)
#define mpd_qdiv_uint __MANGLE_ME(mpd_qdiv_uint)
#define mpd_qdivint __MANGLE_ME(mpd_qdivint)
#define mpd_qdivmod __MANGLE_ME(mpd_qdivmod)
#define mpd_qexp __MANGLE_ME(mpd_qexp)
#define mpd_qexport_u16 __MANGLE_ME(mpd_qexport_u16)
#define mpd_qexport_u32 __MANGLE_ME(mpd_qexport_u32)
#define mpd_qfinalize __MANGLE_ME(mpd_qfinalize)
#define mpd_qfloor __MANGLE_ME(mpd_qfloor)
#define mpd_qfma __MANGLE_ME(mpd_qfma)
#define mpd_qformat __MANGLE_ME(mpd_qformat)
#define mpd_qformat_spec __MANGLE_ME(mpd_qformat_spec)
#define mpd_qget_i32 __MANGLE_ME(mpd_qget_i32)
#define mpd_qget_i64 __MANGLE_ME(mpd_qget_i64)
#define mpd_qget_ssize __MANGLE_ME(mpd_qget_ssize)
#define mpd_qget_u32 __MANGLE_ME(mpd_qget_u32)
#define mpd_qget_u64 __MANGLE_ME(mpd_qget_u64)
#define mpd_qget_uint __MANGLE_ME(mpd_qget_uint)
#define mpd_qimport_u16 __MANGLE_ME(mpd_qimport_u16)
#define mpd_qimport_u32 __MANGLE_ME(mpd_qimport_u32)
#define mpd_qinvert __MANGLE_ME(mpd_qinvert)
#define mpd_qinvroot __MANGLE_ME(mpd_qinvroot)
#define mpd_qln __MANGLE_ME(mpd_qln)
#define mpd_qln10 __MANGLE_ME(mpd_qln10)
#define mpd_qlog10 __MANGLE_ME(mpd_qlog10)
#define mpd_qlogb __MANGLE_ME(mpd_qlogb)
#define mpd_qmax __MANGLE_ME(mpd_qmax)
#define mpd_qmax_mag __MANGLE_ME(mpd_qmax_mag)
#define mpd_qmaxcoeff __MANGLE_ME(mpd_qmaxcoeff)
#define mpd_qmin __MANGLE_ME(mpd_qmin)
#define mpd_qmin_mag __MANGLE_ME(mpd_qmin_mag)
#define mpd_qminus __MANGLE_ME(mpd_qminus)
#define mpd_qmul __MANGLE_ME(mpd_qmul)
#define mpd_qmul_i32 __MANGLE_ME(mpd_qmul_i32)
#define mpd_qmul_i64 __MANGLE_ME(mpd_qmul_i64)
#define mpd_qmul_ssize __MANGLE_ME(mpd_qmul_ssize)
#define mpd_qmul_u32 __MANGLE_ME(mpd_qmul_u32)
#define mpd_qmul_u64 __MANGLE_ME(mpd_qmul_u64)
#define mpd_qmul_uint __MANGLE_ME(mpd_qmul_uint)
#define mpd_qncopy __MANGLE_ME(mpd_qncopy)
#define mpd_qnew __MANGLE_ME(mpd_qnew)
#define mpd_qnew_size __MANGLE_ME(mpd_qnew_size)
#define mpd_qnext_minus __MANGLE_ME(mpd_qnext_minus)
#define mpd_qnext_plus __MANGLE_ME(mpd_qnext_plus)
#define mpd_qnext_toward __MANGLE_ME(mpd_qnext_toward)
#define mpd_qor __MANGLE_ME(mpd_qor)
#define mpd_qplus __MANGLE_ME(mpd_qplus)
#define mpd_qpow __MANGLE_ME(mpd_qpow)
#define mpd_qpowmod __MANGLE_ME(mpd_qpowmod)
#define mpd_qquantize __MANGLE_ME(mpd_qquantize)
#define mpd_qreduce __MANGLE_ME(mpd_qreduce)
#define mpd_qrem __MANGLE_ME(mpd_qrem)
#define mpd_qrem_near __MANGLE_ME(mpd_qrem_near)
#define mpd_qrescale __MANGLE_ME(mpd_qrescale)
#define mpd_qrescale_fmt __MANGLE_ME(mpd_qrescale_fmt)
#define mpd_qresize __MANGLE_ME(mpd_qresize)
#define mpd_qresize_zero __MANGLE_ME(mpd_qresize_zero)
#define mpd_qrotate __MANGLE_ME(mpd_qrotate)
#define mpd_qround_to_int __MANGLE_ME(mpd_qround_to_int)
#define mpd_qround_to_intx __MANGLE_ME(mpd_qround_to_intx)
#define mpd_qscaleb __MANGLE_ME(mpd_qscaleb)
#define mpd_qset_i32 __MANGLE_ME(mpd_qset_i32)
#define mpd_qset_i64 __MANGLE_ME(mpd_qset_i64)
#define mpd_qset_ssize __MANGLE_ME(mpd_qset_ssize)
#define mpd_qset_string __MANGLE_ME(mpd_qset_string)
#define mpd_qset_u32 __MANGLE_ME(mpd_qset_u32)
#define mpd_qset_u64 __MANGLE_ME(mpd_qset_u64)
#define mpd_qset_uint __MANGLE_ME(mpd_qset_uint)
#define mpd_qsetclamp __MANGLE_ME(mpd_qsetclamp)
#define mpd_qsetcr __MANGLE_ME(mpd_qsetcr)
#define mpd_qsetemax __MANGLE_ME(mpd_qsetemax)
#define mpd_qsetemin __MANGLE_ME(mpd_qsetemin)
#define mpd_qsetprec __MANGLE_ME(mpd_qsetprec)
#define mpd_qsetround __MANGLE_ME(mpd_qsetround)
#define mpd_qsetstatus __MANGLE_ME(mpd_qsetstatus)
#define mpd_qsettraps __MANGLE_ME(mpd_qsettraps)
#define mpd_qshift __MANGLE_ME(mpd_qshift)
#define mpd_qshiftl __MANGLE_ME(mpd_qshiftl)
#define mpd_qshiftn __MANGLE_ME(mpd_qshiftn)
#define mpd_qshiftr __MANGLE_ME(mpd_qshiftr)
#define mpd_qshiftr_inplace __MANGLE_ME(mpd_qshiftr_inplace)
#define mpd_qsqrt __MANGLE_ME(mpd_qsqrt)
#define mpd_qsset_i32 __MANGLE_ME(mpd_qsset_i32)
#define mpd_qsset_i64 __MANGLE_ME(mpd_qsset_i64)
#define mpd_qsset_ssize __MANGLE_ME(mpd_qsset_ssize)
#define mpd_qsset_u32 __MANGLE_ME(mpd_qsset_u32)
#define mpd_qsset_u64 __MANGLE_ME(mpd_qsset_u64)
#define mpd_qsset_uint __MANGLE_ME(mpd_qsset_uint)
#define mpd_qsub __MANGLE_ME(mpd_qsub)
#define mpd_qsub_i32 __MANGLE_ME(mpd_qsub_i32)
#define mpd_qsub_i64 __MANGLE_ME(mpd_qsub_i64)
#define mpd_qsub_ssize __MANGLE_ME(mpd_qsub_ssize)
#define mpd_qsub_u32 __MANGLE_ME(mpd_qsub_u32)
#define mpd_qsub_u64 __MANGLE_ME(mpd_qsub_u64)
#define mpd_qsub_uint __MANGLE_ME(mpd_qsub_uint)
#define mpd_qtrunc __MANGLE_ME(mpd_qtrunc)
#define mpd_quantize __MANGLE_ME(mpd_quantize)
#define mpd_qxor __MANGLE_ME(mpd_qxor)
#define mpd_radix __MANGLE_ME(mpd_radix)
#define mpd_realloc __MANGLE_ME(mpd_realloc)
#define mpd_realloc_dyn __MANGLE_ME(mpd_realloc_dyn)
#define mpd_reduce __MANGLE_ME(mpd_reduce)
#define mpd_rem __MANGLE_ME(mpd_rem)
#define mpd_rem_near __MANGLE_ME(mpd_rem_near)
#define mpd_rescale __MANGLE_ME(mpd_rescale)
#define mpd_resize __MANGLE_ME(mpd_resize)
#define mpd_resize_zero __MANGLE_ME(mpd_resize_zero)
#define mpd_rotate __MANGLE_ME(mpd_rotate)
#define mpd_round_to_int __MANGLE_ME(mpd_round_to_int)
#define mpd_round_to_intx __MANGLE_ME(mpd_round_to_intx)
#define mpd_same_quantum __MANGLE_ME(mpd_same_quantum)
#define mpd_scaleb __MANGLE_ME(mpd_scaleb)
#define mpd_set_const_data __MANGLE_ME(mpd_set_const_data)
#define mpd_set_dynamic __MANGLE_ME(mpd_set_dynamic)
#define mpd_set_dynamic_data __MANGLE_ME(mpd_set_dynamic_data)
#define mpd_set_flags __MANGLE_ME(mpd_set_flags)
#define mpd_set_i32 __MANGLE_ME(mpd_set_i32)
#define mpd_set_i64 __MANGLE_ME(mpd_set_i64)
#define mpd_set_infinity __MANGLE_ME(mpd_set_infinity)
#define mpd_set_negative __MANGLE_ME(mpd_set_negative)
#define mpd_set_positive __MANGLE_ME(mpd_set_positive)
#define mpd_set_qnan __MANGLE_ME(mpd_set_qnan)
#define mpd_set_shared_data __MANGLE_ME(mpd_set_shared_data)
#define mpd_set_sign __MANGLE_ME(mpd_set_sign)
#define mpd_set_snan __MANGLE_ME(mpd_set_snan)
#define mpd_set_ssize __MANGLE_ME(mpd_set_ssize)
#define mpd_set_static __MANGLE_ME(mpd_set_static)
#define mpd_set_static_data __MANGLE_ME(mpd_set_static_data)
#define mpd_set_string __MANGLE_ME(mpd_set_string)
#define mpd_set_u32 __MANGLE_ME(mpd_set_u32)
#define mpd_set_u64 __MANGLE_ME(mpd_set_u64)
#define mpd_set_uint __MANGLE_ME(mpd_set_uint)
#define mpd_setdigits __MANGLE_ME(mpd_setdigits)
#define mpd_seterror __MANGLE_ME(mpd_seterror)
#define mpd_setminalloc __MANGLE_ME(mpd_setminalloc)
#define mpd_setspecial __MANGLE_ME(mpd_setspecial)
#define mpd_sh_alloc __MANGLE_ME(mpd_sh_alloc)
#define mpd_shift __MANGLE_ME(mpd_shift)
#define mpd_shiftl __MANGLE_ME(mpd_shiftl)
#define mpd_shiftn __MANGLE_ME(mpd_shiftn)
#define mpd_shiftr __MANGLE_ME(mpd_shiftr)
#define mpd_sign __MANGLE_ME(mpd_sign)
#define mpd_signcpy __MANGLE_ME(mpd_signcpy)
#define mpd_sizeinbase __MANGLE_ME(mpd_sizeinbase)
#define mpd_snprint_flags __MANGLE_ME(mpd_snprint_flags)
#define mpd_sqrt __MANGLE_ME(mpd_sqrt)
#define mpd_sset_i32 __MANGLE_ME(mpd_sset_i32)
#define mpd_sset_i64 __MANGLE_ME(mpd_sset_i64)
#define mpd_sset_ssize __MANGLE_ME(mpd_sset_ssize)
#define mpd_sset_u32 __MANGLE_ME(mpd_sset_u32)
#define mpd_sset_u64 __MANGLE_ME(mpd_sset_u64)
#define mpd_sset_uint __MANGLE_ME(mpd_sset_uint)
#define mpd_sub __MANGLE_ME(mpd_sub)
#define mpd_sub_i32 __MANGLE_ME(mpd_sub_i32)
#define mpd_sub_i64 __MANGLE_ME(mpd_sub_i64)
#define mpd_sub_ssize __MANGLE_ME(mpd_sub_ssize)
#define mpd_sub_u32 __MANGLE_ME(mpd_sub_u32)
#define mpd_sub_u64 __MANGLE_ME(mpd_sub_u64)
#define mpd_sub_uint __MANGLE_ME(mpd_sub_uint)
#define mpd_switch_to_dyn __MANGLE_ME(mpd_switch_to_dyn)
#define mpd_switch_to_dyn_zero __MANGLE_ME(mpd_switch_to_dyn_zero)
#define mpd_to_eng __MANGLE_ME(mpd_to_eng)
#define mpd_to_eng_size __MANGLE_ME(mpd_to_eng_size)
#define mpd_to_sci __MANGLE_ME(mpd_to_sci)
#define mpd_to_sci_size __MANGLE_ME(mpd_to_sci_size)
#define mpd_trail_zeros __MANGLE_ME(mpd_trail_zeros)
#define mpd_trunc __MANGLE_ME(mpd_trunc)
#define mpd_uint_zero __MANGLE_ME(mpd_uint_zero)
#define mpd_validate_lconv __MANGLE_ME(mpd_validate_lconv)
#define mpd_version __MANGLE_ME(mpd_version)
#define mpd_word_digits __MANGLE_ME(mpd_word_digits)
#define mpd_xor __MANGLE_ME(mpd_xor)
#define mpd_zerocoeff __MANGLE_ME(mpd_zerocoeff)
#define six_step_fnt __MANGLE_ME(six_step_fnt)
#define std_fnt __MANGLE_ME(std_fnt)
#define std_inv_fnt __MANGLE_ME(std_inv_fnt)
#define std_trans __MANGLE_ME(std_trans)
#define transpose_pow2 __MANGLE_ME(transpose_pow2)
#endif


/******************************************************************************/
/*                                  Version                                   */
/******************************************************************************/

#define MPD_MAJOR_VERSION 2
#define MPD_MINOR_VERSION 4
#define MPD_MICRO_VERSION 0

#define MPD_VERSION "2.4.0"

const char *mpd_version(void);


/******************************************************************************/
/*                              Configuration                                 */
/******************************************************************************/

/* ABI: 32-bit */
#ifdef CONFIG_64
  #error "cannot use CONFIG_64 with 32-bit header."
#endif

#ifndef CONFIG_32
  #define CONFIG_32
#endif


/* BEGIN CONFIG_64 */
#if defined(CONFIG_64)
/* types for modular and base arithmetic */
#define MPD_UINT_MAX UINT64_MAX
#define MPD_BITS_PER_UINT 64
typedef uint64_t mpd_uint_t;  /* unsigned mod type */

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

/* type for exp, digits, len, prec */
#define MPD_SSIZE_MAX INT64_MAX
#define MPD_SSIZE_MIN INT64_MIN
typedef int64_t mpd_ssize_t;
#define _mpd_strtossize strtoll

/* decimal arithmetic */
#define MPD_RADIX 10000000000000000000ULL  /* 10**19 */
#define MPD_RDIGITS 19
#define MPD_MAX_POW10 19
#define MPD_EXPDIGITS 19  /* MPD_EXPDIGITS <= MPD_RDIGITS+1 */

#define MPD_MAXTRANSFORM_2N 4294967296ULL      /* 2**32 */
#define MPD_MAX_PREC 999999999999999999LL
#define MPD_MAX_PREC_LOG2 64
#define MPD_ELIMIT  1000000000000000000LL
#define MPD_MAX_EMAX   999999999999999999LL    /* ELIMIT-1 */
#define MPD_MIN_EMIN  (-999999999999999999LL)  /* -EMAX */
#define MPD_MIN_ETINY (MPD_MIN_EMIN-(MPD_MAX_PREC-1))
#define MPD_EXP_INF 2000000000000000001LL
#define MPD_EXP_CLAMP (-4000000000000000001LL)
#define MPD_MAXIMPORT 105263157894736842L /* ceil((2*MPD_MAX_PREC)/MPD_RDIGITS) */

/* conversion specifiers */
#define PRI_mpd_uint_t PRIu64
#define PRI_mpd_ssize_t PRIi64
/* END CONFIG_64 */


/* BEGIN CONFIG_32 */
#elif defined(CONFIG_32)
/* types for modular and base arithmetic */
#define MPD_UINT_MAX UINT32_MAX
#define MPD_BITS_PER_UINT 32
typedef uint32_t mpd_uint_t;  /* unsigned mod type */

#ifndef LEGACY_COMPILER
#define MPD_UUINT_MAX UINT64_MAX
typedef uint64_t mpd_uuint_t; /* double width unsigned mod type */
#endif

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

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

/* decimal arithmetic */
#define MPD_RADIX 1000000000UL  /* 10**9 */
#define MPD_RDIGITS 9
#define MPD_MAX_POW10 9
#define MPD_EXPDIGITS 10 /* MPD_EXPDIGITS <= MPD_RDIGITS+1 */

#define MPD_MAXTRANSFORM_2N 33554432UL /* 2**25 */
#define MPD_MAX_PREC 425000000L
#define MPD_MAX_PREC_LOG2 32
#define MPD_ELIMIT 425000001L
#define MPD_MAX_EMAX 425000000L        /* ELIMIT-1 */
#define MPD_MIN_EMIN (-425000000L)     /* -EMAX */
#define MPD_MIN_ETINY (MPD_MIN_EMIN-(MPD_MAX_PREC-1))
#define MPD_EXP_INF 1000000001L      /* allows for emax=999999999 in the tests */
#define MPD_EXP_CLAMP (-2000000001L) /* allows for emin=-999999999 in the tests */
#define MPD_MAXIMPORT 94444445L      /* ceil((2*MPD_MAX_PREC)/MPD_RDIGITS) */

/* conversion specifiers */
#define PRI_mpd_uint_t PRIu32
#define PRI_mpd_ssize_t PRIi32
/* END CONFIG_32 */

#else
  #error "define CONFIG_64 or CONFIG_32"
#endif
/* END CONFIG */


#if MPD_SIZE_MAX != MPD_UINT_MAX
  #error "unsupported platform: need mpd_size_t == mpd_uint_t"
#endif


/******************************************************************************/
/*                                Context                                     */
/******************************************************************************/

enum {
    MPD_ROUND_UP,          /* round away from 0               */
    MPD_ROUND_DOWN,        /* round toward 0 (truncate)       */
    MPD_ROUND_CEILING,     /* round toward +infinity          */
    MPD_ROUND_FLOOR,       /* round toward -infinity          */
    MPD_ROUND_HALF_UP,     /* 0.5 is rounded up               */
    MPD_ROUND_HALF_DOWN,   /* 0.5 is rounded down             */
    MPD_ROUND_HALF_EVEN,   /* 0.5 is rounded to even          */
    MPD_ROUND_05UP,        /* round zero or five away from 0  */
    MPD_ROUND_TRUNC,       /* truncate, but set infinity      */
    MPD_ROUND_GUARD
};

enum { MPD_CLAMP_DEFAULT, MPD_CLAMP_IEEE_754, MPD_CLAMP_GUARD };

extern const char *mpd_round_string[MPD_ROUND_GUARD];
extern const char *mpd_clamp_string[MPD_CLAMP_GUARD];


typedef struct mpd_context_t {
    mpd_ssize_t prec;   /* precision */
    mpd_ssize_t emax;   /* max positive exp */
    mpd_ssize_t emin;   /* min negative exp */
    uint32_t traps;     /* status events that should be trapped */
    uint32_t status;    /* status flags */
    uint32_t newtrap;   /* set by mpd_addstatus_raise() */
    int      round;     /* rounding mode */
    int      clamp;     /* clamp mode */
    int      allcr;     /* all functions correctly rounded */
} mpd_context_t;


/* Status flags */
#define MPD_Clamped             0x00000001U
#define MPD_Conversion_syntax   0x00000002U
#define MPD_Division_by_zero    0x00000004U
#define MPD_Division_impossible 0x00000008U
#define MPD_Division_undefined  0x00000010U
#define MPD_Fpu_error           0x00000020U
#define MPD_Inexact             0x00000040U
#define MPD_Invalid_context     0x00000080U
#define MPD_Invalid_operation   0x00000100U
#define MPD_Malloc_error        0x00000200U
#define MPD_Not_implemented     0x00000400U
#define MPD_Overflow            0x00000800U
#define MPD_Rounded             0x00001000U
#define MPD_Subnormal           0x00002000U
#define MPD_Underflow           0x00004000U
#define MPD_Max_status         (0x00008000U-1U)

/* Conditions that result in an IEEE 754 exception */
#define MPD_IEEE_Invalid_operation (MPD_Conversion_syntax |   \
                                    MPD_Division_impossible | \
                                    MPD_Division_undefined |  \
                                    MPD_Fpu_error |           \
                                    MPD_Invalid_context |     \
                                    MPD_Invalid_operation |   \
                                    MPD_Malloc_error)         \

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

/* Default traps */
#define MPD_Traps (MPD_IEEE_Invalid_operation | \
                   MPD_Division_by_zero |       \
                   MPD_Overflow |               \
                   MPD_Underflow)

/* Official name */
#define MPD_Insufficient_storage MPD_Malloc_error

/* IEEE 754 interchange format contexts */
#define MPD_IEEE_CONTEXT_MAX_BITS 512 /* 16*(log2(MPD_MAX_EMAX / 3)-3) */
#define MPD_DECIMAL32 32
#define MPD_DECIMAL64 64
#define MPD_DECIMAL128 128


#define MPD_MINALLOC_MIN 2
#define MPD_MINALLOC_MAX 64
extern mpd_ssize_t MPD_MINALLOC;
extern void (* mpd_traphandler)(mpd_context_t *);
void mpd_dflt_traphandler(mpd_context_t *);

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

void mpd_maxcontext(mpd_context_t *ctx);
void mpd_defaultcontext(mpd_context_t *ctx);
void mpd_basiccontext(mpd_context_t *ctx);
int mpd_ieee_context(mpd_context_t *ctx, int bits);

mpd_ssize_t mpd_getprec(const mpd_context_t *ctx);
mpd_ssize_t mpd_getemax(const mpd_context_t *ctx);
mpd_ssize_t mpd_getemin(const mpd_context_t *ctx);
int mpd_getround(const mpd_context_t *ctx);
uint32_t mpd_gettraps(const mpd_context_t *ctx);
uint32_t mpd_getstatus(const mpd_context_t *ctx);
int mpd_getclamp(const mpd_context_t *ctx);
int mpd_getcr(const mpd_context_t *ctx);

int mpd_qsetprec(mpd_context_t *ctx, mpd_ssize_t prec);
int mpd_qsetemax(mpd_context_t *ctx, mpd_ssize_t emax);
int mpd_qsetemin(mpd_context_t *ctx, mpd_ssize_t emin);
int mpd_qsetround(mpd_context_t *ctx, int newround);
int mpd_qsettraps(mpd_context_t *ctx, uint32_t flags);
int mpd_qsetstatus(mpd_context_t *ctx, uint32_t flags);
int mpd_qsetclamp(mpd_context_t *ctx, int c);
int mpd_qsetcr(mpd_context_t *ctx, int c);
void mpd_addstatus_raise(mpd_context_t *ctx, uint32_t flags);


/******************************************************************************/
/*                           Decimal Arithmetic                               */
/******************************************************************************/

/* mpd_t flags */
#define MPD_POS                 ((uint8_t)0)
#define MPD_NEG                 ((uint8_t)1)
#define MPD_INF                 ((uint8_t)2)
#define MPD_NAN                 ((uint8_t)4)
#define MPD_SNAN                ((uint8_t)8)
#define MPD_SPECIAL (MPD_INF|MPD_NAN|MPD_SNAN)
#define MPD_STATIC              ((uint8_t)16)
#define MPD_STATIC_DATA         ((uint8_t)32)
#define MPD_SHARED_DATA         ((uint8_t)64)
#define MPD_CONST_DATA          ((uint8_t)128)
#define MPD_DATAFLAGS (MPD_STATIC_DATA|MPD_SHARED_DATA|MPD_CONST_DATA)

/* mpd_t */
typedef struct mpd_t {
    uint8_t flags;
    mpd_ssize_t exp;
    mpd_ssize_t digits;
    mpd_ssize_t len;
    mpd_ssize_t alloc;
    mpd_uint_t *data;
} mpd_t;


typedef unsigned char uchar;


/******************************************************************************/
/*                       Quiet, thread-safe functions                         */
/******************************************************************************/

/* format specification */
typedef struct mpd_spec_t {
    mpd_ssize_t min_width; /* minimum field width */
    mpd_ssize_t prec;      /* fraction digits or significant digits */
    char type;             /* conversion specifier */
    char align;            /* alignment */
    char sign;             /* sign printing/alignment */
    char fill[5];          /* fill character */
    const char *dot;       /* decimal point */
    const char *sep;       /* thousands separator */
    const char *grouping;  /* grouping of digits */
} mpd_spec_t;

/* output to a string */
char *mpd_to_sci(const mpd_t *dec, int fmt);
char *mpd_to_eng(const mpd_t *dec, int fmt);
mpd_ssize_t mpd_to_sci_size(char **res, const mpd_t *dec, int fmt);
mpd_ssize_t mpd_to_eng_size(char **res, const mpd_t *dec, int fmt);
int mpd_validate_lconv(mpd_spec_t *spec);
int mpd_parse_fmt_str(mpd_spec_t *spec, const char *fmt, int caps);
char *mpd_qformat_spec(const mpd_t *dec, const mpd_spec_t *spec, const mpd_context_t *ctx, uint32_t *status);
char *mpd_qformat(const mpd_t *dec, const char *fmt, const mpd_context_t *ctx, uint32_t *status);

#define MPD_NUM_FLAGS 15
#define MPD_MAX_FLAG_STRING 208
#define MPD_MAX_FLAG_LIST (MPD_MAX_FLAG_STRING+18)
#define MPD_MAX_SIGNAL_LIST 121
int mpd_snprint_flags(char *dest, int nmemb, uint32_t flags);
int mpd_lsnprint_flags(char *dest, int nmemb, uint32_t flags, const char *flag_string[]);
int mpd_lsnprint_signals(char *dest, int nmemb, uint32_t flags, const char *signal_string[]);

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

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

/* set to NaN with error flags */
void mpd_seterror(mpd_t *result, uint32_t flags, uint32_t *status);
/* set a special with sign and type */
void mpd_setspecial(mpd_t *dec, uint8_t sign, uint8_t type);
/* set coefficient to zero or all nines */
void mpd_zerocoeff(mpd_t *result);
void mpd_qmaxcoeff(mpd_t *result, const mpd_context_t *ctx, uint32_t *status);

/* quietly assign a C integer type to an mpd_t */
void mpd_qset_ssize(mpd_t *result, mpd_ssize_t a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qset_i32(mpd_t *result, int32_t a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qset_uint(mpd_t *result, mpd_uint_t a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qset_u32(mpd_t *result, uint32_t a, const mpd_context_t *ctx, uint32_t *status);
#ifndef LEGACY_COMPILER
void mpd_qset_i64(mpd_t *result, int64_t a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qset_u64(mpd_t *result, uint64_t a, const mpd_context_t *ctx, uint32_t *status);
#endif

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

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

int32_t mpd_qget_i32(const mpd_t *dec, uint32_t *status);
uint32_t mpd_qget_u32(const mpd_t *dec, uint32_t *status);
#ifndef LEGACY_COMPILER
int64_t mpd_qget_i64(const mpd_t *dec, uint32_t *status);
uint64_t mpd_qget_u64(const mpd_t *dec, uint32_t *status);
#endif

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

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

int mpd_qcopy(mpd_t *result, const mpd_t *a,  uint32_t *status);
mpd_t *mpd_qncopy(const mpd_t *a);
int mpd_qcopy_abs(mpd_t *result, const mpd_t *a, uint32_t *status);
int mpd_qcopy_negate(mpd_t *result, const mpd_t *a, uint32_t *status);
int mpd_qcopy_sign(mpd_t *result, const mpd_t *a, const mpd_t *b, uint32_t *status);

void mpd_qand(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qinvert(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qlogb(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qor(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qscaleb(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qxor(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
int mpd_same_quantum(const mpd_t *a, const mpd_t *b);

void mpd_qrotate(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
int mpd_qshiftl(mpd_t *result, const mpd_t *a, mpd_ssize_t n, uint32_t *status);
mpd_uint_t mpd_qshiftr(mpd_t *result, const mpd_t *a, mpd_ssize_t n, uint32_t *status);
mpd_uint_t mpd_qshiftr_inplace(mpd_t *result, mpd_ssize_t n);
void mpd_qshift(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qshiftn(mpd_t *result, const mpd_t *a, mpd_ssize_t n, const mpd_context_t *ctx, uint32_t *status);

int mpd_qcmp(const mpd_t *a, const mpd_t *b, uint32_t *status);
int mpd_qcompare(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
int mpd_qcompare_signal(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
int mpd_cmp_total(const mpd_t *a, const mpd_t *b);
int mpd_cmp_total_mag(const mpd_t *a, const mpd_t *b);
int mpd_compare_total(mpd_t *result, const mpd_t *a, const mpd_t *b);
int mpd_compare_total_mag(mpd_t *result, const mpd_t *a, const mpd_t *b);

void mpd_qround_to_intx(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qround_to_int(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qtrunc(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qfloor(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qceil(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);

void mpd_qabs(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmax(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmax_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmin(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmin_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qminus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qplus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qnext_minus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qnext_plus(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qnext_toward(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qquantize(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qrescale(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, const mpd_context_t *ctx, uint32_t *status);
void mpd_qrescale_fmt(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, const mpd_context_t *ctx, uint32_t *status);
void mpd_qreduce(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qadd(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qadd_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qadd_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qadd_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qadd_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsub(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsub_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsub_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsub_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsub_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmul(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmul_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmul_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmul_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmul_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
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);
void mpd_qdiv(mpd_t *q, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qdiv_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qdiv_i32(mpd_t *result, const mpd_t *a, int32_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qdiv_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qdiv_u32(mpd_t *result, const mpd_t *a, uint32_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qdivint(mpd_t *q, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qrem(mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qrem_near(mpd_t *r, const mpd_t *a, const mpd_t *b, const mpd_context_t *ctx, uint32_t *status);
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);
void mpd_qpow(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_context_t *ctx, uint32_t *status);
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);
void mpd_qexp(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qln10(mpd_t *result, mpd_ssize_t prec, uint32_t *status);
void mpd_qln(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qlog10(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsqrt(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qinvroot(mpd_t *result, const mpd_t *a, const mpd_context_t *ctx, uint32_t *status);

#ifndef LEGACY_COMPILER
void mpd_qadd_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qadd_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsub_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsub_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmul_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qmul_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qdiv_i64(mpd_t *result, const mpd_t *a, int64_t b, const mpd_context_t *ctx, uint32_t *status);
void mpd_qdiv_u64(mpd_t *result, const mpd_t *a, uint64_t b, const mpd_context_t *ctx, uint32_t *status);
#endif


size_t mpd_sizeinbase(const mpd_t *a, uint32_t base);
void mpd_qimport_u16(mpd_t *result, const uint16_t *srcdata, size_t srclen,
                     uint8_t srcsign, uint32_t srcbase,
                     const mpd_context_t *ctx, uint32_t *status);
void mpd_qimport_u32(mpd_t *result, const uint32_t *srcdata, size_t srclen,
                     uint8_t srcsign, uint32_t srcbase,
                     const mpd_context_t *ctx, uint32_t *status);
size_t mpd_qexport_u16(uint16_t **rdata, size_t rlen, uint32_t base,
                       const mpd_t *src, uint32_t *status);
size_t mpd_qexport_u32(uint32_t **rdata, size_t rlen, uint32_t base,
                       const mpd_t *src, uint32_t *status);


/******************************************************************************/
/*                           Signalling functions                             */
/******************************************************************************/

char *mpd_format(const mpd_t *dec, const char *fmt, mpd_context_t *ctx);
void mpd_import_u16(mpd_t *result, const uint16_t *srcdata, size_t srclen, uint8_t srcsign, uint32_t base, mpd_context_t *ctx);
void mpd_import_u32(mpd_t *result, const uint32_t *srcdata, size_t srclen, uint8_t srcsign, uint32_t base, mpd_context_t *ctx);
size_t mpd_export_u16(uint16_t **rdata, size_t rlen, uint32_t base, const mpd_t *src, mpd_context_t *ctx);
size_t mpd_export_u32(uint32_t **rdata, size_t rlen, uint32_t base, const mpd_t *src, mpd_context_t *ctx);
void mpd_finalize(mpd_t *result, mpd_context_t *ctx);
int mpd_check_nan(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
int mpd_check_nans(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_set_string(mpd_t *result, const char *s, mpd_context_t *ctx);
void mpd_maxcoeff(mpd_t *result, mpd_context_t *ctx);
void mpd_sset_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx);
void mpd_sset_i32(mpd_t *result, int32_t a, mpd_context_t *ctx);
void mpd_sset_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx);
void mpd_sset_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx);
void mpd_set_ssize(mpd_t *result, mpd_ssize_t a, mpd_context_t *ctx);
void mpd_set_i32(mpd_t *result, int32_t a, mpd_context_t *ctx);
void mpd_set_uint(mpd_t *result, mpd_uint_t a, mpd_context_t *ctx);
void mpd_set_u32(mpd_t *result, uint32_t a, mpd_context_t *ctx);
#ifndef LEGACY_COMPILER
void mpd_set_i64(mpd_t *result, int64_t a, mpd_context_t *ctx);
void mpd_set_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx);
#endif
mpd_ssize_t mpd_get_ssize(const mpd_t *a, mpd_context_t *ctx);
mpd_uint_t mpd_get_uint(const mpd_t *a, mpd_context_t *ctx);
mpd_uint_t mpd_abs_uint(const mpd_t *a, mpd_context_t *ctx);
int32_t mpd_get_i32(const mpd_t *a, mpd_context_t *ctx);
uint32_t mpd_get_u32(const mpd_t *a, mpd_context_t *ctx);
#ifndef LEGACY_COMPILER
int64_t mpd_get_i64(const mpd_t *a, mpd_context_t *ctx);
uint64_t mpd_get_u64(const mpd_t *a, mpd_context_t *ctx);
#endif
void mpd_and(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_copy(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_canonical(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_copy_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_copy_negate(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_copy_sign(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_invert(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_logb(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_or(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_rotate(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_scaleb(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_shiftl(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx);
mpd_uint_t mpd_shiftr(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx);
void mpd_shiftn(mpd_t *result, const mpd_t *a, mpd_ssize_t n, mpd_context_t *ctx);
void mpd_shift(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_xor(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_abs(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
int mpd_cmp(const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
int mpd_compare(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
int mpd_compare_signal(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_add(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_add_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
void mpd_add_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
void mpd_add_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
void mpd_add_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
void mpd_sub(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_sub_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
void mpd_sub_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
void mpd_sub_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
void mpd_sub_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
void mpd_div(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_div_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
void mpd_div_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
void mpd_div_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
void mpd_div_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
void mpd_divmod(mpd_t *q, mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_divint(mpd_t *q, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_exp(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_fma(mpd_t *result, const mpd_t *a, const mpd_t *b, const mpd_t *c, mpd_context_t *ctx);
void mpd_ln(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_log10(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_max(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_max_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_min(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_min_mag(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_mul(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_mul_ssize(mpd_t *result, const mpd_t *a, mpd_ssize_t b, mpd_context_t *ctx);
void mpd_mul_i32(mpd_t *result, const mpd_t *a, int32_t b, mpd_context_t *ctx);
void mpd_mul_uint(mpd_t *result, const mpd_t *a, mpd_uint_t b, mpd_context_t *ctx);
void mpd_mul_u32(mpd_t *result, const mpd_t *a, uint32_t b, mpd_context_t *ctx);
void mpd_next_minus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_next_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_next_toward(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_plus(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_pow(mpd_t *result, const mpd_t *base, const mpd_t *exp, mpd_context_t *ctx);
void mpd_powmod(mpd_t *result, const mpd_t *base, const mpd_t *exp, const mpd_t *mod, mpd_context_t *ctx);
void mpd_quantize(mpd_t *result, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_rescale(mpd_t *result, const mpd_t *a, mpd_ssize_t exp, mpd_context_t *ctx);
void mpd_reduce(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_rem(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_rem_near(mpd_t *r, const mpd_t *a, const mpd_t *b, mpd_context_t *ctx);
void mpd_round_to_intx(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_round_to_int(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_trunc(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_floor(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_ceil(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_sqrt(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);
void mpd_invroot(mpd_t *result, const mpd_t *a, mpd_context_t *ctx);

#ifndef LEGACY_COMPILER
void mpd_add_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
void mpd_add_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
void mpd_sub_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
void mpd_sub_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
void mpd_div_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
void mpd_div_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
void mpd_mul_i64(mpd_t *result, const mpd_t *a, int64_t b, mpd_context_t *ctx);
void mpd_mul_u64(mpd_t *result, const mpd_t *a, uint64_t b, mpd_context_t *ctx);
#endif


/******************************************************************************/
/*                          Configuration specific                            */
/******************************************************************************/

#ifdef CONFIG_64
void mpd_qsset_i64(mpd_t *result, int64_t a, const mpd_context_t *ctx, uint32_t *status);
void mpd_qsset_u64(mpd_t *result, uint64_t a, const mpd_context_t *ctx, uint32_t *status);
void mpd_sset_i64(mpd_t *result, int64_t a, mpd_context_t *ctx);
void mpd_sset_u64(mpd_t *result, uint64_t a, mpd_context_t *ctx);
#endif


/******************************************************************************/
/*                       Get attributes of a decimal                          */
/******************************************************************************/

mpd_ssize_t mpd_adjexp(const mpd_t *dec);
mpd_ssize_t mpd_etiny(const mpd_context_t *ctx);
mpd_ssize_t mpd_etop(const mpd_context_t *ctx);
mpd_uint_t mpd_msword(const mpd_t *dec);
int mpd_word_digits(mpd_uint_t word);
/* most significant digit of a word */
mpd_uint_t mpd_msd(mpd_uint_t word);
/* least significant digit of a word */
mpd_uint_t mpd_lsd(mpd_uint_t word);
/* coefficient size needed to store 'digits' */
mpd_ssize_t mpd_digits_to_size(mpd_ssize_t digits);
/* number of digits in the exponent, undefined for MPD_SSIZE_MIN */
int mpd_exp_digits(mpd_ssize_t exp);
int mpd_iscanonical(const mpd_t *dec UNUSED);
int mpd_isfinite(const mpd_t *dec);
int mpd_isinfinite(const mpd_t *dec);
int mpd_isinteger(const mpd_t *dec);
int mpd_isnan(const mpd_t *dec);
int mpd_isnegative(const mpd_t *dec);
int mpd_ispositive(const mpd_t *dec);
int mpd_isqnan(const mpd_t *dec);
int mpd_issigned(const mpd_t *dec);
int mpd_issnan(const mpd_t *dec);
int mpd_isspecial(const mpd_t *dec);
int mpd_iszero(const mpd_t *dec);
/* undefined for special numbers */
int mpd_iszerocoeff(const mpd_t *dec);
int mpd_isnormal(const mpd_t *dec, const mpd_context_t *ctx);
int mpd_issubnormal(const mpd_t *dec, const mpd_context_t *ctx);
/* odd word */
int mpd_isoddword(mpd_uint_t word);
/* odd coefficient */
int mpd_isoddcoeff(const mpd_t *dec);
/* odd decimal, only defined for integers */
int mpd_isodd(const mpd_t *dec);
/* even decimal, only defined for integers */
int mpd_iseven(const mpd_t *dec);
/* 0 if dec is positive, 1 if dec is negative */
uint8_t mpd_sign(const mpd_t *dec);
/* 1 if dec is positive, -1 if dec is negative */
int mpd_arith_sign(const mpd_t *dec);
long mpd_radix(void);
int mpd_isdynamic(const mpd_t *dec);
int mpd_isstatic(const mpd_t *dec);
int mpd_isdynamic_data(const mpd_t *dec);
int mpd_isstatic_data(const mpd_t *dec);
int mpd_isshared_data(const mpd_t *dec);
int mpd_isconst_data(const mpd_t *dec);
mpd_ssize_t mpd_trail_zeros(const mpd_t *dec);


/******************************************************************************/
/*                       Set attributes of a decimal                          */
/******************************************************************************/

/* set number of decimal digits in the coefficient */
void mpd_setdigits(mpd_t *result);
void mpd_set_sign(mpd_t *result, uint8_t sign);
/* copy sign from another decimal */
void mpd_signcpy(mpd_t *result, const mpd_t *a);
void mpd_set_infinity(mpd_t *result);
void mpd_set_qnan(mpd_t *result);
void mpd_set_snan(mpd_t *result);
void mpd_set_negative(mpd_t *result);
void mpd_set_positive(mpd_t *result);
void mpd_set_dynamic(mpd_t *result);
void mpd_set_static(mpd_t *result);
void mpd_set_dynamic_data(mpd_t *result);
void mpd_set_static_data(mpd_t *result);
void mpd_set_shared_data(mpd_t *result);
void mpd_set_const_data(mpd_t *result);
void mpd_clear_flags(mpd_t *result);
void mpd_set_flags(mpd_t *result, uint8_t flags);
void mpd_copy_flags(mpd_t *result, const mpd_t *a);


/******************************************************************************/
/*                              Error Macros                                  */
/******************************************************************************/

#define mpd_err_fatal(...) \
    do {fprintf(stderr, "%s:%d: error: ", __FILE__, __LINE__); \
        fprintf(stderr, __VA_ARGS__);  fputc('\n', stderr);    \
        abort();                                               \
    } while (0)
#define mpd_err_warn(...) \
    do {fprintf(stderr, "%s:%d: warning: ", __FILE__, __LINE__); \
        fprintf(stderr, __VA_ARGS__); fputc('\n', stderr);       \
    } while (0)


/******************************************************************************/
/*                            Memory handling                                 */
/******************************************************************************/

extern void *(* const mpd_mallocfunc)(size_t size);
extern void *(* const mpd_callocfunc)(size_t nmemb, size_t size);
extern void *(* const mpd_reallocfunc)(void *ptr, size_t size);
extern void (* const mpd_free)(void *ptr);

void *mpd_callocfunc_em(size_t nmemb, size_t size);

void *mpd_alloc(mpd_size_t nmemb, mpd_size_t size);
void *mpd_calloc(mpd_size_t nmemb, mpd_size_t size);
void *mpd_realloc(void *ptr, mpd_size_t nmemb, mpd_size_t size, uint8_t *err);
void *mpd_sh_alloc(mpd_size_t struct_size, mpd_size_t nmemb, mpd_size_t size);

mpd_t *mpd_qnew(void);
mpd_t *mpd_new(mpd_context_t *ctx);
mpd_t *mpd_qnew_size(mpd_ssize_t size);
void mpd_del(mpd_t *dec);

void mpd_uint_zero(mpd_uint_t *dest, mpd_size_t len);
int mpd_qresize(mpd_t *result, mpd_ssize_t size, uint32_t *status);
int mpd_qresize_zero(mpd_t *result, mpd_ssize_t size, uint32_t *status);
void mpd_minalloc(mpd_t *result);

int mpd_resize(mpd_t *result, mpd_ssize_t size, mpd_context_t *ctx);
int mpd_resize_zero(mpd_t *result, mpd_ssize_t size, mpd_context_t *ctx);


#ifdef __cplusplus
  #ifdef MPD_CLEAR_STDC_LIMIT_MACROS
    #undef MPD_CLEAR_STDC_LIMIT_MACROS
    #undef __STDC_LIMIT_MACROS
  #endif
} /* END extern "C" */
#endif


#endif /* MPDECIMAL_H */