| 1 | /* _Float128 overrides for building ldbl-128 as _Float128. |
|---|
| 2 | Copyright (C) 2017-2024 Free Software Foundation, Inc. |
|---|
| 3 | This file is part of the GNU C Library. |
|---|
| 4 | |
|---|
| 5 | The GNU C Library is free software; you can redistribute it and/or |
|---|
| 6 | modify it under the terms of the GNU Lesser General Public |
|---|
| 7 | License as published by the Free Software Foundation; either |
|---|
| 8 | version 2.1 of the License, or (at your option) any later version. |
|---|
| 9 | |
|---|
| 10 | The GNU C Library is distributed in the hope that it will be useful, |
|---|
| 11 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
|---|
| 12 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
|---|
| 13 | Lesser General Public License for more details. |
|---|
| 14 | |
|---|
| 15 | You should have received a copy of the GNU Lesser General Public |
|---|
| 16 | License along with the GNU C Library; if not, see |
|---|
| 17 | <https://www.gnu.org/licenses/>. */ |
|---|
| 18 | |
|---|
| 19 | /* This must be included before the function renames below. */ |
|---|
| 20 | #include <gmp.h> |
|---|
| 21 | #include <math.h> |
|---|
| 22 | #undef HUGE_VALL |
|---|
| 23 | #define HUGE_VALL HUGE_VAL_F128 |
|---|
| 24 | #include <math/mul_splitl.h> |
|---|
| 25 | |
|---|
| 26 | /* This must be included before the renames of types and macros from |
|---|
| 27 | it. */ |
|---|
| 28 | #include <ieee754.h> |
|---|
| 29 | |
|---|
| 30 | /* Renames derived from math_private.h. */ |
|---|
| 31 | #include <math_private.h> |
|---|
| 32 | #include <fenv_private.h> |
|---|
| 33 | #include <ieee754_float128.h> |
|---|
| 34 | #define ieee854_long_double_shape_type ieee854_float128_shape_type |
|---|
| 35 | #define ieee854_long_double ieee854_float128 |
|---|
| 36 | |
|---|
| 37 | #undef GET_LDOUBLE_LSW64 |
|---|
| 38 | #undef GET_LDOUBLE_MSW64 |
|---|
| 39 | #undef GET_LDOUBLE_WORDS64 |
|---|
| 40 | #undef SET_LDOUBLE_LSW64 |
|---|
| 41 | #undef SET_LDOUBLE_MSW64 |
|---|
| 42 | #undef SET_LDOUBLE_WORDS64 |
|---|
| 43 | #define GET_LDOUBLE_LSW64(x,y) GET_FLOAT128_LSW64 (x, y) |
|---|
| 44 | #define GET_LDOUBLE_MSW64(x,y) GET_FLOAT128_MSW64 (x, y) |
|---|
| 45 | #define GET_LDOUBLE_WORDS64(x,y,z) GET_FLOAT128_WORDS64 (x, y, z) |
|---|
| 46 | #define SET_LDOUBLE_LSW64(x,y) SET_FLOAT128_LSW64 (x, y) |
|---|
| 47 | #define SET_LDOUBLE_MSW64(x,y) SET_FLOAT128_MSW64 (x, y) |
|---|
| 48 | #define SET_LDOUBLE_WORDS64(x,y,z) SET_FLOAT128_WORDS64 (x, y, z) |
|---|
| 49 | |
|---|
| 50 | #undef IEEE854_LONG_DOUBLE_BIAS |
|---|
| 51 | #define IEEE854_LONG_DOUBLE_BIAS IEEE854_FLOAT128_BIAS |
|---|
| 52 | |
|---|
| 53 | #ifdef SET_RESTORE_ROUNDF128 |
|---|
| 54 | # undef SET_RESTORE_ROUNDL |
|---|
| 55 | # define SET_RESTORE_ROUNDL(RM) SET_RESTORE_ROUNDF128 (RM) |
|---|
| 56 | #endif |
|---|
| 57 | |
|---|
| 58 | #ifdef libc_feholdexcept_setroundf128 |
|---|
| 59 | # undef libc_feholdexcept_setroundl |
|---|
| 60 | # define libc_feholdexcept_setroundl(ENV, RM) \ |
|---|
| 61 | libc_feholdexcept_setroundf128 (ENV, RM) |
|---|
| 62 | #endif |
|---|
| 63 | |
|---|
| 64 | #ifdef libc_feupdateenv_testf128 |
|---|
| 65 | # undef libc_feupdateenv_testl |
|---|
| 66 | # define libc_feupdateenv_testl(ENV, EX) libc_feupdateenv_testf128 (ENV, EX) |
|---|
| 67 | #endif |
|---|
| 68 | |
|---|
| 69 | /* misc macros from the header below. */ |
|---|
| 70 | #include <fix-fp-int-convert-overflow.h> |
|---|
| 71 | #undef FIX_LDBL_LONG_CONVERT_OVERFLOW |
|---|
| 72 | #undef FIX_LDBL_LLONG_CONVERT_OVERFLOW |
|---|
| 73 | #define FIX_LDBL_LONG_CONVERT_OVERFLOW FIX_FLT128_LONG_CONVERT_OVERFLOW |
|---|
| 74 | #define FIX_LDBL_LLONG_CONVERT_OVERFLOW FIX_FLT128_LLONG_CONVERT_OVERFLOW |
|---|
| 75 | |
|---|
| 76 | |
|---|
| 77 | /* float.h constants. */ |
|---|
| 78 | #include <float.h> |
|---|
| 79 | #undef LDBL_DIG |
|---|
| 80 | #undef LDBL_EPSILON |
|---|
| 81 | #undef LDBL_MANT_DIG |
|---|
| 82 | #undef LDBL_MAX |
|---|
| 83 | #undef LDBL_MAX_10_EXP |
|---|
| 84 | #undef LDBL_MAX_EXP |
|---|
| 85 | #undef LDBL_MIN |
|---|
| 86 | #undef LDBL_MIN_10_EXP |
|---|
| 87 | #undef LDBL_MIN_EXP |
|---|
| 88 | #undef LDBL_TRUE_MIN |
|---|
| 89 | #define LDBL_DIG FLT128_DIG |
|---|
| 90 | #define LDBL_EPSILON FLT128_EPSILON |
|---|
| 91 | #define LDBL_MANT_DIG FLT128_MANT_DIG |
|---|
| 92 | #define LDBL_MAX FLT128_MAX |
|---|
| 93 | #define LDBL_MAX_10_EXP FLT128_MAX_10_EXP |
|---|
| 94 | #define LDBL_MAX_EXP FLT128_MAX_EXP |
|---|
| 95 | #define LDBL_MIN FLT128_MIN |
|---|
| 96 | #define LDBL_MIN_10_EXP FLT128_MIN_10_EXP |
|---|
| 97 | #define LDBL_MIN_EXP FLT128_MIN_EXP |
|---|
| 98 | #define LDBL_TRUE_MIN FLT128_TRUE_MIN |
|---|
| 99 | |
|---|
| 100 | |
|---|
| 101 | /* math.h GNU constants. */ |
|---|
| 102 | #undef M_El |
|---|
| 103 | #undef M_LOG2El |
|---|
| 104 | #undef M_LOG10El |
|---|
| 105 | #undef M_LN2l |
|---|
| 106 | #undef M_LN10l |
|---|
| 107 | #undef M_PIl |
|---|
| 108 | #undef M_PI_2l |
|---|
| 109 | #undef M_PI_4l |
|---|
| 110 | #undef M_1_PIl |
|---|
| 111 | #undef M_2_PIl |
|---|
| 112 | #undef M_2_SQRTPIl |
|---|
| 113 | #undef M_SQRT2l |
|---|
| 114 | #undef M_SQRT1_2l |
|---|
| 115 | #define M_El M_Ef128 |
|---|
| 116 | #define M_LOG2El M_LOG2Ef128 |
|---|
| 117 | #define M_LOG10El M_LOG10Ef128 |
|---|
| 118 | #define M_LN2l M_LN2f128 |
|---|
| 119 | #define M_LN10l M_LN10f128 |
|---|
| 120 | #define M_PIl M_PIf128 |
|---|
| 121 | #define M_PI_2l M_PI_2f128 |
|---|
| 122 | #define M_PI_4l M_PI_4f128 |
|---|
| 123 | #define M_1_PIl M_1_PIf128 |
|---|
| 124 | #define M_2_PIl M_2_PIf128 |
|---|
| 125 | #define M_2_SQRTPIl M_2_SQRTPIf128 |
|---|
| 126 | #define M_SQRT2l M_SQRT2f128 |
|---|
| 127 | #define M_SQRT1_2l M_SQRT1_2f128 |
|---|
| 128 | |
|---|
| 129 | |
|---|
| 130 | #include <libm-alias-ldouble.h> |
|---|
| 131 | #include <libm-alias-float128.h> |
|---|
| 132 | #undef libm_alias_ldouble_r |
|---|
| 133 | #define libm_alias_ldouble_r(from, to, r) libm_alias_float128_r (from, to, r) |
|---|
| 134 | |
|---|
| 135 | |
|---|
| 136 | #include <math/math-narrow.h> |
|---|
| 137 | #undef libm_alias_float_ldouble |
|---|
| 138 | #define libm_alias_float_ldouble(func) libm_alias_float32_float128 (func) |
|---|
| 139 | #undef libm_alias_double_ldouble |
|---|
| 140 | #define libm_alias_double_ldouble(func) libm_alias_float64_float128 (func) |
|---|
| 141 | #undef libm_alias_ldouble_narrow |
|---|
| 142 | #define libm_alias_ldouble_narrow(from, to) \ |
|---|
| 143 | libm_alias_float128_narrow (from, to) |
|---|
| 144 | |
|---|
| 145 | #include <math-use-builtins.h> |
|---|
| 146 | #undef USE_NEARBYINTL_BUILTIN |
|---|
| 147 | #define USE_NEARBYINTL_BUILTIN USE_NEARBYINTF128_BUILTIN |
|---|
| 148 | #undef USE_RINTL_BUILTIN |
|---|
| 149 | #define USE_RINTL_BUILTIN USE_RINTF128_BUILTIN |
|---|
| 150 | #undef USE_FLOORL_BUILTIN |
|---|
| 151 | #define USE_FLOORL_BUILTIN USE_FLOORF128_BUILTIN |
|---|
| 152 | #undef USE_CEILL_BUILTIN |
|---|
| 153 | #define USE_CEILL_BUILTIN USE_CEILF128_BUILTIN |
|---|
| 154 | #undef USE_TRUNCL_BUILTIN |
|---|
| 155 | #define USE_TRUNCL_BUILTIN USE_TRUNCF128_BUILTIN |
|---|
| 156 | #undef USE_ROUNDL_BUILTIN |
|---|
| 157 | #define USE_ROUNDL_BUILTIN USE_ROUNDF128_BUILTIN |
|---|
| 158 | #undef USE_ROUNDEVENL_BUILTIN |
|---|
| 159 | #define USE_ROUNDEVENL_BUILTIN USE_ROUNDEVENF128_BUILTIN |
|---|
| 160 | #undef USE_COPYSIGNL_BUILTIN |
|---|
| 161 | #define USE_COPYSIGNL_BUILTIN USE_COPYSIGNF128_BUILTIN |
|---|
| 162 | #undef USE_FMAL_BUILTIN |
|---|
| 163 | #define USE_FMAL_BUILTIN USE_FMAF128_BUILTIN |
|---|
| 164 | #undef USE_LRINTL_BUILTIN |
|---|
| 165 | #define USE_LRINTL_BUILTIN USE_LRINTF128_BUILTIN |
|---|
| 166 | #undef USE_LLRINTL_BUILTIN |
|---|
| 167 | #define USE_LLRINTL_BUILTIN USE_LLRINTF128_BUILTIN |
|---|
| 168 | #undef USE_LOGBL_BUILTIN |
|---|
| 169 | #define USE_LOGBL_BUILTIN USE_LOGBF128_BUILTIN |
|---|
| 170 | |
|---|
| 171 | /* IEEE function renames. */ |
|---|
| 172 | #define __ieee754_acoshl __ieee754_acoshf128 |
|---|
| 173 | #define __ieee754_acosl __ieee754_acosf128 |
|---|
| 174 | #define __ieee754_asinhl __ieee754_asinhf128 |
|---|
| 175 | #define __ieee754_asinl __ieee754_asinf128 |
|---|
| 176 | #define __ieee754_atan2l __ieee754_atan2f128 |
|---|
| 177 | #define __ieee754_atanhl __ieee754_atanhf128 |
|---|
| 178 | #define __ieee754_coshl __ieee754_coshf128 |
|---|
| 179 | #define __ieee754_cosl __ieee754_cosf128 |
|---|
| 180 | #define __ieee754_exp10l __ieee754_exp10f128 |
|---|
| 181 | #define __ieee754_exp2l __ieee754_exp2f128 |
|---|
| 182 | #define __ieee754_expl __ieee754_expf128 |
|---|
| 183 | #define __ieee754_fmodl __ieee754_fmodf128 |
|---|
| 184 | #define __ieee754_gammal_r __ieee754_gammaf128_r |
|---|
| 185 | #define __ieee754_hypotl __ieee754_hypotf128 |
|---|
| 186 | #define __ieee754_ilogbl __ieee754_ilogbf128 |
|---|
| 187 | #define __ieee754_j0l __ieee754_j0f128 |
|---|
| 188 | #define __ieee754_j1l __ieee754_j1f128 |
|---|
| 189 | #define __ieee754_jnl __ieee754_jnf128 |
|---|
| 190 | #define __ieee754_lgammal_r __ieee754_lgammaf128_r |
|---|
| 191 | #define __ieee754_log10l __ieee754_log10f128 |
|---|
| 192 | #define __ieee754_log2l __ieee754_log2f128 |
|---|
| 193 | #define __ieee754_logl __ieee754_logf128 |
|---|
| 194 | #define __ieee754_powl __ieee754_powf128 |
|---|
| 195 | #define __ieee754_rem_pio2l __ieee754_rem_pio2f128 |
|---|
| 196 | #define __ieee754_remainderl __ieee754_remainderf128 |
|---|
| 197 | #define __ieee754_sinhl __ieee754_sinhf128 |
|---|
| 198 | #define __ieee754_sqrtl __ieee754_sqrtf128 |
|---|
| 199 | #define __ieee754_y0l __ieee754_y0f128 |
|---|
| 200 | #define __ieee754_y1l __ieee754_y1f128 |
|---|
| 201 | #define __ieee754_ynl __ieee754_ynf128 |
|---|
| 202 | |
|---|
| 203 | |
|---|
| 204 | /* finite math entry points. */ |
|---|
| 205 | #define __acoshl_finite __acoshf128_finite |
|---|
| 206 | #define __acosl_finite __acosf128_finite |
|---|
| 207 | #define __asinl_finite __asinf128_finite |
|---|
| 208 | #define __atan2l_finite __atan2f128_finite |
|---|
| 209 | #define __atanhl_finite __atanhf128_finite |
|---|
| 210 | #define __coshl_finite __coshf128_finite |
|---|
| 211 | #define __cosl_finite __cosf128_finite |
|---|
| 212 | #define __exp10l_finite __exp10f128_finite |
|---|
| 213 | #define __exp2l_finite __exp2f128_finite |
|---|
| 214 | #define __expl_finite __expf128_finite |
|---|
| 215 | #define __fmodl_finite __fmodf128_finite |
|---|
| 216 | #define __hypotl_finite __hypotf128_finite |
|---|
| 217 | #define __ilogbl_finite __ilogbf128_finite |
|---|
| 218 | #define __j0l_finite __j0f128_finite |
|---|
| 219 | #define __j1l_finite __j1f128_finite |
|---|
| 220 | #define __jnl_finite __jnf128_finite |
|---|
| 221 | #define __lgammal_r_finite __lgammaf128_r_finite |
|---|
| 222 | #define __log10l_finite __log10f128_finite |
|---|
| 223 | #define __log2l_finite __log2f128_finite |
|---|
| 224 | #define __logl_finite __logf128_finite |
|---|
| 225 | #define __powl_finite __powf128_finite |
|---|
| 226 | #define __remainderl_finite __remainderf128_finite |
|---|
| 227 | #define __sinhl_finite __sinhf128_finite |
|---|
| 228 | #define __y0l_finite __y0f128_finite |
|---|
| 229 | #define __y1l_finite __y1f128_finite |
|---|
| 230 | #define __ynl_finite __ynf128_finite |
|---|
| 231 | |
|---|
| 232 | |
|---|
| 233 | /* internal function names. */ |
|---|
| 234 | #define __asinhl __asinhf128 |
|---|
| 235 | #define __atanl __atanf128 |
|---|
| 236 | #define __cbrtl __cbrtf128 |
|---|
| 237 | #define __ceill __ceilf128 |
|---|
| 238 | #define __copysignl __copysignf128 |
|---|
| 239 | #define __cosl __cosf128 |
|---|
| 240 | #define __erfcl __erfcf128 |
|---|
| 241 | #define __erfl __erff128 |
|---|
| 242 | #define __expl __expf128 |
|---|
| 243 | #define __expm1l __expm1f128 |
|---|
| 244 | #define __fabsl __fabsf128 |
|---|
| 245 | #define __fdiml __fdimf128 |
|---|
| 246 | #define __finitel __finitef128 |
|---|
| 247 | #define __floorl __floorf128 |
|---|
| 248 | #define __fmal __fmaf128 |
|---|
| 249 | #define __fmaxl __fmaxf128 |
|---|
| 250 | #define __fminl __fminf128 |
|---|
| 251 | #define __fpclassifyl __fpclassifyf128 |
|---|
| 252 | #define __frexpl __frexpf128 |
|---|
| 253 | #define __gammal_r_finite __gammaf128_r_finite |
|---|
| 254 | #define __getpayloadl __getpayloadf128 |
|---|
| 255 | #define __isinfl __isinff128 |
|---|
| 256 | #define __isnanl __isnanf128 |
|---|
| 257 | #define __issignalingl __issignalingf128 |
|---|
| 258 | #define __ldexpl __ldexpf128 |
|---|
| 259 | #define __llrintl __llrintf128 |
|---|
| 260 | #define __llroundl __llroundf128 |
|---|
| 261 | #define __log1pl __log1pf128 |
|---|
| 262 | #define __logbl __logbf128 |
|---|
| 263 | #define __logl __logf128 |
|---|
| 264 | #define __lrintl __lrintf128 |
|---|
| 265 | #define __lroundl __lroundf128 |
|---|
| 266 | #define __modfl __modff128 |
|---|
| 267 | #define __nearbyintl __nearbyintf128 |
|---|
| 268 | #define __nextafterl __nextafterf128 |
|---|
| 269 | #define __nextdownl __nextdownf128 |
|---|
| 270 | #define __nextupl __nextupf128 |
|---|
| 271 | #define __remquol __remquof128 |
|---|
| 272 | #define __rintl __rintf128 |
|---|
| 273 | #define __roundevenl __roundevenf128 |
|---|
| 274 | #define __roundl __roundf128 |
|---|
| 275 | #define __scalblnl __scalblnf128 |
|---|
| 276 | #define __scalbnl __scalbnf128 |
|---|
| 277 | #define __signbitl __signbitf128 |
|---|
| 278 | #define __sincosl __sincosf128 |
|---|
| 279 | #define __sinl __sinf128 |
|---|
| 280 | #define __sqrtl __sqrtf128 |
|---|
| 281 | #define __tanhl __tanhf128 |
|---|
| 282 | #define __tanl __tanf128 |
|---|
| 283 | #define __totalorderl __totalorderf128 |
|---|
| 284 | #define __totalorder_compatl __totalorder_compatf128 |
|---|
| 285 | #define __totalordermagl __totalordermagf128 |
|---|
| 286 | #define __totalordermag_compatl __totalordermag_compatf128 |
|---|
| 287 | #define __truncl __truncf128 |
|---|
| 288 | #define __x2y2m1l __x2y2m1f128 |
|---|
| 289 | |
|---|
| 290 | #define __faddl __f32addf128 |
|---|
| 291 | #define __daddl __f64addf128 |
|---|
| 292 | #define __fdivl __f32divf128 |
|---|
| 293 | #define __ddivl __f64divf128 |
|---|
| 294 | #define __ffmal __f32fmaf128 |
|---|
| 295 | #define __dfmal __f64fmaf128 |
|---|
| 296 | #define __fmull __f32mulf128 |
|---|
| 297 | #define __dmull __f64mulf128 |
|---|
| 298 | #define __fsqrtl __f32sqrtf128 |
|---|
| 299 | #define __dsqrtl __f64sqrtf128 |
|---|
| 300 | #define __fsubl __f32subf128 |
|---|
| 301 | #define __dsubl __f64subf128 |
|---|
| 302 | |
|---|
| 303 | /* Used on __finite compat alias. */ |
|---|
| 304 | #define __acosl __acosf128 |
|---|
| 305 | #define __acoshl __acoshf128 |
|---|
| 306 | #define __asinl __asinf128 |
|---|
| 307 | #define __atan2l __atan2f128 |
|---|
| 308 | #define __atanhl __atanhf128 |
|---|
| 309 | #define __coshl __coshf128 |
|---|
| 310 | #define __exp10l __exp10f128 |
|---|
| 311 | #define __expl __expf128 |
|---|
| 312 | #define __fmodl __fmodf128 |
|---|
| 313 | #define __gammal_r __gammaf128_r |
|---|
| 314 | #define __hypotl __hypotf128 |
|---|
| 315 | #define __j0l __j0f128 |
|---|
| 316 | #define __j1l __j1f128 |
|---|
| 317 | #define __jnl __jnf128 |
|---|
| 318 | #define __lgammal_r __lgammaf128_r |
|---|
| 319 | #define __log10l __log10f128 |
|---|
| 320 | #define __log2l __log2f128 |
|---|
| 321 | #define __logl __logf128 |
|---|
| 322 | #define __powl __powf128 |
|---|
| 323 | #define __remainderl __remainderf128 |
|---|
| 324 | #define __sinhl __sinhf128 |
|---|
| 325 | #define __y0l __y0f128 |
|---|
| 326 | #define __y1l __y1f128 |
|---|
| 327 | #define __ynl __ynf128 |
|---|
| 328 | |
|---|
| 329 | /* __nexttowardf128 is not _Float128 API. */ |
|---|
| 330 | #define __nexttowardl __nexttowardf128_do_not_use |
|---|
| 331 | #define nexttowardl nexttowardf128_do_not_use |
|---|
| 332 | |
|---|
| 333 | |
|---|
| 334 | /* public entry points. */ |
|---|
| 335 | #define asinhl asinhf128 |
|---|
| 336 | #define atanl atanf128 |
|---|
| 337 | #define cbrtl cbrtf128 |
|---|
| 338 | #define ceill ceilf128 |
|---|
| 339 | #define copysignl copysignf128 |
|---|
| 340 | #define cosl cosf128 |
|---|
| 341 | #define erfcl erfcf128 |
|---|
| 342 | #define erfl erff128 |
|---|
| 343 | #define expl expf128 |
|---|
| 344 | #define expm1l expm1f128 |
|---|
| 345 | #define fabsl fabsf128 |
|---|
| 346 | #define fdiml fdimf128 |
|---|
| 347 | #define finitel finitef128_do_not_use |
|---|
| 348 | #define floorl floorf128 |
|---|
| 349 | #define fmal fmaf128 |
|---|
| 350 | #define fmaxl fmaxf128 |
|---|
| 351 | #define fminl fminf128 |
|---|
| 352 | #define frexpl frexpf128 |
|---|
| 353 | #define getpayloadl getpayloadf128 |
|---|
| 354 | #define isinfl isinff128_do_not_use |
|---|
| 355 | #define isnanl isnanf128_do_not_use |
|---|
| 356 | #define ldexpl ldexpf128 |
|---|
| 357 | #define llrintl llrintf128 |
|---|
| 358 | #define llroundl llroundf128 |
|---|
| 359 | #define log1pl log1pf128 |
|---|
| 360 | #define logbl logbf128 |
|---|
| 361 | #define logl logf128 |
|---|
| 362 | #define lrintl lrintf128 |
|---|
| 363 | #define lroundl lroundf128 |
|---|
| 364 | #define modfl modff128 |
|---|
| 365 | #define nanl nanf128 |
|---|
| 366 | #define nearbyintl nearbyintf128 |
|---|
| 367 | #define nextafterl nextafterf128 |
|---|
| 368 | #define nextdownl nextdownf128 |
|---|
| 369 | #define nextupl nextupf128 |
|---|
| 370 | #define remquol remquof128 |
|---|
| 371 | #define rintl rintf128 |
|---|
| 372 | #define roundevenl roundevenf128 |
|---|
| 373 | #define roundl roundf128 |
|---|
| 374 | #define scalbnl scalbnf128 |
|---|
| 375 | #define sincosl sincosf128 |
|---|
| 376 | #define sinl sinf128 |
|---|
| 377 | #define sqrtl sqrtf128 |
|---|
| 378 | #define tanhl tanhf128 |
|---|
| 379 | #define tanl tanf128 |
|---|
| 380 | #define totalorderl totalorderf128 |
|---|
| 381 | #define totalordermagl totalordermagf128 |
|---|
| 382 | #define truncl truncf128 |
|---|
| 383 | |
|---|
| 384 | |
|---|
| 385 | /* misc internal renames. */ |
|---|
| 386 | #define __builtin_fmal __builtin_fmaf128 |
|---|
| 387 | #define __expl_table __expf128_table |
|---|
| 388 | #define __gamma_productl __gamma_productf128 |
|---|
| 389 | #define __kernel_cosl __kernel_cosf128 |
|---|
| 390 | #define __kernel_rem_pio2l __kernel_rem_pio2f128 |
|---|
| 391 | #define __kernel_sincosl __kernel_sincosf128 |
|---|
| 392 | #define __kernel_sinl __kernel_sinf128 |
|---|
| 393 | #define __kernel_tanl __kernel_tanf128 |
|---|
| 394 | #define __lgamma_negl __lgamma_negf128 |
|---|
| 395 | #define __lgamma_productl __lgamma_productf128 |
|---|
| 396 | #define __mpn_extract_float128 |
|---|
| 397 | #define __sincosl_table __sincosf128_table |
|---|
| 398 | #define mul_splitl mul_splitf128 |
|---|
| 399 | |
|---|
| 400 | /* Builtin renames. */ |
|---|
| 401 | #define __builtin_copysignl __builtin_copysignf128 |
|---|
| 402 | #define __builtin_signbitl __builtin_signbit |
|---|
| 403 | #define __builtin_nearbyintl __builtin_nearbyintf128 |
|---|
| 404 | #define __builtin_rintl __builtin_rintf128 |
|---|
| 405 | #define __builtin_floorl __builtin_floorf128 |
|---|
| 406 | #define __builtin_ceill __builtin_ceilf128 |
|---|
| 407 | #define __builtin_truncl __builtin_truncf128 |
|---|
| 408 | #define __builtin_roundl __builtin_roundf128 |
|---|
| 409 | #define __builtin_roundevenl __builtin_roundevenf128 |
|---|
| 410 | #define __builtin_lrintl __builtin_lrintf128 |
|---|
| 411 | #define __builtin_llrintl __builtin_llrintf128 |
|---|
| 412 | #define __builtin_copysignl __builtin_copysignf128 |
|---|
| 413 | #define __builtin_fabsl __builtin_fabsf128 |
|---|
| 414 | |
|---|
| 415 | /* Get the constant suffix from bits/floatn-compat.h. */ |
|---|
| 416 | #define L(x) __f128 (x) |
|---|
| 417 | |
|---|
| 418 | static inline void |
|---|
| 419 | mul_splitf128 (_Float128 *hi, _Float128 *lo, _Float128 x, _Float128 y) |
|---|
| 420 | { |
|---|
| 421 | #ifdef __FP_FAST_FMAF128 |
|---|
| 422 | /* Fast built-in fused multiply-add. */ |
|---|
| 423 | *hi = x * y; |
|---|
| 424 | *lo = __builtin_fmal (x, y, -*hi); |
|---|
| 425 | #else |
|---|
| 426 | /* Apply Dekker's algorithm. */ |
|---|
| 427 | *hi = x * y; |
|---|
| 428 | # define C ((1LL << (FLT128_MANT_DIG + 1) / 2) + 1) |
|---|
| 429 | _Float128 x1 = x * C; |
|---|
| 430 | _Float128 y1 = y * C; |
|---|
| 431 | # undef C |
|---|
| 432 | x1 = (x - x1) + x1; |
|---|
| 433 | y1 = (y - y1) + y1; |
|---|
| 434 | _Float128 x2 = x - x1; |
|---|
| 435 | _Float128 y2 = y - y1; |
|---|
| 436 | *lo = (((x1 * y1 - *hi) + x1 * y2) + x2 * y1) + x2 * y2; |
|---|
| 437 | #endif |
|---|
| 438 | } |
|---|
| 439 | |
|---|
