1/* GCC Quad-Precision Math Library
2 Copyright (C) 2010, 2011 Free Software Foundation, Inc.
3 Written by Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>
4
5This file is part of the libquadmath library.
6Libquadmath is free software; you can redistribute it and/or
7modify it under the terms of the GNU Library General Public
8License as published by the Free Software Foundation; either
9version 2 of the License, or (at your option) any later version.
10
11Libquadmath is distributed in the hope that it will be useful,
12but WITHOUT ANY WARRANTY; without even the implied warranty of
13MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14Library General Public License for more details.
15
16You should have received a copy of the GNU Library General Public
17License along with libquadmath; see the file COPYING.LIB. If
18not, write to the Free Software Foundation, Inc., 51 Franklin Street - Fifth Floor,
19Boston, MA 02110-1301, USA. */
20
21#ifndef QUADMATH_H
22#define QUADMATH_H
23
24#include <stdlib.h>
25
26#ifdef __cplusplus
27extern "C" {
28#endif
29
30/* Define the complex type corresponding to __float128
31 ("_Complex __float128" is not allowed) */
32#if (!defined(_ARCH_PPC)) || defined(__LONG_DOUBLE_IEEE128__)
33typedef _Complex float __attribute__((mode(TC))) __complex128;
34#else
35typedef _Complex float __attribute__((mode(KC))) __complex128;
36#endif
37
38#ifdef __cplusplus
39# define __quadmath_throw throw ()
40# define __quadmath_nth(fct) fct throw ()
41#else
42# define __quadmath_throw __attribute__((__nothrow__))
43# define __quadmath_nth(fct) __attribute__((__nothrow__)) fct
44#endif
45
46/* Prototypes for real functions */
47extern __float128 acosq (__float128) __quadmath_throw;
48extern __float128 acoshq (__float128) __quadmath_throw;
49extern __float128 asinq (__float128) __quadmath_throw;
50extern __float128 asinhq (__float128) __quadmath_throw;
51extern __float128 atanq (__float128) __quadmath_throw;
52extern __float128 atanhq (__float128) __quadmath_throw;
53extern __float128 atan2q (__float128, __float128) __quadmath_throw;
54extern __float128 cbrtq (__float128) __quadmath_throw;
55extern __float128 ceilq (__float128) __quadmath_throw;
56extern __float128 copysignq (__float128, __float128) __quadmath_throw;
57extern __float128 coshq (__float128) __quadmath_throw;
58extern __float128 cosq (__float128) __quadmath_throw;
59extern __float128 erfq (__float128) __quadmath_throw;
60extern __float128 erfcq (__float128) __quadmath_throw;
61extern __float128 expq (__float128) __quadmath_throw;
62extern __float128 expm1q (__float128) __quadmath_throw;
63extern __float128 fabsq (__float128) __quadmath_throw;
64extern __float128 fdimq (__float128, __float128) __quadmath_throw;
65extern int finiteq (__float128) __quadmath_throw;
66extern __float128 floorq (__float128) __quadmath_throw;
67extern __float128 fmaq (__float128, __float128, __float128) __quadmath_throw;
68extern __float128 fmaxq (__float128, __float128) __quadmath_throw;
69extern __float128 fminq (__float128, __float128) __quadmath_throw;
70extern __float128 fmodq (__float128, __float128) __quadmath_throw;
71extern __float128 frexpq (__float128, int *) __quadmath_throw;
72extern __float128 hypotq (__float128, __float128) __quadmath_throw;
73extern int isinfq (__float128) __quadmath_throw;
74extern int ilogbq (__float128) __quadmath_throw;
75extern int isnanq (__float128) __quadmath_throw;
76extern __float128 j0q (__float128) __quadmath_throw;
77extern __float128 j1q (__float128) __quadmath_throw;
78extern __float128 jnq (int, __float128) __quadmath_throw;
79extern __float128 ldexpq (__float128, int) __quadmath_throw;
80extern __float128 lgammaq (__float128) __quadmath_throw;
81extern long long int llrintq (__float128) __quadmath_throw;
82extern long long int llroundq (__float128) __quadmath_throw;
83extern __float128 logbq (__float128) __quadmath_throw;
84extern __float128 logq (__float128) __quadmath_throw;
85extern __float128 log10q (__float128) __quadmath_throw;
86extern __float128 log2q (__float128) __quadmath_throw;
87extern __float128 log1pq (__float128) __quadmath_throw;
88extern long int lrintq (__float128) __quadmath_throw;
89extern long int lroundq (__float128) __quadmath_throw;
90extern __float128 modfq (__float128, __float128 *) __quadmath_throw;
91extern __float128 nanq (const char *) __quadmath_throw;
92extern __float128 nearbyintq (__float128) __quadmath_throw;
93extern __float128 nextafterq (__float128, __float128) __quadmath_throw;
94extern __float128 powq (__float128, __float128) __quadmath_throw;
95extern __float128 remainderq (__float128, __float128) __quadmath_throw;
96extern __float128 remquoq (__float128, __float128, int *) __quadmath_throw;
97extern __float128 rintq (__float128) __quadmath_throw;
98extern __float128 roundq (__float128) __quadmath_throw;
99extern __float128 scalblnq (__float128, long int) __quadmath_throw;
100extern __float128 scalbnq (__float128, int) __quadmath_throw;
101extern int signbitq (__float128) __quadmath_throw;
102extern void sincosq (__float128, __float128 *, __float128 *) __quadmath_throw;
103extern __float128 sinhq (__float128) __quadmath_throw;
104extern __float128 sinq (__float128) __quadmath_throw;
105extern __float128 sqrtq (__float128) __quadmath_throw;
106extern __float128 tanq (__float128) __quadmath_throw;
107extern __float128 tanhq (__float128) __quadmath_throw;
108extern __float128 tgammaq (__float128) __quadmath_throw;
109extern __float128 truncq (__float128) __quadmath_throw;
110extern __float128 y0q (__float128) __quadmath_throw;
111extern __float128 y1q (__float128) __quadmath_throw;
112extern __float128 ynq (int, __float128) __quadmath_throw;
113
114
115/* Prototypes for complex functions */
116extern __float128 cabsq (__complex128) __quadmath_throw;
117extern __float128 cargq (__complex128) __quadmath_throw;
118extern __float128 cimagq (__complex128) __quadmath_throw;
119extern __float128 crealq (__complex128) __quadmath_throw;
120extern __complex128 cacosq (__complex128) __quadmath_throw;
121extern __complex128 cacoshq (__complex128) __quadmath_throw;
122extern __complex128 casinq (__complex128) __quadmath_throw;
123extern __complex128 casinhq (__complex128) __quadmath_throw;
124extern __complex128 catanq (__complex128) __quadmath_throw;
125extern __complex128 catanhq (__complex128) __quadmath_throw;
126extern __complex128 ccosq (__complex128) __quadmath_throw;
127extern __complex128 ccoshq (__complex128) __quadmath_throw;
128extern __complex128 cexpq (__complex128) __quadmath_throw;
129extern __complex128 cexpiq (__float128) __quadmath_throw;
130extern __complex128 clogq (__complex128) __quadmath_throw;
131extern __complex128 clog10q (__complex128) __quadmath_throw;
132extern __complex128 conjq (__complex128) __quadmath_throw;
133extern __complex128 cpowq (__complex128, __complex128) __quadmath_throw;
134extern __complex128 cprojq (__complex128) __quadmath_throw;
135extern __complex128 csinq (__complex128) __quadmath_throw;
136extern __complex128 csinhq (__complex128) __quadmath_throw;
137extern __complex128 csqrtq (__complex128) __quadmath_throw;
138extern __complex128 ctanq (__complex128) __quadmath_throw;
139extern __complex128 ctanhq (__complex128) __quadmath_throw;
140
141
142/* Prototypes for string <-> __float128 conversion functions */
143extern __float128 strtoflt128 (const char *, char **) __quadmath_throw;
144extern int quadmath_snprintf (char *str, size_t size,
145 const char *format, ...) __quadmath_throw;
146
147
148/* Macros */
149#define FLT128_MAX 1.18973149535723176508575932662800702e4932Q
150#define FLT128_MIN 3.36210314311209350626267781732175260e-4932Q
151#define FLT128_EPSILON 1.92592994438723585305597794258492732e-34Q
152#define FLT128_DENORM_MIN 6.475175119438025110924438958227646552e-4966Q
153#define FLT128_MANT_DIG 113
154#define FLT128_MIN_EXP (-16381)
155#define FLT128_MAX_EXP 16384
156#define FLT128_DIG 33
157#define FLT128_MIN_10_EXP (-4931)
158#define FLT128_MAX_10_EXP 4932
159
160
161#define HUGE_VALQ __builtin_huge_valq()
162/* The following alternative is valid, but brings the warning:
163 (floating constant exceeds range of ‘__float128’) */
164/* #define HUGE_VALQ (__extension__ 0x1.0p32767Q) */
165
166#define M_Eq 2.7182818284590452353602874713526625Q /* e */
167#define M_LOG2Eq 1.4426950408889634073599246810018921Q /* log_2 e */
168#define M_LOG10Eq 0.4342944819032518276511289189166051Q /* log_10 e */
169#define M_LN2q 0.6931471805599453094172321214581766Q /* log_e 2 */
170#define M_LN10q 2.3025850929940456840179914546843642Q /* log_e 10 */
171#define M_PIq 3.1415926535897932384626433832795029Q /* pi */
172#define M_PI_2q 1.5707963267948966192313216916397514Q /* pi/2 */
173#define M_PI_4q 0.7853981633974483096156608458198757Q /* pi/4 */
174#define M_1_PIq 0.3183098861837906715377675267450287Q /* 1/pi */
175#define M_2_PIq 0.6366197723675813430755350534900574Q /* 2/pi */
176#define M_2_SQRTPIq 1.1283791670955125738961589031215452Q /* 2/sqrt(pi) */
177#define M_SQRT2q 1.4142135623730950488016887242096981Q /* sqrt(2) */
178#define M_SQRT1_2q 0.7071067811865475244008443621048490Q /* 1/sqrt(2) */
179
180#define __quadmath_extern_inline \
181 extern inline __attribute__ ((__gnu_inline__))
182
183__quadmath_extern_inline __float128
184__quadmath_nth (cimagq (__complex128 __z))
185{
186 return __imag__ __z;
187}
188
189__quadmath_extern_inline __float128
190__quadmath_nth (crealq (__complex128 __z))
191{
192 return __real__ __z;
193}
194
195__quadmath_extern_inline __complex128
196__quadmath_nth (conjq (__complex128 __z))
197{
198 return __extension__ ~__z;
199}
200
201#ifdef __cplusplus
202}
203#endif
204
205#endif
206