Use sparse triangular solvers for sparse triangular solves. Fixes #13792

JuliaLang · Oct 27, 2015 · 6cce310 · 6cce310
1 parent 5287930
commit 6cce310
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/base/sparse/linalg.jl b/base/sparse/linalg.jl
@@ -235,6 +235,9 @@ function bwdTriSolve!(A::SparseMatrixCSC, B::AbstractVecOrMat)
    return B
 end
 
+A_ldiv_B!{T,Ti}(L::LowerTriangular{T,SparseMatrixCSC{T,Ti}}, B::StridedVecOrMat) = fwdTriSolve!(L.data, B)
+A_ldiv_B!{T,Ti}(L::UpperTriangular{T,SparseMatrixCSC{T,Ti}}, B::StridedVecOrMat) = bwdTriSolve!(L.data, B)
+
 ## triu, tril
 
 function triu{Tv,Ti}(S::SparseMatrixCSC{Tv,Ti}, k::Integer=0)