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