invoke_simd.hpp

Go to the documentation of this file.

//==------ invoke_simd.hpp - SYCL invoke_simd extension --*- C++ -*---------==//
//
// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
// See https://llvm.org/LICENSE.txt for license information.
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
//
// ===--------------------------------------------------------------------=== //
// Implemenation of the sycl_ext_oneapi_invoke_simd extension.
// https://github.com/intel/llvm/blob/sycl/sycl/doc/extensions/proposed/sycl_ext_oneapi_invoke_simd.asciidoc
// ===--------------------------------------------------------------------=== //
 
#pragma once
 
#include <sycl/ext/oneapi/experimental/detail/invoke_simd_types.hpp>
#include <sycl/ext/oneapi/experimental/uniform.hpp>
 
#include <sycl/detail/boost/mp11.hpp>
#include <sycl/sub_group.hpp>
 
#include <functional>
 
// TODOs:
// * (a) TODO enforce constness of a functor/lambda's () operator
// * (b) TODO support lambdas and functors in BE
 
template <bool IsFunc, class SpmdRet, class HelperFunc,
          class... UserSimdFuncAndSpmdArgs, class = std::enable_if_t<!IsFunc>>
__DPCPP_SYCL_EXTERNAL __regcall SpmdRet
__builtin_invoke_simd(HelperFunc helper, const void *obj,
                      UserSimdFuncAndSpmdArgs... args)
#ifdef __SYCL_DEVICE_ONLY__
    ;
#else
{
  // __builtin_invoke_simd is not supported on the host device yet
  throw sycl::exception(sycl::errc::feature_not_supported,
                        "__builtin_invoke_simd is not supported on host");
}
#endif // __SYCL_DEVICE_ONLY__
 
template <bool IsFunc, class SpmdRet, class HelperFunc,
          class... UserSimdFuncAndSpmdArgs, class = std::enable_if_t<IsFunc>>
__DPCPP_SYCL_EXTERNAL __regcall SpmdRet
__builtin_invoke_simd(HelperFunc helper, UserSimdFuncAndSpmdArgs... args)
#ifdef __SYCL_DEVICE_ONLY__
    ;
#else
{
  // __builtin_invoke_simd is not supported on the host device yet
  throw sycl::exception(sycl::errc::feature_not_supported,
                        "__builtin_invoke_simd is not supported on host");
}
#endif // __SYCL_DEVICE_ONLY__
 
namespace sycl {
inline namespace _V1 {
 
namespace ext::oneapi::experimental {
 
// --- Helpers
namespace detail {
 
namespace __MP11_NS = sycl::detail::boost::mp11;
 
// This structure performs the SPMD-to-SIMD parameter type conversion as defined
// by the spec.
template <class T, int N, class = void> struct spmd2simd;
// * `uniform<T>` converts to `T`
template <class T, int N> struct spmd2simd<uniform<T>, N> {
  using type = T;
};
// * tuple of types converts to tuple of converted tuple element types.
template <class... T, int N> struct spmd2simd<std::tuple<T...>, N> {
  using type = std::tuple<typename spmd2simd<T, N>::type...>;
};
// * arithmetic type converts to a simd vector with this element type and the
//   width equal to caller's subgroup size and passed as the `N` template
//   argument.
template <class T, int N>
struct spmd2simd<T, N, std::enable_if_t<std::is_arithmetic_v<T>>> {
  using type = simd<T, N>;
};
 
// * bool converts to `simd_mask` with a user specified element type.
// Arbitrarily use unsigned char for the element type for subgroup size
// deduction and rely on the implicit conversion operator for the the actual
// user type.
template <int N> struct spmd2simd<bool, N> {
  using type = simd_mask<unsigned char, N>;
};
 
// This structure performs the SIMD-to-SPMD return type conversion as defined
// by the spec.
template <class, class = void> struct simd2spmd;
// * `uniform<T>` stays the same
template <class T> struct simd2spmd<uniform<T>> {
  using type = uniform<T>;
};
// * `simd<T, N>` converts to T
template <class T, int N> struct simd2spmd<simd<T, N>> {
  using type = T;
};
// * tuple of types converts to tuple of converted tuple element types.
template <class... T> struct simd2spmd<std::tuple<T...>> {
  using type = std::tuple<typename simd2spmd<T>::type...>;
};
// * arithmetic type T converts to `uniform<T>`
template <class T>
struct simd2spmd<T, std::enable_if_t<std::is_arithmetic_v<T>>> {
  using type = uniform<T>;
};
 
// * `simd_mask` converts to bool
template <class T, int N> struct simd2spmd<simd_mask<T, N>> {
  using type = bool;
};
 
template <> struct simd2spmd<void> { using type = void; };
 
// Determine number of elements in a simd type.
template <class T> struct simd_size {
  static constexpr int value = 1; // 1 element in any type by default
};
 
// * Specialization for the simd type.
template <class T, int N> struct simd_size<simd<T, N>> {
  static constexpr int value = N;
};
 
// Check if given type is uniform.
template <class T> struct is_uniform_type : std::false_type {};
template <class T> struct is_uniform_type<uniform<T>> : std::true_type {
  using type = T;
};
 
// Check if given type is simd or simd_mask.
template <class T> struct is_simd_or_mask_type : std::false_type {};
template <class T, int N>
struct is_simd_or_mask_type<simd<T, N>> : std::true_type {};
template <class T, int N>
struct is_simd_or_mask_type<simd_mask<T, N>> : std::true_type {};
 
// Checks if all the types in the parameter pack are uniform<T>.
template <class... SpmdArgs> struct all_uniform_types {
  constexpr operator bool() {
    using TypeList = __MP11_NS::mp_list<SpmdArgs...>;
    return __MP11_NS::mp_all_of<TypeList, is_uniform_type>::value;
  }
};
 
// "Unwraps" a value of the `uniform` type (used before passing to SPMD
//  arguments to the __builtin_invoke_simd):
// - the case when there is nothing to unwrap
template <typename T> struct unwrap_uniform {
  static auto impl(T val) { return val; }
  using type = T;
};
 
// - the real unwrapping case
template <typename T> struct unwrap_uniform<uniform<T>> {
  static T impl(uniform<T> val) { return val; }
  using type = T;
};
 
// Verify the callee return type matches the subgroup size as is required by the
// spec. For example: simd<int, 8> foo(simd<int,16>); The return type vector
// length (8) does not match the subgroup size (16).
template <auto SgSize, typename SimdRet>
constexpr void verify_return_type_matches_sg_size() {
  if constexpr (is_simd_or_mask_type<SimdRet>::value) {
    constexpr auto RetVecLength = SimdRet::size();
    static_assert(RetVecLength == SgSize,
                  "invoke_simd callee return type vector length must match "
                  "kernel subgroup size");
  }
}
 
// Deduces subgroup size of the caller based on given SIMD callable and
// corresponding SPMD arguments it is being invoke with via invoke_simd.
// Basically, for each supported subgroup size, this meta-function finds out if
// the callable can be invoked by C++ rules given the SPMD arguments transformed
// as prescribed by the spec assuming this subgroup size. One and only one
// subgroup size should conform.
template <class SimdCallable, class... SpmdArgs> struct sg_size {
  template <class N>
  using IsInvocableSgSize = __MP11_NS::mp_bool<std::is_invocable_v<
      SimdCallable, typename spmd2simd<SpmdArgs, N::value>::type...>>;
 
  __DPCPP_SYCL_EXTERNAL constexpr operator int() {
    using SupportedSgSizes = __MP11_NS::mp_list_c<int, 1, 2, 4, 8, 16, 32>;
    using InvocableSgSizes =
        __MP11_NS::mp_copy_if<SupportedSgSizes, IsInvocableSgSize>;
    constexpr auto found_invoke_simd_target =
        __MP11_NS::mp_empty<InvocableSgSizes>::value != 1;
    if constexpr (found_invoke_simd_target) {
      static_assert((__MP11_NS::mp_size<InvocableSgSizes>::value == 1) &&
                    "multiple invoke_simd targets found");
      return __MP11_NS::mp_front<InvocableSgSizes>::value;
    }
    static_assert(
        found_invoke_simd_target,
        "No callable invoke_simd target found. Confirm the "
        "invoke_simd invocation argument types are convertible to the "
        "invoke_simd target argument types");
  }
};
 
// Determine the return type of a SIMD callable.
template <int N, class SimdCallable, class... SpmdArgs>
using SimdRetType =
    std::invoke_result_t<SimdCallable,
                         typename spmd2simd<SpmdArgs, N>::type...>;
// Determine the return type of an invoke_simd based on the return type of a
// SIMD callable.
template <int N, class SimdCallable, class... SpmdArgs>
using SpmdRetType =
    typename simd2spmd<SimdRetType<N, SimdCallable, SpmdArgs...>>::type;
 
template <class SimdCallable, class... SpmdArgs>
static constexpr int get_sg_size() {
  if constexpr (all_uniform_types<SpmdArgs...>()) {
    using SimdRet = std::invoke_result_t<SimdCallable, SpmdArgs...>;
 
    if constexpr (is_simd_or_mask_type<SimdRet>::value) {
      return simd_size<SimdRet>::value;
    } else {
      // fully uniform function - subgroup size does not matter
      return 0;
    }
  } else {
    return sg_size<SimdCallable, SpmdArgs...>();
  }
}
 
// This function is a wrapper around a call to a functor with field or a lambda
// with captures. Note __regcall - this is needed for efficient argument
// forwarding.
template <int N, class Callable, class... T>
[[intel::device_indirectly_callable]] __DPCPP_SYCL_EXTERNAL __regcall detail::
    SimdRetType<N, Callable, T...>
    simd_obj_call_helper(const void *obj_ptr,
                         typename detail::spmd2simd<T, N>::type... simd_args) {
  auto f =
      *reinterpret_cast<const std::remove_reference_t<Callable> *>(obj_ptr);
  return f(simd_args...);
}
 
// This function is a wrapper around a call to a function.
template <int N, class Callable, class... T>
[[intel::device_indirectly_callable]] __DPCPP_SYCL_EXTERNAL __regcall detail::
    SimdRetType<N, Callable, T...>
    simd_func_call_helper(Callable f,
                          typename detail::spmd2simd<T, N>::type... simd_args) {
  return f(simd_args...);
}
 
#ifdef _GLIBCXX_RELEASE
#if _GLIBCXX_RELEASE < 10
#define __INVOKE_SIMD_USE_STD_IS_FUNCTION_WA
#endif // _GLIBCXX_RELEASE < 10
#endif // _GLIBCXX_RELEASE
 
#ifdef __INVOKE_SIMD_USE_STD_IS_FUNCTION_WA
// TODO This is a workaround for libstdc++ version 9 buggy behavior which
// returns false in the code below. Version 10 works fine. Once required
// minimum libstdc++ version is bumped to 10, this w/a should be removed.
//   template <class F> bool foo(F &&f) {
//     return std::is_function_v<std::remove_reference_t<F>>;
//   }
// where F is a function type with __regcall.
template <class F> struct is_regcall_function_ptr_or_ref : std::false_type {};
 
template <class Ret, class... Args>
struct is_regcall_function_ptr_or_ref<Ret(__regcall &)(Args...)>
    : std::true_type {};
 
template <class Ret, class... Args>
struct is_regcall_function_ptr_or_ref<Ret(__regcall *)(Args...)>
    : std::true_type {};
 
template <class Ret, class... Args>
struct is_regcall_function_ptr_or_ref<Ret(__regcall *&)(Args...)>
    : std::true_type {};
 
template <class F>
static constexpr bool is_regcall_function_ptr_or_ref_v =
    is_regcall_function_ptr_or_ref<F>::value;
#endif // __INVOKE_SIMD_USE_STD_IS_FUNCTION_WA
 
template <class Callable>
static constexpr bool is_function_ptr_or_ref_v =
    std::is_function_v<std::remove_pointer_t<std::remove_reference_t<Callable>>>
#ifdef __INVOKE_SIMD_USE_STD_IS_FUNCTION_WA
    || is_regcall_function_ptr_or_ref_v<Callable>
#endif // __INVOKE_SIMD_USE_STD_IS_FUNCTION_WA
    ;
 
template <typename Callable> struct remove_ref_from_func_ptr_ref_type {
  using type = Callable;
};
 
template <typename Ret, typename... Args>
struct remove_ref_from_func_ptr_ref_type<Ret (*&)(Args...)> {
  using type = Ret (*)(Args...);
};
 
template <typename Ret, typename... Args>
struct remove_ref_from_func_ptr_ref_type<Ret(__regcall *&)(Args...)> {
  using type = Ret(__regcall *)(Args...);
};
 
template <typename T>
using remove_ref_from_func_ptr_ref_type_t =
    typename remove_ref_from_func_ptr_ref_type<T>::type;
 
template <typename T> struct strip_regcall_from_function_ptr;
 
template <typename Ret, typename... Args>
struct strip_regcall_from_function_ptr<Ret (*)(Args...)> {
  using type = Ret (*)(Args...);
};
 
template <typename Ret, typename... Args>
struct strip_regcall_from_function_ptr<Ret(__regcall *)(Args...)> {
  using type = Ret (*)(Args...);
};
 
template <typename T>
using strip_regcall_from_function_ptr_t =
    typename strip_regcall_from_function_ptr<T>::type;
 
template <typename T> struct is_non_trivially_copyable_uniform {
  static constexpr bool value =
      is_uniform_type<T>::value &&
      !std::is_trivially_copyable_v<typename unwrap_uniform<T>::type>;
};
 
template <> struct is_non_trivially_copyable_uniform<void> {
  static constexpr bool value = false;
};
 
template <typename T>
inline constexpr bool is_non_trivially_copyable_uniform_v =
    is_non_trivially_copyable_uniform<T>::value;
 
template <typename Ret, typename... Args>
constexpr bool has_ref_arg(Ret (*)(Args...)) {
  return (... || std::is_reference_v<Args>);
}
 
template <typename Ret, typename... Args>
constexpr bool has_ref_ret(Ret (*)(Args...)) {
  return std::is_reference_v<Ret>;
}
 
template <typename Ret, typename... Args>
constexpr bool has_non_uniform_struct_ret(Ret (*)(Args...)) {
  return std::is_class_v<Ret> && !is_simd_or_mask_type<Ret>::value &&
         !is_uniform_type<Ret>::value;
}
 
template <typename Ret, typename... Args>
constexpr bool has_non_trivially_copyable_uniform_ret(Ret (*)(Args...)) {
  return is_non_trivially_copyable_uniform_v<Ret>;
}
 
template <class Callable> constexpr void verify_callable() {
  if constexpr (is_function_ptr_or_ref_v<Callable>) {
    using RemoveRef =
        remove_ref_from_func_ptr_ref_type_t<std::remove_reference_t<Callable>>;
    using FuncPtrType =
        std::conditional_t<std::is_pointer_v<RemoveRef>, RemoveRef,
                           std::add_pointer_t<RemoveRef>>;
    using FuncPtrNoCC = strip_regcall_from_function_ptr_t<FuncPtrType>;
    constexpr FuncPtrNoCC obj = {};
    constexpr bool callable_has_ref_ret = has_ref_ret(obj);
    static_assert(
        !callable_has_ref_ret,
        "invoke_simd does not support callables returning references");
    constexpr bool callable_has_ref_arg = has_ref_arg(obj);
    static_assert(
        !callable_has_ref_arg,
        "invoke_simd does not support callables with reference arguments");
#ifndef __INVOKE_SIMD_ENABLE_STRUCTS
    constexpr bool callable_has_non_uniform_struct_ret =
        has_non_uniform_struct_ret(obj);
    static_assert(!callable_has_non_uniform_struct_ret,
                  "invoke_simd does not support callables returning "
                  "non-uniform structures");
#endif
#ifdef __SYCL_DEVICE_ONLY__
    constexpr bool callable_has_uniform_non_trivially_copyable_ret =
        has_non_trivially_copyable_uniform_ret(obj);
    static_assert(!callable_has_uniform_non_trivially_copyable_ret,
                  "invoke_simd does not support callables returning uniforms "
                  "that are not trivially copyable");
#endif
  }
}
 
template <class... Ts>
constexpr void verify_no_uniform_non_trivially_copyable_args() {
#ifdef __SYCL_DEVICE_ONLY__
  constexpr bool has_non_trivially_copyable_uniform_arg =
      (... || is_non_trivially_copyable_uniform_v<Ts>);
  static_assert(!has_non_trivially_copyable_uniform_arg,
                "Uniform arguments must be trivially copyable");
#endif
}
 
template <class... Ts> constexpr void verify_no_non_uniform_struct_args() {
#if defined(__SYCL_DEVICE_ONLY__) && !defined(__INVOKE_SIMD_ENABLE_STRUCTS)
  constexpr bool has_non_uniform_struct_arg =
      (... || (std::is_class_v<Ts> && !is_simd_or_mask_type<Ts>::value &&
               !is_uniform_type<Ts>::value));
  static_assert(!has_non_uniform_struct_arg,
                "Structure arguments must be uniform");
#endif
}
 
template <class Callable, class... Ts>
constexpr void verify_valid_args_and_ret() {
  verify_no_uniform_non_trivially_copyable_args<Ts...>();
 
  verify_no_non_uniform_struct_args<Ts...>();
 
  verify_callable<Callable>();
}
 
} // namespace detail
 
// --- The main API
 
// TODO works only for functions and pointers to functions now,
// enable for lambda functions and functors.
template <class Callable, class... T>
__attribute__((always_inline)) auto invoke_simd(sycl::sub_group sg,
                                                Callable &&f, T... args) {
  // If the invoke_simd call site is fully uniform, then it does not matter
  // what the subgroup size is and arguments don't need widening and return
  // value does not need shrinking by this library or SPMD compiler, so 0
  // is fine in this case.
  detail::verify_valid_args_and_ret<Callable, T...>();
  constexpr int N = detail::get_sg_size<Callable, T...>();
  using RetSpmd = detail::SpmdRetType<N, Callable, T...>;
  detail::verify_return_type_matches_sg_size<
      N, detail::SimdRetType<N, Callable, T...>>();
  constexpr bool is_function = detail::is_function_ptr_or_ref_v<Callable>;
 
  if constexpr (is_function) {
    // The variables typed as pointer to a function become lvalue-reference
    // when passed to invoke_simd() as universal pointers. That creates an
    // additional indirection, which is resolved automatically by the compiler
    // for the caller side of __builtin_invoke_simd, but which must be resolved
    // manually during the creation of simd_func_call_helper.
    // The class remove_ref_from_func_ptr_ref_type is used removes that
    // unwanted indirection.
    return __builtin_invoke_simd<true /*function*/, RetSpmd>(
        detail::simd_func_call_helper<
            N, detail::remove_ref_from_func_ptr_ref_type_t<Callable>, T...>,
        f, detail::unwrap_uniform<T>::impl(args)...);
  } else {
    // TODO support functors and lambdas which are handled in this branch.
    // The limiting factor for now is that the LLVMIR data flow analysis
    // implemented in LowerInvokeSimd.cpp which, finds actual invoke_simd
    // target function, can't handle this case yet.
    return __builtin_invoke_simd<false /*functor/lambda*/, RetSpmd>(
        detail::simd_obj_call_helper<N, Callable, T...>, &f,
        detail::unwrap_uniform<T>::impl(args)...);
  }
// TODO Temporary macro and assert to enable API compilation testing.
// LowerInvokeSimd.cpp does not support this case yet.
#ifndef __INVOKE_SIMD_ENABLE_ALL_CALLABLES
  static_assert(is_function &&
                "invoke_simd does not support functors or lambdas yet");
#endif // __INVOKE_SIMD_ENABLE_ALL_CALLABLES
}
 
} // namespace ext::oneapi::experimental
} // namespace _V1
} // namespace sycl