From d44f24c47cc0d03135fabb7c764ab23e9d5d1152 Mon Sep 17 00:00:00 2001
From: akr <akr@b2dd03c8-39d4-4d8f-98ff-823fe69b080e>
Date: Sat, 9 Feb 2008 14:45:43 +0000
Subject: * missing/tgamma.c (tgamma): use lgamma_r if available.

git-svn-id: svn+ssh://ci.ruby-lang.org/ruby/trunk@15420 b2dd03c8-39d4-4d8f-98ff-823fe69b080e
---
 ChangeLog        |  4 ++++
 missing/tgamma.c | 31 ++++++++++++++++++++++++++++++-
 2 files changed, 34 insertions(+), 1 deletion(-)

diff --git a/ChangeLog b/ChangeLog
index 079640fff9..2a6aee983e 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,7 @@
+Sat Feb  9 23:44:29 2008  Tanaka Akira  <akr@fsij.org>
+
+	* missing/tgamma.c (tgamma): use lgamma_r if available.
+
 Sat Feb  9 23:22:52 2008  Nobuyoshi Nakada  <nobu@ruby-lang.org>
 
 	* ext/bigdecimal/extconf.rb: simplified the condition.
diff --git a/missing/tgamma.c b/missing/tgamma.c
index 4ace1bcc5f..91f83fcb82 100644
--- a/missing/tgamma.c
+++ b/missing/tgamma.c
@@ -10,6 +10,35 @@ reference - Haruhiko Okumura: C-gengo niyoru saishin algorithm jiten
     gamma.c -- Gamma function
 ***********************************************************/
 #include <math.h>
+#include <errno.h>
+#include "ruby/config.h"
+
+#ifdef HAVE_LGAMMA_R
+
+double tgamma(double x)
+{
+    int sign;
+    double d;
+    if (x == 0.0) { /* Pole Error */
+        errno = ERANGE;
+        return 1/x < 0 ? -HUGE_VAL : HUGE_VAL;
+    }
+    if (x < 0) {
+        int sign;
+        static double zero = 0.0;
+        double i, f;
+        f = modf(-x, &i);
+        if (f == 0.0) { /* Domain Error */
+            errno = EDOM;
+            return zero/zero;
+        }
+    }
+    d = lgamma_r(x, &sign);
+    return sign * exp(d);
+}
+
+#else
+
 #include <errno.h>
 #define PI      3.14159265358979324  /* $\pi$ */
 #define LOG_2PI 1.83787706640934548  /* $\log 2\pi$ */
@@ -61,4 +90,4 @@ double tgamma(double x)  /* Gamma function */
     }
     return exp(loggamma(x));
 }
-
+#endif
-- 
cgit v1.2.3