From faf42780192f65ff0cd15842f79579e9786ce4b5 Mon Sep 17 00:00:00 2001
From: Daniel Richtmann <daniel.richtmann@physik.uni-regensburg.de>
Date: Mon, 29 Jan 2018 17:21:42 +0100
Subject: [PATCH] Use 2 passes of GS in coarse operator construction

---
 lib/algorithms/CoarsenedMatrix.h | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/lib/algorithms/CoarsenedMatrix.h b/lib/algorithms/CoarsenedMatrix.h
index c76cdd14..59025276 100644
--- a/lib/algorithms/CoarsenedMatrix.h
+++ b/lib/algorithms/CoarsenedMatrix.h
@@ -334,10 +334,9 @@ namespace Grid {
 
       CoarseVector iProj(Grid()); 
       CoarseVector oProj(Grid()); 
-      CoarseScalar InnerProd(Grid()); 
 
       // Orthogonalise the subblocks over the basis
-      blockOrthogonalise(InnerProd,Subspace.subspace);
+      Subspace.Orthogonalise();
 
       // Compute the matrix elements of linop between this orthonormal
       // set of vectors.