xref: /dflybsd-src/contrib/mpc/src/set_x.c (revision 04febcfb30580676d3e95f58a16c5137ee478b32)
1*d30dc8cbSJohn Marino /* mpc_set_x -- Set the real part of a complex number
2*d30dc8cbSJohn Marino    (imaginary part equals +0 regardless of rounding mode).
3*d30dc8cbSJohn Marino 
4*d30dc8cbSJohn Marino Copyright (C) 2008, 2009, 2010, 2011 INRIA
5*d30dc8cbSJohn Marino 
6*d30dc8cbSJohn Marino This file is part of GNU MPC.
7*d30dc8cbSJohn Marino 
8*d30dc8cbSJohn Marino GNU MPC is free software; you can redistribute it and/or modify it under
9*d30dc8cbSJohn Marino the terms of the GNU Lesser General Public License as published by the
10*d30dc8cbSJohn Marino Free Software Foundation; either version 3 of the License, or (at your
11*d30dc8cbSJohn Marino option) any later version.
12*d30dc8cbSJohn Marino 
13*d30dc8cbSJohn Marino GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY
14*d30dc8cbSJohn Marino WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
15*d30dc8cbSJohn Marino FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
16*d30dc8cbSJohn Marino more details.
17*d30dc8cbSJohn Marino 
18*d30dc8cbSJohn Marino You should have received a copy of the GNU Lesser General Public License
19*d30dc8cbSJohn Marino along with this program. If not, see http://www.gnu.org/licenses/ .
20*d30dc8cbSJohn Marino */
21*d30dc8cbSJohn Marino 
22*d30dc8cbSJohn Marino #include "config.h"
23*d30dc8cbSJohn Marino 
24*d30dc8cbSJohn Marino #ifdef HAVE_INTTYPES_H
25*d30dc8cbSJohn Marino # include <inttypes.h> /* for intmax_t */
26*d30dc8cbSJohn Marino #else
27*d30dc8cbSJohn Marino # ifdef HAVE_STDINT_H
28*d30dc8cbSJohn Marino #  include <stdint.h>
29*d30dc8cbSJohn Marino # endif
30*d30dc8cbSJohn Marino #endif
31*d30dc8cbSJohn Marino 
32*d30dc8cbSJohn Marino #ifdef HAVE_COMPLEX_H
33*d30dc8cbSJohn Marino # include <complex.h>
34*d30dc8cbSJohn Marino #endif
35*d30dc8cbSJohn Marino 
36*d30dc8cbSJohn Marino #include "mpc-impl.h"
37*d30dc8cbSJohn Marino 
38*d30dc8cbSJohn Marino #define MPC_SET_X(real_t, z, real_value, rnd)     \
39*d30dc8cbSJohn Marino   {                                                                     \
40*d30dc8cbSJohn Marino     int _inex_re, _inex_im;                                             \
41*d30dc8cbSJohn Marino     _inex_re = (mpfr_set_ ## real_t) (mpc_realref (z), (real_value), MPC_RND_RE (rnd)); \
42*d30dc8cbSJohn Marino     _inex_im = mpfr_set_ui (mpc_imagref (z), 0, MPC_RND_IM (rnd)); \
43*d30dc8cbSJohn Marino     return MPC_INEX (_inex_re, _inex_im);                               \
44*d30dc8cbSJohn Marino   }
45*d30dc8cbSJohn Marino 
46*d30dc8cbSJohn Marino int
mpc_set_fr(mpc_ptr a,mpfr_srcptr b,mpc_rnd_t rnd)47*d30dc8cbSJohn Marino mpc_set_fr (mpc_ptr a, mpfr_srcptr b, mpc_rnd_t rnd)
48*d30dc8cbSJohn Marino    MPC_SET_X (fr, a, b, rnd)
49*d30dc8cbSJohn Marino 
50*d30dc8cbSJohn Marino int
51*d30dc8cbSJohn Marino mpc_set_d (mpc_ptr a, double b, mpc_rnd_t rnd)
52*d30dc8cbSJohn Marino    MPC_SET_X (d, a, b, rnd)
53*d30dc8cbSJohn Marino 
54*d30dc8cbSJohn Marino int
55*d30dc8cbSJohn Marino mpc_set_ld (mpc_ptr a, long double b, mpc_rnd_t rnd)
56*d30dc8cbSJohn Marino    MPC_SET_X (ld, a, b, rnd)
57*d30dc8cbSJohn Marino 
58*d30dc8cbSJohn Marino int
59*d30dc8cbSJohn Marino mpc_set_ui (mpc_ptr a, unsigned long int b, mpc_rnd_t rnd)
60*d30dc8cbSJohn Marino    MPC_SET_X (ui, a, b, rnd)
61*d30dc8cbSJohn Marino 
62*d30dc8cbSJohn Marino int
63*d30dc8cbSJohn Marino mpc_set_si (mpc_ptr a, long int b, mpc_rnd_t rnd)
64*d30dc8cbSJohn Marino    MPC_SET_X (si, a, b, rnd)
65*d30dc8cbSJohn Marino 
66*d30dc8cbSJohn Marino int
67*d30dc8cbSJohn Marino mpc_set_z (mpc_ptr a, mpz_srcptr b, mpc_rnd_t rnd)
68*d30dc8cbSJohn Marino    MPC_SET_X (z, a, b, rnd)
69*d30dc8cbSJohn Marino 
70*d30dc8cbSJohn Marino int
71*d30dc8cbSJohn Marino mpc_set_q (mpc_ptr a, mpq_srcptr b, mpc_rnd_t rnd)
72*d30dc8cbSJohn Marino    MPC_SET_X (q, a, b, rnd)
73*d30dc8cbSJohn Marino 
74*d30dc8cbSJohn Marino int
75*d30dc8cbSJohn Marino mpc_set_f (mpc_ptr a, mpf_srcptr b, mpc_rnd_t rnd)
76*d30dc8cbSJohn Marino    MPC_SET_X (f, a, b, rnd)
77*d30dc8cbSJohn Marino 
78*d30dc8cbSJohn Marino #ifdef _MPC_H_HAVE_INTMAX_T
79*d30dc8cbSJohn Marino int
80*d30dc8cbSJohn Marino mpc_set_uj (mpc_ptr a, uintmax_t b, mpc_rnd_t rnd)
81*d30dc8cbSJohn Marino    MPC_SET_X (uj, a, b, rnd)
82*d30dc8cbSJohn Marino 
83*d30dc8cbSJohn Marino int
84*d30dc8cbSJohn Marino mpc_set_sj (mpc_ptr a, intmax_t b, mpc_rnd_t rnd)
85*d30dc8cbSJohn Marino    MPC_SET_X (sj, a, b, rnd)
86*d30dc8cbSJohn Marino #endif
87*d30dc8cbSJohn Marino 
88*d30dc8cbSJohn Marino #ifdef HAVE_COMPLEX_H
89*d30dc8cbSJohn Marino int
90*d30dc8cbSJohn Marino mpc_set_dc (mpc_ptr a, double _Complex b, mpc_rnd_t rnd) {
91*d30dc8cbSJohn Marino    return mpc_set_d_d (a, creal (b), cimag (b), rnd);
92*d30dc8cbSJohn Marino }
93*d30dc8cbSJohn Marino 
94*d30dc8cbSJohn Marino int
mpc_set_ldc(mpc_ptr a,long double _Complex b,mpc_rnd_t rnd)95*d30dc8cbSJohn Marino mpc_set_ldc (mpc_ptr a, long double _Complex b, mpc_rnd_t rnd) {
96*d30dc8cbSJohn Marino    return mpc_set_ld_ld (a, creall (b), cimagl (b), rnd);
97*d30dc8cbSJohn Marino }
98*d30dc8cbSJohn Marino #endif
99*d30dc8cbSJohn Marino 
100*d30dc8cbSJohn Marino void
mpc_set_nan(mpc_ptr a)101*d30dc8cbSJohn Marino mpc_set_nan (mpc_ptr a) {
102*d30dc8cbSJohn Marino    mpfr_set_nan (mpc_realref (a));
103*d30dc8cbSJohn Marino    mpfr_set_nan (mpc_imagref (a));
104*d30dc8cbSJohn Marino }
105