/////////////////////////////////////////////////////////////// // Copyright 2013 John Maddock. Distributed under the Boost // Software License, Version 1.0. (See accompanying file // LICENSE_1_0.txt or copy at https://www.boost.org/LICENSE_1_0.txt // // Comparison operators for cpp_int_backend: // #ifndef BOOST_MP_DETAIL_BITSCAN_HPP #define BOOST_MP_DETAIL_BITSCAN_HPP #include #include #if (defined(BOOST_MSVC) || (defined(__clang__) && defined(__c2__)) || (defined(BOOST_INTEL) && defined(_MSC_VER))) && (defined(_M_IX86) || defined(_M_X64)) #include #endif namespace boost { namespace multiprecision { namespace detail { template inline BOOST_MP_CXX14_CONSTEXPR unsigned find_lsb_default(Unsigned mask) { unsigned result = 0; while (!(mask & 1u)) { mask >>= 1; ++result; } return result; } template inline BOOST_MP_CXX14_CONSTEXPR unsigned find_msb_default(Unsigned mask) { unsigned index = 0; while (mask) { ++index; mask >>= 1; } return --index; } template inline BOOST_MP_CXX14_CONSTEXPR unsigned find_lsb(Unsigned mask, const std::integral_constant&) { return find_lsb_default(mask); } template inline BOOST_MP_CXX14_CONSTEXPR unsigned find_msb(Unsigned mask, const std::integral_constant&) { return find_msb_default(mask); } #if (defined(BOOST_MSVC) || (defined(__clang__) && defined(__c2__)) || (defined(BOOST_INTEL) && defined(_MSC_VER))) && (defined(_M_IX86) || defined(_M_X64)) #pragma intrinsic(_BitScanForward, _BitScanReverse) BOOST_FORCEINLINE unsigned find_lsb(unsigned long mask, const std::integral_constant&) { unsigned long result; _BitScanForward(&result, mask); return result; } BOOST_FORCEINLINE unsigned find_msb(unsigned long mask, const std::integral_constant&) { unsigned long result; _BitScanReverse(&result, mask); return result; } #ifdef _M_X64 #pragma intrinsic(_BitScanForward64, _BitScanReverse64) BOOST_FORCEINLINE unsigned find_lsb(unsigned __int64 mask, const std::integral_constant&) { unsigned long result; _BitScanForward64(&result, mask); return result; } template BOOST_FORCEINLINE unsigned find_msb(Unsigned mask, const std::integral_constant&) { unsigned long result; _BitScanReverse64(&result, mask); return result; } #endif template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_lsb(Unsigned mask) { using ui_type = typename boost::multiprecision::detail::make_unsigned::type; using tag_type = typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned long), std::integral_constant, #ifdef _M_X64 typename std::conditional< sizeof(Unsigned) <= sizeof(__int64), std::integral_constant, std::integral_constant >::type #else std::integral_constant #endif >::type; #ifndef BOOST_MP_NO_CONSTEXPR_DETECTION if (BOOST_MP_IS_CONST_EVALUATED(mask)) { return find_lsb_default(mask); } else #endif return find_lsb(static_cast(mask), tag_type()); } template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_msb(Unsigned mask) { using ui_type = typename boost::multiprecision::detail::make_unsigned::type; using tag_type = typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned long), std::integral_constant, #ifdef _M_X64 typename std::conditional< sizeof(Unsigned) <= sizeof(__int64), std::integral_constant, std::integral_constant >::type #else std::integral_constant #endif >::type; #ifndef BOOST_MP_NO_CONSTEXPR_DETECTION if (BOOST_MP_IS_CONST_EVALUATED(mask)) { return find_msb_default(mask); } else #endif return find_msb(static_cast(mask), tag_type()); } #elif defined(BOOST_GCC) || defined(__clang__) || (defined(BOOST_INTEL) && defined(__GNUC__)) BOOST_FORCEINLINE unsigned find_lsb(unsigned mask, std::integral_constant const&) { return __builtin_ctz(mask); } BOOST_FORCEINLINE unsigned find_lsb(unsigned long mask, std::integral_constant const&) { return __builtin_ctzl(mask); } BOOST_FORCEINLINE unsigned find_lsb(boost::ulong_long_type mask, std::integral_constant const&) { return __builtin_ctzll(mask); } BOOST_FORCEINLINE unsigned find_msb(unsigned mask, std::integral_constant const&) { return sizeof(unsigned) * CHAR_BIT - 1 - __builtin_clz(mask); } BOOST_FORCEINLINE unsigned find_msb(unsigned long mask, std::integral_constant const&) { return sizeof(unsigned long) * CHAR_BIT - 1 - __builtin_clzl(mask); } BOOST_FORCEINLINE unsigned find_msb(boost::ulong_long_type mask, std::integral_constant const&) { return sizeof(boost::ulong_long_type) * CHAR_BIT - 1 - __builtin_clzll(mask); } #ifdef BOOST_HAS_INT128 __extension__ using uint128_type = unsigned __int128; BOOST_FORCEINLINE unsigned find_msb(uint128_type mask, std::integral_constant const&) { union { uint128_type v; std::uint64_t sv[2]; } val; val.v = mask; #if BOOST_ENDIAN_LITTLE_BYTE if (val.sv[1]) return find_msb(val.sv[1], std::integral_constant()) + 64; return find_msb(val.sv[0], std::integral_constant()); #else if (val.sv[0]) return find_msb(val.sv[0], std::integral_constant()) + 64; return find_msb(val.sv[1], std::integral_constant()); #endif } BOOST_FORCEINLINE unsigned find_lsb(uint128_type mask, std::integral_constant const&) { union { uint128_type v; std::uint64_t sv[2]; } val; val.v = mask; #if BOOST_ENDIAN_LITTLE_BYTE if (val.sv[0] == 0) return find_lsb(val.sv[1], std::integral_constant()) + 64; return find_lsb(val.sv[0], std::integral_constant()); #else if (val.sv[1] == 0) return find_lsb(val.sv[0], std::integral_constant()) + 64; return find_lsb(val.sv[1], std::integral_constant()); #endif } #endif template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_lsb(Unsigned mask) { using ui_type = typename boost::multiprecision::detail::make_unsigned::type; using tag_type = typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned), std::integral_constant, typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned long), std::integral_constant, typename std::conditional< sizeof(Unsigned) <= sizeof(boost::ulong_long_type), std::integral_constant, std::integral_constant >::type>::type>::type; #ifndef BOOST_MP_NO_CONSTEXPR_DETECTION if (BOOST_MP_IS_CONST_EVALUATED(mask)) { return find_lsb_default(mask); } else #endif return find_lsb(static_cast(mask), tag_type()); } template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_msb(Unsigned mask) { using ui_type = typename boost::multiprecision::detail::make_unsigned::type; using tag_type = typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned), std::integral_constant, typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned long), std::integral_constant, typename std::conditional< sizeof(Unsigned) <= sizeof(boost::ulong_long_type), std::integral_constant, std::integral_constant >::type>::type>::type; #ifndef BOOST_MP_NO_CONSTEXPR_DETECTION if (BOOST_MP_IS_CONST_EVALUATED(mask)) { return find_msb_default(mask); } else #endif return find_msb(static_cast(mask), tag_type()); } #elif defined(BOOST_INTEL) BOOST_FORCEINLINE unsigned find_lsb(unsigned mask, std::integral_constant const&) { return _bit_scan_forward(mask); } BOOST_FORCEINLINE unsigned find_msb(unsigned mask, std::integral_constant const&) { return _bit_scan_reverse(mask); } template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_lsb(Unsigned mask) { using ui_type = typename boost::multiprecision::detail::make_unsigned::type; using tag_type = typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned), std::integral_constant, std::integral_constant >::type; #ifndef BOOST_MP_NO_CONSTEXPR_DETECTION if (BOOST_MP_IS_CONST_EVALUATED(mask)) { return find_lsb_default(mask); } else #endif return find_lsb(static_cast(mask), tag_type()); } template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_msb(Unsigned mask) { using ui_type = typename boost::multiprecision::detail::make_unsigned::type; using tag_type = typename std::conditional< sizeof(Unsigned) <= sizeof(unsigned), std::integral_constant, std::integral_constant >::type; #ifndef BOOST_MP_NO_CONSTEXPR_DETECTION if (BOOST_MP_IS_CONST_EVALUATED(mask)) { return find_msb_default(mask); } else #endif return find_msb(static_cast(mask), tag_type()); } #else template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_lsb(Unsigned mask) { return find_lsb(mask, std::integral_constant()); } template BOOST_FORCEINLINE BOOST_MP_CXX14_CONSTEXPR unsigned find_msb(Unsigned mask) { return find_msb(mask, std::integral_constant()); } #endif }}} // namespace boost::multiprecision::detail #endif