We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent aa106e2 commit d14d835Copy full SHA for d14d835
src/abstractblocksparsearray/abstractblocksparsematrix.jl
@@ -67,7 +67,7 @@ function svd!(
67
@assert is_block_permutation_matrix(A) "Cannot keep sparsity: use `svd` to convert to `BlockedMatrix"
68
U, S, Vt = _allocate_svd_output(A, full, alg)
69
for bI in eachblockstoredindex(A)
70
- bUSV = svd!(A[bI]; full, alg)
+ bUSV = svd!(@view!(A[bI]); full, alg)
71
brow, bcol = Tuple(bI)
72
U[brow, bcol] = bUSV.U
73
S[bcol] = bUSV.S
0 commit comments