1 | //===-- Unittests for e^x - 1 ---------------------------------------------===// |
2 | // |
3 | // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. |
4 | // See https://llvm.org/LICENSE.txt for license information. |
5 | // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception |
6 | // |
7 | //===----------------------------------------------------------------------===// |
8 | |
9 | #include "hdr/math_macros.h" |
10 | #include "src/__support/FPUtil/FPBits.h" |
11 | #include "src/__support/libc_errno.h" |
12 | #include "src/math/expm1.h" |
13 | #include "test/UnitTest/FPMatcher.h" |
14 | #include "test/UnitTest/Test.h" |
15 | |
16 | #include <stdint.h> |
17 | |
18 | using LlvmLibcExpm1Test = LIBC_NAMESPACE::testing::FPTest<double>; |
19 | |
20 | TEST_F(LlvmLibcExpm1Test, SpecialNumbers) { |
21 | EXPECT_FP_EQ_WITH_EXCEPTION(aNaN, LIBC_NAMESPACE::expm1(sNaN), FE_INVALID); |
22 | EXPECT_MATH_ERRNO(0); |
23 | |
24 | EXPECT_FP_EQ(aNaN, LIBC_NAMESPACE::expm1(aNaN)); |
25 | EXPECT_FP_EQ(inf, LIBC_NAMESPACE::expm1(inf)); |
26 | EXPECT_FP_EQ_ALL_ROUNDING(-1.0, LIBC_NAMESPACE::expm1(neg_inf)); |
27 | EXPECT_FP_EQ_WITH_EXCEPTION(inf, LIBC_NAMESPACE::expm1(0x1.0p20), |
28 | FE_OVERFLOW); |
29 | EXPECT_FP_EQ_ALL_ROUNDING(zero, LIBC_NAMESPACE::expm1(zero)); |
30 | EXPECT_FP_EQ_ALL_ROUNDING(neg_zero, LIBC_NAMESPACE::expm1(neg_zero)); |
31 | // |x| < 2^-53, expm1(x) = x |
32 | EXPECT_FP_EQ(-0x1.23456789abcdep-55, |
33 | LIBC_NAMESPACE::expm1(-0x1.23456789abcdep-55)); |
34 | EXPECT_FP_EQ(0x1.23456789abcdep-55, |
35 | LIBC_NAMESPACE::expm1(0x1.23456789abcdep-55)); |
36 | // log(2^-54) |
37 | EXPECT_FP_EQ(-1.0, LIBC_NAMESPACE::expm1(-0x1.2b708872320e2p5)); |
38 | } |
39 | |
40 | #ifdef LIBC_TEST_FTZ_DAZ |
41 | |
42 | using namespace LIBC_NAMESPACE::testing; |
43 | |
44 | TEST_F(LlvmLibcExpm1Test, FTZMode) { |
45 | ModifyMXCSR mxcsr(FTZ); |
46 | |
47 | EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::expm1(min_denormal)); |
48 | EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::expm1(max_denormal)); |
49 | } |
50 | |
51 | TEST_F(LlvmLibcExpm1Test, DAZMode) { |
52 | ModifyMXCSR mxcsr(DAZ); |
53 | |
54 | EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::expm1(min_denormal)); |
55 | EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::expm1(max_denormal)); |
56 | } |
57 | |
58 | TEST_F(LlvmLibcExpm1Test, FTZDAZMode) { |
59 | ModifyMXCSR mxcsr(FTZ | DAZ); |
60 | |
61 | EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::expm1(min_denormal)); |
62 | EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::expm1(max_denormal)); |
63 | } |
64 | |
65 | #endif |
66 | |