blob: 553bcda6a65b0369564a931142ca86d6651b34c4 [file] [log] [blame]
Austin Schuh189376f2018-12-20 22:11:15 +11001// 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_HALF_H
9#define EIGEN_SPECIALFUNCTIONS_HALF_H
10
11namespace Eigen {
12namespace numext {
13
14#if EIGEN_HAS_C99_MATH
15template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half lgamma(const Eigen::half& a) {
16 return Eigen::half(Eigen::numext::lgamma(static_cast<float>(a)));
17}
18template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half digamma(const Eigen::half& a) {
19 return Eigen::half(Eigen::numext::digamma(static_cast<float>(a)));
20}
21template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half zeta(const Eigen::half& x, const Eigen::half& q) {
22 return Eigen::half(Eigen::numext::zeta(static_cast<float>(x), static_cast<float>(q)));
23}
24template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half polygamma(const Eigen::half& n, const Eigen::half& x) {
25 return Eigen::half(Eigen::numext::polygamma(static_cast<float>(n), static_cast<float>(x)));
26}
27template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half erf(const Eigen::half& a) {
28 return Eigen::half(Eigen::numext::erf(static_cast<float>(a)));
29}
30template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half erfc(const Eigen::half& a) {
31 return Eigen::half(Eigen::numext::erfc(static_cast<float>(a)));
32}
33template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half igamma(const Eigen::half& a, const Eigen::half& x) {
34 return Eigen::half(Eigen::numext::igamma(static_cast<float>(a), static_cast<float>(x)));
35}
36template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half igammac(const Eigen::half& a, const Eigen::half& x) {
37 return Eigen::half(Eigen::numext::igammac(static_cast<float>(a), static_cast<float>(x)));
38}
39template<> EIGEN_STRONG_INLINE EIGEN_DEVICE_FUNC Eigen::half betainc(const Eigen::half& a, const Eigen::half& b, const Eigen::half& x) {
40 return Eigen::half(Eigen::numext::betainc(static_cast<float>(a), static_cast<float>(b), static_cast<float>(x)));
41}
42#endif
43
44} // end namespace numext
45} // end namespace Eigen
46
47#endif // EIGEN_SPECIALFUNCTIONS_HALF_H