From e98fb284fbb880e8a4d59856057f0cbdfa6b721c Mon Sep 17 00:00:00 2001 From: Sascha Timme Date: Fri, 8 May 2020 14:47:32 +0200 Subject: [PATCH] Fix Polynomials deprecation and allow v1 --- Project.toml | 2 +- src/math/elementary/trig.jl | 19 +++++++++---------- 2 files changed, 10 insertions(+), 11 deletions(-) diff --git a/Project.toml b/Project.toml index debb8b38..2880cf15 100644 --- a/Project.toml +++ b/Project.toml @@ -20,7 +20,7 @@ SpecialFunctions = "276daf66-3868-5448-9aa4-cd146d93841b" [compat] GenericSVD = "0" GenericSchur = "0" -Polynomials = "0" +Polynomials = "0.8, 1" Quadmath = "0" Requires = "0, 1.0" SpecialFunctions = "0" diff --git a/src/math/elementary/trig.jl b/src/math/elementary/trig.jl index fc957c42..29c996a4 100644 --- a/src/math/elementary/trig.jl +++ b/src/math/elementary/trig.jl @@ -233,7 +233,7 @@ function tan(x::Double64) isnan(x) && return x isinf(x) && throw(DomainError("tan(x) only defined for finite x")) abs(HI(x)) >= 0.36815538909255385 && return Double64(tan(Float128(x))) # (15/128 * pi) - + abs(mod1pi(x-Double64(pi)/2)) <= eps(one(DoubleFloat{Float64})) && return DoubleFloat{Float64}(Inf) iszero(x) && return zero(typeof(x)) !isfinite(x) && return nan(typeof(x)) @@ -287,8 +287,8 @@ const tan0qrtrpi_denomcoeffs = [ Double64(1.5859442637424446e-9, 6.30689666197131e-26) ]; -const tan0qrtrpi_numerpoly = Poly(tan0qrtrpi_numercoeffs); -const tan0qrtrpi_denompoly = Poly(tan0qrtrpi_denomcoeffs); +const tan0qrtrpi_numerpoly = Polynomial(tan0qrtrpi_numercoeffs); +const tan0qrtrpi_denompoly = Polynomial(tan0qrtrpi_denomcoeffs); function tan0qrtrpi(x::Double64) numer = polyval(tan0qrtrpi_numerpoly, x) @@ -298,7 +298,7 @@ end function csc(x::DoubleFloat{T}) where {T<:IEEEFloat} - isnan(x) && return x + isnan(x) && return x isinf(x) && throw(DomainError("csc(x) only defined for finite x")) return inv(sin(x)) end @@ -320,8 +320,8 @@ function sinpi(x::DoubleFloat{T}) where {T<:IEEEFloat} isnan(x) && return x isinf(x) && throw(DomainError("sinpi(x) only defined for finite x")) return DoubleFloat{T}(sinpi(Quadmath.Float128(x))) - #= - y = Double64(x) + #= + y = Double64(x) hi,lo = mul322(pi_1o1_t64, HILO(y)) y = Double64(hi, lo) z = sin(y) @@ -333,8 +333,8 @@ function cospi(x::DoubleFloat{T}) where {T<:IEEEFloat} isnan(x) && return x isinf(x) && throw(DomainError("cospi(x) only defined for finite x")) return DoubleFloat{T}(cospi(Quadmath.Float128(x))) - #= - y = Double64(x) + #= + y = Double64(x) hi,lo = mul322(pi_1o1_t64, HILO(y)) y = Double64(hi, lo) z = cos(y) @@ -347,11 +347,10 @@ function tanpi(x::DoubleFloat{T}) where {T<:IEEEFloat} isinf(x) && throw(DomainError("tanpi(x) only defined for finite x")) return sinpi(x)/cospi(x) #= - y = Double64(x) + y = Double64(x) hi,lo = mul322(pi_1o1_t64, HILO(y)) y = Double64(hi, lo) z = tan(y) return DoubleFloat{T}(z) =# end -