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