generalize to sparse rhs

Author: dkarrasch

src/solvers/cholmod.jl

@@ -15,7 +15,7 @@ import Base: (*), convert, copy, eltype, getindex, getproperty, show, size,
 using Base: require_one_based_indexing
 
 using LinearAlgebra
-using LinearAlgebra: RealHermSymComplexHerm, AdjOrTrans, AdjOrTransStridedMat
+using LinearAlgebra: RealHermSymComplexHerm, AdjOrTrans, AdjOrTransAbsMat
 import LinearAlgebra: (\), AdjointFactorization,
                  cholesky, cholesky!, det, diag, ishermitian, isposdef,
                  issuccess, issymmetric, ldiv!, ldlt, ldlt!, logdet,
@@ -1928,7 +1928,7 @@ function \(adjL::AdjointFactorization{<:VTypes,<:Factor}, B::StridedMatrix)
     L = adjL.parent
     return Matrix(solve(CHOLMOD_A, L, Dense(B)))
 end
-(\)(adjL::AdjointFactorization{<:VTypes,<:Factor}, B::LinearAlgebra.AdjOrTransStridedMat) = adjL \ copy(B)
+(\)(adjL::AdjointFactorization{<:VTypes,<:Factor}, B::AdjOrTransAbsMat) = adjL \ copy(B)
 
 const RealHermSymComplexHermSSL{Ti, Tr} = Union{
     Symmetric{Tr, SparseMatrixCSC{Tr, Ti}},

test/cholmod.jl

@@ -1026,9 +1026,12 @@ end
 
     B = rand(n, 2)
     Bt = Matrix(B')
+    Bts = sparse(B')
 
     F = factorize(A')
     @test F \ B ≈ F \ Bt'
+    @test F \ B ≈ F \ Bts'
+    @test issparse(F \ Bts')
 end
 
 end # for Tv ∈ (Float32, Float64)