diff --git a/src/Constants.jl b/src/Constants.jl index 9a91fcd..a412d5d 100644 --- a/src/Constants.jl +++ b/src/Constants.jl @@ -5,4 +5,3 @@ const ZETA3_F64 = 1.2020569031595943 # zeta(3) with Float64 precision const ZETA4_F64 = 1.0823232337111382 # zeta(4) with Float64 precision const ZETA5_F64 = 1.0369277551433699 # zeta(5) with Float64 precision const ZETA6_F64 = 1.0173430619844491 # zeta(6) with Float64 precision -const EULERGAMMA_F64 = 0.57721566490153286 # Euler gamma constant with Float64 precision diff --git a/src/Harmonic.jl b/src/Harmonic.jl index af541e6..fd6cf72 100644 --- a/src/Harmonic.jl +++ b/src/Harmonic.jl @@ -5,7 +5,7 @@ function harmonic(n::Integer, ::Type{T})::T where T elseif n == 1 one(T) elseif issimplefloat(T) && n >= 20 - convert(T, EULERGAMMA_F64 + digamma(n + 1)) + convert(T, Base.MathConstants.eulergamma + digamma(n + 1)) else one(T) + sum(x -> inv(T(x)), 2:n) end