Update the type of the temporary buffers

JuliaSmoothOptimizers · Jan 5, 2025 · bcc7be2 · bcc7be2
1 parent 25af95b
commit bcc7be2
Showing 1 changed file with 8 additions and 8 deletions.
diff --git 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, :)