-
Notifications
You must be signed in to change notification settings - Fork 30
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge branch 'master' into implement-clip
- Loading branch information
Showing
12 changed files
with
1,609 additions
and
5 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
172 changes: 172 additions & 0 deletions
172
dpctl/tensor/libtensor/include/kernels/elementwise_functions/cbrt.hpp
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,172 @@ | ||
| //=== cbrt.hpp - Unary function CBRT ------ *-C++-*--/===// | ||
| // | ||
| // Data Parallel Control (dpctl) | ||
| // | ||
| // Copyright 2020-2023 Intel Corporation | ||
| // | ||
| // Licensed under the Apache License, Version 2.0 (the "License"); | ||
| // you may not use this file except in compliance with the License. | ||
| // You may obtain a copy of the License at | ||
| // | ||
| // http://www.apache.org/licenses/LICENSE-2.0 | ||
| // | ||
| // Unless required by applicable law or agreed to in writing, software | ||
| // distributed under the License is distributed on an "AS IS" BASIS, | ||
| // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. | ||
| // See the License for the specific language governing permissions and | ||
| // limitations under the License. | ||
| // | ||
| //===---------------------------------------------------------------------===// | ||
| /// | ||
| /// \file | ||
| /// This file defines kernels for elementwise evaluation of CBRT(x) | ||
| /// function that compute a square root. | ||
| //===---------------------------------------------------------------------===// | ||
|
|
||
| #pragma once | ||
| #include <CL/sycl.hpp> | ||
| #include <cmath> | ||
| #include <cstddef> | ||
| #include <cstdint> | ||
| #include <type_traits> | ||
|
|
||
| #include "kernels/elementwise_functions/common.hpp" | ||
|
|
||
| #include "utils/offset_utils.hpp" | ||
| #include "utils/type_dispatch.hpp" | ||
| #include "utils/type_utils.hpp" | ||
| #include <pybind11/pybind11.h> | ||
|
|
||
| namespace dpctl | ||
| { | ||
| namespace tensor | ||
| { | ||
| namespace kernels | ||
| { | ||
| namespace cbrt | ||
| { | ||
|
|
||
| namespace py = pybind11; | ||
| namespace td_ns = dpctl::tensor::type_dispatch; | ||
|
|
||
| template <typename argT, typename resT> struct CbrtFunctor | ||
| { | ||
|
|
||
| // is function constant for given argT | ||
| using is_constant = typename std::false_type; | ||
| // constant value, if constant | ||
| // constexpr resT constant_value = resT{}; | ||
| // is function defined for sycl::vec | ||
| using supports_vec = typename std::false_type; | ||
| // do both argTy and resTy support sugroup store/load operation | ||
| using supports_sg_loadstore = typename std::true_type; | ||
|
|
||
| resT operator()(const argT &in) const | ||
| { | ||
| return sycl::cbrt(in); | ||
| } | ||
| }; | ||
|
|
||
| template <typename argTy, | ||
| typename resTy = argTy, | ||
| unsigned int vec_sz = 4, | ||
| unsigned int n_vecs = 2> | ||
| using CbrtContigFunctor = elementwise_common:: | ||
| UnaryContigFunctor<argTy, resTy, CbrtFunctor<argTy, resTy>, vec_sz, n_vecs>; | ||
|
|
||
| template <typename argTy, typename resTy, typename IndexerT> | ||
| using CbrtStridedFunctor = elementwise_common:: | ||
| UnaryStridedFunctor<argTy, resTy, IndexerT, CbrtFunctor<argTy, resTy>>; | ||
|
|
||
| template <typename T> struct CbrtOutputType | ||
| { | ||
| using value_type = typename std::disjunction< // disjunction is C++17 | ||
| // feature, supported by DPC++ | ||
| td_ns::TypeMapResultEntry<T, sycl::half, sycl::half>, | ||
| td_ns::TypeMapResultEntry<T, float, float>, | ||
| td_ns::TypeMapResultEntry<T, double, double>, | ||
| td_ns::DefaultResultEntry<void>>::result_type; | ||
| }; | ||
|
|
||
| template <typename T1, typename T2, unsigned int vec_sz, unsigned int n_vecs> | ||
| class cbrt_contig_kernel; | ||
|
|
||
| template <typename argTy> | ||
| sycl::event cbrt_contig_impl(sycl::queue &exec_q, | ||
| size_t nelems, | ||
| const char *arg_p, | ||
| char *res_p, | ||
| const std::vector<sycl::event> &depends = {}) | ||
| { | ||
| return elementwise_common::unary_contig_impl< | ||
| argTy, CbrtOutputType, CbrtContigFunctor, cbrt_contig_kernel>( | ||
| exec_q, nelems, arg_p, res_p, depends); | ||
| } | ||
|
|
||
| template <typename fnT, typename T> struct CbrtContigFactory | ||
| { | ||
| fnT get() | ||
| { | ||
| if constexpr (std::is_same_v<typename CbrtOutputType<T>::value_type, | ||
| void>) { | ||
| fnT fn = nullptr; | ||
| return fn; | ||
| } | ||
| else { | ||
| fnT fn = cbrt_contig_impl<T>; | ||
| return fn; | ||
| } | ||
| } | ||
| }; | ||
|
|
||
| template <typename fnT, typename T> struct CbrtTypeMapFactory | ||
| { | ||
| /*! @brief get typeid for output type of std::cbrt(T x) */ | ||
| std::enable_if_t<std::is_same<fnT, int>::value, int> get() | ||
| { | ||
| using rT = typename CbrtOutputType<T>::value_type; | ||
| return td_ns::GetTypeid<rT>{}.get(); | ||
| } | ||
| }; | ||
|
|
||
| template <typename T1, typename T2, typename T3> class cbrt_strided_kernel; | ||
|
|
||
| template <typename argTy> | ||
| sycl::event | ||
| cbrt_strided_impl(sycl::queue &exec_q, | ||
| size_t nelems, | ||
| int nd, | ||
| const py::ssize_t *shape_and_strides, | ||
| const char *arg_p, | ||
| py::ssize_t arg_offset, | ||
| char *res_p, | ||
| py::ssize_t res_offset, | ||
| const std::vector<sycl::event> &depends, | ||
| const std::vector<sycl::event> &additional_depends) | ||
| { | ||
| return elementwise_common::unary_strided_impl< | ||
| argTy, CbrtOutputType, CbrtStridedFunctor, cbrt_strided_kernel>( | ||
| exec_q, nelems, nd, shape_and_strides, arg_p, arg_offset, res_p, | ||
| res_offset, depends, additional_depends); | ||
| } | ||
|
|
||
| template <typename fnT, typename T> struct CbrtStridedFactory | ||
| { | ||
| fnT get() | ||
| { | ||
| if constexpr (std::is_same_v<typename CbrtOutputType<T>::value_type, | ||
| void>) { | ||
| fnT fn = nullptr; | ||
| return fn; | ||
| } | ||
| else { | ||
| fnT fn = cbrt_strided_impl<T>; | ||
| return fn; | ||
| } | ||
| } | ||
| }; | ||
|
|
||
| } // namespace cbrt | ||
| } // namespace kernels | ||
| } // namespace tensor | ||
| } // namespace dpctl |
Oops, something went wrong.