diff options
Diffstat (limited to 'mpc/src/mpc-impl.h')
| -rw-r--r-- | mpc/src/mpc-impl.h | 194 |
1 files changed, 194 insertions, 0 deletions
diff --git a/mpc/src/mpc-impl.h b/mpc/src/mpc-impl.h new file mode 100644 index 0000000000..b2aaa90eb3 --- /dev/null +++ b/mpc/src/mpc-impl.h @@ -0,0 +1,194 @@ +/* mpc-impl.h -- Internal include file for mpc. + +Copyright (C) 2002, 2004, 2005, 2008, 2009, 2010, 2011, 2012 INRIA + +This file is part of GNU MPC. + +GNU MPC is free software; you can redistribute it and/or modify it under +the terms of the GNU Lesser General Public License as published by the +Free Software Foundation; either version 3 of the License, or (at your +option) any later version. + +GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY +WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS +FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for +more details. + +You should have received a copy of the GNU Lesser General Public License +along with this program. If not, see http://www.gnu.org/licenses/ . +*/ + +#ifndef __MPC_IMPL_H +#define __MPC_IMPL_H +#define __MPC_LIBRARY_BUILD + /* to indicate we are inside the library build */ + +#include "config.h" +#ifdef HAVE_STDLIB_H +#include <stdlib.h> +#endif +#include "mpc.h" + +/* + * Miscellaneous useful macros + */ + +#define MPC_MIN(h,i) ((h) < (i) ? (h) : (i)) +#define MPC_MAX(h,i) ((h) > (i) ? (h) : (i)) + +/* Safe absolute value (to avoid possible integer overflow) */ +/* type is the target (unsigned) type (copied from mpfr-impl.h) */ +#ifdef SAFE_ABS +#undef SAFE_ABS +#endif +#define SAFE_ABS(type,x) ((x) >= 0 ? (type)(x) : -(type)(x)) + + +/* + * MPFR constants and macros + */ + +#ifndef BITS_PER_MP_LIMB +#define BITS_PER_MP_LIMB mp_bits_per_limb +#endif + +#define MPFR_SIGNBIT(x) (mpfr_signbit (x) ? -1 : 1) +#define MPC_MPFR_SIGN(x) (mpfr_zero_p (x) ? 0 : MPFR_SIGNBIT (x)) + /* should be called MPFR_SIGN, but this is taken in mpfr.h */ +#define MPFR_CHANGE_SIGN(x) mpfr_neg(x,x,GMP_RNDN) +#define MPFR_COPYSIGN(x,y,z,rnd) (mpfr_nan_p (z) ? \ + mpfr_setsign (x, y, 0, rnd) : \ + mpfr_copysign (x, y, z, rnd)) + /* work around spurious signs in nan */ +#define MPFR_ADD_ONE_ULP(x) mpfr_add_one_ulp (x, GMP_RNDN) +#define MPFR_SUB_ONE_ULP(x) mpfr_sub_one_ulp (x, GMP_RNDN) + /* drop unused rounding mode from macroes */ +#define MPFR_SWAP(a,b) do { mpfr_srcptr tmp; tmp = a; a = b; b = tmp; } while (0) + + +/* + * Macro implementing rounding away from zero, to ease compatibility with + * mpfr < 3. f is the complete function call with a rounding mode of + * MPFR_RNDA, rop the name of the variable containing the result; it is + * already contained in f, but needs to be repeated so that the macro can + * modify the variable. + * Usage: replace each call to a function such as + * mpfr_add (rop, a, b, MPFR_RNDA) + * by + * ROUND_AWAY (mpfr_add (rop, a, b, MPFR_RNDA), rop) +*/ +#if MPFR_VERSION_MAJOR < 3 + /* round towards zero, add 1 ulp if not exact */ +#define MPFR_RNDA GMP_RNDZ +#define ROUND_AWAY(f,rop) \ + ((f) ? MPFR_ADD_ONE_ULP (rop), MPFR_SIGNBIT (rop) : 0) +#else +#define ROUND_AWAY(f,rop) \ + (f) +#endif /* mpfr < 3 */ + +#if MPFR_VERSION_MAJOR < 3 +/* declare missing functions, defined in get_version.c */ +__MPC_DECLSPEC void mpfr_set_zero (mpfr_ptr, int); +__MPC_DECLSPEC int mpfr_regular_p (mpfr_srcptr); +#endif /* mpfr < 3 */ + + +/* + * MPC macros + */ + +#define MPC_PREC_RE(x) (mpfr_get_prec(mpc_realref(x))) +#define MPC_PREC_IM(x) (mpfr_get_prec(mpc_imagref(x))) +#define MPC_MAX_PREC(x) MPC_MAX(MPC_PREC_RE(x), MPC_PREC_IM(x)) + +#define INV_RND(r) \ + (((r) == GMP_RNDU) ? GMP_RNDD : (((r) == GMP_RNDD) ? GMP_RNDU : (r))) + +#define mpc_inf_p(z) (mpfr_inf_p(mpc_realref(z))||mpfr_inf_p(mpc_imagref(z))) + /* Convention in C99 (G.3): z is regarded as an infinity if at least one of + its parts is infinite */ +#define mpc_zero_p(z) (mpfr_zero_p(mpc_realref(z))&&mpfr_zero_p(mpc_imagref(z))) + /* Convention in C99 (G.3): z is regarded as a zero if each of its parts is + a zero */ +#define mpc_fin_p(z) (mpfr_number_p(mpc_realref(z))&&mpfr_number_p(mpc_imagref(z))) + /* Convention in C99 (G.3): z is regarded as finite if both its parts are */ +#define mpc_nan_p(z) ((mpfr_nan_p(mpc_realref(z)) && !mpfr_inf_p(mpc_imagref(z))) || (mpfr_nan_p(mpc_imagref(z)) && !mpfr_inf_p(mpc_realref(z)))) + /* Consider as NaN all other numbers containing at least one NaN */ + + +/* + * ASSERT macros + */ + +#ifdef NDEBUG +#define MPC_ASSERT(expr) \ + do { \ + } while (0) +#else +#define MPC_ASSERT(expr) \ + do { \ + if (!(expr)) \ + { \ + fprintf (stderr, "%s:%d: MPC assertion failed: %s\n", \ + __FILE__, __LINE__, #expr); \ + abort(); \ + } \ + } while (0) +#endif + + +/* + * Debug macros + */ + +#define MPC_OUT(x) \ +do { \ + printf (#x "[%lu,%lu]=", (unsigned long int) MPC_PREC_RE (x), \ + (unsigned long int) MPC_PREC_IM (x)); \ + mpc_out_str (stdout, 2, 0, x, MPC_RNDNN); \ + printf ("\n"); \ +} while (0) + +#define MPFR_OUT(x) \ +do { \ + printf (#x "[%lu]=", (unsigned long int) mpfr_get_prec (x)); \ + mpfr_out_str (stdout, 2, 0, x, GMP_RNDN); \ + printf ("\n"); \ +} while (0) + + +/* + * Constants + */ + +#ifndef MUL_KARATSUBA_THRESHOLD +#define MUL_KARATSUBA_THRESHOLD 23 +#endif + + +/* + * Define internal functions + */ + +#if defined (__cplusplus) +extern "C" { +#endif + + +__MPC_DECLSPEC int mpc_mul_naive (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t); +__MPC_DECLSPEC int mpc_mul_karatsuba (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t); +__MPC_DECLSPEC int mpc_fma_naive (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t); +__MPC_DECLSPEC int mpc_pow_usi (mpc_ptr, mpc_srcptr, unsigned long, int, mpc_rnd_t); +__MPC_DECLSPEC char* mpc_alloc_str (size_t); +__MPC_DECLSPEC char* mpc_realloc_str (char*, size_t, size_t); +__MPC_DECLSPEC void mpc_free_str (char*); +__MPC_DECLSPEC mpfr_prec_t mpc_ceil_log2 (mpfr_prec_t); +__MPC_DECLSPEC int set_pi_over_2 (mpfr_ptr, int, mpfr_rnd_t); + +#if defined (__cplusplus) +} +#endif + + +#endif |