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