Add test case from a recent glibc bug.
[gnulib.git] / tests / test-log1p-ieee.h
blob28f793f62a7361d60d656826b30d306e0f67c60e
1 /* Test of log1p*() function family.
2 Copyright (C) 2012-2018 Free Software Foundation, Inc.
4 This program is free software: you can redistribute it and/or modify
5 it under the terms of the GNU General Public License as published by
6 the Free Software Foundation; either version 3 of the License, or
7 (at your option) any later version.
9 This program is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 GNU General Public License for more details.
14 You should have received a copy of the GNU General Public License
15 along with this program. If not, see <https://www.gnu.org/licenses/>. */
17 static void
18 test_function (void)
20 /* [MX] shaded specification in POSIX. */
22 /* Negative arguments. */
23 ASSERT (ISNAN (LOG1P (- L_(1.00001))));
24 ASSERT (ISNAN (LOG1P (- L_(2.0))));
25 ASSERT (ISNAN (LOG1P (- L_(1000000.0))));
26 ASSERT (ISNAN (LOG1P (- INFINITY)));
28 /* NaN. */
29 ASSERT (ISNAN (LOG1P (NAN)));
31 /* Only rational value. */
33 DOUBLE z = LOG1P (L_(0.0));
34 ASSERT (z == L_(0.0));
35 ASSERT (!signbit (z));
38 DOUBLE z = LOG1P (MINUS_ZERO);
39 ASSERT (z == L_(0.0));
40 ASSERT (!!signbit (z) == !!signbit (MINUS_ZERO));
43 /* Infinity. */
44 ASSERT (LOG1P (INFINITY) == INFINITY);