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
18using LlvmLibcExpm1Test = LIBC_NAMESPACE::testing::FPTest<double>;
19
20TEST_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
42using namespace LIBC_NAMESPACE::testing;
43
44TEST_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
51TEST_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
58TEST_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

source code of libc/test/src/math/smoke/expm1_test.cpp