Merge remote-tracking branch 'origin/master'

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.2.1"
+version = "1.2.2"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"

src/extras/misc.jl

@@ -62,3 +62,10 @@ end
 function Base.:(-)(m::Matrix{DoubleFloat{T}}, x::T) where {T}
     return m .- x
 end
+
+# for getting fast floatmin2, used in Givens rotations in LinearAlgebra
+# these values were computed with the existing code saved
+# floatmin2(::Type{T}) where {T} = (twopar = 2one(T); twopar^trunc(Integer,log(floatmin(T)/eps(T))/log(twopar)/twopar))
+LinearAlgebra.floatmin2(::Type{Double64}) = Double64(reinterpret(Float64, 0x2350000000000000), 0.0)
+LinearAlgebra.floatmin2(::Type{Double32}) = Double32(reinterpret(Float32, 0x2c000000), 0.0f0)
+LinearAlgebra.floatmin2(::Type{Double16}) = Double16(Float16(8.0), Float16(0.0))

test/specialvalues.jl

@@ -61,3 +61,8 @@ end
     @test isnan(asech(T(NaN)))
     @test isnan(acoth(T(NaN)))
 end
+    
+@testset "floatmin2 $T" for T in (Double16, Double32, Double64)
+    trueval = (twopar = 2one(T); twopar^trunc(Integer,log(floatmin(T)/eps(T))/log(twopar)/twopar))
+    @test LinearAlgebra.floatmin2(T) == trueval 
+end