From 738278785652a2040d4d29fe10a3131403c6407d Mon Sep 17 00:00:00 2001
From: Daniel Richtmann <daniel.richtmann@physik.uni-regensburg.de>
Date: Wed, 8 Nov 2017 14:23:55 +0100
Subject: [PATCH] Some minor changes

---
 lib/algorithms/iterative/GeneralisedMinimalResidual.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/lib/algorithms/iterative/GeneralisedMinimalResidual.h b/lib/algorithms/iterative/GeneralisedMinimalResidual.h
index 9851987c..62a1699f 100644
--- a/lib/algorithms/iterative/GeneralisedMinimalResidual.h
+++ b/lib/algorithms/iterative/GeneralisedMinimalResidual.h
@@ -247,8 +247,8 @@ class GeneralisedMinimalResidual : public OperatorFunction<Field> {
     for (int i = iter; i >= 0; i--) {
       y[i] = gamma[i];
       for (int k = i + 1; k <= iter; k++)
-        y[i] -= H(k, i) * y[k];
-      y[i] /= H(i, i);
+        y[i] = y[i] - H(k, i) * y[k];
+      y[i] = y[i] / H(i, i);
     }
 
     // TODO: Use axpys or similar for these