https://github.com/Sukumarsawant updated https://github.com/llvm/llvm-project/pull/207739
>From c654b5d3c49df5f3c7d0f8f117c572f37123ce46 Mon Sep 17 00:00:00 2001 From: Sukumarsawant <[email protected]> Date: Mon, 6 Jul 2026 19:42:15 +0530 Subject: [PATCH 1/4] 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 ([email protected]): // 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 <[email protected]> Date: Mon, 6 Jul 2026 20:14:26 +0530 Subject: [PATCH 2/4] 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 <[email protected]> Date: Mon, 6 Jul 2026 20:24:46 +0530 Subject: [PATCH 3/4] 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); >From eded35c87fbda6c30af54d8592d036d087438c09 Mon Sep 17 00:00:00 2001 From: Sukumarsawant <[email protected]> Date: Mon, 6 Jul 2026 20:32:14 +0530 Subject: [PATCH 4/4] nit --- libc/test/src/math/SqrtTest.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/libc/test/src/math/SqrtTest.h b/libc/test/src/math/SqrtTest.h index d59f5309a085c..6ce552e571845 100644 --- a/libc/test/src/math/SqrtTest.h +++ b/libc/test/src/math/SqrtTest.h @@ -46,7 +46,7 @@ class SqrtTest : public LIBC_NAMESPACE::testing::FEnvSafeTest { for (StorageType i = 0, v = 0; i <= COUNT; ++i, v += STEP) { FPBits x_bits(v); InType x = x_bits.get_val(); - if (x_bits.is_nan() || (x < 0)) + if (x_bits.is_nan() || (x < InType(0))) continue; ASSERT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Sqrt, x, func(x), 0.5); } _______________________________________________ llvm-branch-commits mailing list [email protected] https://lists.llvm.org/cgi-bin/mailman/listinfo/llvm-branch-commits
