diff --git a/src/block_minres.jl b/src/block_minres.jl index f9f222186..004472003 100644 --- a/src/block_minres.jl +++ b/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