1/* mpc.h -- Include file for mpc.
2
3Copyright (C) 2002, 2003, 2004, 2005, 2007, 2008, 2009, 2010, 2011, 2012, 2016, 2017, 2018, 2020 INRIA
4
5This file is part of GNU MPC.
6
7GNU MPC is free software; you can redistribute it and/or modify it under
8the terms of the GNU Lesser General Public License as published by the
9Free Software Foundation; either version 3 of the License, or (at your
10option) any later version.
11
12GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY
13WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
14FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
15more details.
16
17You should have received a copy of the GNU Lesser General Public License
18along with this program. If not, see http://www.gnu.org/licenses/ .
19*/
20
21#ifndef __MPC_H
22#define __MPC_H
23
24#include "gmp.h"
25#include "mpfr.h"
26
27/* Define MPC version number */
28#define MPC_VERSION_MAJOR 1
29#define MPC_VERSION_MINOR 2
30#define MPC_VERSION_PATCHLEVEL 1
31#define MPC_VERSION_STRING "1.2.1"
32
33/* Macros dealing with MPC VERSION */
34#define MPC_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c))
35#define MPC_VERSION \
36 MPC_VERSION_NUM(MPC_VERSION_MAJOR,MPC_VERSION_MINOR,MPC_VERSION_PATCHLEVEL)
37
38/* Check if stdint.h/inttypes.h is included */
39#if defined (INTMAX_C) && defined (UINTMAX_C)
40#define _MPC_H_HAVE_INTMAX_T 1
41#endif
42
43/* Return values */
44
45/* Transform negative to 2, positive to 1, leave 0 unchanged.
46 Warning: since inex is evaluated two times, we should avoid
47 MPC_INEX(mpc_mul (...), mpc_mul (...)) */
48#define MPC_INEX_POS(inex) (((inex) < 0) ? 2 : ((inex) == 0) ? 0 : 1)
49/* Transform 2 to negative, 1 to positive, leave 0 unchanged */
50#define MPC_INEX_NEG(inex) (((inex) == 2) ? -1 : ((inex) == 0) ? 0 : 1)
51
52/* The global inexact flag is made of (real flag) + 4 * (imaginary flag), where
53 each of the real and imaginary inexact flag are:
54 0 when the result is exact (no rounding error)
55 1 when the result is larger than the exact value
56 2 when the result is smaller than the exact value */
57#define MPC_INEX(inex_re, inex_im) \
58 (MPC_INEX_POS(inex_re) | (MPC_INEX_POS(inex_im) << 2))
59#define MPC_INEX_RE(inex) MPC_INEX_NEG((inex) & 3)
60#define MPC_INEX_IM(inex) MPC_INEX_NEG((inex) >> 2)
61
62/* For functions computing two results, the return value is
63 inexact1+16*inexact2, which is 0 iif both results are exact. */
64#define MPC_INEX12(inex1, inex2) (inex1 | (inex2 << 4))
65#define MPC_INEX1(inex) (inex & 15)
66#define MPC_INEX2(inex) (inex >> 4)
67
68/* Definition of rounding modes */
69
70/* a complex rounding mode is just a pair of two real rounding modes
71 we reserve four bits for a real rounding mode. */
72typedef int mpc_rnd_t;
73
74#define MPC_RND(r1,r2) (((int)(r1)) + ((int)(r2) << 4))
75#define MPC_RND_RE(x) ((mpfr_rnd_t)((x) & 0x0F))
76#define MPC_RND_IM(x) ((mpfr_rnd_t)((x) >> 4))
77
78#define MPC_RNDNN MPC_RND (MPFR_RNDN,MPFR_RNDN)
79#define MPC_RNDNZ MPC_RND (MPFR_RNDN,MPFR_RNDZ)
80#define MPC_RNDNU MPC_RND (MPFR_RNDN,MPFR_RNDU)
81#define MPC_RNDND MPC_RND (MPFR_RNDN,MPFR_RNDD)
82
83#define MPC_RNDZN MPC_RND (MPFR_RNDZ,MPFR_RNDN)
84#define MPC_RNDZZ MPC_RND (MPFR_RNDZ,MPFR_RNDZ)
85#define MPC_RNDZU MPC_RND (MPFR_RNDZ,MPFR_RNDU)
86#define MPC_RNDZD MPC_RND (MPFR_RNDZ,MPFR_RNDD)
87
88#define MPC_RNDUN MPC_RND (MPFR_RNDU,MPFR_RNDN)
89#define MPC_RNDUZ MPC_RND (MPFR_RNDU,MPFR_RNDZ)
90#define MPC_RNDUU MPC_RND (MPFR_RNDU,MPFR_RNDU)
91#define MPC_RNDUD MPC_RND (MPFR_RNDU,MPFR_RNDD)
92
93#define MPC_RNDDN MPC_RND (MPFR_RNDD,MPFR_RNDN)
94#define MPC_RNDDZ MPC_RND (MPFR_RNDD,MPFR_RNDZ)
95#define MPC_RNDDU MPC_RND (MPFR_RNDD,MPFR_RNDU)
96#define MPC_RNDDD MPC_RND (MPFR_RNDD,MPFR_RNDD)
97
98
99/* Definitions of types and their semantics */
100
101typedef struct {
102 mpfr_t re;
103 mpfr_t im;
104}
105__mpc_struct;
106
107typedef __mpc_struct mpc_t[1];
108typedef __mpc_struct *mpc_ptr;
109typedef const __mpc_struct *mpc_srcptr;
110
111/* Support for WINDOWS DLL, see
112 http://lists.gforge.inria.fr/pipermail/mpc-discuss/2011-November/000990.html;
113 when building the DLL, export symbols, otherwise behave as GMP */
114#if defined (__MPC_LIBRARY_BUILD) && __GMP_LIBGMP_DLL
115#define __MPC_DECLSPEC __GMP_DECLSPEC_EXPORT
116#else
117#define __MPC_DECLSPEC __GMP_DECLSPEC
118#endif
119
120#if defined (__cplusplus)
121extern "C" {
122#endif
123
124__MPC_DECLSPEC int mpc_add (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t);
125__MPC_DECLSPEC int mpc_add_fr (mpc_ptr, mpc_srcptr, mpfr_srcptr, mpc_rnd_t);
126__MPC_DECLSPEC int mpc_add_si (mpc_ptr, mpc_srcptr, long int, mpc_rnd_t);
127__MPC_DECLSPEC int mpc_add_ui (mpc_ptr, mpc_srcptr, unsigned long int, mpc_rnd_t);
128__MPC_DECLSPEC int mpc_sub (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t);
129__MPC_DECLSPEC int mpc_sub_fr (mpc_ptr, mpc_srcptr, mpfr_srcptr, mpc_rnd_t);
130__MPC_DECLSPEC int mpc_fr_sub (mpc_ptr, mpfr_srcptr, mpc_srcptr, mpc_rnd_t);
131__MPC_DECLSPEC int mpc_sub_ui (mpc_ptr, mpc_srcptr, unsigned long int, mpc_rnd_t);
132__MPC_DECLSPEC int mpc_ui_ui_sub (mpc_ptr, unsigned long int, unsigned long int, mpc_srcptr, mpc_rnd_t);
133__MPC_DECLSPEC int mpc_mul (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t);
134__MPC_DECLSPEC int mpc_mul_fr (mpc_ptr, mpc_srcptr, mpfr_srcptr, mpc_rnd_t);
135__MPC_DECLSPEC int mpc_mul_ui (mpc_ptr, mpc_srcptr, unsigned long int, mpc_rnd_t);
136__MPC_DECLSPEC int mpc_mul_si (mpc_ptr, mpc_srcptr, long int, mpc_rnd_t);
137__MPC_DECLSPEC int mpc_mul_i (mpc_ptr, mpc_srcptr, int, mpc_rnd_t);
138__MPC_DECLSPEC int mpc_sqr (mpc_ptr, mpc_srcptr, mpc_rnd_t);
139__MPC_DECLSPEC int mpc_div (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t);
140__MPC_DECLSPEC int mpc_pow (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t);
141__MPC_DECLSPEC int mpc_pow_fr (mpc_ptr, mpc_srcptr, mpfr_srcptr, mpc_rnd_t);
142__MPC_DECLSPEC int mpc_pow_ld (mpc_ptr, mpc_srcptr, long double, mpc_rnd_t);
143__MPC_DECLSPEC int mpc_pow_d (mpc_ptr, mpc_srcptr, double, mpc_rnd_t);
144__MPC_DECLSPEC int mpc_pow_si (mpc_ptr, mpc_srcptr, long, mpc_rnd_t);
145__MPC_DECLSPEC int mpc_pow_ui (mpc_ptr, mpc_srcptr, unsigned long, mpc_rnd_t);
146__MPC_DECLSPEC int mpc_pow_z (mpc_ptr, mpc_srcptr, mpz_srcptr, mpc_rnd_t);
147__MPC_DECLSPEC int mpc_div_fr (mpc_ptr, mpc_srcptr, mpfr_srcptr, mpc_rnd_t);
148__MPC_DECLSPEC int mpc_fr_div (mpc_ptr, mpfr_srcptr, mpc_srcptr, mpc_rnd_t);
149__MPC_DECLSPEC int mpc_div_ui (mpc_ptr, mpc_srcptr, unsigned long int, mpc_rnd_t);
150__MPC_DECLSPEC int mpc_ui_div (mpc_ptr, unsigned long int, mpc_srcptr, mpc_rnd_t);
151__MPC_DECLSPEC int mpc_div_2ui (mpc_ptr, mpc_srcptr, unsigned long int, mpc_rnd_t);
152__MPC_DECLSPEC int mpc_mul_2ui (mpc_ptr, mpc_srcptr, unsigned long int, mpc_rnd_t);
153__MPC_DECLSPEC int mpc_div_2si (mpc_ptr, mpc_srcptr, long int, mpc_rnd_t);
154__MPC_DECLSPEC int mpc_mul_2si (mpc_ptr, mpc_srcptr, long int, mpc_rnd_t);
155__MPC_DECLSPEC int mpc_conj (mpc_ptr, mpc_srcptr, mpc_rnd_t);
156__MPC_DECLSPEC int mpc_neg (mpc_ptr, mpc_srcptr, mpc_rnd_t);
157__MPC_DECLSPEC int mpc_sum (mpc_ptr, const mpc_ptr *, unsigned long, mpc_rnd_t);
158__MPC_DECLSPEC int mpc_dot (mpc_ptr, const mpc_ptr *, const mpc_ptr *, unsigned long, mpc_rnd_t);
159__MPC_DECLSPEC int mpc_norm (mpfr_ptr, mpc_srcptr, mpfr_rnd_t);
160__MPC_DECLSPEC int mpc_abs (mpfr_ptr, mpc_srcptr, mpfr_rnd_t);
161__MPC_DECLSPEC int mpc_sqrt (mpc_ptr, mpc_srcptr, mpc_rnd_t);
162__MPC_DECLSPEC int mpc_set (mpc_ptr, mpc_srcptr, mpc_rnd_t);
163__MPC_DECLSPEC int mpc_set_d (mpc_ptr, double, mpc_rnd_t);
164__MPC_DECLSPEC int mpc_set_d_d (mpc_ptr, double, double, mpc_rnd_t);
165__MPC_DECLSPEC int mpc_set_ld (mpc_ptr, long double, mpc_rnd_t);
166__MPC_DECLSPEC int mpc_set_ld_ld (mpc_ptr, long double, long double, mpc_rnd_t);
167__MPC_DECLSPEC int mpc_set_f (mpc_ptr, mpf_srcptr, mpc_rnd_t);
168__MPC_DECLSPEC int mpc_set_f_f (mpc_ptr, mpf_srcptr, mpf_srcptr, mpc_rnd_t);
169__MPC_DECLSPEC int mpc_set_fr (mpc_ptr, mpfr_srcptr, mpc_rnd_t);
170__MPC_DECLSPEC int mpc_set_fr_fr (mpc_ptr, mpfr_srcptr, mpfr_srcptr, mpc_rnd_t);
171__MPC_DECLSPEC int mpc_set_q (mpc_ptr, mpq_srcptr, mpc_rnd_t);
172__MPC_DECLSPEC int mpc_set_q_q (mpc_ptr, mpq_srcptr, mpq_srcptr, mpc_rnd_t);
173__MPC_DECLSPEC int mpc_set_si (mpc_ptr, long int, mpc_rnd_t);
174__MPC_DECLSPEC int mpc_set_si_si (mpc_ptr, long int, long int, mpc_rnd_t);
175__MPC_DECLSPEC int mpc_set_ui (mpc_ptr, unsigned long int, mpc_rnd_t);
176__MPC_DECLSPEC int mpc_set_ui_ui (mpc_ptr, unsigned long int, unsigned long int, mpc_rnd_t);
177__MPC_DECLSPEC int mpc_set_z (mpc_ptr, mpz_srcptr, mpc_rnd_t);
178__MPC_DECLSPEC int mpc_set_z_z (mpc_ptr, mpz_srcptr, mpz_srcptr, mpc_rnd_t);
179__MPC_DECLSPEC void mpc_swap (mpc_ptr, mpc_ptr);
180__MPC_DECLSPEC int mpc_fma (mpc_ptr, mpc_srcptr, mpc_srcptr, mpc_srcptr, mpc_rnd_t);
181
182__MPC_DECLSPEC void mpc_set_nan (mpc_ptr);
183
184__MPC_DECLSPEC int mpc_real (mpfr_ptr, mpc_srcptr, mpfr_rnd_t);
185__MPC_DECLSPEC int mpc_imag (mpfr_ptr, mpc_srcptr, mpfr_rnd_t);
186__MPC_DECLSPEC int mpc_arg (mpfr_ptr, mpc_srcptr, mpfr_rnd_t);
187__MPC_DECLSPEC int mpc_proj (mpc_ptr, mpc_srcptr, mpc_rnd_t);
188__MPC_DECLSPEC int mpc_cmp (mpc_srcptr, mpc_srcptr);
189__MPC_DECLSPEC int mpc_cmp_si_si (mpc_srcptr, long int, long int);
190__MPC_DECLSPEC int mpc_cmp_abs (mpc_srcptr, mpc_srcptr);
191__MPC_DECLSPEC int mpc_exp (mpc_ptr, mpc_srcptr, mpc_rnd_t);
192__MPC_DECLSPEC int mpc_log (mpc_ptr, mpc_srcptr, mpc_rnd_t);
193__MPC_DECLSPEC int mpc_log10 (mpc_ptr, mpc_srcptr, mpc_rnd_t);
194__MPC_DECLSPEC int mpc_sin (mpc_ptr, mpc_srcptr, mpc_rnd_t);
195__MPC_DECLSPEC int mpc_cos (mpc_ptr, mpc_srcptr, mpc_rnd_t);
196__MPC_DECLSPEC int mpc_sin_cos (mpc_ptr, mpc_ptr, mpc_srcptr, mpc_rnd_t, mpc_rnd_t);
197__MPC_DECLSPEC int mpc_tan (mpc_ptr, mpc_srcptr, mpc_rnd_t);
198__MPC_DECLSPEC int mpc_sinh (mpc_ptr, mpc_srcptr, mpc_rnd_t);
199__MPC_DECLSPEC int mpc_cosh (mpc_ptr, mpc_srcptr, mpc_rnd_t);
200__MPC_DECLSPEC int mpc_tanh (mpc_ptr, mpc_srcptr, mpc_rnd_t);
201__MPC_DECLSPEC int mpc_asin (mpc_ptr, mpc_srcptr, mpc_rnd_t);
202__MPC_DECLSPEC int mpc_acos (mpc_ptr, mpc_srcptr, mpc_rnd_t);
203__MPC_DECLSPEC int mpc_atan (mpc_ptr, mpc_srcptr, mpc_rnd_t);
204__MPC_DECLSPEC int mpc_asinh (mpc_ptr, mpc_srcptr, mpc_rnd_t);
205__MPC_DECLSPEC int mpc_acosh (mpc_ptr, mpc_srcptr, mpc_rnd_t);
206__MPC_DECLSPEC int mpc_atanh (mpc_ptr, mpc_srcptr, mpc_rnd_t);
207__MPC_DECLSPEC int mpc_rootofunity (mpc_ptr, unsigned long int, unsigned long int, mpc_rnd_t);
208__MPC_DECLSPEC void mpc_clear (mpc_ptr);
209__MPC_DECLSPEC int mpc_urandom (mpc_ptr, gmp_randstate_t);
210__MPC_DECLSPEC void mpc_init2 (mpc_ptr, mpfr_prec_t);
211__MPC_DECLSPEC void mpc_init3 (mpc_ptr, mpfr_prec_t, mpfr_prec_t);
212__MPC_DECLSPEC mpfr_prec_t mpc_get_prec (mpc_srcptr x);
213__MPC_DECLSPEC void mpc_get_prec2 (mpfr_prec_t *pr, mpfr_prec_t *pi, mpc_srcptr x);
214__MPC_DECLSPEC void mpc_set_prec (mpc_ptr, mpfr_prec_t);
215__MPC_DECLSPEC const char * mpc_get_version (void);
216
217__MPC_DECLSPEC int mpc_strtoc (mpc_ptr, const char *, char **, int, mpc_rnd_t);
218__MPC_DECLSPEC int mpc_set_str (mpc_ptr, const char *, int, mpc_rnd_t);
219__MPC_DECLSPEC char * mpc_get_str (int, size_t, mpc_srcptr, mpc_rnd_t);
220__MPC_DECLSPEC void mpc_free_str (char *);
221
222/* declare certain functions only if appropriate headers have been included */
223#ifdef _MPC_H_HAVE_INTMAX_T
224__MPC_DECLSPEC int mpc_set_sj (mpc_ptr, intmax_t, mpc_rnd_t);
225__MPC_DECLSPEC int mpc_set_uj (mpc_ptr, uintmax_t, mpc_rnd_t);
226__MPC_DECLSPEC int mpc_set_sj_sj (mpc_ptr, intmax_t, intmax_t, mpc_rnd_t);
227__MPC_DECLSPEC int mpc_set_uj_uj (mpc_ptr, uintmax_t, uintmax_t, mpc_rnd_t);
228#endif
229
230#ifdef _Complex_I
231__MPC_DECLSPEC int mpc_set_dc (mpc_ptr, double _Complex, mpc_rnd_t);
232__MPC_DECLSPEC int mpc_set_ldc (mpc_ptr, long double _Complex, mpc_rnd_t);
233__MPC_DECLSPEC double _Complex mpc_get_dc (mpc_srcptr, mpc_rnd_t);
234__MPC_DECLSPEC long double _Complex mpc_get_ldc (mpc_srcptr, mpc_rnd_t);
235#endif
236
237#ifdef _GMP_H_HAVE_FILE
238__MPC_DECLSPEC int mpc_inp_str (mpc_ptr, FILE *, size_t *, int, mpc_rnd_t);
239__MPC_DECLSPEC size_t mpc_out_str (FILE *, int, size_t, mpc_srcptr, mpc_rnd_t);
240#endif
241
242#if defined (__cplusplus)
243}
244#endif
245
246#define mpc_realref(x) ((x)->re)
247#define mpc_imagref(x) ((x)->im)
248
249#define mpc_cmp_si(x, y) \
250 ( mpc_cmp_si_si ((x), (y), 0l) )
251#define mpc_ui_sub(x, y, z, r) mpc_ui_ui_sub (x, y, 0ul, z, r)
252
253/*
254 Define a fake mpfr_set_fr so that, for instance, mpc_set_fr_z would
255 be defined as follows:
256 mpc_set_fr_z (mpc_t rop, mpfr_t x, mpz_t y, mpc_rnd_t rnd)
257 MPC_SET_X_Y (fr, z, rop, x, y, rnd)
258*/
259#ifndef mpfr_set_fr
260#define mpfr_set_fr mpfr_set
261#endif
262#define MPC_SET_X_Y(real_t, imag_t, z, real_value, imag_value, rnd) \
263 { \
264 int _inex_re, _inex_im; \
265 _inex_re = (mpfr_set_ ## real_t) (mpc_realref (z), (real_value), MPC_RND_RE (rnd)); \
266 _inex_im = (mpfr_set_ ## imag_t) (mpc_imagref (z), (imag_value), MPC_RND_IM (rnd)); \
267 return MPC_INEX (_inex_re, _inex_im); \
268 }
269
270#endif /* ifndef __MPC_H */
271

source code of include/mpc.h