native_simd_utility.hpp

// Copyright Take Vos 2022, 2023.
// 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)
 
#pragma once
 
#include "../utility/utility.hpp"
#include "../macros.hpp"
#include <type_traits>
 
#ifdef HI_HAS_SSE
#include <xmmintrin.h>
#endif
#ifdef HI_HAS_SSE2
#include <emmintrin.h>
#endif
#ifdef HI_HAS_SSE3
#include <pmmintrin.h>
#endif
#ifdef HI_HAS_SSSE3
#include <tmmintrin.h>
#endif
#ifdef HI_HAS_SSE4_1
#include <smmintrin.h>
#endif
#ifdef HI_HAS_SSE4_2
#include <nmmintrin.h>
#endif
#ifdef HI_HAS_AVX
#include <immintrin.h>
#endif
#ifdef HI_HAS_AVX512F
#include <immintrin.h>
#endif
 
 
 
namespace hi { inline namespace v1 {
 
#ifdef HI_HAS_SSE
enum class native_rounding_mode : int {
    nearest = _MM_FROUND_TO_NEAREST_INT | _MM_FROUND_NO_EXC,
    negative_infinite = _MM_FROUND_TO_NEG_INF | _MM_FROUND_NO_EXC,
    positive_infinite = _MM_FROUND_TO_POS_INF | _MM_FROUND_NO_EXC,
    truncate = _MM_FROUND_TO_ZERO | _MM_FROUND_NO_EXC,
 
    current = _MM_FROUND_CUR_DIRECTION
};
#endif
 
namespace detail {
 
template<fixed_string SourceElements, size_t NumElements>
[[nodiscard]] constexpr size_t native_swizzle_to_packed_indices() noexcept
{
    constexpr auto max_elements = sizeof(size_t) == 8 ? 16 : 8;
    static_assert(NumElements > 1 and NumElements <= max_elements and std::has_single_bit(NumElements));
    static_assert(SourceElements.size() <= NumElements);
 
    constexpr auto shift = std::bit_width(NumElements - 1);
 
    auto r = 0_uz;
 
    auto i = NumElements;
 
    // Unspecified elements will retain their original position.
    for (; i != SourceElements.size(); --i) {
        r <<= shift;
        r |= i - 1;
    }
 
    for (; i != 0; --i) {
        hilet c = SourceElements[i - 1];
 
        r <<= shift;
        if (c >= '0' and c <= '9') {
            // Numeric elements will retain their original position.
            r |= i - 1;
 
        } else if (c >= 'a' and c <= 'p') {
            r |= char_cast<size_t>(c - 'a');
 
        } else if (c >= 'x' and c <= 'z') {
            r |= char_cast<size_t>(c - 'x');
 
        } else if (c == 'w') {
            r |= 3;
 
        } else {
            hi_no_default();
        }
    }
 
    return r;
}
 
template<fixed_string SourceElements, size_t NumElements, char Value>
[[nodiscard]] constexpr size_t native_swizzle_to_mask() noexcept
{
    static_assert(NumElements > 1 and NumElements <= sizeof(size_t) * CHAR_BIT and std::has_single_bit(NumElements));
    static_assert(SourceElements.size() <= NumElements);
 
    auto r = 0_uz;
 
    auto i = NumElements;
 
    // Unspecified elements will be treated as '0'.
    for (; i != SourceElements.size(); --i) {
        r <<= 1;
        r |= wide_cast<size_t>(Value == '0');
    }
 
    for (; i != 0; --i) {
        hilet c = SourceElements[i - 1];
 
        r <<= 1;
        r |= wide_cast<size_t>(c == Value);
    }
 
    return r;
}
 
} // namespace detail
 
template<typename T, size_t N>
struct native_simd;
 
}} // namespace hi::v1
 
template<typename T, size_t N>
struct std::equal_to<::hi::native_simd<T, N>> {
    [[nodiscard]] constexpr bool operator()(::hi::native_simd<T, N> const& lhs, ::hi::native_simd<T, N> const& rhs) noexcept
    {
        return equal(lhs, rhs);
    }
};
 
namespace testing::internal {
 
// Add equality operator to Google-test internal namespace so that ASSERT_EQ() work.
template<typename T, size_t N>
inline bool operator==(::hi::native_simd<T, N> lhs, ::hi::native_simd<T, N> rhs) noexcept
{
    return std::equal_to<::hi::native_simd<T, N>>{}(lhs, rhs);
}
 
// Add equality operator to Google-test internal namespace so that ASSERT_NE() work.
template<typename T, size_t N>
inline bool operator!=(::hi::native_simd<T, N> lhs, ::hi::native_simd<T, N> rhs) noexcept
{
    return not std::equal_to<::hi::native_simd<T, N>>{}(lhs, rhs);
}
 
} // namespace testing::internal