[llvm-branch-commits] [libc] [libc] make sqrtf128 use emulated float128 (PR #207739)
via llvm-branch-commits
llvm-branch-commits at lists.llvm.org
Mon Jul 6 07:55:07 PDT 2026
https://github.com/Sukumarsawant updated https://github.com/llvm/llvm-project/pull/207739
>From c654b5d3c49df5f3c7d0f8f117c572f37123ce46 Mon Sep 17 00:00:00 2001
From: Sukumarsawant <sawantsukumar at gmail.com>
Date: Mon, 6 Jul 2026 19:42:15 +0530
Subject: [PATCH 1/3] make sqrtf128 use emulated float128
---
libc/shared/math/sqrtf128.h | 9 +++-----
libc/src/__support/math/CMakeLists.txt | 2 +-
libc/src/__support/math/sqrtf128.h | 15 ++++++-------
libc/src/math/generic/CMakeLists.txt | 1 +
libc/src/math/generic/sqrtf128.cpp | 5 ++++-
libc/src/math/sqrtf128.h | 5 ++++-
libc/test/shared/shared_math_test.cpp | 2 +-
libc/test/src/math/CMakeLists.txt | 1 +
libc/test/src/math/smoke/CMakeLists.txt | 1 +
libc/test/src/math/smoke/sqrtf128_test.cpp | 25 ++++++++++++----------
libc/test/src/math/sqrtf128_test.cpp | 9 ++++----
11 files changed, 41 insertions(+), 34 deletions(-)
diff --git a/libc/shared/math/sqrtf128.h b/libc/shared/math/sqrtf128.h
index 609b8921ab36b..ed0c1c5c63906 100644
--- a/libc/shared/math/sqrtf128.h
+++ b/libc/shared/math/sqrtf128.h
@@ -9,13 +9,12 @@
#ifndef LLVM_LIBC_SHARED_MATH_SQRTF128_H
#define LLVM_LIBC_SHARED_MATH_SQRTF128_H
-#include "include/llvm-libc-types/float128.h"
-
-#ifdef LIBC_TYPES_HAS_FLOAT128
-
#include "shared/libc_common.h"
+#include "src/__support/FPUtil/float128.h"
#include "src/__support/math/sqrtf128.h"
+using LIBC_NAMESPACE::fputil::Float128;
+
namespace LIBC_NAMESPACE_DECL {
namespace shared {
@@ -24,6 +23,4 @@ using math::sqrtf128;
} // namespace shared
} // namespace LIBC_NAMESPACE_DECL
-#endif // LIBC_TYPES_HAS_FLOAT128
-
#endif // LLVM_LIBC_SHARED_MATH_SQRTF128_H
diff --git a/libc/src/__support/math/CMakeLists.txt b/libc/src/__support/math/CMakeLists.txt
index a2699355a18e8..78b4f55cca47e 100644
--- a/libc/src/__support/math/CMakeLists.txt
+++ b/libc/src/__support/math/CMakeLists.txt
@@ -5394,7 +5394,7 @@ add_header_library(
libc.src.__support.common
libc.src.__support.macros.optimization
libc.src.__support.uint128
- libc.include.llvm-libc-types.float128
+ libc.src.__support.FPUtil.float128
)
add_header_library(
diff --git a/libc/src/__support/math/sqrtf128.h b/libc/src/__support/math/sqrtf128.h
index 14dcbec1acbb6..9710ffe92eef0 100644
--- a/libc/src/__support/math/sqrtf128.h
+++ b/libc/src/__support/math/sqrtf128.h
@@ -9,18 +9,17 @@
#ifndef LLVM_LIBC_SRC___SUPPORT_MATH_SQRTF128_H
#define LLVM_LIBC_SRC___SUPPORT_MATH_SQRTF128_H
-#include "include/llvm-libc-types/float128.h"
-
-#ifdef LIBC_TYPES_HAS_FLOAT128
-
#include "src/__support/CPP/bit.h"
#include "src/__support/FPUtil/FEnvImpl.h"
#include "src/__support/FPUtil/FPBits.h"
+#include "src/__support/FPUtil/float128.h"
#include "src/__support/FPUtil/rounding_mode.h"
#include "src/__support/common.h"
#include "src/__support/macros/optimization.h"
#include "src/__support/uint128.h"
+using LIBC_NAMESPACE::fputil::Float128;
+
// Compute sqrtf128 with correct rounding for all rounding modes using integer
// arithmetic by Alexei Sibidanov (sibid at uvic.ca):
// https://github.com/sibidanov/llvm-project/tree/as_sqrt_v2
@@ -280,9 +279,9 @@ LIBC_INLINE constexpr uint64_t rsqrt_approx(uint64_t m) {
} // namespace sqrtf128_internal
-LIBC_INLINE float128 sqrtf128(float128 x) {
+LIBC_INLINE Float128 sqrtf128(Float128 x) {
using namespace sqrtf128_internal;
- using FPBits = fputil::FPBits<float128>;
+ using FPBits = fputil::FPBits<Float128>;
// Get rounding mode.
uint32_t rm = fputil::get_round();
@@ -441,12 +440,10 @@ LIBC_INLINE float128 sqrtf128(float128 x) {
// if(frac) fputil::raise_except_if_required(FE_INEXACT);
v += static_cast<UInt128>(e2) << FPBits::FRACTION_LEN; // place exponent
- return cpp::bit_cast<float128>(v);
+ return cpp::bit_cast<Float128>(v);
}
} // namespace math
} // namespace LIBC_NAMESPACE_DECL
-#endif // LIBC_TYPES_HAS_FLOAT128
-
#endif // LLVM_LIBC_SRC___SUPPORT_MATH_SQRTF128_H
diff --git a/libc/src/math/generic/CMakeLists.txt b/libc/src/math/generic/CMakeLists.txt
index 6769f4fdd4dcb..2a24329b93a9c 100644
--- a/libc/src/math/generic/CMakeLists.txt
+++ b/libc/src/math/generic/CMakeLists.txt
@@ -2685,6 +2685,7 @@ add_entrypoint_object(
../sqrtf128.h
DEPENDS
libc.src.__support.math.sqrtf128
+ libc.src.__support.FPUtil.float128
)
add_entrypoint_object(
diff --git a/libc/src/math/generic/sqrtf128.cpp b/libc/src/math/generic/sqrtf128.cpp
index 47e4da97581ba..9134324ace67b 100644
--- a/libc/src/math/generic/sqrtf128.cpp
+++ b/libc/src/math/generic/sqrtf128.cpp
@@ -7,11 +7,14 @@
//===----------------------------------------------------------------------===//
#include "src/math/sqrtf128.h"
+#include "src/__support/FPUtil/float128.h"
#include "src/__support/math/sqrtf128.h"
+using LIBC_NAMESPACE::fputil::Float128;
+
namespace LIBC_NAMESPACE_DECL {
-LLVM_LIBC_FUNCTION(float128, sqrtf128, (float128 x)) {
+LLVM_LIBC_FUNCTION(Float128, sqrtf128, (Float128 x)) {
return math::sqrtf128(x);
}
diff --git a/libc/src/math/sqrtf128.h b/libc/src/math/sqrtf128.h
index 81f4c12eff90d..80604826c087c 100644
--- a/libc/src/math/sqrtf128.h
+++ b/libc/src/math/sqrtf128.h
@@ -9,12 +9,15 @@
#ifndef LLVM_LIBC_SRC_MATH_SQRTF128_H
#define LLVM_LIBC_SRC_MATH_SQRTF128_H
+#include "src/__support/FPUtil/float128.h"
#include "src/__support/macros/config.h"
#include "src/__support/macros/properties/types.h"
+using LIBC_NAMESPACE::fputil::Float128;
+
namespace LIBC_NAMESPACE_DECL {
-float128 sqrtf128(float128 x);
+Float128 sqrtf128(Float128 x);
} // namespace LIBC_NAMESPACE_DECL
diff --git a/libc/test/shared/shared_math_test.cpp b/libc/test/shared/shared_math_test.cpp
index c71e44f71b23a..ac8fd3f511d50 100644
--- a/libc/test/shared/shared_math_test.cpp
+++ b/libc/test/shared/shared_math_test.cpp
@@ -583,6 +583,7 @@ using LIBC_NAMESPACE::fputil::Float128;
TEST(LlvmLibcSharedMathTest, AllEmuFloat128) {
EXPECT_FP_EQ(Float128(0.0), LIBC_NAMESPACE::shared::ceilf128(Float128(0.0)));
+ EXPECT_FP_EQ(Float128(1.0), LIBC_NAMESPACE::shared::sqrtf128(Float128(1.0)));
}
#ifdef LIBC_TYPES_HAS_FLOAT128
@@ -608,7 +609,6 @@ TEST(LlvmLibcSharedMathTest, AllFloat128) {
EXPECT_FP_EQ(float128(0.0), LIBC_NAMESPACE::shared::logbf128(float128(1.0)));
EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::shared::dfmaf128(
float128(0.0), float128(0.0), float128(0.0)));
- EXPECT_FP_EQ(float128(1.0), LIBC_NAMESPACE::shared::sqrtf128(float128(1.0)));
EXPECT_FP_EQ(0.0, LIBC_NAMESPACE::shared::dsqrtf128(float128(0.0)));
EXPECT_EQ(0L, LIBC_NAMESPACE::shared::llogbf128(float128(1.0)));
diff --git a/libc/test/src/math/CMakeLists.txt b/libc/test/src/math/CMakeLists.txt
index c4944081bb2a2..a07ba41bafb82 100644
--- a/libc/test/src/math/CMakeLists.txt
+++ b/libc/test/src/math/CMakeLists.txt
@@ -1790,6 +1790,7 @@ add_fp_unittest(
SqrtTest.h
DEPENDS
libc.src.math.sqrtf128
+ libc.src.__support.FPUtil.float128
)
add_fp_unittest(
diff --git a/libc/test/src/math/smoke/CMakeLists.txt b/libc/test/src/math/smoke/CMakeLists.txt
index 7876270493e11..e5b49cb0bbe39 100644
--- a/libc/test/src/math/smoke/CMakeLists.txt
+++ b/libc/test/src/math/smoke/CMakeLists.txt
@@ -3600,6 +3600,7 @@ add_fp_unittest(
SqrtTest.h
DEPENDS
libc.src.math.sqrtf128
+ libc.src.__support.FPUtil.float128
)
add_fp_unittest(
diff --git a/libc/test/src/math/smoke/sqrtf128_test.cpp b/libc/test/src/math/smoke/sqrtf128_test.cpp
index 2fca74842d313..5dbc063b79946 100644
--- a/libc/test/src/math/smoke/sqrtf128_test.cpp
+++ b/libc/test/src/math/smoke/sqrtf128_test.cpp
@@ -8,14 +8,17 @@
#include "SqrtTest.h"
+#include "src/__support/FPUtil/float128.h"
#include "src/__support/uint128.h"
#include "src/math/sqrtf128.h"
-LIST_SQRT_TESTS(float128, LIBC_NAMESPACE::sqrtf128);
+using LIBC_NAMESPACE::fputil::Float128;
+
+LIST_SQRT_TESTS(Float128, LIBC_NAMESPACE::sqrtf128);
TEST_F(LlvmLibcSqrtTest, HardToRound) {
using LIBC_NAMESPACE::fputil::testing::RoundingMode;
- using FPBits = LIBC_NAMESPACE::fputil::FPBits<float128>;
+ using FPBits = LIBC_NAMESPACE::fputil::FPBits<Float128>;
// Since there is no exact half cases for square root I encode the
// round direction in the sign of the result. E.g. if the number is
@@ -23,7 +26,7 @@ TEST_F(LlvmLibcSqrtTest, HardToRound) {
// (the absolute value). Thus I can test not only hard to round
// cases for the round to nearest mode but also the directional
// modes.
- float128 HARD_TO_ROUND[][2] = {
+ Float128 HARD_TO_ROUND[][2] = {
{0x0.000000dee2f5b6a26c8f07f05442p-16382q,
-0x1.ddbd8763a617cff753e2a31083p-8204q},
{0x0.000000c86d174c5ad8ae54a548e7p-16382q,
@@ -85,9 +88,9 @@ TEST_F(LlvmLibcSqrtTest, HardToRound) {
0x1.ffffffffffffffffffffffffffffp+0q},
};
- auto rnd = [](float128 x, RoundingMode rm) -> float128 {
+ auto rnd = [](Float128 x, RoundingMode rm) -> Float128 {
bool is_neg = x < 0;
- float128 y = is_neg ? -x : x;
+ Float128 y = is_neg ? -x : x;
FPBits ybits(y);
if (is_neg &&
@@ -107,7 +110,7 @@ TEST_F(LlvmLibcSqrtTest, HardToRound) {
}
// Exact results for subnormal arguments
- float128 EXACT_SUBNORMAL[][2] = {
+ Float128 EXACT_SUBNORMAL[][2] = {
{0x0.0000000000000000000000000001p-16382q, 0x1p-8247q},
{0x0.0000000000000000000000000004p-16382q, 0x1p-8246q},
{0x0.0000000000001000000000000000p-16382q, 0x1p-8217q},
@@ -121,17 +124,17 @@ TEST_F(LlvmLibcSqrtTest, HardToRound) {
// Check exact cases starting from small numbers
for (unsigned k = 1; k < 100 * 100; ++k) {
unsigned k2 = k * k;
- float128 x = static_cast<float128>(k2);
- float128 y = static_cast<float128>(k);
+ Float128 x = static_cast<Float128>(k2);
+ Float128 y = static_cast<Float128>(k);
EXPECT_FP_EQ_ALL_ROUNDING(y, LIBC_NAMESPACE::sqrtf128(x));
};
// Then from the largest number.
uint64_t k0 = 101904826760412362ULL;
for (uint64_t k = k0; k > k0 - 10000; --k) {
- float128 k_f128 = static_cast<float128>(k);
- float128 x = k_f128 * k_f128;
- float128 y = static_cast<float128>(k);
+ Float128 k_f128 = static_cast<Float128>(k);
+ Float128 x = k_f128 * k_f128;
+ Float128 y = static_cast<Float128>(k);
EXPECT_FP_EQ_ALL_ROUNDING(y, LIBC_NAMESPACE::sqrtf128(x));
}
}
diff --git a/libc/test/src/math/sqrtf128_test.cpp b/libc/test/src/math/sqrtf128_test.cpp
index 25229f834d33c..4f9ea74d1e204 100644
--- a/libc/test/src/math/sqrtf128_test.cpp
+++ b/libc/test/src/math/sqrtf128_test.cpp
@@ -7,15 +7,16 @@
//===----------------------------------------------------------------------===//
#include "SqrtTest.h"
-
+#include "src/__support/FPUtil/float128.h"
+#include "src/__support/integer_literals.h"
#include "src/math/sqrtf128.h"
-#include "src/__support/integer_literals.h"
+using LIBC_NAMESPACE::fputil::Float128;
-LIST_SQRT_TESTS(float128, LIBC_NAMESPACE::sqrtf128)
+LIST_SQRT_TESTS(Float128, LIBC_NAMESPACE::sqrtf128)
TEST_F(LlvmLibcSqrtTest, SpecialInputs) {
- constexpr float128 INPUTS[] = {
+ constexpr Float128 INPUTS[] = {
0x0.000000dee2f5b6a26c8f07f05442p-16382q,
0x0.000000c86d174c5ad8ae54a548e7p-16382q,
0x0.000020ab15cfe0b8e488e128f535p-16382q,
>From 485884adcbf66b63205ad2f1283487de5954f31e Mon Sep 17 00:00:00 2001
From: Sukumarsawant <sawantsukumar at gmail.com>
Date: Mon, 6 Jul 2026 20:14:26 +0530
Subject: [PATCH 2/3] chore: casting to f128
---
libc/test/src/math/smoke/sqrtf128_test.cpp | 133 +++++++++++----------
libc/test/src/math/sqrtf128_test.cpp | 36 +++---
2 files changed, 87 insertions(+), 82 deletions(-)
diff --git a/libc/test/src/math/smoke/sqrtf128_test.cpp b/libc/test/src/math/smoke/sqrtf128_test.cpp
index 5dbc063b79946..8eafb9dc7ca6d 100644
--- a/libc/test/src/math/smoke/sqrtf128_test.cpp
+++ b/libc/test/src/math/smoke/sqrtf128_test.cpp
@@ -27,65 +27,65 @@ TEST_F(LlvmLibcSqrtTest, HardToRound) {
// cases for the round to nearest mode but also the directional
// modes.
Float128 HARD_TO_ROUND[][2] = {
- {0x0.000000dee2f5b6a26c8f07f05442p-16382q,
- -0x1.ddbd8763a617cff753e2a31083p-8204q},
- {0x0.000000c86d174c5ad8ae54a548e7p-16382q,
- 0x1.c507bb538940719890851ec1ca88p-8204q},
- {0x0.000020ab15cfe0b8e488e128f535p-16382q,
- -0x1.6dccb402560213bc0d62d62e910bp-8201q},
- {0x0.0000219e97732a9970f2511989bap-16382q,
- 0x1.73163d28be706f4b5052791e28a5p-8201q},
- {0x0.000026e477546ae99ef57066f9fdp-16382q,
- -0x1.8f20dd0d0c570a23ea59bc2bf009p-8201q},
- {0x0.00002d0f88d27a496b3e533f5067p-16382q,
- 0x1.ad9d4abe9f047225a7352bcc52c1p-8201q},
- {0x1.0000000000000000000000000001p+0q, 0x1p+0q},
- {0x1.0000000000000000000000000002p+0q,
- -0x1.0000000000000000000000000001p+0q},
- {0x1.0000000000000000000000000003p+0q,
- 0x1.0000000000000000000000000001p+0q},
- {0x1.0000000000000000000000000005p+0q,
- 0x1.0000000000000000000000000002p+0q},
- {0x1.0000000000000000000000000006p+0q,
- -0x1.0000000000000000000000000003p+0q},
- {0x1.1d4c381cbf3a0aa15b9aee344892p+0q,
- 0x1.0e408c3fadc5e64b449c63673f4bp+0q},
- {0x1.2af17a4ae6f93d11310c49c11b59p+0q,
- -0x1.14a3bdf0ea5231f12d421a5dbe33p+0q},
- {0x1.96f893bf29fb91e0fbe19a46d0c8p+0q,
- 0x1.42c6bf6202e66f2295807dee44d9p+0q},
- {0x1.97fb3839925b66804c429289cce8p+0q,
- -0x1.432d4049ac1c85a241f333d326e9p+0q},
- {0x1.be1d900eaeb1533f0f19cc15c7e6p+0q,
- 0x1.51f1715154da44f3bf11f3d96c2dp+0q},
- {0x1.c4f5074269525063a26051a0ad27p+0q,
- 0x1.54864e9b1daa4d9135ff00663366p+0q},
- {0x1.035cb5f298a801dc4be9b1f8cd97p+1q,
- -0x1.6c688775bffcb3f507ba11d0abb9p+0q},
- {0x1.274be02380427e709beab4dedeb4p+1q,
- -0x1.84d5763281f2318422392e506b1cp+0q},
- {0x1.64e797cfdbaa3f7e2f33279dbc6p+1q,
- 0x1.ab79b164e255b26eca00ff99cc99p+0q},
- {0x1.693a741358c9dac44a570a7e9f6cp+1q,
- 0x1.ae0e8eaeab25bb0c40ee0c2693d3p+0q},
- {0x1.8275db3fc4d822596047adcb71b9p+1q,
- -0x1.bcd2bfb653e37a5dbe0ccc2cd917p+0q},
- {0x1.83280bb98c4a7b88bd6f535899d9p+1q,
- 0x1.bd39409dfd1990dd6a7f8211bb27p+0q},
- {0x1.d78d8352b48608b510bfd5c75315p+1q,
- -0x1.eb5c420f15adce0ed2bde5a241cep+0q},
- {0x1.e3e4774f564b526edff84ce46668p+1q,
- 0x1.f1bf73c0523a19b4bb639c98c0b5p+0q},
- {0x1.fffffffffffffffffffffffffffap+1q,
- -0x1.fffffffffffffffffffffffffffdp+0q},
- {0x1.fffffffffffffffffffffffffffbp+1q,
- 0x1.fffffffffffffffffffffffffffdp+0q},
- {0x1.fffffffffffffffffffffffffffdp+1q,
- 0x1.fffffffffffffffffffffffffffep+0q},
- {0x1.fffffffffffffffffffffffffffep+1q,
- -0x1.ffffffffffffffffffffffffffffp+0q},
- {0x1.ffffffffffffffffffffffffffffp+1q,
- 0x1.ffffffffffffffffffffffffffffp+0q},
+ {Float128(0x0.000000dee2f5b6a26c8f07f05442p-16382q),
+ Float128(-0x1.ddbd8763a617cff753e2a31083p-8204q)},
+ {Float128(0x0.000000c86d174c5ad8ae54a548e7p-16382q),
+ Float128(0x1.c507bb538940719890851ec1ca88p-8204q)},
+ {Float128(0x0.000020ab15cfe0b8e488e128f535p-16382q),
+ Float128(-0x1.6dccb402560213bc0d62d62e910bp-8201q)},
+ {Float128(0x0.0000219e97732a9970f2511989bap-16382q),
+ Float128(0x1.73163d28be706f4b5052791e28a5p-8201q)},
+ {Float128(0x0.000026e477546ae99ef57066f9fdp-16382q),
+ Float128(-0x1.8f20dd0d0c570a23ea59bc2bf009p-8201q)},
+ {Float128(0x0.00002d0f88d27a496b3e533f5067p-16382q),
+ Float128(0x1.ad9d4abe9f047225a7352bcc52c1p-8201q)},
+ {Float128(0x1.0000000000000000000000000001p+0q), Float128(0x1p+0q)},
+ {Float128(0x1.0000000000000000000000000002p+0q),
+ Float128(-0x1.0000000000000000000000000001p+0q)},
+ {Float128(0x1.0000000000000000000000000003p+0q),
+ Float128(0x1.0000000000000000000000000001p+0q)},
+ {Float128(0x1.0000000000000000000000000005p+0q),
+ Float128(0x1.0000000000000000000000000002p+0q)},
+ {Float128(0x1.0000000000000000000000000006p+0q),
+ Float128(-0x1.0000000000000000000000000003p+0q)},
+ {Float128(0x1.1d4c381cbf3a0aa15b9aee344892p+0q),
+ Float128(0x1.0e408c3fadc5e64b449c63673f4bp+0q)},
+ {Float128(0x1.2af17a4ae6f93d11310c49c11b59p+0q),
+ Float128(-0x1.14a3bdf0ea5231f12d421a5dbe33p+0q)},
+ {Float128(0x1.96f893bf29fb91e0fbe19a46d0c8p+0q),
+ Float128(0x1.42c6bf6202e66f2295807dee44d9p+0q)},
+ {Float128(0x1.97fb3839925b66804c429289cce8p+0q),
+ Float128(-0x1.432d4049ac1c85a241f333d326e9p+0q)},
+ {Float128(0x1.be1d900eaeb1533f0f19cc15c7e6p+0q),
+ Float128(0x1.51f1715154da44f3bf11f3d96c2dp+0q)},
+ {Float128(0x1.c4f5074269525063a26051a0ad27p+0q),
+ Float128(0x1.54864e9b1daa4d9135ff00663366p+0q)},
+ {Float128(0x1.035cb5f298a801dc4be9b1f8cd97p+1q),
+ Float128(-0x1.6c688775bffcb3f507ba11d0abb9p+0q)},
+ {Float128(0x1.274be02380427e709beab4dedeb4p+1q),
+ Float128(-0x1.84d5763281f2318422392e506b1cp+0q)},
+ {Float128(0x1.64e797cfdbaa3f7e2f33279dbc6p+1q),
+ Float128(0x1.ab79b164e255b26eca00ff99cc99p+0q)},
+ {Float128(0x1.693a741358c9dac44a570a7e9f6cp+1q),
+ Float128(0x1.ae0e8eaeab25bb0c40ee0c2693d3p+0q)},
+ {Float128(0x1.8275db3fc4d822596047adcb71b9p+1q),
+ Float128(-0x1.bcd2bfb653e37a5dbe0ccc2cd917p+0q)},
+ {Float128(0x1.83280bb98c4a7b88bd6f535899d9p+1q),
+ Float128(0x1.bd39409dfd1990dd6a7f8211bb27p+0q)},
+ {Float128(0x1.d78d8352b48608b510bfd5c75315p+1q),
+ Float128(-0x1.eb5c420f15adce0ed2bde5a241cep+0q)},
+ {Float128(0x1.e3e4774f564b526edff84ce46668p+1q),
+ Float128(0x1.f1bf73c0523a19b4bb639c98c0b5p+0q)},
+ {Float128(0x1.fffffffffffffffffffffffffffap+1q),
+ Float128(-0x1.fffffffffffffffffffffffffffdp+0q)},
+ {Float128(0x1.fffffffffffffffffffffffffffbp+1q),
+ Float128(0x1.fffffffffffffffffffffffffffdp+0q)},
+ {Float128(0x1.fffffffffffffffffffffffffffdp+1q),
+ Float128(0x1.fffffffffffffffffffffffffffep+0q)},
+ {Float128(0x1.fffffffffffffffffffffffffffep+1q),
+ Float128(-0x1.ffffffffffffffffffffffffffffp+0q)},
+ {Float128(0x1.ffffffffffffffffffffffffffffp+1q),
+ Float128(0x1.ffffffffffffffffffffffffffffp+0q)},
};
auto rnd = [](Float128 x, RoundingMode rm) -> Float128 {
@@ -111,11 +111,16 @@ TEST_F(LlvmLibcSqrtTest, HardToRound) {
// Exact results for subnormal arguments
Float128 EXACT_SUBNORMAL[][2] = {
- {0x0.0000000000000000000000000001p-16382q, 0x1p-8247q},
- {0x0.0000000000000000000000000004p-16382q, 0x1p-8246q},
- {0x0.0000000000001000000000000000p-16382q, 0x1p-8217q},
- {0x0.0000000000010000000000000000p-16382q, 0x1p-8215q},
- {0x0.0000000000100000000000000000p-16382q, 0x1p-8213q},
+ {Float128(0x0.0000000000000000000000000001p-16382q),
+ Float128(0x1p-8247q)},
+ {Float128(0x0.0000000000000000000000000004p-16382q),
+ Float128(0x1p-8246q)},
+ {Float128(0x0.0000000000001000000000000000p-16382q),
+ Float128(0x1p-8217q)},
+ {Float128(0x0.0000000000010000000000000000p-16382q),
+ Float128(0x1p-8215q)},
+ {Float128(0x0.0000000000100000000000000000p-16382q),
+ Float128(0x1p-8213q)},
};
for (auto t : EXACT_SUBNORMAL)
diff --git a/libc/test/src/math/sqrtf128_test.cpp b/libc/test/src/math/sqrtf128_test.cpp
index 4f9ea74d1e204..764403d32d1be 100644
--- a/libc/test/src/math/sqrtf128_test.cpp
+++ b/libc/test/src/math/sqrtf128_test.cpp
@@ -17,24 +17,24 @@ LIST_SQRT_TESTS(Float128, LIBC_NAMESPACE::sqrtf128)
TEST_F(LlvmLibcSqrtTest, SpecialInputs) {
constexpr Float128 INPUTS[] = {
- 0x0.000000dee2f5b6a26c8f07f05442p-16382q,
- 0x0.000000c86d174c5ad8ae54a548e7p-16382q,
- 0x0.000020ab15cfe0b8e488e128f535p-16382q,
- 0x0.0000219e97732a9970f2511989bap-16382q,
- 0x0.000026e477546ae99ef57066f9fdp-16382q,
- 0x0.00002d0f88d27a496b3e533f5067p-16382q,
- 0x1.0000000000000000000000000001p+0q,
- 0x1.0000000000000000000000000003p+0q,
- 0x1.0000000000000000000000000005p+0q,
- 0x1.2af17a4ae6f93d11310c49c11b59p+0q,
- 0x1.c4f5074269525063a26051a0ad27p+0q,
- 0x1.035cb5f298a801dc4be9b1f8cd97p+1q,
- 0x1.274be02380427e709beab4dedeb4p+1q,
- 0x1.64e797cfdbaa3f7e2f33279dbc6p+1q,
- 0x1.d78d8352b48608b510bfd5c75315p+1q,
- 0x1.fffffffffffffffffffffffffffbp+1q,
- 0x1.fffffffffffffffffffffffffffdp+1q,
- 0x1.ffffffffffffffffffffffffffffp+1q,
+ Float128(0x0.000000dee2f5b6a26c8f07f05442p-16382q),
+ Float128(0x0.000000c86d174c5ad8ae54a548e7p-16382q),
+ Float128(0x0.000020ab15cfe0b8e488e128f535p-16382q),
+ Float128(0x0.0000219e97732a9970f2511989bap-16382q),
+ Float128(0x0.000026e477546ae99ef57066f9fdp-16382q),
+ Float128(0x0.00002d0f88d27a496b3e533f5067p-16382q),
+ Float128(0x1.0000000000000000000000000001p+0q),
+ Float128(0x1.0000000000000000000000000003p+0q),
+ Float128(0x1.0000000000000000000000000005p+0q),
+ Float128(0x1.2af17a4ae6f93d11310c49c11b59p+0q),
+ Float128(0x1.c4f5074269525063a26051a0ad27p+0q),
+ Float128(0x1.035cb5f298a801dc4be9b1f8cd97p+1q),
+ Float128(0x1.274be02380427e709beab4dedeb4p+1q),
+ Float128(0x1.64e797cfdbaa3f7e2f33279dbc6p+1q),
+ Float128(0x1.d78d8352b48608b510bfd5c75315p+1q),
+ Float128(0x1.fffffffffffffffffffffffffffbp+1q),
+ Float128(0x1.fffffffffffffffffffffffffffdp+1q),
+ Float128(0x1.ffffffffffffffffffffffffffffp+1q),
};
for (auto input : INPUTS) {
>From b048b26dbb150aa7118d8692410ada422ec284ca Mon Sep 17 00:00:00 2001
From: Sukumarsawant <sawantsukumar at gmail.com>
Date: Mon, 6 Jul 2026 20:24:46 +0530
Subject: [PATCH 3/3] nit
---
libc/test/src/math/smoke/sqrtf128_test.cpp | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/libc/test/src/math/smoke/sqrtf128_test.cpp b/libc/test/src/math/smoke/sqrtf128_test.cpp
index 8eafb9dc7ca6d..54299c78098b1 100644
--- a/libc/test/src/math/smoke/sqrtf128_test.cpp
+++ b/libc/test/src/math/smoke/sqrtf128_test.cpp
@@ -89,7 +89,7 @@ TEST_F(LlvmLibcSqrtTest, HardToRound) {
};
auto rnd = [](Float128 x, RoundingMode rm) -> Float128 {
- bool is_neg = x < 0;
+ bool is_neg = x < Float128(0);
Float128 y = is_neg ? -x : x;
FPBits ybits(y);
More information about the llvm-branch-commits
mailing list