From 7d25d9d717ae98685f27a5da4b9ab533fca28213 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Mateus=20Ara=C3=BAjo?= Date: Thu, 4 Apr 2024 08:00:50 +0200 Subject: [PATCH] Update src/specfun.jl Co-authored-by: RalphAS --- src/specfun.jl | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/src/specfun.jl b/src/specfun.jl index 7592bdb..c7b1bf0 100644 --- a/src/specfun.jl +++ b/src/specfun.jl @@ -28,6 +28,10 @@ lgamma(x::Float128) = Float128(@ccall(libquadmath.lgammaq(x::Cfloat128)::Cfloat128)) function logabsgamma(x::Float128) - sign = x >= 0 ? 1 : 2*mod(ceil(Int64,x),2)-1 + if isfinite(x) + sign = x >= 0 ? 1 : 2*mod(ceil(Int64,x),2)-1 + else + sign = 1 + end return lgamma(x), sign end