diff options
Diffstat (limited to 'mpc/src/set_x.c')
-rw-r--r-- | mpc/src/set_x.c | 104 |
1 files changed, 104 insertions, 0 deletions
diff --git a/mpc/src/set_x.c b/mpc/src/set_x.c new file mode 100644 index 0000000000..94ec12d96b --- /dev/null +++ b/mpc/src/set_x.c @@ -0,0 +1,104 @@ +/* mpc_set_x -- Set the real part of a complex number + (imaginary part equals +0 regardless of rounding mode). + +Copyright (C) 2008, 2009, 2010, 2011 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/ . +*/ + +#include "config.h" + +#ifdef HAVE_INTTYPES_H +# include <inttypes.h> /* for intmax_t */ +#else +# ifdef HAVE_STDINT_H +# include <stdint.h> +# endif +#endif + +#ifdef HAVE_COMPLEX_H +# include <complex.h> +#endif + +#include "mpc-impl.h" + +#define MPC_SET_X(real_t, z, real_value, rnd) \ + { \ + int _inex_re, _inex_im; \ + _inex_re = (mpfr_set_ ## real_t) (mpc_realref (z), (real_value), MPC_RND_RE (rnd)); \ + _inex_im = mpfr_set_ui (mpc_imagref (z), 0, MPC_RND_IM (rnd)); \ + return MPC_INEX (_inex_re, _inex_im); \ + } + +int +mpc_set_fr (mpc_ptr a, mpfr_srcptr b, mpc_rnd_t rnd) + MPC_SET_X (fr, a, b, rnd) + +int +mpc_set_d (mpc_ptr a, double b, mpc_rnd_t rnd) + MPC_SET_X (d, a, b, rnd) + +int +mpc_set_ld (mpc_ptr a, long double b, mpc_rnd_t rnd) + MPC_SET_X (ld, a, b, rnd) + +int +mpc_set_ui (mpc_ptr a, unsigned long int b, mpc_rnd_t rnd) + MPC_SET_X (ui, a, b, rnd) + +int +mpc_set_si (mpc_ptr a, long int b, mpc_rnd_t rnd) + MPC_SET_X (si, a, b, rnd) + +int +mpc_set_z (mpc_ptr a, mpz_srcptr b, mpc_rnd_t rnd) + MPC_SET_X (z, a, b, rnd) + +int +mpc_set_q (mpc_ptr a, mpq_srcptr b, mpc_rnd_t rnd) + MPC_SET_X (q, a, b, rnd) + +int +mpc_set_f (mpc_ptr a, mpf_srcptr b, mpc_rnd_t rnd) + MPC_SET_X (f, a, b, rnd) + +#ifdef _MPC_H_HAVE_INTMAX_T +int +mpc_set_uj (mpc_ptr a, uintmax_t b, mpc_rnd_t rnd) + MPC_SET_X (uj, a, b, rnd) + +int +mpc_set_sj (mpc_ptr a, intmax_t b, mpc_rnd_t rnd) + MPC_SET_X (sj, a, b, rnd) +#endif + +#ifdef HAVE_COMPLEX_H +int +mpc_set_dc (mpc_ptr a, double _Complex b, mpc_rnd_t rnd) { + return mpc_set_d_d (a, creal (b), cimag (b), rnd); +} + +int +mpc_set_ldc (mpc_ptr a, long double _Complex b, mpc_rnd_t rnd) { + return mpc_set_ld_ld (a, creall (b), cimagl (b), rnd); +} +#endif + +void +mpc_set_nan (mpc_ptr a) { + mpfr_set_nan (mpc_realref (a)); + mpfr_set_nan (mpc_imagref (a)); +} |