simd_sve.h
// Simd SVE specific implementations -*- C++ -*-
// Copyright The GNU Toolchain Authors.
//
// This file is part of the GNU ISO C++ Library. This library is free
// software; you can redistribute it and/or modify it under the
// terms of the GNU General Public License as published by the
// Free Software Foundation; either version 3, or (at your option)
// any later version.
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
// Under Section 7 of GPL version 3, you are granted additional
// permissions described in the GCC Runtime Library Exception, version
// 3.1, as published by the Free Software Foundation.
// You should have received a copy of the GNU General Public License and
// a copy of the GCC Runtime Library Exception along with this program;
// see the files COPYING3 and COPYING.RUNTIME respectively. If not, see
// <http://www.gnu.org/licenses/>.
#ifndef _GLIBCXX_EXPERIMENTAL_SIMD_SVE_H_
#define _GLIBCXX_EXPERIMENTAL_SIMD_SVE_H_
#if __cplusplus >= 201703L
#if !_GLIBCXX_SIMD_HAVE_SVE
#error "simd_sve.h may only be included when SVE on ARM is available"
#endif
_GLIBCXX_SIMD_BEGIN_NAMESPACE
// Helper function mapping to sve supported types
template <typename _Tp>
constexpr auto
__get_sve_value_type()
{
if constexpr (is_integral_v<_Tp>)
{
if constexpr (is_signed_v<_Tp>)
{
if constexpr (sizeof(_Tp) == 1)
return int8_t{};
else if constexpr (sizeof(_Tp) == 2)
return int16_t{};
else if constexpr (sizeof(_Tp) == 4)
return int32_t{};
else if constexpr (sizeof(_Tp) == 8)
return int64_t{};
else
return _Tp{};
}
else
{
if constexpr (sizeof(_Tp) == 1)
return uint8_t{};
else if constexpr (sizeof(_Tp) == 2)
return uint16_t{};
else if constexpr (sizeof(_Tp) == 4)
return uint32_t{};
else if constexpr (sizeof(_Tp) == 8)
return uint64_t{};
else
return _Tp{};
}
}
else
{
if constexpr (is_floating_point_v<_Tp>)
{
if constexpr (sizeof(_Tp) == 4)
return float32_t{};
else if constexpr (sizeof(_Tp) == 8)
return float64_t{};
else
return _Tp{};
}
}
}
template <typename _Tp>
using __get_sve_value_type_t = decltype(__get_sve_value_type<_Tp>());
typedef svbool_t __sve_bool_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
template <typename _Tp, size_t _Np>
struct __sve_vector_type;
template <typename _Tp, size_t _Np>
using __sve_vector_type_t = typename __sve_vector_type<_Tp, _Np>::type;
template <size_t _Np>
struct __sve_vector_type<int8_t, _Np>
{
typedef svint8_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(int8_t __dup)
{ return svdup_s8(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b8(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<uint8_t, _Np>
{
typedef svuint8_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(uint8_t __dup)
{ return svdup_u8(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b8(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<int16_t, _Np>
{
typedef svint16_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(int16_t __dup)
{ return svdup_s16(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b16(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<uint16_t, _Np>
{
typedef svuint16_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(uint16_t __dup)
{ return svdup_u16(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b16(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<int32_t, _Np>
{
typedef svint32_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(int32_t __dup)
{ return svdup_s32(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b32(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<uint32_t, _Np>
{
typedef svuint32_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(uint32_t __dup)
{ return svdup_u32(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b32(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<int64_t, _Np>
{
typedef svint64_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(int64_t __dup)
{ return svdup_s64(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b64(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<uint64_t, _Np>
{
typedef svuint64_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(uint64_t __dup)
{ return svdup_u64(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b64(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<float, _Np>
{
typedef svfloat32_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(float __dup)
{ return svdup_f32(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b32(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <size_t _Np>
struct __sve_vector_type<double, _Np>
{
typedef svfloat64_t __sve_vlst_type __attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static __sve_vlst_type
__sve_broadcast(double __dup)
{ return svdup_f64(__dup); }
inline static __sve_bool_type
__sve_active_mask()
{ return svwhilelt_b64(size_t(0), _Np); };
using type = __sve_vlst_type;
};
template <typename _Tp, size_t _Np>
struct __sve_vector_type
: __sve_vector_type<__get_sve_value_type_t<_Tp>, _Np>
{};
template <size_t _Size>
struct __sve_mask_type
{
static_assert((_Size & (_Size - 1)) != 0, "This trait may only be used for non-power-of-2 "
"sizes. Power-of-2 sizes must be specialized.");
using type = typename __sve_mask_type<std::__bit_ceil(_Size)>::type;
};
template <size_t _Size>
using __sve_mask_type_t = typename __sve_mask_type<_Size>::type;
template <>
struct __sve_mask_type<1>
{
using type = __sve_bool_type;
using __sve_mask_uint_type = uint8_t;
typedef svuint8_t __sve_mask_vector_type
__attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static auto
__sve_mask_active_count(type __active_mask, type __pred)
{ return svcntp_b8(__active_mask, __pred); }
inline static type
__sve_mask_first_true()
{ return svptrue_pat_b8(SV_VL1); }
inline static type
__sve_mask_next_true(type __active_mask, type __pred)
{ return svpnext_b8(__active_mask, __pred); }
inline static bool
__sve_mask_get(type __active_mask, size_t __i)
{ return __sve_mask_vector_type(svdup_u8_z(__active_mask, 1))[__i] != 0;}
inline static const __sve_mask_vector_type __index0123 = svindex_u8(0, 1);
};
template <>
struct __sve_mask_type<2>
{
using type = __sve_bool_type;
using __sve_mask_uint_type = uint16_t;
typedef svuint16_t __sve_mask_vector_type
__attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static auto
__sve_mask_active_count(type __active_mask, type __pred)
{ return svcntp_b16(__active_mask, __pred); }
inline static type
__sve_mask_first_true()
{ return svptrue_pat_b16(SV_VL1); }
inline static type
__sve_mask_next_true(type __active_mask, type __pred)
{ return svpnext_b16(__active_mask, __pred); }
inline static bool
__sve_mask_get(type __active_mask, size_t __i)
{ return __sve_mask_vector_type(svdup_u16_z(__active_mask, 1))[__i] != 0;}
inline static const __sve_mask_vector_type __index0123 = svindex_u16(0, 1);
};
template <>
struct __sve_mask_type<4>
{
using type = __sve_bool_type;
using __sve_mask_uint_type = uint32_t;
typedef svuint32_t __sve_mask_vector_type
__attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static auto
__sve_mask_active_count(type __active_mask, type __pred)
{ return svcntp_b32(__active_mask, __pred); }
inline static type
__sve_mask_first_true()
{ return svptrue_pat_b32(SV_VL1); }
inline static type
__sve_mask_next_true(type __active_mask, type __pred)
{ return svpnext_b32(__active_mask, __pred); }
inline static bool
__sve_mask_get(type __active_mask, size_t __i)
{ return __sve_mask_vector_type(svdup_u32_z(__active_mask, 1))[__i] != 0;}
inline static const __sve_mask_vector_type __index0123 = svindex_u32(0, 1);
};
template <>
struct __sve_mask_type<8>
{
using type = __sve_bool_type;
using __sve_mask_uint_type = uint64_t;
typedef svuint64_t __sve_mask_vector_type
__attribute__((arm_sve_vector_bits(__ARM_FEATURE_SVE_BITS)));
inline static auto
__sve_mask_active_count(type __active_mask, type __pred)
{ return svcntp_b64(__active_mask, __pred); }
inline static type
__sve_mask_first_true()
{ return svptrue_pat_b64(SV_VL1); }
inline static type
__sve_mask_next_true(type __active_mask, type __pred)
{ return svpnext_b64(__active_mask, __pred); }
inline static bool
__sve_mask_get(type __active_mask, size_t __i)
{ return __sve_mask_vector_type(svdup_u64_z(__active_mask, 1))[__i] != 0;}
inline static const __sve_mask_vector_type __index0123 = svindex_u64(0, 1);
};
template <typename _To, typename _From>
_GLIBCXX_SIMD_INTRINSIC constexpr auto
__sve_reinterpret_cast(_From __v)
{
if constexpr (std::is_same_v<_To, int32_t>)
return svreinterpret_s32(__v);
else if constexpr (std::is_same_v<_To, int64_t>)
return svreinterpret_s64(__v);
else if constexpr (std::is_same_v<_To, float32_t>)
return svreinterpret_f32(__v);
else if constexpr (std::is_same_v<_To, float64_t>)
return svreinterpret_f64(__v);
else
__assert_unreachable<_To>(); // add more cases if needed.
}
template <typename _Tp, size_t _Width>
struct _SveSimdWrapper
{
static_assert(__is_vectorizable_v<_Tp>);
static_assert(_Width >= 2); // 1 doesn't make sense, use _Tp directly then
using _BuiltinType = __sve_vector_type_t<_Tp, _Width>;
using value_type = _Tp;
static inline constexpr size_t _S_full_size = sizeof(_BuiltinType) / sizeof(value_type);
static inline constexpr int _S_size = _Width;
static inline constexpr bool _S_is_partial = _S_full_size != _S_size;
_BuiltinType _M_data;
_GLIBCXX_SIMD_INTRINSIC constexpr _SveSimdWrapper<_Tp, _S_full_size>
__as_full_vector() const
{ return _M_data; }
_GLIBCXX_SIMD_INTRINSIC constexpr
_SveSimdWrapper(initializer_list<_Tp> __init)
: _M_data(__generate_from_n_evaluations<_Width, _BuiltinType>(
[&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
return __init.begin()[__i.value];
}))
{}
_GLIBCXX_SIMD_INTRINSIC constexpr
_SveSimdWrapper() = default;
_GLIBCXX_SIMD_INTRINSIC constexpr
_SveSimdWrapper(const _SveSimdWrapper&) = default;
_GLIBCXX_SIMD_INTRINSIC constexpr
_SveSimdWrapper(_SveSimdWrapper&&) = default;
_GLIBCXX_SIMD_INTRINSIC constexpr _SveSimdWrapper&
operator=(const _SveSimdWrapper&) = default;
_GLIBCXX_SIMD_INTRINSIC constexpr _SveSimdWrapper&
operator=(_SveSimdWrapper&&) = default;
_GLIBCXX_SIMD_INTRINSIC constexpr
_SveSimdWrapper(__sve_vector_type_t<_Tp, _Width> __x)
: _M_data(__x)
{}
template <typename... _As, typename = enable_if_t<((is_same_v<simd_abi::scalar, _As> && ...)
&& sizeof...(_As) <= _Width)>>
_GLIBCXX_SIMD_INTRINSIC constexpr
operator _SimdTuple<_Tp, _As...>() const
{
return __generate_from_n_evaluations<sizeof...(_As), _SimdTuple<_Tp, _As...>>(
[&](auto __i) constexpr _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
return _M_data[int(__i)];
});
}
_GLIBCXX_SIMD_INTRINSIC constexpr
operator const _BuiltinType&() const
{ return _M_data; }
_GLIBCXX_SIMD_INTRINSIC constexpr
operator _BuiltinType&()
{ return _M_data; }
_GLIBCXX_SIMD_INTRINSIC constexpr _Tp
operator[](size_t __i) const
{ return _M_data[__i]; }
template <size_t __i>
_GLIBCXX_SIMD_INTRINSIC constexpr _Tp
operator[](_SizeConstant<__i>) const
{ return _M_data[__i]; }
_GLIBCXX_SIMD_INTRINSIC constexpr void
_M_set(size_t __i, _Tp __x)
{
_M_data[__i] = __x;
}
_GLIBCXX_SIMD_INTRINSIC constexpr bool
_M_is_constprop() const
{ return false; }
_GLIBCXX_SIMD_INTRINSIC constexpr bool
_M_is_constprop_none_of() const
{ return false; }
_GLIBCXX_SIMD_INTRINSIC constexpr bool
_M_is_constprop_all_of() const
{ return false; }
};
template <size_t _Bits, size_t _Width>
struct _SveMaskWrapper
{
using _BuiltinSveMaskType = __sve_mask_type<_Bits>;
using _BuiltinSveVectorType = __sve_vector_type<__int_with_sizeof_t<_Bits>, _Width>;
using _BuiltinType = typename _BuiltinSveMaskType::type;
using value_type = bool;
static constexpr size_t _S_full_size = sizeof(_BuiltinType);
_GLIBCXX_SIMD_INTRINSIC constexpr _SveMaskWrapper<_Bits, _S_full_size>
__as_full_vector() const
{ return _M_data; }
_GLIBCXX_SIMD_INTRINSIC constexpr
_SveMaskWrapper() = default;
_GLIBCXX_SIMD_INTRINSIC constexpr
_SveMaskWrapper(_BuiltinType __k)
: _M_data(__k)
{};
_GLIBCXX_SIMD_INTRINSIC
operator const _BuiltinType&() const
{ return _M_data; }
_GLIBCXX_SIMD_INTRINSIC
operator _BuiltinType&()
{ return _M_data; }
_GLIBCXX_SIMD_INTRINSIC _BuiltinType
__intrin() const
{ return _M_data; }
_GLIBCXX_SIMD_INTRINSIC constexpr value_type
operator[](size_t __i) const
{
return _BuiltinSveMaskType::__sve_mask_get(_M_data, __i);
}
template <size_t __i>
_GLIBCXX_SIMD_INTRINSIC constexpr value_type
operator[](_SizeConstant<__i>) const
{
return _BuiltinSveMaskType::__sve_mask_get(_M_data, __i);
}
_GLIBCXX_SIMD_INTRINSIC constexpr void
_M_set(size_t __i, value_type __x)
{
_BuiltinType __index
= svcmpeq(_BuiltinSveVectorType::__sve_active_mask(), _BuiltinSveMaskType::__index0123,
typename _BuiltinSveMaskType::__sve_mask_uint_type(__i));
if (__x)
_M_data = svorr_z(_BuiltinSveVectorType::__sve_active_mask(), _M_data, __index);
else
_M_data = svbic_z(_BuiltinSveVectorType::__sve_active_mask(), _M_data, __index);
}
_GLIBCXX_SIMD_INTRINSIC constexpr bool
_M_is_constprop() const
{ return false; }
_GLIBCXX_SIMD_INTRINSIC constexpr bool
_M_is_constprop_none_of() const
{ return false; }
_GLIBCXX_SIMD_INTRINSIC constexpr bool
_M_is_constprop_all_of() const
{ return false; }
_BuiltinType _M_data;
};
struct _CommonImplSve;
template <typename _Abi, typename = __detail::__odr_helper>
struct _SimdImplSve;
template <typename _Abi, typename = __detail::__odr_helper>
struct _MaskImplSve;
template <int _UsedBytes, int>
struct simd_abi::_SveAbi
{
template <typename _Tp>
static constexpr size_t _S_size = _UsedBytes / sizeof(_Tp);
struct _IsValidAbiTag
: __bool_constant<(_UsedBytes > 1)>
{};
template <typename _Tp>
struct _IsValidSizeFor
: __bool_constant<(_UsedBytes / sizeof(_Tp) > 1 && _UsedBytes % sizeof(_Tp) == 0
&& _UsedBytes <= __sve_vectorized_size_bytes)>
{};
template <typename _Tp>
struct _IsValid
: conjunction<_IsValidAbiTag, __bool_constant<__have_sve>,
__bool_constant<(__vectorized_sizeof<_Tp>() > sizeof(_Tp))>,
_IsValidSizeFor<_Tp>>
{};
template <typename _Tp>
static constexpr bool _S_is_valid_v = _IsValid<_Tp>::value;
using _CommonImpl = _CommonImplSve;
using _SimdImpl = _SimdImplSve<_SveAbi<_UsedBytes>>;
using _MaskImpl = _MaskImplSve<_SveAbi<_UsedBytes>>;
template <typename _Tp>
using _MaskMember = _SveMaskWrapper<sizeof(_Tp), _S_size<_Tp>>;
template <typename _Tp, bool = _S_is_valid_v<_Tp>>
struct __traits : _InvalidTraits
{};
template <typename _Tp>
struct __traits<_Tp, true>
{
using _IsValid = true_type;
using _SimdImpl = _SimdImplSve<_SveAbi<_UsedBytes>>;
using _MaskImpl = _MaskImplSve<_SveAbi<_UsedBytes>>;
using _SimdMember = _SveSimdWrapper<_Tp, _S_size<_Tp>>; // sve vector type
using _MaskMember = _SveMaskWrapper<sizeof(_Tp), _S_size<_Tp>>; // sve mask type
static constexpr size_t _S_simd_align = alignof(_SimdMember);
static constexpr size_t _S_mask_align = alignof(_MaskMember);
static constexpr size_t _S_full_size = _SimdMember::_S_full_size;
static constexpr bool _S_is_partial = _SimdMember::_S_is_partial;
struct _SimdBase
{
_GLIBCXX_SIMD_ALWAYS_INLINE explicit
operator __sve_vector_type_t<_Tp, _S_size<_Tp>>() const
{ return __data(*static_cast<const simd<_Tp, _SveAbi<_UsedBytes>>*>(this)); }
};
class _SimdCastType
{
using _Ap = __sve_vector_type_t<_Tp, _S_size<_Tp>>;
_SimdMember _M_data;
public:
_GLIBCXX_SIMD_ALWAYS_INLINE constexpr
_SimdCastType(_Ap __a)
: _M_data(__a)
{}
_GLIBCXX_SIMD_ALWAYS_INLINE constexpr
operator _SimdMember() const
{ return _M_data; }
};
struct _MaskBase
{
_GLIBCXX_SIMD_ALWAYS_INLINE explicit
operator __sve_mask_type_t<sizeof(_Tp)>() const
{
return __data(*static_cast<const simd_mask<_Tp, _SveAbi<_UsedBytes>>*>(this));
}
};
class _MaskCastType
{
using _Ap = __sve_mask_type_t<sizeof(_Tp)>;
_Ap _M_data;
public:
_GLIBCXX_SIMD_ALWAYS_INLINE constexpr
_MaskCastType(_Ap __a)
: _M_data(__a)
{}
_GLIBCXX_SIMD_ALWAYS_INLINE constexpr
operator _MaskMember() const
{ return _M_data; }
};
};
template <typename _Tp>
static constexpr size_t _S_full_size = __traits<_Tp>::_S_full_size;
template <typename _Tp>
static constexpr bool _S_is_partial = __traits<_Tp>::_S_is_partial;
};
template <typename _Tp, size_t _Np>
using __sve_mask = __sve_mask_type<sizeof(_Tp)>;
struct _CommonImplSve
{
// _S_converts_via_decomposition
// This lists all cases where a __vector_convert needs to fall back to
// conversion of individual scalars (i.e. decompose the input vector into
// scalars, convert, compose output vector). In those cases, _S_masked_load &
// _S_masked_store prefer to use the _S_bit_iteration implementation.
template <typename _From, typename _To, size_t _ToSize>
static inline constexpr bool __converts_via_decomposition_v = sizeof(_From) != sizeof(_To);
template <typename _Tp, typename _Up, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_load(const _Up* __p, _SveMaskWrapper<sizeof(_Tp), _Np> __k)
{
using _STp = __get_sve_value_type_t<_Tp>;
using _SUp = __get_sve_value_type_t<_Up>;
using _V = __sve_vector_type_t<_Tp, _Np>;
const _SUp* __up = reinterpret_cast<const _SUp*>(__p);
if constexpr (std::is_same_v<_Tp, _Up>)
return _V(svld1(__k._M_data, __up));
if constexpr (std::is_integral_v<_Tp> && std::is_integral_v<_Up>
&& (sizeof(_Tp) > sizeof(_Up)))
{
if constexpr (std::is_same_v<_SUp, int8_t>)
{
if constexpr (std::is_same_v<_STp, int16_t>)
return _V(svld1sb_s16(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint16_t>)
return _V(svld1sb_u16(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, int32_t>)
return _V(svld1sb_s32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint32_t>)
return _V(svld1sb_u32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, int64_t>)
return _V(svld1sb_s64(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint64_t>)
return _V(svld1sb_u64(__k._M_data, __up));
}
if constexpr (std::is_same_v<_SUp, uint8_t>)
{
if constexpr (std::is_same_v<_STp, int16_t>)
return _V(svld1ub_s16(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint16_t>)
return _V(svld1ub_u16(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, int32_t>)
return _V(svld1ub_s32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint32_t>)
return _V(svld1ub_u32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, int64_t>)
return _V(svld1ub_s64(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint64_t>)
return _V(svld1ub_u64(__k._M_data, __up));
}
if constexpr (std::is_same_v<_SUp, int16_t>)
{
if constexpr (std::is_same_v<_STp, int32_t>)
return _V(svld1sh_s32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint32_t>)
return _V(svld1sh_u32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, int64_t>)
return _V(svld1sh_s64(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint64_t>)
return _V(svld1sh_u64(__k._M_data, __up));
}
if constexpr (std::is_same_v<_SUp, uint16_t>)
{
if constexpr (std::is_same_v<_STp, int32_t>)
return _V(svld1uh_s32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint32_t>)
return _V(svld1uh_u32(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, int64_t>)
return _V(svld1uh_s64(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint64_t>)
return _V(svld1uh_u64(__k._M_data, __up));
}
if constexpr (std::is_same_v<_SUp, int32_t>)
{
if constexpr (std::is_same_v<_STp, int64_t>)
return _V(svld1sw_s64(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint64_t>)
return _V(svld1sw_u64(__k._M_data, __up));
}
if constexpr (std::is_same_v<_SUp, uint32_t>)
{
if constexpr (std::is_same_v<_STp, int64_t>)
return _V(svld1uw_s64(__k._M_data, __up));
if constexpr (std::is_same_v<_STp, uint64_t>)
return _V(svld1uw_u64(__k._M_data, __up));
}
}
return __generate_from_n_evaluations<_Np, __sve_vector_type_t<_Tp, _Np>>(
[&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
return __k[__i] ? static_cast<_Tp>(__p[__i]) : _Tp{};
});
}
template <typename _Tp, typename _Up, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_store(_Up* __p, _SveSimdWrapper<_Tp, _Np> __x, _SveMaskWrapper<sizeof(_Tp), _Np> __k)
{
using _SUp = __get_sve_value_type_t<_Up>;
using _STp = __get_sve_value_type_t<_Tp>;
_SUp* __up = reinterpret_cast<_SUp*>(__p);
if constexpr (std::is_same_v<_Tp, _Up>)
return svst1(__k._M_data, __up, __x);
if constexpr (std::is_integral_v<_Tp> && std::is_integral_v<_Up>
&& (sizeof(_Tp) > sizeof(_Up)))
{
if constexpr (std::is_same_v<_SUp, int8_t> && std::is_signed_v<_STp>)
return svst1b(__k._M_data, __up, __x);
if constexpr (std::is_same_v<_SUp, uint8_t> && std::is_unsigned_v<_STp>)
return svst1b(__k._M_data, __up, __x);
if constexpr (std::is_same_v<_SUp, int16_t> && std::is_signed_v<_STp>)
return svst1h(__k._M_data, __up, __x);
if constexpr (std::is_same_v<_SUp, uint16_t> && std::is_unsigned_v<_STp>)
return svst1h(__k._M_data, __up, __x);
if constexpr (std::is_same_v<_SUp, int32_t> && std::is_signed_v<_STp>)
return svst1w(__k._M_data, __up, __x);
if constexpr (std::is_same_v<_SUp, uint32_t> && std::is_unsigned_v<_STp>)
return svst1w(__k._M_data, __up, __x);
}
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
if (__k[__i])
__p[__i] = static_cast<_Up>(__x[__i]);
});
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_blend(_SveMaskWrapper<sizeof(_Tp), _Np> __k, _SveSimdWrapper<_Tp, _Np> __at0,
_SveSimdWrapper<_Tp, _Np> __at1)
{ return svsel(__k._M_data, __at1._M_data, __at0._M_data); }
template <size_t _Np, bool _Sanitized>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_store_bool_array(_BitMask<_Np, _Sanitized> __x, bool* __mem)
{
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
__mem[__i] = __x[__i];
});
}
};
template <typename _Abi, typename>
struct _SimdImplSve
{
template <typename _Tp>
using _MaskMember = typename _Abi::template _MaskMember<_Tp>;
template <typename _Tp>
using _SimdMember = typename _Abi::template __traits<_Tp>::_SimdMember;
using _CommonImpl = typename _Abi::_CommonImpl;
using _SuperImpl = typename _Abi::_SimdImpl;
using _MaskImpl = typename _Abi::_MaskImpl;
template <typename _Tp>
static constexpr size_t _S_full_size = _Abi::template _S_full_size<_Tp>;
template <typename _Tp>
static constexpr size_t _S_size = _Abi::template _S_size<_Tp>;
template <typename _Tp>
using _TypeTag = _Tp*;
using abi_type = _Abi;
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr auto
_S_broadcast(_Tp __x) noexcept
{
return __sve_vector_type<_Tp, __sve_vectorized_size_bytes / sizeof(_Tp)>
::__sve_broadcast(__x);
}
template <typename _Fp, typename _Tp>
inline static constexpr _SimdMember<_Tp>
_S_generator(_Fp&& __gen, _TypeTag<_Tp>)
{
constexpr size_t _Np = _S_size<_Tp>;
_SveSimdWrapper<_Tp, _Np> __ret;
__execute_n_times<_S_size<_Tp>>(
[&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA { __ret._M_set(__i, __gen(__i)); });
return __ret;
}
template <typename _Tp, typename _Up>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SimdMember<_Tp>
_S_load(const _Up* __mem, _TypeTag<_Tp>) noexcept
{
constexpr size_t _Np = _S_size<_Tp>;
_SimdMember<_Tp> __ret = _CommonImpl::template _S_load<_Tp, _Up, _Np>(
__mem, _SveMaskWrapper<sizeof(_Tp), _Np>{
__sve_vector_type<_Tp, _Np>::__sve_active_mask()});
return __ret;
}
template <typename _Tp, size_t _Np, typename _Up>
static constexpr inline _SveSimdWrapper<_Tp, _Np>
_S_masked_load(_SveSimdWrapper<_Tp, _Np> __merge, _MaskMember<_Tp> __k, const _Up* __mem)
noexcept
{
__sve_vector_type_t<_Tp, _Np> __v
= _CommonImpl::template _S_load<_Tp, _Up, _Np>(__mem, __k);
__sve_vector_type_t<_Tp, _Np> __ret = svsel(__k._M_data, __v, __merge._M_data);
return __ret;
}
template <typename _Tp, typename _Up>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_store(_SimdMember<_Tp> __v, _Up* __mem, _TypeTag<_Tp>) noexcept
{
constexpr size_t _Np = _S_size<_Tp>;
_CommonImpl::template _S_store<_Tp, _Up, _Np>(
__mem, __v, __sve_vector_type<_Tp, _Np>::__sve_active_mask());
}
template <typename _Tp, typename _Up, size_t _Np>
static constexpr inline void
_S_masked_store(const _SveSimdWrapper<_Tp, _Np> __v, _Up* __mem,
const _SveMaskWrapper<sizeof(_Tp), _Np> __k) noexcept
{ _CommonImpl::template _S_store<_Tp, _Up, _Np>(__mem, __v, __k); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_negate(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{
return svcmpeq(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data,
__sve_vector_type<_Tp, _Np>::__sve_broadcast(_Tp{}));
}
template <typename _Tp, typename _BinaryOperation>
_GLIBCXX_SIMD_INTRINSIC static constexpr _Tp
_S_reduce(simd<_Tp, _Abi> __x, _BinaryOperation&& __binary_op)
{
auto __x_data = __x._M_data;
constexpr size_t _Np = simd_size_v<_Tp, _Abi>;
using __sve_vec_t = __sve_vector_type_t<_Tp, _Np>;
std::size_t __i = __x.size();
for (; (__i % 2) != 1; __i /= 2)
{
__x_data = __binary_op(simd<_Tp, _Abi>(
__private_init, _SveSimdWrapper<_Tp, _Np>(
__sve_vec_t(svuzp1(__x_data, __x_data)))),
simd<_Tp, _Abi>(
__private_init, _SveSimdWrapper<_Tp, _Np>(
__sve_vec_t(svuzp2(__x_data, __x_data))))
)._M_data;
}
_Tp __res = __x_data[0];
for (size_t __ri = 1; __ri != __i; __ri++)
__res = __binary_op(__x_data[__ri], __res);
return __res;
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _Tp
_S_reduce(simd<_Tp, _Abi> __x, plus<>)
{
return svaddv(__sve_vector_type<_Tp, _S_size<_Tp>>::__sve_active_mask(), __x._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _Tp
_S_reduce(simd<_Tp, _Abi> __x, bit_and<>)
{
return svandv(__sve_vector_type<_Tp, _S_size<_Tp>>::__sve_active_mask(), __x._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _Tp
_S_reduce(simd<_Tp, _Abi> __x, bit_or<>)
{
return svorv(__sve_vector_type<_Tp, _S_size<_Tp>>::__sve_active_mask(), __x._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _Tp
_S_reduce(simd<_Tp, _Abi> __x, bit_xor<>)
{
return sveorv(__sve_vector_type<_Tp, _S_size<_Tp>>::__sve_active_mask(), __x._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _Tp
_S_reduce(simd<_Tp, _Abi> __x, __detail::_Maximum())
{
return svmaxv(__sve_vector_type<_Tp, _S_size<_Tp>>::__sve_active_mask(), __x._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _Tp
_S_reduce(simd<_Tp, _Abi> __x, __detail::_Minimum())
{
return svminv(__sve_vector_type<_Tp, _S_size<_Tp>>::__sve_active_mask(), __x._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_NORMAL_MATH _GLIBCXX_SIMD_INTRINSIC static constexpr
__sve_vector_type_t<_Tp, _Np>
_S_min(_SveSimdWrapper<_Tp, _Np> __a, _SveSimdWrapper<_Tp, _Np> __b)
{
return svmin_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __a._M_data, __b._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_NORMAL_MATH _GLIBCXX_SIMD_INTRINSIC static constexpr
__sve_vector_type_t<_Tp, _Np>
_S_max(_SveSimdWrapper<_Tp, _Np> __a, _SveSimdWrapper<_Tp, _Np> __b)
{
return svmax_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __a._M_data, __b._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_NORMAL_MATH _GLIBCXX_SIMD_INTRINSIC static constexpr
pair<_SveSimdWrapper<_Tp, _Np>, _SveSimdWrapper<_Tp, _Np>>
_S_minmax(_SveSimdWrapper<_Tp, _Np> __a, _SveSimdWrapper<_Tp, _Np> __b)
{
return {
svmin_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __a._M_data, __b._M_data),
svmax_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __a._M_data, __b._M_data)
};
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_complement(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{
if constexpr (is_floating_point_v<_Tp>)
{
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
return __sve_reinterpret_cast<_Tp>(
svnot_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__sve_reinterpret_cast<_Ip>(__x)));
}
else
return svnot_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveSimdWrapper<_Tp, _Np>
_S_unary_minus(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{
return svmul_x(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data,
static_cast<_Tp>(-1));
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_plus(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{ return __x._M_data + __y._M_data; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_minus(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{ return __x._M_data - __y._M_data; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_multiplies(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{ return __x._M_data * __y._M_data; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_divides(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
__sve_vector_type_t<_Tp, _Np> __y_padded
= svsel(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__y._M_data, __sve_vector_type<_Tp, _Np>::__sve_broadcast(1));
return __x._M_data / __y_padded;
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_modulus(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
__sve_vector_type_t<_Tp, _Np> __y_padded
= svsel(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__y._M_data, __sve_vector_type<_Tp, _Np>::__sve_broadcast(1));
return __x._M_data % __y_padded;
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_bit_and(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
if constexpr (is_floating_point_v<_Tp>)
{
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
return __sve_reinterpret_cast<_Tp>(
svand_x(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__sve_reinterpret_cast<_Ip>(__x), __sve_reinterpret_cast<_Ip>(__y)));
}
else
return svand_x(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_bit_or(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
if constexpr (is_floating_point_v<_Tp>)
{
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
return __sve_reinterpret_cast<_Tp>(
svorr_x(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__sve_reinterpret_cast<_Ip>(__x), __sve_reinterpret_cast<_Ip>(__y)));
}
else
return svorr_x(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_bit_xor(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
if constexpr (is_floating_point_v<_Tp>)
{
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
return __sve_reinterpret_cast<_Tp>(
sveor_x(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__sve_reinterpret_cast<_Ip>(__x), __sve_reinterpret_cast<_Ip>(__y)));
}
else
return sveor_x(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static __sve_vector_type_t<_Tp, _Np>
_S_bit_shift_left(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{ return __x._M_data << __y._M_data; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static __sve_vector_type_t<_Tp, _Np>
_S_bit_shift_right(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{ return __x._M_data >> __y._M_data; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_bit_shift_left(_SveSimdWrapper<_Tp, _Np> __x, int __y)
{ return __x._M_data << __y; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr __sve_vector_type_t<_Tp, _Np>
_S_bit_shift_right(_SveSimdWrapper<_Tp, _Np> __x, int __y)
{ return __x._M_data >> __y; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_increment(_SveSimdWrapper<_Tp, _Np>& __x)
{ __x = __x._M_data + 1; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_decrement(_SveSimdWrapper<_Tp, _Np>& __x)
{ __x = __x._M_data - 1; }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_equal_to(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
return svcmpeq(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_not_equal_to(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
return svcmpne(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_less(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
return svcmplt(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_less_equal(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
return svcmple(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data);
}
// simd.math
#define _GLIBCXX_SIMD_MATH_FALLBACK(__name) \
template <typename _Tp, size_t _Np, typename... _More> \
static _SveSimdWrapper<_Tp, _Np> _S_##__name(const _SveSimdWrapper<_Tp, _Np>& __x, \
const _More&... __more) \
{ \
_SveSimdWrapper<_Tp, _Np> __r; \
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA { \
__r._M_set(__i, __name(__x[__i], __more[__i]...)); \
}); \
return __r; \
}
#define _GLIBCXX_SIMD_MATH_FALLBACK_FIXEDRET(_RetTp, __name) \
template <typename _Tp, typename... _More> \
static auto _S_##__name(const _Tp& __x, const _More&... __more) \
{ \
return __fixed_size_storage_t<_RetTp, _Tp::_S_size>::_S_generate( \
[&](auto __meta) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA { \
return __meta._S_generator( \
[&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA { \
return __name(__x[__meta._S_offset + __i], \
__more[__meta._S_offset + __i]...); \
}, static_cast<_RetTp*>(nullptr)); \
}); \
}
_GLIBCXX_SIMD_MATH_FALLBACK(acos)
_GLIBCXX_SIMD_MATH_FALLBACK(asin)
_GLIBCXX_SIMD_MATH_FALLBACK(atan)
_GLIBCXX_SIMD_MATH_FALLBACK(atan2)
_GLIBCXX_SIMD_MATH_FALLBACK(cos)
_GLIBCXX_SIMD_MATH_FALLBACK(sin)
_GLIBCXX_SIMD_MATH_FALLBACK(tan)
_GLIBCXX_SIMD_MATH_FALLBACK(acosh)
_GLIBCXX_SIMD_MATH_FALLBACK(asinh)
_GLIBCXX_SIMD_MATH_FALLBACK(atanh)
_GLIBCXX_SIMD_MATH_FALLBACK(cosh)
_GLIBCXX_SIMD_MATH_FALLBACK(sinh)
_GLIBCXX_SIMD_MATH_FALLBACK(tanh)
_GLIBCXX_SIMD_MATH_FALLBACK(exp)
_GLIBCXX_SIMD_MATH_FALLBACK(exp2)
_GLIBCXX_SIMD_MATH_FALLBACK(expm1)
_GLIBCXX_SIMD_MATH_FALLBACK_FIXEDRET(int, ilogb)
_GLIBCXX_SIMD_MATH_FALLBACK(log)
_GLIBCXX_SIMD_MATH_FALLBACK(log10)
_GLIBCXX_SIMD_MATH_FALLBACK(log1p)
_GLIBCXX_SIMD_MATH_FALLBACK(log2)
_GLIBCXX_SIMD_MATH_FALLBACK(logb)
// modf implemented in simd_math.h
_GLIBCXX_SIMD_MATH_FALLBACK(scalbn)
_GLIBCXX_SIMD_MATH_FALLBACK(scalbln)
_GLIBCXX_SIMD_MATH_FALLBACK(cbrt)
_GLIBCXX_SIMD_MATH_FALLBACK(pow)
_GLIBCXX_SIMD_MATH_FALLBACK(erf)
_GLIBCXX_SIMD_MATH_FALLBACK(erfc)
_GLIBCXX_SIMD_MATH_FALLBACK(lgamma)
_GLIBCXX_SIMD_MATH_FALLBACK(tgamma)
_GLIBCXX_SIMD_MATH_FALLBACK_FIXEDRET(long, lrint)
_GLIBCXX_SIMD_MATH_FALLBACK_FIXEDRET(long long, llrint)
_GLIBCXX_SIMD_MATH_FALLBACK_FIXEDRET(long, lround)
_GLIBCXX_SIMD_MATH_FALLBACK_FIXEDRET(long long, llround)
_GLIBCXX_SIMD_MATH_FALLBACK(fmod)
_GLIBCXX_SIMD_MATH_FALLBACK(remainder)
template <typename _Tp, size_t _Np>
static _SveSimdWrapper<_Tp, _Np>
_S_remquo(const _SveSimdWrapper<_Tp, _Np> __x, const _SveSimdWrapper<_Tp, _Np> __y,
__fixed_size_storage_t<int, _Np>* __z)
{
_SveSimdWrapper<_Tp, _Np> __r{};
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
int __tmp;
__r._M_set(__i, remquo(__x[__i], __y[__i], &__tmp));
__z->_M_set(__i, __tmp);
});
return __r;
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static __fixed_size_storage_t<int, _Np>
_S_fpclassify(_SveSimdWrapper<_Tp, _Np> __x)
{
__fixed_size_storage_t<int, _Np> __r{};
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
__r._M_set(__i, std::fpclassify(__x[__i]));
});
return __r;
}
// copysign in simd_math.h
_GLIBCXX_SIMD_MATH_FALLBACK(nextafter)
_GLIBCXX_SIMD_MATH_FALLBACK(fdim)
#undef _GLIBCXX_SIMD_MATH_FALLBACK
#undef _GLIBCXX_SIMD_MATH_FALLBACK_FIXEDRET
template <typename _Tp, size_t _Np, typename _Op>
static constexpr _MaskMember<_Tp>
__fp_cmp(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y, _Op __op)
{
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
using _VI = __sve_vector_type_t<_Ip, _Np>;
using _WI = _SveSimdWrapper<_Ip, _Np>;
const _WI __fmv = __sve_vector_type<_Ip, _Np>::__sve_broadcast(__finite_max_v<_Ip>);
const _WI __zerov = __sve_vector_type<_Ip, _Np>::__sve_broadcast(0);
const _WI __xn = _VI(__sve_reinterpret_cast<_Ip>(__x));
const _WI __yn = _VI(__sve_reinterpret_cast<_Ip>(__y));
const _WI __xp
= svsel(_S_less(__xn, __zerov), _S_unary_minus(_WI(_S_bit_and(__xn, __fmv))), __xn);
const _WI __yp
= svsel(_S_less(__yn, __zerov), _S_unary_minus(_WI(_S_bit_and(__yn, __fmv))), __yn);
return svbic_z(__sve_vector_type<_Ip, _Np>::__sve_active_mask(), __op(__xp, __yp)._M_data,
_SuperImpl::_S_isunordered(__x, __y)._M_data);
}
template <typename _Tp, size_t _Np>
static constexpr _MaskMember<_Tp>
_S_isgreater(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y) noexcept
{ return __fp_cmp(__x, __y, [](auto __xp, auto __yp) { return _S_less(__yp, __xp); }); }
template <typename _Tp, size_t _Np>
static constexpr _MaskMember<_Tp>
_S_isgreaterequal(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y) noexcept
{ return __fp_cmp(__x, __y, [](auto __xp, auto __yp) { return _S_less_equal(__yp, __xp); }); }
template <typename _Tp, size_t _Np>
static constexpr _MaskMember<_Tp>
_S_isless(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y) noexcept
{ return __fp_cmp(__x, __y, [](auto __xp, auto __yp) { return _S_less(__xp, __yp); }); }
template <typename _Tp, size_t _Np>
static constexpr _MaskMember<_Tp>
_S_islessequal(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y) noexcept
{ return __fp_cmp(__x, __y, [](auto __xp, auto __yp) { return _S_less_equal(__xp, __yp); }); }
template <typename _Tp, size_t _Np>
static constexpr _MaskMember<_Tp>
_S_islessgreater(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y) noexcept
{
return svbic_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(),
_SuperImpl::_S_not_equal_to(__x, __y)._M_data,
_SuperImpl::_S_isunordered(__x, __y)._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_abs(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svabs_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_fabs(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svabs_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_sqrt(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svsqrt_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_ldexp(_SveSimdWrapper<_Tp, _Np> __x, __fixed_size_storage_t<int, _Np> __y) noexcept
{
auto __sve_register = __y.first;
if constexpr (std::is_same_v<_Tp, float>)
return svscale_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data,
__sve_register._M_data);
else
{
__sve_vector_type_t<int64_t, _Np> __sve_d_register = svunpklo(__sve_register);
return svscale_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data,
__sve_d_register);
}
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_fma(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y,
_SveSimdWrapper<_Tp, _Np> __z)
{
return svmad_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data,
__z._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_fmax(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
return svmaxnm_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_fmin(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
return svminnm_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _MaskMember<_Tp>
_S_isfinite([[maybe_unused]] _SveSimdWrapper<_Tp, _Np> __x)
{
#if __FINITE_MATH_ONLY__
return __sve_vector_type_t<_Tp, _Np>::__sve_all_true_mask();
#else
// if all exponent bits are set, __x is either inf or NaN
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
const __sve_vector_type_t<_Ip, _Np> __absn = __sve_reinterpret_cast<_Ip>(_S_abs(__x));
const __sve_vector_type_t<_Ip, _Np> __maxn
= __sve_reinterpret_cast<_Ip>(
__sve_vector_type<_Tp, _Np>::__sve_broadcast(__finite_max_v<_Tp>));
return _S_less_equal(_SveSimdWrapper<_Ip, _Np>{__absn}, _SveSimdWrapper<_Ip, _Np>{__maxn});
#endif
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _MaskMember<_Tp>
_S_isinf([[maybe_unused]] _SveSimdWrapper<_Tp, _Np> __x)
{
#if __FINITE_MATH_ONLY__
return {}; // false
#else
return _S_equal_to<_Tp, _Np>(_S_abs(__x), _S_broadcast(__infinity_v<_Tp>));
#endif
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _MaskMember<_Tp>
_S_isnan([[maybe_unused]] _SveSimdWrapper<_Tp, _Np> __x)
{
#if __FINITE_MATH_ONLY__
return {}; // false
#else
return svcmpuo(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __x._M_data);
#endif
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _MaskMember<_Tp>
_S_isnormal(_SveSimdWrapper<_Tp, _Np> __x)
{
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
using _V = __sve_vector_type_t<_Ip, _Np>;
using _VW = _SveSimdWrapper<_Ip, _Np>;
const _V __absn = __sve_reinterpret_cast<_Ip>(_S_abs(__x));
const _V __minn = __sve_reinterpret_cast<_Ip>(
__sve_vector_type<_Tp, _Np>::__sve_broadcast(__norm_min_v<_Tp>));
#if __FINITE_MATH_ONLY__
return _S_greater_equal(_VW{__absn}, _VW{__minn});
#else
const _V __maxn = __sve_reinterpret_cast<_Ip>(
__sve_vector_type<_Tp, _Np>::__sve_broadcast(__finite_max_v<_Tp>));
return _MaskImpl::_S_bit_and(_S_less_equal(_VW{__minn}, _VW{__absn}),
_S_less_equal(_VW{__absn}, _VW{__maxn}));
#endif
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _MaskMember<_Tp>
_S_signbit(_SveSimdWrapper<_Tp, _Np> __x)
{
using _Ip = __get_sve_value_type_t<__int_for_sizeof_t<_Tp>>;
using _V = __sve_vector_type_t<_Ip, _Np>;
using _VW = _SveSimdWrapper<_Ip, _Np>;
const _V __xn = __sve_reinterpret_cast<_Ip>(__x);
const _V __zeron = __sve_vector_type<_Ip, _Np>::__sve_broadcast(0);
return _S_less(_VW{__xn}, _VW{__zeron});
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _MaskMember<_Tp>
_S_isunordered(_SveSimdWrapper<_Tp, _Np> __x, _SveSimdWrapper<_Tp, _Np> __y)
{
return svcmpuo(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data, __y._M_data);
}
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_nearbyint(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svrinti_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_rint(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return _SuperImpl::_S_nearbyint(__x); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_trunc(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svrintz_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_round(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svrinta_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_floor(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svrintm_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static _SveSimdWrapper<_Tp, _Np>
_S_ceil(_SveSimdWrapper<_Tp, _Np> __x) noexcept
{ return svrintp_z(__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __x._M_data); }
template <typename _Tp, size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_masked_assign(_SveMaskWrapper<_Bits, _Np> __k, _SveSimdWrapper<_Tp, _Np>& __lhs,
__type_identity_t<_SveSimdWrapper<_Tp, _Np>> __rhs)
{ __lhs = _CommonImpl::_S_blend(__k, __lhs, __rhs); }
template <typename _Tp, size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_masked_assign(_SveMaskWrapper<_Bits, _Np> __k, _SveSimdWrapper<_Tp, _Np>& __lhs,
__type_identity_t<_Tp> __rhs)
{ __lhs = _CommonImpl::_S_blend(__k, __lhs, __data(simd<_Tp, _Abi>(__rhs))); }
template <typename _Op, typename _Tp, size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_masked_cassign(const _SveMaskWrapper<_Bits, _Np> __k, _SveSimdWrapper<_Tp, _Np>& __lhs,
const __type_identity_t<_SveSimdWrapper<_Tp, _Np>> __rhs, _Op __op)
{
__lhs = _CommonImpl::_S_blend(__k, __lhs,
_SveSimdWrapper<_Tp, _Np>(__op(_SuperImpl{}, __lhs, __rhs)));
}
template <typename _Op, typename _Tp, size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_masked_cassign(const _SveMaskWrapper<_Bits, _Np> __k, _SveSimdWrapper<_Tp, _Np>& __lhs,
const __type_identity_t<_Tp> __rhs, _Op __op)
{ _S_masked_cassign(__k, __lhs, _S_broadcast(__rhs), __op); }
template <typename _Tp, size_t _Np, typename _Up>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_set(_SveSimdWrapper<_Tp, _Np>& __v, int __i, _Up&& __x) noexcept
{ __v._M_set(__i, static_cast<_Up&&>(__x)); }
template <template <typename> class _Op, typename _Tp, size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveSimdWrapper<_Tp, _Np>
_S_masked_unary(const _SveMaskWrapper<_Bits, _Np> __k, const _SveSimdWrapper<_Tp, _Np> __v)
{
auto __vv = simd<_Tp, _Abi>{__private_init, __v};
_Op<decltype(__vv)> __op;
return _CommonImpl::_S_blend(__k, __v, __data(__op(__vv)));
}
};
template <typename _Abi, typename>
struct _MaskImplSve
{
template <typename _Tp>
using _MaskMember = typename _Abi::template _MaskMember<_Tp>;
template <typename _Tp>
using _TypeTag = _Tp*;
template <typename _Tp>
static constexpr size_t _S_size = simd_size_v<_Tp, _Abi>;
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_broadcast(bool __x)
{
constexpr size_t _Np = simd_size_v<_Tp, _Abi>;
__sve_bool_type __tr = __sve_vector_type<_Tp, _Np>::__sve_active_mask();
__sve_bool_type __fl = svpfalse_b();
return __x ? __tr : __fl;
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_load(const bool* __mem)
{
constexpr size_t _Np = simd_size_v<_Tp, _Abi>;
const uint8_t* __p = reinterpret_cast<const uint8_t*>(__mem);
__sve_bool_type __u8_active_mask = __sve_vector_type<uint8_t, _Np>::__sve_active_mask();
__sve_vector_type_t<uint8_t, _Np> __u8_vec_mask_load = svld1(__u8_active_mask, __p);
__sve_bool_type __u8_mask = svcmpne(__u8_active_mask, __u8_vec_mask_load, 0);
__sve_bool_type __tp_mask = __u8_mask;
for (size_t __up_size = 1; __up_size != sizeof(_Tp); __up_size *= 2)
{
__tp_mask = svunpklo(__tp_mask);
}
_SveMaskWrapper<sizeof(_Tp), simd_size_v<_Tp, _Abi>> __r{__tp_mask};
return __r;
}
template <size_t _Bits, size_t _Np>
static inline _SveMaskWrapper<_Bits, _Np>
_S_masked_load(_SveMaskWrapper<_Bits, _Np> __merge, _SveMaskWrapper<_Bits, _Np> __mask,
const bool* __mem) noexcept
{
_SveMaskWrapper<_Bits, _Np> __r;
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
if (__mask[__i])
__r._M_set(__i, __mem[__i]);
else
__r._M_set(__i, __merge[__i]);
});
return __r;
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_store(_SveMaskWrapper<_Bits, _Np> __v, bool* __mem) noexcept
{
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
__mem[__i] = __v[__i];
});
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr void
_S_masked_store(const _SveMaskWrapper<_Bits, _Np> __v, bool* __mem,
const _SveMaskWrapper<_Bits, _Np> __k) noexcept
{
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
if (__k[__i])
__mem[__i] = __v[__i];
});
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SanitizedBitMask<_Np>
_S_to_bits(_SveMaskWrapper<_Bits, _Np> __x)
{
_ULLong __r = 0;
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
__r |= _ULLong(__x[__i]) << __i;
});
return __r;
}
template <size_t _Np, typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static _MaskMember<_Tp>
_S_from_bitmask(_SanitizedBitMask<_Np> __bits, _TypeTag<_Tp>)
{
_SveMaskWrapper<sizeof(_Tp), _Np> __r;
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
__r._M_set(__i, __bits[__i]);
});
return __r;
}
template <typename _Tp, typename _Up, typename _UAbi>
_GLIBCXX_SIMD_INTRINSIC static constexpr auto
_S_convert(simd_mask<_Up, _UAbi> __x)
{
using _R = _SveMaskWrapper<sizeof(_Tp), simd_size_v<_Tp, _Abi>>;
if constexpr (__is_scalar_abi<_UAbi>())
{
_R __r{__sve_bool_type(svpfalse())};
__r._M_set(0, __data(__x));
return __r;
}
if constexpr (__is_sve_abi<_UAbi>())
{
if constexpr (sizeof(_Up) == sizeof(_Tp))
return __data(__x);
if constexpr (sizeof(_Up) < sizeof(_Tp))
{
__sve_bool_type __xmdata = __data(__x)._M_data;
__sve_bool_type __r = __xmdata;
for (size_t __up_size = sizeof(_Up); __up_size != sizeof(_Tp); __up_size *= 2)
{
__r = svunpklo(__r);
}
return _R{__r};
}
else
{
_R __r{__sve_bool_type(svpfalse())};
constexpr size_t __min_size
= std::min(simd_size_v<_Tp, _Abi>, simd_mask<_Up, _UAbi>::size());
__execute_n_times<__min_size>(
[&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA { __r._M_set(__i, __x[__i]); });
return __r;
}
}
if constexpr (__is_neon_abi<_UAbi>())
{
_R __r{__sve_bool_type(svpfalse())};
constexpr size_t __min_size
= std::min(simd_size_v<_Tp, _Abi>, simd_mask<_Up, _UAbi>::size());
__execute_n_times<__min_size>(
[&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA { __r._M_set(__i, __x[__i]); });
return __r;
}
if constexpr (__is_fixed_size_abi<_UAbi>())
{
return _S_convert<_Tp>(__data(__x));
}
return _R{};
}
template <typename _Tp, size_t _Np, bool _Sanitized>
_GLIBCXX_SIMD_INTRINSIC static constexpr _MaskMember<_Tp>
_S_convert(_BitMask<_Np, _Sanitized> __x)
{
_MaskMember<_Tp> __r{};
__execute_n_times<_Np>([&](auto __i) _GLIBCXX_SIMD_ALWAYS_INLINE_LAMBDA {
__r._M_set(__i, __x[__i]);
});
return __r;
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveMaskWrapper<_Bits, _Np>
_S_logical_and(const _SveMaskWrapper<_Bits, _Np>& __x, const _SveMaskWrapper<_Bits, _Np>& __y)
{
return svand_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveMaskWrapper<_Bits, _Np>
_S_logical_or(const _SveMaskWrapper<_Bits, _Np>& __x, const _SveMaskWrapper<_Bits, _Np>& __y)
{
return svorr_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveMaskWrapper<_Bits, _Np>
_S_bit_not(const _SveMaskWrapper<_Bits, _Np>& __x)
{
return svnot_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__x._M_data);
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveMaskWrapper<_Bits, _Np>
_S_bit_and(const _SveMaskWrapper<_Bits, _Np>& __x, const _SveMaskWrapper<_Bits, _Np>& __y)
{
return svand_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveMaskWrapper<_Bits, _Np>
_S_bit_or(const _SveMaskWrapper<_Bits, _Np>& __x, const _SveMaskWrapper<_Bits, _Np>& __y)
{
return svorr_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static constexpr _SveMaskWrapper<_Bits, _Np>
_S_bit_xor(const _SveMaskWrapper<_Bits, _Np>& __x, const _SveMaskWrapper<_Bits, _Np>& __y)
{
return sveor_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__x._M_data, __y._M_data);
}
template <size_t _Bits, size_t _Np>
static constexpr void
_S_set(_SveMaskWrapper<_Bits, _Np>& __k, int __i, bool __x) noexcept
{
auto __index = svcmpeq(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__sve_mask_type<_Bits>::__index0123,
typename __sve_mask_type<_Bits>::__sve_mask_uint_type(__i));
if (__x)
__k._M_data = svorr_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__k._M_data, __index);
else
__k._M_data = svbic_z(_SveMaskWrapper<_Bits, _Np>::_BuiltinSveVectorType::__sve_active_mask(),
__k._M_data, __index);
}
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static void
_S_masked_assign(_SveMaskWrapper<_Bits, _Np> __k, _SveMaskWrapper<_Bits, _Np>& __lhs,
_SveMaskWrapper<_Bits, _Np> __rhs)
{ __lhs._M_data = svsel(__k._M_data, __rhs._M_data, __lhs._M_data); }
template <size_t _Bits, size_t _Np>
_GLIBCXX_SIMD_INTRINSIC static void
_S_masked_assign(_SveMaskWrapper<_Bits, _Np> __k, _SveMaskWrapper<_Bits, _Np>& __lhs,
bool __rhs)
{
__lhs._M_data
= svsel(__k._M_data, _S_broadcast<__int_with_sizeof_t<_Bits>>(__rhs), __lhs._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static int
_S_popcount(simd_mask<_Tp, _Abi> __k)
{
constexpr size_t _Np = simd_size_v<_Tp, _Abi>;
return __sve_mask_type<sizeof(_Tp)>::__sve_mask_active_count(
__sve_vector_type<_Tp, _Np>::__sve_active_mask(), __k._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static bool
_S_all_of(simd_mask<_Tp, _Abi> __k)
{ return _S_popcount(__k) == simd_size_v<_Tp, _Abi>; }
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static bool
_S_any_of(simd_mask<_Tp, _Abi> __k)
{
return svptest_any(__sve_vector_type<_Tp, simd_size_v<_Tp, _Abi>>::__sve_active_mask(),
__k._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static bool
_S_none_of(simd_mask<_Tp, _Abi> __k)
{
return !svptest_any(__sve_vector_type<_Tp, simd_size_v<_Tp, _Abi>>::__sve_active_mask(),
__k._M_data);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static bool
_S_some_of(simd_mask<_Tp, _Abi> __k)
{
int __msk_count = _S_popcount(__k);
return (__msk_count > 0) && (__msk_count < (int) simd_size_v<_Tp, _Abi>);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static int
_S_find_first_set(simd_mask<_Tp, _Abi> __k)
{
return svclastb(svpfirst(__k._M_data, svpfalse()),
-1, __sve_mask_type<sizeof(_Tp)>::__index0123);
}
template <typename _Tp>
_GLIBCXX_SIMD_INTRINSIC static int
_S_find_last_set(simd_mask<_Tp, _Abi> __k)
{ return svclastb(__k._M_data, -1, __sve_mask_type<sizeof(_Tp)>::__index0123); }
};
_GLIBCXX_SIMD_END_NAMESPACE
#endif // __cplusplus >= 201703L
#endif // _GLIBCXX_EXPERIMENTAL_SIMD_SVE_H_
// vim: sw=2 noet ts=8 sts=2 tw=100
std::true_type
__bool_constant< true > true_type
The type used as a compile-time boolean with true value.
Definition type_traits:111
std::enable_if_t
typename enable_if< _Cond, _Tp >::type enable_if_t
Alias template for enable_if.
Definition type_traits:2696
std::begin
_Tp * begin(valarray< _Tp > &__va) noexcept
Return an iterator pointing to the first element of the valarray.
Definition valarray:1227
std::min
constexpr const _Tp & min(const _Tp &, const _Tp &)
This does what you think it does.
Definition stl_algobase.h:233
std::auto_ptr_ref
Definition auto_ptr.h:49