// (C) Copyright Matt Borland 2021. // Use, modification and distribution are subject to the // Boost Software License, Version 1.0. (See accompanying file // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt) // // Constepxr implementation of abs (see c.math.abs secion 26.8.2 of the ISO standard) #ifndef BOOST_MATH_CCMATH_ABS #define BOOST_MATH_CCMATH_ABS #include #include #include #include #include #include #include namespace boost::math::ccmath { namespace detail { template constexpr T abs_impl(T x) noexcept { if (boost::math::ccmath::isnan(x)) { return std::numeric_limits::quiet_NaN(); } else if (x == static_cast(-0)) { return static_cast(0); } if constexpr (std::is_integral_v) { BOOST_MATH_ASSERT(x != (std::numeric_limits::min)()); } return x >= 0 ? x : -x; } } // Namespace detail template , bool> = true> constexpr T abs(T x) noexcept { if(BOOST_MATH_IS_CONSTANT_EVALUATED(x)) { return detail::abs_impl(x); } else { using std::abs; return abs(x); } } // If abs() is called with an argument of type X for which is_unsigned_v is true and if X // cannot be converted to int by integral promotion (7.3.7), the program is ill-formed. template , bool> = true> constexpr T abs(T x) noexcept { if constexpr (std::is_convertible_v) { return detail::abs_impl(static_cast(x)); } else { static_assert(sizeof(T) == 0, "Taking the absolute value of an unsigned value not covertible to int is UB."); return T(0); // Unreachable, but suppresses warnings } } constexpr long int labs(long int j) noexcept { return boost::math::ccmath::abs(j); } constexpr long long int llabs(long long int j) noexcept { return boost::math::ccmath::abs(j); } } // Namespaces #endif // BOOST_MATH_CCMATH_ABS