Added preliminary support for mpfr_lngamma and tests.

TODO: add support for negative numbers with small exponent.
Re-enable the generic tests when this is done.


git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@4413 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 12 insertions, 0 deletions

diff --git a/tests/tlngamma.c b/tests/tlngamma.c
index 0c7f81351..30dc9d848 100644
--- a/tests/tlngamma.c
+++ b/tests/tlngamma.c
@@ -208,6 +208,18 @@ special (void)
       exit (1);
     }
 
+  mpfr_set_prec (x, 256);
+  mpfr_set_prec (y, 32);
+  mpfr_set_si_2exp (x, -1, 200, GMP_RNDN);
+  mpfr_sub_ui (x, x, 1, GMP_RNDN);
+  mpfr_div_2ui (x, x, 1, GMP_RNDN);
+  mpfr_lngamma (y, x, GMP_RNDN);
+  if (!mpfr_nan_p (y))
+    {
+      printf ("Error for lngamma(-2^199-0.5)\n");
+      exit (1);
+    }
+
   mpfr_clear (x);
   mpfr_clear (y);
 }