From eb0d82c77863b5e52758f012da09c3ef6480a7a9 Mon Sep 17 00:00:00 2001 From: Bruno Haible Date: Thu, 18 Jan 2024 19:46:53 +0100 Subject: [PATCH] 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. --- ChangeLog | 6 ++++++ m4/log10l.m4 | 24 ++++++++++++++++-------- 2 files changed, 22 insertions(+), 8 deletions(-) diff --git a/ChangeLog b/ChangeLog index ea0c6970f6..37924abdb8 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,9 @@ +2024-01-18 Bruno Haible + + 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 doc: Update regarding *l functions on NetBSD 9.3. diff --git a/m4/log10l.m4 b/m4/log10l.m4 index aefa2df029..ee352b11f6 100644 --- a/m4/log10l.m4 +++ b/m4/log10l.m4 @@ -1,4 +1,4 @@ -# 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, @@ -115,6 +115,12 @@ long double gy; 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; @@ -126,19 +132,21 @@ int main (int argc, char *argv[]) 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; } ]])], -- 2.39.5