1 // This file is part of Eigen, a lightweight C++ template library
2 // for linear algebra.
3 //
4 // This Source Code Form is subject to the terms of the Mozilla
5 // Public License v. 2.0. If a copy of the MPL was not distributed
6 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
7
8 #ifndef EIGEN_SPECIALFUNCTIONS_BFLOAT16_H
9 #define EIGEN_SPECIALFUNCTIONS_BFLOAT16_H
10
11 namespace Eigen {
12 namespace numext {
13
14 #if EIGEN_HAS_C99_MATH
lgamma(const Eigen::bfloat16 & a)15 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 lgamma(const Eigen::bfloat16& a) {
16 return Eigen::bfloat16(Eigen::numext::lgamma(static_cast<float>(a)));
17 }
digamma(const Eigen::bfloat16 & a)18 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 digamma(const Eigen::bfloat16& a) {
19 return Eigen::bfloat16(Eigen::numext::digamma(static_cast<float>(a)));
20 }
zeta(const Eigen::bfloat16 & x,const Eigen::bfloat16 & q)21 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 zeta(const Eigen::bfloat16& x, const Eigen::bfloat16& q) {
22 return Eigen::bfloat16(Eigen::numext::zeta(static_cast<float>(x), static_cast<float>(q)));
23 }
polygamma(const Eigen::bfloat16 & n,const Eigen::bfloat16 & x)24 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 polygamma(const Eigen::bfloat16& n, const Eigen::bfloat16& x) {
25 return Eigen::bfloat16(Eigen::numext::polygamma(static_cast<float>(n), static_cast<float>(x)));
26 }
erf(const Eigen::bfloat16 & a)27 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 erf(const Eigen::bfloat16& a) {
28 return Eigen::bfloat16(Eigen::numext::erf(static_cast<float>(a)));
29 }
erfc(const Eigen::bfloat16 & a)30 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 erfc(const Eigen::bfloat16& a) {
31 return Eigen::bfloat16(Eigen::numext::erfc(static_cast<float>(a)));
32 }
ndtri(const Eigen::bfloat16 & a)33 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 ndtri(const Eigen::bfloat16& a) {
34 return Eigen::bfloat16(Eigen::numext::ndtri(static_cast<float>(a)));
35 }
igamma(const Eigen::bfloat16 & a,const Eigen::bfloat16 & x)36 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 igamma(const Eigen::bfloat16& a, const Eigen::bfloat16& x) {
37 return Eigen::bfloat16(Eigen::numext::igamma(static_cast<float>(a), static_cast<float>(x)));
38 }
39 template <>
igamma_der_a(const Eigen::bfloat16 & a,const Eigen::bfloat16 & x)40 EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 igamma_der_a(const Eigen::bfloat16& a, const Eigen::bfloat16& x) {
41 return Eigen::bfloat16(Eigen::numext::igamma_der_a(static_cast<float>(a), static_cast<float>(x)));
42 }
43 template <>
gamma_sample_der_alpha(const Eigen::bfloat16 & alpha,const Eigen::bfloat16 & sample)44 EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 gamma_sample_der_alpha(const Eigen::bfloat16& alpha, const Eigen::bfloat16& sample) {
45 return Eigen::bfloat16(Eigen::numext::gamma_sample_der_alpha(static_cast<float>(alpha), static_cast<float>(sample)));
46 }
igammac(const Eigen::bfloat16 & a,const Eigen::bfloat16 & x)47 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 igammac(const Eigen::bfloat16& a, const Eigen::bfloat16& x) {
48 return Eigen::bfloat16(Eigen::numext::igammac(static_cast<float>(a), static_cast<float>(x)));
49 }
betainc(const Eigen::bfloat16 & a,const Eigen::bfloat16 & b,const Eigen::bfloat16 & x)50 template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::bfloat16 betainc(const Eigen::bfloat16& a, const Eigen::bfloat16& b, const Eigen::bfloat16& x) {
51 return Eigen::bfloat16(Eigen::numext::betainc(static_cast<float>(a), static_cast<float>(b), static_cast<float>(x)));
52 }
53 #endif
54
55 } // end namespace numext
56 } // end namespace Eigen
57
58 #endif // EIGEN_SPECIALFUNCTIONS_BFLOAT16_H
59