| // This file is part of Eigen, a lightweight C++ template library |
| // for linear algebra. |
| // |
| // Copyright (C) 2014 Benoit Steiner <benoit.steiner.goog@gmail.com> |
| // |
| // This Source Code Form is subject to the terms of the Mozilla |
| // Public License v. 2.0. If a copy of the MPL was not distributed |
| // with this file, You can obtain one at http://mozilla.org/MPL/2.0/. |
| |
| #ifndef EIGEN_CUDA_SPECIALFUNCTIONS_H |
| #define EIGEN_CUDA_SPECIALFUNCTIONS_H |
| |
| namespace Eigen { |
| |
| namespace internal { |
| |
| // Make sure this is only available when targeting a GPU: we don't want to |
| // introduce conflicts between these packet_traits definitions and the ones |
| // we'll use on the host side (SSE, AVX, ...) |
| #if defined(__CUDACC__) && defined(EIGEN_USE_GPU) |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 plgamma<float4>(const float4& a) { |
| return make_float4(lgammaf(a.x), lgammaf(a.y), lgammaf(a.z), lgammaf(a.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 |
| plgamma<double2>(const double2& a) { |
| using numext::lgamma; |
| return make_double2(lgamma(a.x), lgamma(a.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 pdigamma<float4>(const float4& a) { |
| using numext::digamma; |
| return make_float4(digamma(a.x), digamma(a.y), digamma(a.z), digamma(a.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 |
| pdigamma<double2>(const double2& a) { |
| using numext::digamma; |
| return make_double2(digamma(a.x), digamma(a.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 pzeta<float4>(const float4& x, |
| const float4& q) { |
| using numext::zeta; |
| return make_float4(zeta(x.x, q.x), zeta(x.y, q.y), zeta(x.z, q.z), |
| zeta(x.w, q.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 pzeta<double2>(const double2& x, |
| const double2& q) { |
| using numext::zeta; |
| return make_double2(zeta(x.x, q.x), zeta(x.y, q.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 |
| ppolygamma<float4>(const float4& n, const float4& x) { |
| using numext::polygamma; |
| return make_float4(polygamma(n.x, x.x), polygamma(n.y, x.y), |
| polygamma(n.z, x.z), polygamma(n.w, x.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 |
| ppolygamma<double2>(const double2& n, const double2& x) { |
| using numext::polygamma; |
| return make_double2(polygamma(n.x, x.x), polygamma(n.y, x.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 perf<float4>(const float4& a) { |
| return make_float4(erff(a.x), erff(a.y), erff(a.z), erff(a.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 perf<double2>(const double2& a) { |
| using numext::erf; |
| return make_double2(erf(a.x), erf(a.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 perfc<float4>(const float4& a) { |
| using numext::erfc; |
| return make_float4(erfc(a.x), erfc(a.y), erfc(a.z), erfc(a.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 perfc<double2>(const double2& a) { |
| using numext::erfc; |
| return make_double2(erfc(a.x), erfc(a.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 pigamma<float4>(const float4& a, |
| const float4& x) { |
| using numext::igamma; |
| return make_float4(igamma(a.x, x.x), igamma(a.y, x.y), igamma(a.z, x.z), |
| igamma(a.w, x.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 |
| pigamma<double2>(const double2& a, const double2& x) { |
| using numext::igamma; |
| return make_double2(igamma(a.x, x.x), igamma(a.y, x.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 pigammac<float4>(const float4& a, |
| const float4& x) { |
| using numext::igammac; |
| return make_float4(igammac(a.x, x.x), igammac(a.y, x.y), igammac(a.z, x.z), |
| igammac(a.w, x.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 |
| pigammac<double2>(const double2& a, const double2& x) { |
| using numext::igammac; |
| return make_double2(igammac(a.x, x.x), igammac(a.y, x.y)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE float4 pbetainc<float4>(const float4& a, |
| const float4& b, |
| const float4& x) { |
| using numext::betainc; |
| return make_float4(betainc(a.x, b.x, x.x), betainc(a.y, b.y, x.y), |
| betainc(a.z, b.z, x.z), betainc(a.w, b.w, x.w)); |
| } |
| |
| template <> |
| EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE double2 |
| pbetainc<double2>(const double2& a, const double2& b, const double2& x) { |
| using numext::betainc; |
| return make_double2(betainc(a.x, b.x, x.x), betainc(a.y, b.y, x.y)); |
| } |
| |
| #endif |
| |
| } // end namespace internal |
| |
| } // end namespace Eigen |
| |
| #endif // EIGEN_CUDA_SPECIALFUNCTIONS_H |