Create new sparse matrix for R factor in sparse least square solve (#…

…30289) instead of using a view to avoid slow fallback in back substitution. Fixes #30288 (cherry picked from commit 87c5f36)
JuliaLang · Feb 20, 2020 · c8887c4 · c8887c4
1 parent 33f02ce
commit c8887c4
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/stdlib/SuiteSparse/src/spqr.jl b/stdlib/SuiteSparse/src/spqr.jl
@@ -391,7 +391,7 @@ function _ldiv_basic(F::QRSparse, B::StridedVecOrMat)
     X[rnk + 1:end, :] .= 0
 
     # Solve R*X = B
-    LinearAlgebra.ldiv!(UpperTriangular(view(F.R, Base.OneTo(rnk), Base.OneTo(rnk))),
+    LinearAlgebra.ldiv!(UpperTriangular(F.R[Base.OneTo(rnk), Base.OneTo(rnk)]),
                         view(X0, Base.OneTo(rnk), :))
 
     # Apply right permutation and extract solution from X