Merge pull request #462 from JuliaArrays/ChrisRackauckas-patch-2

Fix SparseMatrixCSC with alternative Integer types

Author: ChrisRackauckas

ext/ArrayInterfaceSparseArraysExt.jl

@@ -14,25 +14,25 @@ function findstructralnz(x::SparseMatrixCSC)
     (rowind, colind)
 end
 
-function bunchkaufman_instance(A::SparseMatrixCSC)
-    bunchkaufman(sparse(similar(A, 1, 1)), check = false)
+function bunchkaufman_instance(A::SparseMatrixCSC{Tv, Ti}) where {Tv, Ti}
+    bunchkaufman(SparseMatrixCSC{Tv, Ti}(similar(A, 1, 1)), check = false)
 end
 
-function cholesky_instance(A::Union{SparseMatrixCSC,Symmetric{<:Number,<:SparseMatrixCSC}}, pivot = DEFAULT_CHOLESKY_PIVOT)
-    cholesky(sparse(similar(A, 1, 1)), check = false)
+function cholesky_instance(A::Union{SparseMatrixCSC{Tv, Ti},Symmetric{<:Number,<:SparseMatrixCSC{Tv, Ti}}}, pivot = DEFAULT_CHOLESKY_PIVOT) where {Tv, Ti}
+    cholesky(SparseMatrixCSC{Tv, Ti}(similar(A, 1, 1)), check = false)
 end
 
-function ldlt_instance(A::SparseMatrixCSC)
-    ldlt(sparse(similar(A, 1, 1)), check=false)
+function ldlt_instance(A::SparseMatrixCSC{Tv, Ti}) where {Tv, Ti}
+    ldlt(SparseMatrixCSC{Tv, Ti}(similar(A, 1, 1)), check=false)
 end
 
 # Could be optimized but this should work for any real case.
-function lu_instance(jac_prototype::SparseMatrixCSC, pivot = DEFAULT_CHOLESKY_PIVOT)
-    lu(sparse(rand(1,1)))
+function lu_instance(jac_prototype::SparseMatrixCSC{Tv, Ti}, pivot = DEFAULT_CHOLESKY_PIVOT) where {Tv, Ti}
+    lu(SparseMatrixCSC{Tv, Ti}(rand(1,1)))
 end
 
-function qr_instance(jac_prototype::SparseMatrixCSC, pivot = DEFAULT_CHOLESKY_PIVOT)
-    qr(sparse(rand(1,1)))
+function qr_instance(jac_prototype::SparseMatrixCSC{Tv, Ti}, pivot = DEFAULT_CHOLESKY_PIVOT) where {Tv, Ti}
+    qr(SparseMatrixCSC{Tv, Ti}(rand(1,1)))
 end
 
 end