summaryrefslogtreecommitdiff
path: root/lib/libm/__rem_pio2l.c
diff options
context:
space:
mode:
authorKacper <kacper@mail.openlinux.dev>2025-12-25 19:24:38 +0100
committerKacper <kacper@mail.openlinux.dev>2025-12-25 20:35:03 +0100
commita984eb367c032dbe2577f01238c3d1268526be70 (patch)
tree437fef40379b2758b129ccea39df3570fa2d145e /lib/libm/__rem_pio2l.c
parent8834571b202cf4dc9c649cfb096c213b6ecf1566 (diff)
Clang-tidy fixes
Diffstat (limited to 'lib/libm/__rem_pio2l.c')
-rw-r--r--lib/libm/__rem_pio2l.c66
1 files changed, 27 insertions, 39 deletions
diff --git a/lib/libm/__rem_pio2l.c b/lib/libm/__rem_pio2l.c
index f5d28dde..7c980ba4 100644
--- a/lib/libm/__rem_pio2l.c
+++ b/lib/libm/__rem_pio2l.c
@@ -27,9 +27,7 @@ static const long double toint = 1.5 / LDBL_EPSILON;
#if LDBL_MANT_DIG == 64
/* u ~< 0x1p25*pi/2 */
-#define SMALL(u) \
- ((((u).i.se & 0x7fffU) << 16 | (u).i.m >> 48) < \
- ((0x3fff + 25) << 16 | 0x921f >> 1 | 0x8000))
+#define SMALL(u) ((((u).i.se & 0x7fffU) << 16 | (u).i.m >> 48) < ((0x3fff + 25) << 16 | 0x921f >> 1 | 0x8000))
#define QUOBITS(x) ((uint32_t)(int32_t)(x) & 0x7fffffff)
#define ROUND1 22
#define ROUND2 61
@@ -46,47 +44,37 @@ static const long double toint = 1.5 / LDBL_EPSILON;
*/
static const double pio2_1 = 1.57079632679597125389e+00, /* 0x3FF921FB,
0x54444000 */
- pio2_2 = -1.07463465549783099519e-12, /* -0x12e7b967674000.0p-92 */
- pio2_3 = 6.36831716351370313614e-25; /* 0x18a2e037074000.0p-133 */
-static const long double
- pio4 = 0x1.921fb54442d1846ap-1L,
- invpio2 = 6.36619772367581343076e-01L, /* 0xa2f9836e4e44152a.0p-64
- */
- pio2_1t = -1.07463465549719416346e-12L, /* -0x973dcb3b399d747f.0p-103 */
- pio2_2t = 6.36831716351095013979e-25L, /* 0xc51701b839a25205.0p-144 */
- pio2_3t = -2.75299651904407171810e-37L; /* -0xbb5bf6c7ddd660ce.0p-185 */
+ pio2_2 = -1.07463465549783099519e-12, /* -0x12e7b967674000.0p-92 */
+ pio2_3 = 6.36831716351370313614e-25; /* 0x18a2e037074000.0p-133 */
+static const long double pio4 = 0x1.921fb54442d1846ap-1L,
+ invpio2 = 6.36619772367581343076e-01L, /* 0xa2f9836e4e44152a.0p-64
+ */
+ pio2_1t = -1.07463465549719416346e-12L, /* -0x973dcb3b399d747f.0p-103 */
+ pio2_2t = 6.36831716351095013979e-25L, /* 0xc51701b839a25205.0p-144 */
+ pio2_3t = -2.75299651904407171810e-37L; /* -0xbb5bf6c7ddd660ce.0p-185 */
#elif LDBL_MANT_DIG == 113
/* u ~< 0x1p45*pi/2 */
-#define SMALL(u) \
- (((u.i.se & 0x7fffU) << 16 | u.i.top) < ((0x3fff + 45) << 16 | 0x921f))
+#define SMALL(u) (((u.i.se & 0x7fffU) << 16 | u.i.top) < ((0x3fff + 45) << 16 | 0x921f))
#define QUOBITS(x) ((uint32_t)(int64_t)x & 0x7fffffff)
#define ROUND1 51
#define ROUND2 119
#define NX 5
#define NY 3
-static const long double
- pio4 = 0x1.921fb54442d18469898cc51701b8p-1L,
- invpio2 =
- 6.3661977236758134307553505349005747e-01L, /* 0x145f306dc9c882a53f84eafa3ea6a.0p-113
- */
- pio2_1 =
- 1.5707963267948966192292994253909555e+00L, /* 0x1921fb54442d18469800000000000.0p-112
- */
- pio2_1t =
- 2.0222662487959507323996846200947577e-21L, /* 0x13198a2e03707344a4093822299f3.0p-181
- */
- pio2_2 =
- 2.0222662487959507323994779168837751e-21L, /* 0x13198a2e03707344a400000000000.0p-181
- */
- pio2_2t =
- 2.0670321098263988236496903051604844e-43L, /* 0x127044533e63a0105df531d89cd91.0p-254
- */
- pio2_3 =
- 2.0670321098263988236499468110329591e-43L, /* 0x127044533e63a0105e00000000000.0p-254
- */
- pio2_3t =
- -2.5650587247459238361625433492959285e-65L; /* -0x159c4ec64ddaeb5f78671cbfb2210.0p-327
- */
+static const long double pio4 = 0x1.921fb54442d18469898cc51701b8p-1L,
+ invpio2 = 6.3661977236758134307553505349005747e-01L, /* 0x145f306dc9c882a53f84eafa3ea6a.0p-113
+ */
+ pio2_1 = 1.5707963267948966192292994253909555e+00L, /* 0x1921fb54442d18469800000000000.0p-112
+ */
+ pio2_1t = 2.0222662487959507323996846200947577e-21L, /* 0x13198a2e03707344a4093822299f3.0p-181
+ */
+ pio2_2 = 2.0222662487959507323994779168837751e-21L, /* 0x13198a2e03707344a400000000000.0p-181
+ */
+ pio2_2t = 2.0670321098263988236496903051604844e-43L, /* 0x127044533e63a0105df531d89cd91.0p-254
+ */
+ pio2_3 = 2.0670321098263988236499468110329591e-43L, /* 0x127044533e63a0105e00000000000.0p-254
+ */
+ pio2_3t = -2.5650587247459238361625433492959285e-65L; /* -0x159c4ec64ddaeb5f78671cbfb2210.0p-327
+ */
#endif
int __rem_pio2l(long double x, long double *y)
@@ -131,8 +119,8 @@ int __rem_pio2l(long double x, long double *y)
ey = u.i.se & 0x7fff;
if (ex - ey > ROUND2) { /* 3rd iteration, good to
180/316 bits */
- t = r; /* will cover all possible cases (not
- verified for ld128) */
+ t = r; /* will cover all possible cases (not
+ verified for ld128) */
w = fn * pio2_3;
r = t - w;
w = fn * pio2_3t - ((t - r) - w);