[Scipy-svn] r4146 - trunk/scipy/special/tests

[scipy-svn@scip...]

Author: cookedm
Date: 2008-04-17 16:11:39 -0500 (Thu, 17 Apr 2008)
New Revision: 4146

Modified:
   trunk/scipy/special/tests/test_basic.py
Log:
Add test case for #299: gammaincinv gives incorrect answers


Modified: trunk/scipy/special/tests/test_basic.py
===================================================================
--- trunk/scipy/special/tests/test_basic.py	2008-04-17 21:05:42 UTC (rev 4145)
+++ trunk/scipy/special/tests/test_basic.py	2008-04-17 21:11:39 UTC (rev 4146)
@@ -1051,6 +1051,12 @@
         y = gammaincinv(.4,.4)
         x = gammainc(.4,y)
         assert_almost_equal(x,0.4,1)
+        y = gammainc(10, 0.05)
+        x = gammaincinv(10, 2.5715803516000736e-20)
+        assert_almost_equal(0.05, x, decimal=10)
+        assert_almost_equal(y, 2.5715803516000736e-20, decimal=10)
+        x = gammaincinv(50, 8.20754777388471303050299243573393e-18)
+        assert_almost_equal(11.0, x, decimal=10)
 
     def test_rgamma(self):
         rgam = rgamma(8)