Merge branch 'factorizations' of https://github.com/lkdvos/BlockSparseArrays.jl into factorizations

Author: lkdvos

src/abstractblocksparsearray/abstractblocksparsematrix.jl

@@ -49,9 +49,10 @@ function _allocate_svd_output(A::AbstractBlockSparseMatrix, full::Bool, ::Algori
     slengths[col] = min(brows[row], bcols[col])
   end
 
-  U = similar(A, axes(A, 1), blockedrange(slengths))
-  S = similar(A, real(eltype(A)), blockedrange(slengths))
-  Vt = similar(A, blockedrange(slengths), axes(A, 2))
+  s_axis = blockedrange(slengths)
+  U = similar(A, axes(A, 1), s_axis)
+  S = similar(A, real(eltype(A)), s_axis)
+  Vt = similar(A, s_axis, axes(A, 2))
 
   return U, S, Vt
 end
@@ -67,10 +68,10 @@ function svd!(
   U, S, Vt = _allocate_svd_output(A, full, alg)
   for bI in eachblockstoredindex(A)
     bUSV = svd!(A[bI]; full, alg)
-    brow, bcol = Int.(Tuple(bI))
-    U[Block(brow, bcol)] = bUSV.U
-    S[Block(bcol)] = bUSV.S
-    Vt[Block(bcol, bcol)] = bUSV.Vt
+    brow, bcol = Tuple(bI)
+    U[brow, bcol] = bUSV.U
+    S[bcol] = bUSV.S
+    Vt[bcol, bcol] = bUSV.Vt
   end
   return SVD(U, S, Vt)
 end