[libcxx-commits] [libcxx] [libc++] simplify the midpoint function (PR #81717)

via libcxx-commits libcxx-commits at lists.llvm.org
Wed Feb 14 01:35:03 PST 2024


https://github.com/rilysh created https://github.com/llvm/llvm-project/pull/81717

Right now we've a nested ternary for the midpoint function, but this can be simplified a bit more, using if..else. This also slightly increases the readability of that function.

>From 6f38fc2c896a14db13ae440d48ec3dc114ae9faf Mon Sep 17 00:00:00 2001
From: rilysh <nightquick at proton.me>
Date: Wed, 14 Feb 2024 14:57:11 +0530
Subject: [PATCH] [libc++] simplify the midpoint function

Right now we've a nested ternary for the midpoint function, but this
can be simplified a bit more, using if..else. This also slightly
increases the readability of that function.
---
 libcxx/include/__numeric/midpoint.h | 18 ++++++++++--------
 1 file changed, 10 insertions(+), 8 deletions(-)

diff --git a/libcxx/include/__numeric/midpoint.h b/libcxx/include/__numeric/midpoint.h
index 5d715c21d8eaca..4232610ce27be4 100644
--- a/libcxx/include/__numeric/midpoint.h
+++ b/libcxx/include/__numeric/midpoint.h
@@ -67,14 +67,16 @@ template <class _Fp>
 _LIBCPP_HIDE_FROM_ABI constexpr enable_if_t<is_floating_point_v<_Fp>, _Fp> midpoint(_Fp __a, _Fp __b) noexcept {
   constexpr _Fp __lo = numeric_limits<_Fp>::min() * 2;
   constexpr _Fp __hi = numeric_limits<_Fp>::max() / 2;
-  return std::__fp_abs(__a) <= __hi && std::__fp_abs(__b) <= __hi
-           ? // typical case: overflow is impossible
-             (__a + __b) / 2
-           :                                             // always correctly rounded
-             std::__fp_abs(__a) < __lo ? __a + __b / 2 : // not safe to halve a
-                 std::__fp_abs(__b) < __lo ? __a / 2 + __b
-                                           : // not safe to halve b
-                 __a / 2 + __b / 2;          // otherwise correctly rounded
+
+  // typical case: overflow is impossible
+  if (std::__fp_abs(__a) <= __hi && std::__fp_abs(__b) <= __hi)
+    return (__a + __b) / 2; // always correctly rounded
+  else if (std::__fp_abs(__a) < __lo)
+    return __a + __b / 2; // not safe to halve a
+  else if (std::__fp_abs(__b) < __lo)
+    return __a / 2 + __b; // not safe to halve b
+  else
+    return __a / 2 + __b / 2; // otherwise correctly rounded
 }
 
 #endif // _LIBCPP_STD_VER >= 20



More information about the libcxx-commits mailing list