Merge pull request #121 from JuliaMath/infnorm

Infnorm

Author: JeffreySarnoff

Project.toml

@@ -4,7 +4,7 @@ keywords = ["math", "floating-point", "doubledouble", "extended-precision", "acc
 license = "MIT"
 desc = "extended precision floating point types using value pairs"
 repo = "https://github.com/JuliaMath/DoubleFloats.jl.git"
-version = "1.1.15"
+version = "1.1.16"
 
 [deps]
 GenericSVD = "01680d73-4ee2-5a08-a1aa-533608c188bb"

src/math/linearalgebra/support.jl

@@ -8,7 +8,7 @@ end
     return r === c
 end
 
-function norm(v::Array{DoubleFloat{T},N}, p::Real=2.0) where {N, T<:IEEEFloat}
+function LinearAlgebra.norm(v::Array{DoubleFloat{T},N}, p::Real=2.0) where {N, T<:IEEEFloat}
     isempty(v) && return zero(DoubleFloat{T})
 
     if isinf(p)
@@ -22,7 +22,7 @@ function norm(v::Array{DoubleFloat{T},N}, p::Real=2.0) where {N, T<:IEEEFloat}
     end    
 end
 
-function norm(v::Array{Complex{DoubleFloat{T}},N}, p::Real=2.0) where {N, T<:IEEEFloat}
+function LinearAlgebra.norm(v::Array{Complex{DoubleFloat{T}},N}, p::Real=2.0) where {N, T<:IEEEFloat}
     isempty(v) && return zero(DoubleFloat{T})
 
     if isinf(p)

src/math/prearith/prearith.jl

@@ -100,17 +100,21 @@ ulp(::Type{D}) where {T<:AbstractFloat, D<:DoubleFloat{T}} = D(posulp(poslulp(on
 
 
 function Base.Math.nextfloat(x::DoubleFloat{T}) where {T<:IEEEFloat}
-    !isfinite(x) && return(x)
-    if !iszero(LO(x))
+    if !isfinite(x)
+       isnan(x) && return x
+       signbit(x) ? floatmin(DoubleFloat{T}) : x
+    elseif !iszero(LO(x))
         DoubleFloat{T}(HI(x)) + nextfloat(LO(x))
     else
         DoubleFloat{T}(HI(x), ulp(ulp(HI(x))))
     end
 end
 
 function Base.Math.prevfloat(x::DoubleFloat{T}) where {T<:IEEEFloat}
-    !isfinite(x) && return(x)
-    if !iszero(LO(x))
+    if !isfinite(x)
+        isnan(x) && return x
+        signbit(x) ? x : floatmax(DoubleFloat{T})
+    elseif !iszero(LO(x))
         DoubleFloat{T}(HI(x)) + prevfloat(LO(x))
     else
         DoubleFloat{T}(HI(x), -ulp(ulp(HI(x))))