mirror of
https://git.savannah.gnu.org/git/gnulib.git
synced 2026-10-24 00:48:56 +03:00
f1d55e9f84
Done through sed -e 's/SIZEOF \([(][^()]*[)]\)/countof \1/g' * lib/boot-time.c: Include <stdcountof.h>. * lib/readutmp.c: Include <stdcountof.h>. * lib/boot-time-aux.h (get_linux_boot_time_fallback, get_openbsd_boot_time, get_windows_boot_time): Use countof. * lib/duplocale.c: Include <stdcountof.h>. (SIZEOF): Remove macro. (duplocale): Use countof instead. * lib/fatal-signal.c: Include <stdcountof.h>. (SIZEOF): Remove macro. (num_fatal_signals, actions_allocated): Use countof instead. * lib/iconv_open.c: Include <stdcountof.h>. (SIZEOF): Remove macro. (rpl_iconv_open): Use countof instead. * lib/term-style-control.c: Include <stdcountof.h>. (SIZEOF): Remove macro. (num_job_control_signals): Use countof instead. * lib/uniconv/u16-conv-to-enc.c: Include <stdcountof.h>. (SIZEOF): Remove macro. * lib/uniconv/u32-conv-to-enc.c: Include <stdcountof.h>. (SIZEOF): Remove macro. * lib/uniconv/u-conv-to-enc.h (FUNC): Use countof instead. * lib/uniconv/u16-strconv-to-enc.c: Include <stdcountof.h>. (SIZEOF): Remove macro. * lib/uniconv/u32-strconv-to-enc.c: Include <stdcountof.h>. (SIZEOF): Remove macro. * lib/uniconv/u-strconv-to-enc.h (FUNC): Use countof instead. * lib/uniname/uniname.c: Include <stdcountof.h>. (SIZEOF): Remove macro. (unicode_name_word, unicode_name_word_lookup, unicode_code_to_index, unicode_index_to_code, unicode_character_name, unicode_name_character): Use countof instead. * lib/uniwidth/width.c: Include <stdcountof.h>. (SIZEOF): Remove macro. (uc_width): Use countof instead. * lib/wait-process.c: Include <stdcountof.h>. (SIZEOF): Remove macro. (slaves_allocated): Use countof instead. * lib/gen-uni-tables.c (countof): Renamed from SIZEOF. * modules/boot-time (Depends-on): Add stdcountof-h. * modules/readutmp (Depends-on): Likewise. * modules/duplocale (Depends-on): Likewise. * modules/fatal-signal (Depends-on): Likewise. * modules/iconv_open (Depends-on): Likewise. * modules/term-style-control (Depends-on): Likewise. * modules/uniconv/u16-conv-to-enc (Depends-on): Likewise. * modules/uniconv/u32-conv-to-enc (Depends-on): Likewise. * modules/uniconv/u16-strconv-to-enc (Depends-on): Likewise. * modules/uniconv/u32-strconv-to-enc (Depends-on): Likewise. * modules/uniname/uniname (Depends-on): Likewise. * modules/uniwidth/width (Depends-on): Likewise. * modules/wait-process (Depends-on): Likewise. * tests/macros.h (SIZEOF): Remove macro. * tests/**/*.[hc]: Use countof instead of SIZEOF. Include <stdcountof.h> as needed. * modules/**/*-tests (Depends-on): Add stdcountof-h if needed.
89 lines
2.9 KiB
C
89 lines
2.9 KiB
C
/* Test of exp*() function family.
|
|
Copyright (C) 2012-2026 Free Software Foundation, Inc.
|
|
|
|
This program is free software: you can redistribute it and/or modify
|
|
it under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
This program is distributed in the hope that it will be useful,
|
|
but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
GNU General Public License for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with this program. If not, see <https://www.gnu.org/licenses/>. */
|
|
|
|
static void
|
|
test_function (void)
|
|
{
|
|
const DOUBLE TWO_MANT_DIG =
|
|
/* Assume MANT_DIG <= 5 * 31.
|
|
Use the identity
|
|
n = floor(n/5) + floor((n+1)/5) + ... + floor((n+4)/5). */
|
|
(DOUBLE) (1U << ((MANT_DIG - 1) / 5))
|
|
* (DOUBLE) (1U << ((MANT_DIG - 1 + 1) / 5))
|
|
* (DOUBLE) (1U << ((MANT_DIG - 1 + 2) / 5))
|
|
* (DOUBLE) (1U << ((MANT_DIG - 1 + 3) / 5))
|
|
* (DOUBLE) (1U << ((MANT_DIG - 1 + 4) / 5));
|
|
|
|
/* Randomized tests. */
|
|
{
|
|
/* Error bound, in ulps. */
|
|
const DOUBLE err_bound =
|
|
(sizeof (DOUBLE) > sizeof (double) ?
|
|
#if defined __i386__ && defined __FreeBSD__
|
|
/* On FreeBSD/x86 6.4, the 'long double' type really has only 53 bits of
|
|
precision in the compiler but 64 bits of precision at runtime. See
|
|
<https://lists.gnu.org/r/bug-gnulib/2008-07/msg00063.html>.
|
|
The compiler has truncated all 'long double' literals in expl.c to
|
|
53 bits of precision. */
|
|
L_(1350.0)
|
|
#else
|
|
L_(4.0)
|
|
#endif
|
|
: L_(2.0));
|
|
|
|
for (int i = 0; i < countof (RANDOM); i++)
|
|
{
|
|
DOUBLE x = L_(16.0) * RANDOM[i]; /* 0.0 <= x <= 16.0 */
|
|
DOUBLE y = EXP (x);
|
|
DOUBLE z = EXP (- x);
|
|
DOUBLE err = y * z - L_(1.0);
|
|
ASSERT (y > L_(0.0));
|
|
ASSERT (z > L_(0.0));
|
|
ASSERT (err > - err_bound / TWO_MANT_DIG
|
|
&& err < err_bound / TWO_MANT_DIG);
|
|
}
|
|
}
|
|
|
|
{
|
|
/* Error bound, in ulps. */
|
|
const DOUBLE err_bound =
|
|
(sizeof (DOUBLE) > sizeof (double) ?
|
|
#if defined __i386__ && defined __NetBSD__
|
|
L_(25000.0)
|
|
#elif defined __i386__ && defined __FreeBSD__
|
|
L_(2400.0)
|
|
#else
|
|
L_(24.0)
|
|
#endif
|
|
: L_(24.0));
|
|
|
|
for (int i = 0; i < countof (RANDOM) / 5; i++)
|
|
for (int j = 0; j < countof (RANDOM) / 5; j++)
|
|
{
|
|
DOUBLE x = L_(32.0) * RANDOM[i] - L_(16.0); /* -16.0 <= x <= 16.0 */
|
|
DOUBLE y = L_(32.0) * RANDOM[j] - L_(16.0); /* -16.0 <= y <= 16.0 */
|
|
DOUBLE z = - x - y;
|
|
/* Approximately x + y + z = 0. */
|
|
DOUBLE err = EXP (x) * EXP (y) * EXP (z) - L_(1.0);
|
|
ASSERT (err > - err_bound / TWO_MANT_DIG
|
|
&& err < err_bound / TWO_MANT_DIG);
|
|
}
|
|
}
|
|
}
|
|
|
|
volatile DOUBLE x;
|
|
DOUBLE y;
|