commit 3cb0a237c195dd470ff5cd7a6e793b74d2d6815d
author Till Hoffmann <tillahoffmann@gmail.com> Fri Apr 01 17:51:39 2016 +0100
committer Till Hoffmann <tillahoffmann@gmail.com> Fri Apr 01 17:51:39 2016 +0100
tree afb4441de31d77de553c662eab8883be68c14ffe
parent 57239f4a8149dbd603ad376e90a0a4574b846710

Fixed suggestions by Eugene Brevdo.

test/array.cpp

diff --git a/test/array.cpp b/test/array.cpp
index 56d1969..8b0a347 100644
--- a/test/array.cpp
+++ b/test/array.cpp
@@ -332,10 +332,22 @@
     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
+    // Check the polygamma against scipy.special.polygamma examples
     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));
+    
+    // Check the polygamma function over a larger range of values
+    VERIFY_IS_APPROX(numext::polygamma(Scalar(17), Scalar(4.7)), RealScalar(293.334565435));
+    VERIFY_IS_APPROX(numext::polygamma(Scalar(31), Scalar(11.8)), RealScalar(0.445487887616));
+    VERIFY_IS_APPROX(numext::polygamma(Scalar(28), Scalar(17.7)), RealScalar(-2.47810300902e-07));
+    VERIFY_IS_APPROX(numext::polygamma(Scalar(8), Scalar(30.2)), RealScalar(-8.29668781082e-09));
+    /* The following tests only pass for doubles because floats cannot handle the large values of
+       the gamma function.
+    VERIFY_IS_APPROX(numext::polygamma(Scalar(42), Scalar(15.8)), RealScalar(-0.434562276666));
+    VERIFY_IS_APPROX(numext::polygamma(Scalar(147), Scalar(54.1)), RealScalar(0.567742190178));
+    VERIFY_IS_APPROX(numext::polygamma(Scalar(170), Scalar(64)), RealScalar(-0.0108615497927));
+    */
 
     {
       // Test various propreties of igamma & igammac.  These are normalized