use copyto! instead of views to construct blocks

The-ReSolver · May 10, 2023 · 7c0f559 · 7c0f559
1 parent 246b818
commit 7c0f559
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/utils.jl b/src/utils.jl
@@ -15,8 +15,7 @@ function split_into_blocks(Q::M, Nf::Int, No::Int) where {M <: AbstractMatrix}
     # loop over the blocks and take views of snapshot matrix
     for nb in 1:Nb
         offset = (nb - 1)*(Nf - No)
-        # TODO: change to copyto! so that the original matrix isn't affected by windowing
-        Q_blocks[nb] = @view(Q[:, (1 + offset):(Nf + offset)])
+        copyto!(Q_blocks[nb], @view(Q[:, (1 + offset):(Nf + offset)]))
     end
 
     return Q_blocks, Nb