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
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
|
/* GCC Quad-Precision Math Library
Copyright (C) 2010 Free Software Foundation, Inc.
Written by Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>
This file is part of the libiberty library.
Libiberty is free software; you can redistribute it and/or
modify it under the terms of the GNU Library General Public
License as published by the Free Software Foundation; either
version 2 of the License, or (at your option) any later version.
Libiberty 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
Library General Public License for more details.
You should have received a copy of the GNU Library General Public
License along with libiberty; see the file COPYING.LIB. If
not, write to the Free Software Foundation, Inc., 51 Franklin Street - Fifth Floor,
Boston, MA 02110-1301, USA. */
#ifndef QUADMATH_H
#define QUADMATH_H
#include <stdlib.h>
// Define the complex type corresponding to __float128
// ("_Complex __float128" is not allowed)
typedef _Complex float __attribute__((mode(TC))) __complex128;
// Prototypes for real functions
extern __float128 acosq (__float128);
extern __float128 acoshq (__float128);
extern __float128 asinq (__float128);
extern __float128 asinhq (__float128);
extern __float128 atanq (__float128);
extern __float128 atanhq (__float128);
extern __float128 atan2q (__float128, __float128);
extern __float128 cbrtq (const __float128);
extern __float128 ceilq (__float128);
extern __float128 copysignq (__float128, __float128);
extern __float128 coshq (__float128);
extern __float128 cosq (__float128);
extern __float128 erfq (__float128);
extern __float128 erfcq (__float128);
extern __float128 expq (__float128);
extern __float128 expm1q (__float128);
extern __float128 fabsq (__float128);
extern int finiteq (const __float128);
extern __float128 floorq (__float128);
extern __float128 fmaq (__float128, __float128, __float128);
extern __float128 fmodq (__float128, __float128);
extern __float128 frexpq (__float128, int *);
extern __float128 hypotq (__float128, __float128);
extern int isinfq (__float128);
extern int isnanq (const __float128);
extern __float128 j0q (__float128);
extern __float128 j1q (__float128);
extern __float128 jnq (int, __float128);
extern __float128 ldexpq (__float128, int);
extern __float128 lgammaq (__float128);
extern long long int llroundq (__float128);
extern __float128 logq (__float128);
extern __float128 log10q (__float128);
extern __float128 log1pq (__float128);
extern long int lroundq (__float128);
extern __float128 modfq (__float128, __float128 *);
extern __float128 nanq (const char *);
extern __float128 nextafterq (__float128, __float128);
extern __float128 powq (__float128, __float128);
extern __float128 remainderq (__float128, __float128);
extern __float128 rintq (__float128);
extern __float128 roundq (__float128);
extern __float128 scalblnq (__float128, long int);
extern __float128 scalbnq (__float128, int);
extern int signbitq (const __float128);
extern void sincosq (__float128, __float128 *, __float128 *);
extern __float128 sinhq (__float128);
extern __float128 sinq (__float128);
extern __float128 sqrtq (const __float128);
extern __float128 tanq (__float128);
extern __float128 tanhq (__float128);
extern __float128 tgammaq (__float128);
extern __float128 truncq (__float128);
extern __float128 y0q (__float128);
extern __float128 y1q (__float128);
extern __float128 ynq (int, __float128);
// Prototypes for complex functions
extern __float128 cabsq (__complex128);
extern __float128 cargq (__complex128);
extern __complex128 cacosq (__complex128);
extern __complex128 cacoshq (__complex128);
extern __complex128 casinq (__complex128);
extern __complex128 casinhq (__complex128);
extern __complex128 catanq (__complex128);
extern __complex128 catanhq (__complex128);
extern __complex128 ccosq (__complex128);
extern __complex128 ccoshq (__complex128);
extern __complex128 cexpq (__complex128);
extern __complex128 cexpiq (__float128);
extern __complex128 clogq (__complex128);
extern __complex128 clog10q (__complex128);
extern __complex128 cpowq (__complex128, __complex128);
extern __complex128 csinq (__complex128);
extern __complex128 csinhq (__complex128);
extern __complex128 csqrtq (__complex128);
extern __complex128 ctanq (__complex128);
extern __complex128 ctanhq (__complex128);
// Prototypes for our I/O functions
extern int quadmath_strtopQ (const char *, char **, void *);
extern void quadmath_dtoaq (char *, size_t, size_t, __float128);
// Macros
#define FLT128_MAX 1.18973149535723176508575932662800702e4932Q
#define FLT128_MIN 3.36210314311209350626267781732175260e-4932Q
#define FLT128_EPSILON 1.92592994438723585305597794258492732e-34Q
#define FLT128_DENORM_MIN 6.475175119438025110924438958227646552e-4966Q
#define FLT128_MANT_DIG 113
#define FLT128_MIN_EXP (-16381)
#define FLT128_MAX_EXP 16384
// TODO -- One day, we need to add the following macros:
// FLT128_DIG, FLT128_MIN_10_EXP, FLT128_MAX_10_EXP
#define HUGE_VALQ __builtin_huge_valq()
/* The following alternative is valid, but brings the warning:
(floating constant exceeds range of ‘__float128’) */
/* #define HUGE_VALQ (__extension__ 0x1.0p32767Q) */
#define M_Eq 2.7182818284590452353602874713526625Q /* e */
#define M_LOG2Eq 1.4426950408889634073599246810018921Q /* log_2 e */
#define M_LOG10Eq 0.4342944819032518276511289189166051Q /* log_10 e */
#define M_LN2q 0.6931471805599453094172321214581766Q /* log_e 2 */
#define M_LN10q 2.3025850929940456840179914546843642Q /* log_e 10 */
#define M_PIq 3.1415926535897932384626433832795029Q /* pi */
#define M_PI_2q 1.5707963267948966192313216916397514Q /* pi/2 */
#define M_PI_4q 0.7853981633974483096156608458198757Q /* pi/4 */
#define M_1_PIq 0.3183098861837906715377675267450287Q /* 1/pi */
#define M_2_PIq 0.6366197723675813430755350534900574Q /* 2/pi */
#define M_2_SQRTPIq 1.1283791670955125738961589031215452Q /* 2/sqrt(pi) */
#define M_SQRT2q 1.4142135623730950488016887242096981Q /* sqrt(2) */
#define M_SQRT1_2q 0.7071067811865475244008443621048490Q /* 1/sqrt(2) */
#endif
|
