From bcc7be2d6d5702d3206f624fb45669358094d12c Mon Sep 17 00:00:00 2001 From: Alexis Montoison Date: Mon, 4 Nov 2024 01:39:29 -0600 Subject: [PATCH] Update the type of the temporary buffers --- src/block_minres.jl | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/src/block_minres.jl b/src/block_minres.jl index c9d8b5106..a35613b83 100644 --- a/src/block_minres.jl +++ b/src/block_minres.jl @@ -116,14 +116,14 @@ kwargs_block_minres = (:M, :ldiv, :atol, :rtol, :itmax, :timemax, :verbose, :his R₀ = warm_start ? Q : B # Temporary buffers -- should be stored in the solver - Ψₖ = zeros(p, p) - Ωₖ = zeros(p, p) - Ψₖ₊₁ = zeros(p, p) - Πₖ₋₂ = zeros(p, p) - Γbarₖ₋₁ = zeros(p, p) - Γₖ₋₁ = zeros(p, p) - Λbarₖ = zeros(p, p) - Λₖ = zeros(p, p) + Ψₖ = similar(B, p, p) + Ωₖ = similar(B, p, p) + Ψₖ₊₁ = similar(B, p, p) + Πₖ₋₂ = similar(B, p, p) + Γbarₖ₋₁ = similar(B, p, p) + Γₖ₋₁ = similar(B, p, p) + Λbarₖ = similar(B, p, p) + Λₖ = similar(B, p, p) # Define the blocks D1 and D2 D1 = view(D, 1:p, :)