forked from llvm/llvm-project
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
[libc][math][c23] Add {f,d}mul{l,f128} and f16mul{,f,l,f128} C23 math…
… functions Part of llvm#93566. Fixes llvm#94833.
- Loading branch information
1 parent
861a8ed
commit 5007943
Showing
49 changed files
with
1,242 additions
and
364 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,120 @@ | ||
| //===-- Division of IEEE 754 floating-point numbers -------------*- C++ -*-===// | ||
| // | ||
| // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. | ||
| // See https://llvm.org/LICENSE.txt for license information. | ||
| // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| // | ||
| //===----------------------------------------------------------------------===// | ||
|
|
||
| #ifndef LLVM_LIBC_SRC___SUPPORT_FPUTIL_GENERIC_MUL_H | ||
| #define LLVM_LIBC_SRC___SUPPORT_FPUTIL_GENERIC_MUL_H | ||
|
|
||
| #include "hdr/errno_macros.h" | ||
| #include "hdr/fenv_macros.h" | ||
| #include "src/__support/CPP/bit.h" | ||
| #include "src/__support/CPP/type_traits.h" | ||
| #include "src/__support/FPUtil/BasicOperations.h" | ||
| #include "src/__support/FPUtil/FEnvImpl.h" | ||
| #include "src/__support/FPUtil/FPBits.h" | ||
| #include "src/__support/FPUtil/dyadic_float.h" | ||
| #include "src/__support/macros/attributes.h" | ||
| #include "src/__support/macros/config.h" | ||
| #include "src/__support/macros/optimization.h" | ||
|
|
||
| namespace LIBC_NAMESPACE_DECL { | ||
| namespace fputil::generic { | ||
|
|
||
| template <typename OutType, typename InType> | ||
| LIBC_INLINE cpp::enable_if_t<cpp::is_floating_point_v<OutType> && | ||
| cpp::is_floating_point_v<InType> && | ||
| sizeof(OutType) <= sizeof(InType), | ||
| OutType> | ||
| mul(InType x, InType y) { | ||
| using OutFPBits = FPBits<OutType>; | ||
| using OutStorageType = typename OutFPBits::StorageType; | ||
| using InFPBits = FPBits<InType>; | ||
| using InStorageType = typename InFPBits::StorageType; | ||
| // The product of two p-digit numbers is a 2p-digit number. | ||
| using DyadicFloat = DyadicFloat<cpp::bit_ceil( | ||
| 2 * static_cast<size_t>(InFPBits::FRACTION_LEN))>; | ||
| using DyadicMantissaType = typename DyadicFloat::MantissaType; | ||
|
|
||
| InFPBits x_bits(x); | ||
| InFPBits y_bits(y); | ||
|
|
||
| Sign result_sign = x_bits.sign() == y_bits.sign() ? Sign::POS : Sign::NEG; | ||
|
|
||
| if (LIBC_UNLIKELY(x_bits.is_inf_or_nan() || y_bits.is_inf_or_nan() || | ||
| x_bits.is_zero() || y_bits.is_zero())) { | ||
| if (x_bits.is_nan() || y_bits.is_nan()) { | ||
| if (x_bits.is_signaling_nan() || y_bits.is_signaling_nan()) | ||
| raise_except_if_required(FE_INVALID); | ||
|
|
||
| if (x_bits.is_quiet_nan()) { | ||
| InStorageType x_payload = static_cast<InStorageType>(getpayload(x)); | ||
| if ((x_payload & ~(OutFPBits::FRACTION_MASK >> 1)) == 0) | ||
| return OutFPBits::quiet_nan(x_bits.sign(), | ||
| static_cast<OutStorageType>(x_payload)) | ||
| .get_val(); | ||
| } | ||
|
|
||
| if (y_bits.is_quiet_nan()) { | ||
| InStorageType y_payload = static_cast<InStorageType>(getpayload(y)); | ||
| if ((y_payload & ~(OutFPBits::FRACTION_MASK >> 1)) == 0) | ||
| return OutFPBits::quiet_nan(y_bits.sign(), | ||
| static_cast<OutStorageType>(y_payload)) | ||
| .get_val(); | ||
| } | ||
|
|
||
| return OutFPBits::quiet_nan().get_val(); | ||
| } | ||
|
|
||
| if (x_bits.is_inf()) { | ||
| if (y_bits.is_zero()) { | ||
| set_errno_if_required(EDOM); | ||
| raise_except_if_required(FE_INVALID); | ||
| return OutFPBits::quiet_nan().get_val(); | ||
| } | ||
|
|
||
| return OutFPBits::inf(result_sign).get_val(); | ||
| } | ||
|
|
||
| if (y_bits.is_inf()) { | ||
| if (x_bits.is_zero()) { | ||
| set_errno_if_required(EDOM); | ||
| raise_except_if_required(FE_INVALID); | ||
| return OutFPBits::quiet_nan().get_val(); | ||
| } | ||
|
|
||
| return OutFPBits::inf(result_sign).get_val(); | ||
| } | ||
|
|
||
| // Now either x or y is zero, and the other one is finite. | ||
| return OutFPBits::zero(result_sign).get_val(); | ||
| } | ||
|
|
||
| DyadicFloat xd(x); | ||
| DyadicFloat yd(y); | ||
|
|
||
| // The product is either in the [2, 4) range or the [1, 2) range. We initially | ||
| // set the exponent as if the product were in the [2, 4) range. If the product | ||
| // is in the [1, 2) range, then it has 1 leading zero bit and the DyadicFloat | ||
| // constructor will normalize it and adjust the exponent. | ||
| int result_exp = xd.get_unbiased_exponent() + yd.get_unbiased_exponent() + 1 - | ||
| DyadicMantissaType::BITS + 1; | ||
|
|
||
| // We use a single DyadicFloat type, which can fit the 2p-digit product, for | ||
| // both the normalized inputs and the product. So we need to right-shift the | ||
| // normalized input mantissas before multiplying them. | ||
| xd.shift_right(DyadicMantissaType::BITS / 2); | ||
| yd.shift_right(DyadicMantissaType::BITS / 2); | ||
| DyadicMantissaType product = xd.mantissa * yd.mantissa; | ||
|
|
||
| DyadicFloat result(result_sign, result_exp, product); | ||
| return result.template as<OutType, /*ShouldSignalExceptions=*/true>(); | ||
| } | ||
|
|
||
| } // namespace fputil::generic | ||
| } // namespace LIBC_NAMESPACE_DECL | ||
|
|
||
| #endif // LLVM_LIBC_SRC___SUPPORT_FPUTIL_GENERIC_MUL_H |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,21 @@ | ||
| //===-- Implementation header for dmulf128 ----------------------*- C++ -*-===// | ||
| // | ||
| // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. | ||
| // See https://llvm.org/LICENSE.txt for license information. | ||
| // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| // | ||
| //===----------------------------------------------------------------------===// | ||
|
|
||
| #ifndef LLVM_LIBC_SRC_MATH_DMULF128_H | ||
| #define LLVM_LIBC_SRC_MATH_DMULF128_H | ||
|
|
||
| #include "src/__support/macros/config.h" | ||
| #include "src/__support/macros/properties/types.h" | ||
|
|
||
| namespace LIBC_NAMESPACE_DECL { | ||
|
|
||
| double dmulf128(float128 x, float128 y); | ||
|
|
||
| } // namespace LIBC_NAMESPACE_DECL | ||
|
|
||
| #endif // LLVM_LIBC_SRC_MATH_DMULF128_H |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,20 @@ | ||
| //===-- Implementation header for dmull -------------------------*- C++ -*-===// | ||
| // | ||
| // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. | ||
| // See https://llvm.org/LICENSE.txt for license information. | ||
| // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| // | ||
| //===----------------------------------------------------------------------===// | ||
|
|
||
| #ifndef LLVM_LIBC_SRC_MATH_DMULL_H | ||
| #define LLVM_LIBC_SRC_MATH_DMULL_H | ||
|
|
||
| #include "src/__support/macros/config.h" | ||
|
|
||
| namespace LIBC_NAMESPACE_DECL { | ||
|
|
||
| double dmull(long double x, long double y); | ||
|
|
||
| } // namespace LIBC_NAMESPACE_DECL | ||
|
|
||
| #endif // LLVM_LIBC_SRC_MATH_DMULL_H |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,21 @@ | ||
| //===-- Implementation header for f16mul ------------------------*- C++ -*-===// | ||
| // | ||
| // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. | ||
| // See https://llvm.org/LICENSE.txt for license information. | ||
| // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| // | ||
| //===----------------------------------------------------------------------===// | ||
|
|
||
| #ifndef LLVM_LIBC_SRC_MATH_F16MUL_H | ||
| #define LLVM_LIBC_SRC_MATH_F16MUL_H | ||
|
|
||
| #include "src/__support/macros/config.h" | ||
| #include "src/__support/macros/properties/types.h" | ||
|
|
||
| namespace LIBC_NAMESPACE_DECL { | ||
|
|
||
| float16 f16mul(double x, double y); | ||
|
|
||
| } // namespace LIBC_NAMESPACE_DECL | ||
|
|
||
| #endif // LLVM_LIBC_SRC_MATH_F16MUL_H |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,21 @@ | ||
| //===-- Implementation header for f16mulf -----------------------*- C++ -*-===// | ||
| // | ||
| // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. | ||
| // See https://llvm.org/LICENSE.txt for license information. | ||
| // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception | ||
| // | ||
| //===----------------------------------------------------------------------===// | ||
|
|
||
| #ifndef LLVM_LIBC_SRC_MATH_F16MULF_H | ||
| #define LLVM_LIBC_SRC_MATH_F16MULF_H | ||
|
|
||
| #include "src/__support/macros/config.h" | ||
| #include "src/__support/macros/properties/types.h" | ||
|
|
||
| namespace LIBC_NAMESPACE_DECL { | ||
|
|
||
| float16 f16mulf(float x, float y); | ||
|
|
||
| } // namespace LIBC_NAMESPACE_DECL | ||
|
|
||
| #endif // LLVM_LIBC_SRC_MATH_F16MULF_H |
Oops, something went wrong.