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.