1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
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));
}
|
