diff options
author | Till Hoffmann <tillahoffmann@gmail.com> | 2016-04-01 14:35:21 +0100 |
---|---|---|
committer | Till Hoffmann <tillahoffmann@gmail.com> | 2016-04-01 14:35:21 +0100 |
commit | 57239f4a8149dbd603ad376e90a0a4574b846710 (patch) | |
tree | 4ff7041da6cd7121c19271bfb385c813ea66e0a8 /test | |
parent | dd5d390daf3a3a561a772b64f1b602e5f240bf8b (diff) |
Added polygamma function.
Diffstat (limited to 'test')
-rw-r--r-- | test/array.cpp | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/test/array.cpp b/test/array.cpp index 2f0b0f1b6..56d196923 100644 --- a/test/array.cpp +++ b/test/array.cpp @@ -331,6 +331,11 @@ template<typename ArrayType> void array_real(const ArrayType& m) VERIFY_IS_APPROX(numext::zeta(Scalar(3), Scalar(-2.5)), RealScalar(0.054102025820864097)); VERIFY_IS_EQUAL(numext::zeta(Scalar(1), Scalar(1.2345)), // The second scalar does not matter std::numeric_limits<RealScalar>::infinity()); + + // Check the polygamma against scipy.special.polygamma + VERIFY_IS_APPROX(numext::polygamma(Scalar(1), Scalar(2)), RealScalar(0.644934066848)); + VERIFY_IS_APPROX(numext::polygamma(Scalar(1), Scalar(3)), RealScalar(0.394934066848)); + VERIFY_IS_APPROX(numext::polygamma(Scalar(1), Scalar(25.5)), RealScalar(0.0399946696496)); { // Test various propreties of igamma & igammac. These are normalized |