Merge pull request #406 from JuliaArrays/qr_bigfloat

Fix BigFloat QR instance and test better

Author: ChrisRackauckas

src/ArrayInterface.jl

@@ -576,6 +576,10 @@ function qr_instance(A::Matrix{T}) where {T}
     LinearAlgebra.QRCompactWYQ(zeros(T,0,0),zeros(T,0,0))
 end
 
+function qr_instance(A::Matrix{BigFloat})
+    LinearAlgebra.QR(zeros(BigFloat,0,0),zeros(BigFloat,0))
+end
+
 # Could be optimized but this should work for any real case.
 function qr_instance(jac_prototype::SparseMatrixCSC)
     qr(sparse(rand(1,1)))

test/core.jl

@@ -259,3 +259,18 @@ end
     @test !ArrayInterface.ensures_sorted([])
     @test ArrayInterface.ensures_sorted(1:10)
 end
+
+@testset "linearalgebra instances" begin
+    for A in [rand(2,2), rand(Float32,2,2), rand(BigFloat,2,2)]
+        
+        @test ArrayInterface.lu_instance(A) isa typeof(lu(A))
+        @test ArrayInterface.qr_instance(A) isa typeof(qr(A))
+
+        if !(eltype(A) <: BigFloat)
+            @test ArrayInterface.bunchkaufman_instance(A) isa typeof(bunchkaufman(A' * A))
+            @test ArrayInterface.cholesky_instance(A) isa typeof(cholesky(A' * A))
+            @test ArrayInterface.ldlt_instance(A) isa typeof(ldlt(SymTridiagonal(A' * A)))
+            @test ArrayInterface.svd_instance(A) isa typeof(svd(A))
+        end
+    end
+end