+2024-01-18 Bruno Haible <bruno@clisp.org>
+
+ log10l: Work around bug on NetBSD 9.3/arm64.
+ * m4/log10l.m4 (gl_FUNC_LOG10L_WORKS): Add another test case, taken from
+ the unit test.
+
2024-01-18 Bruno Haible <bruno@clisp.org>
doc: Update regarding <math.h> *l functions on NetBSD 9.3.
-# log10l.m4 serial 14
+# log10l.m4 serial 15
dnl Copyright (C) 2011-2024 Free Software Foundation, Inc.
dnl This file is free software; the Free Software Foundation
dnl gives unlimited permission to copy and/or distribute it,
int main (int argc, char *argv[])
{
long double (* volatile my_log10l) (long double) = argc ? log10l : dummy;
+ const long double TWO_LDBL_MANT_DIG = /* 2^LDBL_MANT_DIG */
+ (long double) (1U << ((LDBL_MANT_DIG - 1) / 5))
+ * (long double) (1U << ((LDBL_MANT_DIG - 1 + 1) / 5))
+ * (long double) (1U << ((LDBL_MANT_DIG - 1 + 2) / 5))
+ * (long double) (1U << ((LDBL_MANT_DIG - 1 + 3) / 5))
+ * (long double) (1U << ((LDBL_MANT_DIG - 1 + 4) / 5));
int result = 0;
/* Dummy call, to trigger the AIX 5.1 bug. */
gx = 0.6L;
if (!(gy + gy == gy))
result |= 1;
}
- /* This test fails on musl 1.2.2/arm64, musl 1.2.2/s390x, NetBSD 9.3. */
+ /* This test fails on musl 1.2.2/arm64, musl 1.2.2/s390x,
+ NetBSD 9.3 except arm64. */
{
- const long double TWO_LDBL_MANT_DIG = /* 2^LDBL_MANT_DIG */
- (long double) (1U << ((LDBL_MANT_DIG - 1) / 5))
- * (long double) (1U << ((LDBL_MANT_DIG - 1 + 1) / 5))
- * (long double) (1U << ((LDBL_MANT_DIG - 1 + 2) / 5))
- * (long double) (1U << ((LDBL_MANT_DIG - 1 + 3) / 5))
- * (long double) (1U << ((LDBL_MANT_DIG - 1 + 4) / 5));
long double x = 7.90097792256024576L;
long double err = (my_log10l (x) + my_log10l (1.0L / x)) * TWO_LDBL_MANT_DIG;
if (!(err >= -100.0L && err <= 100.0L))
result |= 2;
}
+ /* This test fails on NetBSD 9.3/arm64. */
+ {
+ long double x = 13.53398352203022253L;
+ long double err = (my_log10l (x) + my_log10l (1.0L / x)) * TWO_LDBL_MANT_DIG;
+ if (!(err >= -100.0L && err <= 100.0L))
+ result |= 4;
+ }
return result;
}
]])],