From 4c351aea173995f6e65bb85af8511d8ee78ba839 Mon Sep 17 00:00:00 2001 From: akr Date: Sat, 9 Feb 2008 02:10:57 +0000 Subject: * 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 --- missing/lgamma_r.c | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) (limited to 'missing/lgamma_r.c') 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 +#include #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); -- cgit v1.2.3