From 217d4a2f5dcb5f48141f49aa6e874e3c704c7e31 Mon Sep 17 00:00:00 2001 From: Eric Teunis de Boone Date: Thu, 27 Feb 2020 15:33:08 +0100 Subject: [PATCH] Week3: fixed error in the Conjugate Gradient Method --- week3/solvers.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/week3/solvers.py b/week3/solvers.py index fb336c3..1f2f0d8 100644 --- a/week3/solvers.py +++ b/week3/solvers.py @@ -81,7 +81,7 @@ def conjugate_gradient(A, b, eps): r_i = r_f - t * np.dot(A, v_f) s = np.dot(r_i, r_i) / np.dot(r_f, r_f) - v_i = r_i - s*v_f + v_i = r_i + s*v_f # Set r and v vectors for next loop r_f = r_i @@ -101,7 +101,7 @@ def conjugate_gradient(A, b, eps): r_i = r_f - t * np.dot(A, v_f) s = np.dot(r_i, r_i) / np.dot(r_f, r_f) - v_i = r_i - s*v_f + v_i = r_i + s*v_f # Save r and v vectors for next loop r_f = r_i