// This file is part of Eigen, a lightweight C++ template library
// for linear algebra.
//
// Copyright (C) 2016 Gael Guennebaud <gael.guennebaud@inria.fr>
//
// 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_SPECIALFUNCTIONS_PACKETMATH_H
#define EIGEN_SPECIALFUNCTIONS_PACKETMATH_H

namespace Eigen {

namespace internal {

/** \internal \returns the ln(|gamma(\a a)|) (coeff-wise) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_DECLARE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS Packet
plgamma(const Packet &a) {
  using numext::lgamma;
  return lgamma(a);
}

/** \internal \returns the derivative of lgamma, psi(\a a) (coeff-wise) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_DECLARE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS Packet
pdigamma(const Packet &a) {
  using numext::digamma;
  return digamma(a);
}

/** \internal \returns the zeta function of two arguments (coeff-wise) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_DECLARE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS Packet
pzeta(const Packet& x, const Packet& q) {
  using numext::zeta;
  return zeta(x, q);
}

/** \internal \returns the polygamma function (coeff-wise) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_DECLARE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS Packet
ppolygamma(const Packet& n, const Packet& x) {
  using numext::polygamma;
  return polygamma(n, x);
}

/** \internal \returns the erf(\a a) (coeff-wise) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_DECLARE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS Packet
perf(const Packet& a) {
  using numext::erf;
  return erf(a);
}

/** \internal \returns the erfc(\a a) (coeff-wise) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_DECLARE_FUNCTION_ALLOWING_MULTIPLE_DEFINITIONS Packet
perfc(const Packet& a) {
  using numext::erfc;
  return erfc(a);
}

/** \internal \returns the incomplete gamma function igamma(\a a, \a x) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE Packet pigamma(const Packet& a,
                                                     const Packet& x) {
  using numext::igamma;
  return igamma(a, x);
}

/** \internal \returns the complementary incomplete gamma function igammac(\a a,
 * \a x) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE Packet pigammac(const Packet& a,
                                                      const Packet& x) {
  using numext::igammac;
  return igammac(a, x);
}

/** \internal \returns the complementary incomplete gamma function betainc(\a a,
 * \a b, \a x) */
template <typename Packet>
EIGEN_DEVICE_FUNC EIGEN_STRONG_INLINE Packet pbetainc(const Packet& a,
                                                      const Packet& b,
                                                      const Packet& x) {
  using numext::betainc;
  return betainc(a, b, x);
}

}  // end namespace internal

}  // end namespace Eigen

#endif  // EIGEN_SPECIALFUNCTIONS_PACKETMATH_H