Skip to content

Commit e70acb3

Browse files
stevengjdkarrasch
andauthored
Update src/svd.jl
Co-authored-by: Daniel Karrasch <[email protected]>
1 parent 6a5ed42 commit e70acb3

File tree

1 file changed

+3
-1
lines changed

1 file changed

+3
-1
lines changed

src/svd.jl

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -299,7 +299,9 @@ function ldiv!(F::SVD{T}, B::AbstractVecOrMat; atol::Real=0, rtol::Real = (eps(r
299299
if k == 0
300300
B[1:n, :] .= 0
301301
else
302-
mul!(view(B, 1:n, :), view(F.Vt, 1:k, :)', view(F.S, 1:k) .\ (view(F.U, :, 1:k)' * _cut_B(B, 1:m)))
302+
temp = view(F.U, :, 1:k)' * _cut_B(B, 1:m)
303+
ldiv!(Diagonal(view(F.S, 1:k)), temp)
304+
mul!(view(B, 1:n, :), view(F.Vt, 1:k, :)', temp)
303305
end
304306
return B
305307
end

0 commit comments

Comments
 (0)