bug #1232: refactor special functions as a new SpecialFunctions module, currently in unsupported/.

1 files changed, 0 insertions, 105 deletions

diff --git a/Eigen/src/Core/GlobalFunctions.h b/Eigen/src/Core/GlobalFunctions.h
index b9c3ec25b..879a93e6b 100644
--- a/Eigen/src/Core/GlobalFunctions.h
+++ b/Eigen/src/Core/GlobalFunctions.h
@@ -174,111 +174,6 @@ namespace Eigen
   }
 #endif
 
-  /** \cpp11 \returns an expression of the coefficient-wise igamma(\a a, \a x) to the given arrays.
-    *
-    * This function computes the coefficient-wise incomplete gamma function.
-    *
-    * \note This function supports only float and double scalar types in c++11 mode. To support other scalar types,
-    * or float/double in non c++11 mode, the user has to provide implementations of igammac(T,T) for any scalar
-    * type T to be supported.
-    *
-    * \sa Eigen::igammac(), Eigen::lgamma()
-    */
-  template<typename Derived,typename ExponentDerived>
-  inline const Eigen::CwiseBinaryOp<Eigen::internal::scalar_igamma_op<typename Derived::Scalar>, const Derived, const ExponentDerived>
-  igamma(const Eigen::ArrayBase<Derived>& a, const Eigen::ArrayBase<ExponentDerived>& x) 
-  {
-    return Eigen::CwiseBinaryOp<Eigen::internal::scalar_igamma_op<typename Derived::Scalar>, const Derived, const ExponentDerived>(
-      a.derived(),
-      x.derived()
-    );
-  }
-
-  /** \cpp11 \returns an expression of the coefficient-wise igammac(\a a, \a x) to the given arrays.
-    *
-    * This function computes the coefficient-wise complementary incomplete gamma function.
-    *
-    * \note This function supports only float and double scalar types in c++11 mode. To support other scalar types,
-    * or float/double in non c++11 mode, the user has to provide implementations of igammac(T,T) for any scalar
-    * type T to be supported.
-    *
-    * \sa Eigen::igamma(), Eigen::lgamma()
-    */
-  template<typename Derived,typename ExponentDerived>
-  inline const Eigen::CwiseBinaryOp<Eigen::internal::scalar_igammac_op<typename Derived::Scalar>, const Derived, const ExponentDerived>
-  igammac(const Eigen::ArrayBase<Derived>& a, const Eigen::ArrayBase<ExponentDerived>& x) 
-  {
-    return Eigen::CwiseBinaryOp<Eigen::internal::scalar_igammac_op<typename Derived::Scalar>, const Derived, const ExponentDerived>(
-      a.derived(),
-      x.derived()
-    );
-  }
-
-  /** \cpp11 \returns an expression of the coefficient-wise polygamma(\a n, \a x) to the given arrays.
-    *
-    * It returns the \a n -th derivative of the digamma(psi) evaluated at \c x.
-    *
-    * \note This function supports only float and double scalar types in c++11 mode. To support other scalar types,
-    * or float/double in non c++11 mode, the user has to provide implementations of polygamma(T,T) for any scalar
-    * type T to be supported.
-    *
-    * \sa Eigen::digamma()
-    */
-  // * \warning Be careful with the order of the parameters: x.polygamma(n) is equivalent to polygamma(n,x)
-  // * \sa ArrayBase::polygamma()
-  template<typename DerivedN,typename DerivedX>
-  inline const Eigen::CwiseBinaryOp<Eigen::internal::scalar_polygamma_op<typename DerivedX::Scalar>, const DerivedN, const DerivedX>
-  polygamma(const Eigen::ArrayBase<DerivedN>& n, const Eigen::ArrayBase<DerivedX>& x)
-  {
-    return Eigen::CwiseBinaryOp<Eigen::internal::scalar_polygamma_op<typename DerivedX::Scalar>, const DerivedN, const DerivedX>(
-      n.derived(),
-      x.derived()
-    );
-  }
-
-  /** \cpp11 \returns an expression of the coefficient-wise betainc(\a x, \a a, \a b) to the given arrays.
-    *
-    * This function computes the regularized incomplete beta function (integral).
-    *
-    * \note This function supports only float and double scalar types in c++11 mode. To support other scalar types,
-    * or float/double in non c++11 mode, the user has to provide implementations of betainc(T,T,T) for any scalar
-    * type T to be supported.
-    *
-    * \sa Eigen::betainc(), Eigen::lgamma()
-    */
-  template<typename ArgADerived, typename ArgBDerived, typename ArgXDerived>
-  inline const Eigen::CwiseTernaryOp<Eigen::internal::scalar_betainc_op<typename ArgXDerived::Scalar>, const ArgADerived, const ArgBDerived, const ArgXDerived>
-  betainc(const Eigen::ArrayBase<ArgADerived>& a, const Eigen::ArrayBase<ArgBDerived>& b, const Eigen::ArrayBase<ArgXDerived>& x)
-  {
-    return Eigen::CwiseTernaryOp<Eigen::internal::scalar_betainc_op<typename ArgXDerived::Scalar>, const ArgADerived, const ArgBDerived, const ArgXDerived>(
-      a.derived(),
-      b.derived(),
-      x.derived()
-    );
-  }
-
-
-  /** \returns an expression of the coefficient-wise zeta(\a x, \a q) to the given arrays.
-    *
-    * It returns the Riemann zeta function of two arguments \a x and \a q:
-    *
-    * \param x is the exposent, it must be > 1
-    * \param q is the shift, it must be > 0
-    *
-    * \note This function supports only float and double scalar types. To support other scalar types, the user has
-    * to provide implementations of zeta(T,T) for any scalar type T to be supported.
-    *
-    * \sa ArrayBase::zeta()
-    */
-  template<typename DerivedX,typename DerivedQ>
-  inline const Eigen::CwiseBinaryOp<Eigen::internal::scalar_zeta_op<typename DerivedX::Scalar>, const DerivedX, const DerivedQ>
-  zeta(const Eigen::ArrayBase<DerivedX>& x, const Eigen::ArrayBase<DerivedQ>& q)
-  {
-    return Eigen::CwiseBinaryOp<Eigen::internal::scalar_zeta_op<typename DerivedX::Scalar>, const DerivedX, const DerivedQ>(
-      x.derived(),
-      q.derived()
-    );
-  }
 
   namespace internal
   {