| 1 | /* Copyright (C) 1997-2022 Free Software Foundation, Inc. |
| 2 | This file is part of the GNU C Library. |
| 3 | |
| 4 | The GNU C Library is free software; you can redistribute it and/or |
| 5 | modify it under the terms of the GNU Lesser General Public |
| 6 | License as published by the Free Software Foundation; either |
| 7 | version 2.1 of the License, or (at your option) any later version. |
| 8 | |
| 9 | The GNU C Library is distributed in the hope that it will be useful, |
| 10 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
| 11 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| 12 | Lesser General Public License for more details. |
| 13 | |
| 14 | You should have received a copy of the GNU Lesser General Public |
| 15 | License along with the GNU C Library; if not, see |
| 16 | <https://www.gnu.org/licenses/>. */ |
| 17 | |
| 18 | /* |
| 19 | * ISO C99: 7.3 Complex arithmetic <complex.h> |
| 20 | */ |
| 21 | |
| 22 | #ifndef _COMPLEX_H |
| 23 | #define _COMPLEX_H 1 |
| 24 | |
| 25 | #define |
| 26 | #include <bits/libc-header-start.h> |
| 27 | |
| 28 | /* Get general and ISO C99 specific information. */ |
| 29 | #include <bits/mathdef.h> |
| 30 | |
| 31 | /* Gather machine-dependent _FloatN type support. */ |
| 32 | #include <bits/floatn.h> |
| 33 | |
| 34 | __BEGIN_DECLS |
| 35 | |
| 36 | /* We might need to add support for more compilers here. But since ISO |
| 37 | C99 is out hopefully all maintained compilers will soon provide the data |
| 38 | types `float complex' and `double complex'. */ |
| 39 | #if __GNUC_PREREQ (2, 7) && !__GNUC_PREREQ (2, 97) |
| 40 | # define _Complex __complex__ |
| 41 | #endif |
| 42 | |
| 43 | #define complex _Complex |
| 44 | |
| 45 | /* Narrowest imaginary unit. This depends on the floating-point |
| 46 | evaluation method. |
| 47 | XXX This probably has to go into a gcc related file. */ |
| 48 | #define _Complex_I (__extension__ 1.0iF) |
| 49 | |
| 50 | /* Another more descriptive name is `I'. |
| 51 | XXX Once we have the imaginary support switch this to _Imaginary_I. */ |
| 52 | #undef I |
| 53 | #define I _Complex_I |
| 54 | |
| 55 | #if defined __USE_ISOC11 && __GNUC_PREREQ (4, 7) |
| 56 | /* Macros to expand into expression of specified complex type. */ |
| 57 | # define CMPLX(x, y) __builtin_complex ((double) (x), (double) (y)) |
| 58 | # define CMPLXF(x, y) __builtin_complex ((float) (x), (float) (y)) |
| 59 | # define CMPLXL(x, y) __builtin_complex ((long double) (x), (long double) (y)) |
| 60 | #endif |
| 61 | |
| 62 | #if __HAVE_FLOAT16 && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 63 | # define CMPLXF16(x, y) __builtin_complex ((_Float16) (x), (_Float16) (y)) |
| 64 | #endif |
| 65 | |
| 66 | #if __HAVE_FLOAT32 && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 67 | # define CMPLXF32(x, y) __builtin_complex ((_Float32) (x), (_Float32) (y)) |
| 68 | #endif |
| 69 | |
| 70 | #if __HAVE_FLOAT64 && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 71 | # define CMPLXF64(x, y) __builtin_complex ((_Float64) (x), (_Float64) (y)) |
| 72 | #endif |
| 73 | |
| 74 | #if __HAVE_FLOAT128 && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 75 | # define CMPLXF128(x, y) __builtin_complex ((_Float128) (x), (_Float128) (y)) |
| 76 | #endif |
| 77 | |
| 78 | #if __HAVE_FLOAT32X && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 79 | # define CMPLXF32X(x, y) __builtin_complex ((_Float32x) (x), (_Float32x) (y)) |
| 80 | #endif |
| 81 | |
| 82 | #if __HAVE_FLOAT64X && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 83 | # define CMPLXF64X(x, y) __builtin_complex ((_Float64x) (x), (_Float64x) (y)) |
| 84 | #endif |
| 85 | |
| 86 | #if __HAVE_FLOAT128X && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 87 | # define CMPLXF128X(x, y) \ |
| 88 | __builtin_complex ((_Float128x) (x), (_Float128x) (y)) |
| 89 | #endif |
| 90 | |
| 91 | /* The file <bits/cmathcalls.h> contains the prototypes for all the |
| 92 | actual math functions. These macros are used for those prototypes, |
| 93 | so we can easily declare each function as both `name' and `__name', |
| 94 | and can declare the float versions `namef' and `__namef'. */ |
| 95 | |
| 96 | #define __MATHCALL(function, args) \ |
| 97 | __MATHDECL (_Mdouble_complex_,function, args) |
| 98 | #define __MATHDECL_IMPL(type, function, args) \ |
| 99 | __MATHDECL_1(type, function, args); \ |
| 100 | __MATHDECL_1(type, __CONCAT(__,function), args) |
| 101 | #define __MATHDECL(type, function, args) \ |
| 102 | __MATHDECL_IMPL(type, function, args) |
| 103 | #define __MATHDECL_1_IMPL(type, function, args) \ |
| 104 | extern type __MATH_PRECNAME(function) args __THROW |
| 105 | #define __MATHDECL_1(type, function, args) \ |
| 106 | __MATHDECL_1_IMPL(type, function, args) |
| 107 | |
| 108 | #define _Mdouble_ double |
| 109 | #define __MATH_PRECNAME(name) name |
| 110 | #include <bits/cmathcalls.h> |
| 111 | #undef _Mdouble_ |
| 112 | #undef __MATH_PRECNAME |
| 113 | |
| 114 | /* Now the float versions. */ |
| 115 | #define _Mdouble_ float |
| 116 | #define __MATH_PRECNAME(name) name##f |
| 117 | #include <bits/cmathcalls.h> |
| 118 | #undef _Mdouble_ |
| 119 | #undef __MATH_PRECNAME |
| 120 | |
| 121 | /* And the long double versions. It is non-critical to define them |
| 122 | here unconditionally since `long double' is required in ISO C99. */ |
| 123 | #if !(defined __NO_LONG_DOUBLE_MATH && defined _LIBC) \ |
| 124 | || defined __LDBL_COMPAT |
| 125 | # ifdef __LDBL_COMPAT |
| 126 | # undef __MATHDECL_1 |
| 127 | # define __MATHDECL_1(type, function, args) \ |
| 128 | extern type __REDIRECT_NTH(__MATH_PRECNAME(function), args, function) |
| 129 | # elif __LDOUBLE_REDIRECTS_TO_FLOAT128_ABI == 1 |
| 130 | # undef __MATHDECL_1 |
| 131 | # undef __MATHDECL |
| 132 | # define __REDIR_TO(function) \ |
| 133 | __ ## function ## ieee128 |
| 134 | # define __MATHDECL_1(type, function, alias, args) \ |
| 135 | extern type __REDIRECT_NTH(__MATH_PRECNAME(function), args, alias) |
| 136 | #define __MATHDECL(type, function, args) \ |
| 137 | __MATHDECL_1(type, function, __REDIR_TO(function), args); \ |
| 138 | __MATHDECL_1(type, __CONCAT(__,function), __REDIR_TO(function), args) |
| 139 | # endif |
| 140 | |
| 141 | # define _Mdouble_ long double |
| 142 | # define __MATH_PRECNAME(name) name##l |
| 143 | # include <bits/cmathcalls.h> |
| 144 | # if defined __LDBL_COMPAT \ |
| 145 | || __LDOUBLE_REDIRECTS_TO_FLOAT128_ABI == 1 |
| 146 | # undef __REDIR_TO |
| 147 | # undef __MATHDECL_1 |
| 148 | # undef __MATHDECL |
| 149 | #define __MATHDECL(type, function, args) \ |
| 150 | __MATHDECL_IMPL(type, function, args) |
| 151 | # define __MATHDECL_1(type, function, args) \ |
| 152 | __MATHDECL_1_IMPL(type, function, args) |
| 153 | # endif |
| 154 | #endif |
| 155 | #undef _Mdouble_ |
| 156 | #undef __MATH_PRECNAME |
| 157 | |
| 158 | #if (__HAVE_DISTINCT_FLOAT16 || (__HAVE_FLOAT16 && !defined _LIBC)) \ |
| 159 | && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 160 | # undef _Mdouble_complex_ |
| 161 | # define _Mdouble_complex_ __CFLOAT16 |
| 162 | # define _Mdouble_ _Float16 |
| 163 | # define __MATH_PRECNAME(name) name##f16 |
| 164 | # include <bits/cmathcalls.h> |
| 165 | # undef _Mdouble_ |
| 166 | # undef __MATH_PRECNAME |
| 167 | # undef _Mdouble_complex_ |
| 168 | #endif |
| 169 | |
| 170 | #if (__HAVE_DISTINCT_FLOAT32 || (__HAVE_FLOAT32 && !defined _LIBC)) \ |
| 171 | && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 172 | # undef _Mdouble_complex_ |
| 173 | # define _Mdouble_complex_ __CFLOAT32 |
| 174 | # define _Mdouble_ _Float32 |
| 175 | # define __MATH_PRECNAME(name) name##f32 |
| 176 | # include <bits/cmathcalls.h> |
| 177 | # undef _Mdouble_ |
| 178 | # undef __MATH_PRECNAME |
| 179 | # undef _Mdouble_complex_ |
| 180 | #endif |
| 181 | |
| 182 | #if (__HAVE_DISTINCT_FLOAT64 || (__HAVE_FLOAT64 && !defined _LIBC)) \ |
| 183 | && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 184 | # undef _Mdouble_complex_ |
| 185 | # define _Mdouble_complex_ __CFLOAT64 |
| 186 | # define _Mdouble_ _Float64 |
| 187 | # define __MATH_PRECNAME(name) name##f64 |
| 188 | # include <bits/cmathcalls.h> |
| 189 | # undef _Mdouble_ |
| 190 | # undef __MATH_PRECNAME |
| 191 | # undef _Mdouble_complex_ |
| 192 | #endif |
| 193 | |
| 194 | #if (__HAVE_DISTINCT_FLOAT128 || (__HAVE_FLOAT128 && !defined _LIBC)) \ |
| 195 | && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 196 | # undef _Mdouble_complex_ |
| 197 | # define _Mdouble_complex_ __CFLOAT128 |
| 198 | # define _Mdouble_ _Float128 |
| 199 | # define __MATH_PRECNAME(name) name##f128 |
| 200 | # include <bits/cmathcalls.h> |
| 201 | # undef _Mdouble_ |
| 202 | # undef __MATH_PRECNAME |
| 203 | # undef _Mdouble_complex_ |
| 204 | #endif |
| 205 | |
| 206 | #if (__HAVE_DISTINCT_FLOAT32X || (__HAVE_FLOAT32X && !defined _LIBC)) \ |
| 207 | && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 208 | # undef _Mdouble_complex_ |
| 209 | # define _Mdouble_complex_ __CFLOAT32X |
| 210 | # define _Mdouble_ _Float32x |
| 211 | # define __MATH_PRECNAME(name) name##f32x |
| 212 | # include <bits/cmathcalls.h> |
| 213 | # undef _Mdouble_ |
| 214 | # undef __MATH_PRECNAME |
| 215 | # undef _Mdouble_complex_ |
| 216 | #endif |
| 217 | |
| 218 | #if (__HAVE_DISTINCT_FLOAT64X || (__HAVE_FLOAT64X && !defined _LIBC)) \ |
| 219 | && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 220 | # undef _Mdouble_complex_ |
| 221 | # define _Mdouble_complex_ __CFLOAT64X |
| 222 | # define _Mdouble_ _Float64x |
| 223 | # define __MATH_PRECNAME(name) name##f64x |
| 224 | # include <bits/cmathcalls.h> |
| 225 | # undef _Mdouble_ |
| 226 | # undef __MATH_PRECNAME |
| 227 | # undef _Mdouble_complex_ |
| 228 | #endif |
| 229 | |
| 230 | #if (__HAVE_DISTINCT_FLOAT128X || (__HAVE_FLOAT128X && !defined _LIBC)) \ |
| 231 | && __GLIBC_USE (IEC_60559_TYPES_EXT) |
| 232 | # undef _Mdouble_complex_ |
| 233 | # define _Mdouble_complex_ __CFLOAT128X |
| 234 | # define _Mdouble_ _Float128x |
| 235 | # define __MATH_PRECNAME(name) name##f128x |
| 236 | # include <bits/cmathcalls.h> |
| 237 | # undef _Mdouble_ |
| 238 | # undef __MATH_PRECNAME |
| 239 | # undef _Mdouble_complex_ |
| 240 | #endif |
| 241 | |
| 242 | #undef __MATHDECL_1_IMPL |
| 243 | #undef __MATHDECL_1 |
| 244 | #undef __MATHDECL |
| 245 | #undef __MATHCALL |
| 246 | |
| 247 | __END_DECLS |
| 248 | |
| 249 | #endif /* complex.h */ |
| 250 | |