From 0a677239d72ee3367f03742f47a2a49057f19103 Mon Sep 17 00:00:00 2001 From: couriersud Date: Sun, 13 Jan 2019 02:16:08 +0100 Subject: [PATCH] Fix bug in GMRES solver. (nw) --- src/lib/netlist/solver/nld_ms_gmres.h | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/lib/netlist/solver/nld_ms_gmres.h b/src/lib/netlist/solver/nld_ms_gmres.h index 92d8a0fd6d9..2a73d45d3d1 100644 --- a/src/lib/netlist/solver/nld_ms_gmres.h +++ b/src/lib/netlist/solver/nld_ms_gmres.h @@ -288,8 +288,8 @@ unsigned matrix_solver_GMRES_t::solve_ilu_gmres (nl_double (& RE vec_set(mr+1, NL_FCONST(0.0), m_g); m_g[0] = rho; - for (std::size_t i = 0; i < mr; i++) - vec_set(mr + 1, NL_FCONST(0.0), m_ht[i]); + for (std::size_t i = 0; i < mr + 1; i++) + vec_set(mr, NL_FCONST(0.0), m_ht[i]); vec_mult_scalar(n, residual, NL_FCONST(1.0) / rho, m_v[0]);