bob

src/block_minres.jl

@@ -180,7 +180,7 @@ kwargs_block_minres = (:M, :ldiv, :atol, :rtol, :itmax, :timemax, :verbose, :his
       mul!(W, A, V)  # Q ← AVₖ
       for i = 1 : inner_iter
         mul!(Ω, V', W)       # Ωₖ = Vₖᴴ * Q
-        (iter ≥ 2) && mul!(Q, ...)  # Q ← Q - βₖ * Vₖ₋₁ * Ψₖᴴ
+        # (iter ≥ 2) && mul!(Q, ...)  # Q ← Q - βₖ * Vₖ₋₁ * Ψₖᴴ
         mul!(Q, V, R, α, β)  # Q = Q - Vₖ * Ωₖ
       end