1 /*
2 * SPDX-License-Identifier: BSD-3-Clause
3 *
4 * Copyright © 2022 Keith Packard
5 *
6 * Redistribution and use in source and binary forms, with or without
7 * modification, are permitted provided that the following conditions
8 * are met:
9 *
10 * 1. Redistributions of source code must retain the above copyright
11 * notice, this list of conditions and the following disclaimer.
12 *
13 * 2. Redistributions in binary form must reproduce the above
14 * copyright notice, this list of conditions and the following
15 * disclaimer in the documentation and/or other materials provided
16 * with the distribution.
17 *
18 * 3. Neither the name of the copyright holder nor the names of its
19 * contributors may be used to endorse or promote products derived
20 * from this software without specific prior written permission.
21 *
22 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
23 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
24 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
25 * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
26 * COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT,
27 * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
28 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
29 * SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
30 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT,
31 * STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
32 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED
33 * OF THE POSSIBILITY OF SUCH DAMAGE.
34 */
35
36 #include "math_ld.h"
37
38 #if !_HAVE_FAST_FMAL && (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113)
39
40 typedef long double FLOAT_T;
41
42 #define FMA fmal
43 #define NEXTAFTER nextafterl
44 #define LDEXP ldexpl
45 #define FREXP frexpl
46 #define SCALBN scalbnl
47 #define ILOGB ilogbl
48 #define COPYSIGN copysignl
49
50 #if LDBL_MANT_DIG == 64
51 #define SPLIT (0x1p32L + 1.0L)
52 #if __LDBL_MIN_EXP__ == -16382
53 #define FLOAT_DENORM_BIAS 0
54 #endif
55 #endif
56 #if LDBL_MANT_DIG == 113
57 #define SPLIT (0x1p57L + 1.0L)
58 #endif
59
60 #define FLOAT_MANT_DIG __LDBL_MANT_DIG__
61 #define FLOAT_MAX_EXP __LDBL_MAX_EXP__
62 #define FLOAT_MIN __LDBL_MIN__
63
64 static inline int
odd_mant(FLOAT_T x)65 odd_mant(FLOAT_T x)
66 {
67 union IEEEl2bits u;
68
69 u.e = x;
70 return u.bits.manl & 1;
71 }
72
73 static unsigned int
EXPONENT(FLOAT_T x)74 EXPONENT(FLOAT_T x)
75 {
76 union IEEEl2bits u;
77
78 u.e = x;
79 return u.bits.exp;
80 }
81
82 #ifdef PICOLIBC_LONG_DOUBLE_NOEXCEPT
83 #define feraiseexcept(x) ((void) (x))
84 #endif
85
86 #include "fma_inc.h"
87
88 _MATH_ALIAS_d_ddd(fma)
89
90 #endif
91
