From 4c351aea173995f6e65bb85af8511d8ee78ba839 Mon Sep 17 00:00:00 2001
From: akr <akr@b2dd03c8-39d4-4d8f-98ff-823fe69b080e>
Date: Sat, 9 Feb 2008 02:10:57 +0000
Subject: [PATCH] * missing/xlgamma_r.c (lgamma_r): return HUGE_VAL for
 non-positive   integers.

git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@15411 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
---
 ChangeLog          | 5 +++++
 missing/lgamma_r.c | 9 +++++----
 2 files changed, 10 insertions(+), 4 deletions(-)

diff --git a/ChangeLog b/ChangeLog
index 0274117f0a..df2aa45fa3 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,8 @@
+Sat Feb  9 11:09:26 2008  Tanaka Akira  <akr@fsij.org>
+
+	* missing/xlgamma_r.c (lgamma_r): return HUGE_VAL for non-positive
+	  integers.
+
 Sat Feb  9 10:03:07 2008  Tanaka Akira  <akr@fsij.org>
 
 	* string.c (rb_str_new4): copy encoding from orig, instead of shared
diff --git a/missing/lgamma_r.c b/missing/lgamma_r.c
index d4402c67b1..ab6c28294a 100644
--- a/missing/lgamma_r.c
+++ b/missing/lgamma_r.c
@@ -12,6 +12,7 @@ reference - Haruhiko Okumura: C-gengo niyoru saishin algorithm jiten
     gamma.c -- Gamma function
 ***********************************************************/
 #include <math.h>
+#include <errno.h>
 #define PI      3.14159265358979324  /* $\pi$ */
 #define LOG_2PI 1.83787706640934548  /* $\log 2\pi$ */
 #define LOG_PI  1.14472988584940017  /* $\log_e \pi$ */
@@ -47,13 +48,13 @@ loggamma(double x)  /* the natural logarithm of the Gamma function. */
 double
 lgamma_r(double x, int *signp)
 {
-    if (x < 0) {
+    if (x <= 0) {
         double i, f, s;
         f = modf(-x, &i);
-        if (f == 0.0) {
-            static const double zero = 0.0;
+        if (f == 0.0) { /* pole error */
             *signp = 1;
-            return 1.0/zero;
+            errno = ERANGE;
+            return HUGE_VAL;
         }
         *signp = (fmod(i, 2.0) != 0.0) ? 1 : -1;
         s = sin(PI * f);