| 1 | //===-- Mimics llvm/Support/MathExtras.h ------------------------*- C++ -*-===// |
| 2 | // Provides useful math functions. |
| 3 | // |
| 4 | // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. |
| 5 | // See https://llvm.org/LICENSE.txt for license information. |
| 6 | // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception |
| 7 | // |
| 8 | //===----------------------------------------------------------------------===// |
| 9 | |
| 10 | #ifndef LLVM_LIBC_SRC___SUPPORT_MATH_EXTRAS_H |
| 11 | #define |
| 12 | |
| 13 | #include "src/__support/CPP/bit.h" // countl_one, countr_zero |
| 14 | #include "src/__support/CPP/limits.h" // CHAR_BIT, numeric_limits |
| 15 | #include "src/__support/CPP/type_traits.h" // is_unsigned_v, is_constant_evaluated |
| 16 | #include "src/__support/macros/attributes.h" // LIBC_INLINE |
| 17 | #include "src/__support/macros/config.h" |
| 18 | |
| 19 | namespace LIBC_NAMESPACE_DECL { |
| 20 | |
| 21 | // Create a bitmask with the count right-most bits set to 1, and all other bits |
| 22 | // set to 0. Only unsigned types are allowed. |
| 23 | template <typename T, size_t count> |
| 24 | LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, T> |
| 25 | mask_trailing_ones() { |
| 26 | constexpr unsigned T_BITS = CHAR_BIT * sizeof(T); |
| 27 | static_assert(count <= T_BITS && "Invalid bit index" ); |
| 28 | return count == 0 ? 0 : (T(-1) >> (T_BITS - count)); |
| 29 | } |
| 30 | |
| 31 | // Create a bitmask with the count left-most bits set to 1, and all other bits |
| 32 | // set to 0. Only unsigned types are allowed. |
| 33 | template <typename T, size_t count> |
| 34 | LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, T> |
| 35 | mask_leading_ones() { |
| 36 | return T(~mask_trailing_ones<T, CHAR_BIT * sizeof(T) - count>()); |
| 37 | } |
| 38 | |
| 39 | // Create a bitmask with the count right-most bits set to 0, and all other bits |
| 40 | // set to 1. Only unsigned types are allowed. |
| 41 | template <typename T, size_t count> |
| 42 | LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, T> |
| 43 | mask_trailing_zeros() { |
| 44 | return mask_leading_ones<T, CHAR_BIT * sizeof(T) - count>(); |
| 45 | } |
| 46 | |
| 47 | // Create a bitmask with the count left-most bits set to 0, and all other bits |
| 48 | // set to 1. Only unsigned types are allowed. |
| 49 | template <typename T, size_t count> |
| 50 | LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, T> |
| 51 | mask_leading_zeros() { |
| 52 | return mask_trailing_ones<T, CHAR_BIT * sizeof(T) - count>(); |
| 53 | } |
| 54 | |
| 55 | // Returns whether 'a + b' overflows, the result is stored in 'res'. |
| 56 | template <typename T> |
| 57 | [[nodiscard]] LIBC_INLINE constexpr bool add_overflow(T a, T b, T &res) { |
| 58 | return __builtin_add_overflow(a, b, &res); |
| 59 | } |
| 60 | |
| 61 | // Returns whether 'a - b' overflows, the result is stored in 'res'. |
| 62 | template <typename T> |
| 63 | [[nodiscard]] LIBC_INLINE constexpr bool sub_overflow(T a, T b, T &res) { |
| 64 | return __builtin_sub_overflow(a, b, &res); |
| 65 | } |
| 66 | |
| 67 | #define RETURN_IF(TYPE, BUILTIN) \ |
| 68 | if constexpr (cpp::is_same_v<T, TYPE>) \ |
| 69 | return BUILTIN(a, b, carry_in, carry_out); |
| 70 | |
| 71 | // Returns the result of 'a + b' taking into account 'carry_in'. |
| 72 | // The carry out is stored in 'carry_out' it not 'nullptr', dropped otherwise. |
| 73 | // We keep the pass by pointer interface for consistency with the intrinsic. |
| 74 | template <typename T> |
| 75 | [[nodiscard]] LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, T> |
| 76 | add_with_carry(T a, T b, T carry_in, T &carry_out) { |
| 77 | if constexpr (!cpp::is_constant_evaluated()) { |
| 78 | #if __has_builtin(__builtin_addcb) |
| 79 | RETURN_IF(unsigned char, __builtin_addcb) |
| 80 | #elif __has_builtin(__builtin_addcs) |
| 81 | RETURN_IF(unsigned short, __builtin_addcs) |
| 82 | #elif __has_builtin(__builtin_addc) |
| 83 | RETURN_IF(unsigned int, __builtin_addc) |
| 84 | #elif __has_builtin(__builtin_addcl) |
| 85 | RETURN_IF(unsigned long, __builtin_addcl) |
| 86 | #elif __has_builtin(__builtin_addcll) |
| 87 | RETURN_IF(unsigned long long, __builtin_addcll) |
| 88 | #endif |
| 89 | } |
| 90 | T sum = {}; |
| 91 | T carry1 = add_overflow(a, b, sum); |
| 92 | T carry2 = add_overflow(sum, carry_in, sum); |
| 93 | carry_out = carry1 | carry2; |
| 94 | return sum; |
| 95 | } |
| 96 | |
| 97 | // Returns the result of 'a - b' taking into account 'carry_in'. |
| 98 | // The carry out is stored in 'carry_out' it not 'nullptr', dropped otherwise. |
| 99 | // We keep the pass by pointer interface for consistency with the intrinsic. |
| 100 | template <typename T> |
| 101 | [[nodiscard]] LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, T> |
| 102 | sub_with_borrow(T a, T b, T carry_in, T &carry_out) { |
| 103 | if constexpr (!cpp::is_constant_evaluated()) { |
| 104 | #if __has_builtin(__builtin_subcb) |
| 105 | RETURN_IF(unsigned char, __builtin_subcb) |
| 106 | #elif __has_builtin(__builtin_subcs) |
| 107 | RETURN_IF(unsigned short, __builtin_subcs) |
| 108 | #elif __has_builtin(__builtin_subc) |
| 109 | RETURN_IF(unsigned int, __builtin_subc) |
| 110 | #elif __has_builtin(__builtin_subcl) |
| 111 | RETURN_IF(unsigned long, __builtin_subcl) |
| 112 | #elif __has_builtin(__builtin_subcll) |
| 113 | RETURN_IF(unsigned long long, __builtin_subcll) |
| 114 | #endif |
| 115 | } |
| 116 | T sub = {}; |
| 117 | T carry1 = sub_overflow(a, b, sub); |
| 118 | T carry2 = sub_overflow(sub, carry_in, sub); |
| 119 | carry_out = carry1 | carry2; |
| 120 | return sub; |
| 121 | } |
| 122 | |
| 123 | #undef RETURN_IF |
| 124 | |
| 125 | template <typename T> |
| 126 | [[nodiscard]] LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, int> |
| 127 | first_leading_zero(T value) { |
| 128 | return value == cpp::numeric_limits<T>::max() ? 0 |
| 129 | : cpp::countl_one(value) + 1; |
| 130 | } |
| 131 | |
| 132 | template <typename T> |
| 133 | [[nodiscard]] LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, int> |
| 134 | first_leading_one(T value) { |
| 135 | return first_leading_zero(static_cast<T>(~value)); |
| 136 | } |
| 137 | |
| 138 | template <typename T> |
| 139 | [[nodiscard]] LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, int> |
| 140 | first_trailing_zero(T value) { |
| 141 | return value == cpp::numeric_limits<T>::max() |
| 142 | ? 0 |
| 143 | : cpp::countr_zero(static_cast<T>(~value)) + 1; |
| 144 | } |
| 145 | |
| 146 | template <typename T> |
| 147 | [[nodiscard]] LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, int> |
| 148 | first_trailing_one(T value) { |
| 149 | return value == 0 ? 0 : cpp::countr_zero(value) + 1; |
| 150 | } |
| 151 | |
| 152 | template <typename T> |
| 153 | [[nodiscard]] LIBC_INLINE constexpr cpp::enable_if_t<cpp::is_unsigned_v<T>, int> |
| 154 | count_zeros(T value) { |
| 155 | return cpp::popcount<T>(static_cast<T>(~value)); |
| 156 | } |
| 157 | |
| 158 | } // namespace LIBC_NAMESPACE_DECL |
| 159 | |
| 160 | #endif // LLVM_LIBC_SRC___SUPPORT_MATH_EXTRAS_H |
| 161 | |