From 3575278b57259da65b45b65fef2167b98ce2784c Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Wed, 12 May 2021 10:57:04 -0400
Subject: [PATCH] Small change

---
 Grid/qcd/action/fermion/SchurFactoredFermionOperator.h | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/Grid/qcd/action/fermion/SchurFactoredFermionOperator.h b/Grid/qcd/action/fermion/SchurFactoredFermionOperator.h
index ffc345c5..d04ad1f2 100644
--- a/Grid/qcd/action/fermion/SchurFactoredFermionOperator.h
+++ b/Grid/qcd/action/fermion/SchurFactoredFermionOperator.h
@@ -66,12 +66,14 @@ template<class Impl>
 class SchurFactoredFermionOperator : public Impl
 {
   INHERIT_IMPL_TYPES(Impl);
+
+public:
   
   FermionOperator<Impl> & DirichletFermOp;
   FermionOperator<Impl> & FermOp; 
   OperatorFunction<FermionField> &Solver;
   Coordinate Block;
-public:
+
   SchurFactoredFermionOperator(FermionOperator<Impl> & _FermOp,
 			       FermionOperator<Impl> & _DirichletFermOp,
 			       OperatorFunction<FermionField> &_Solver,
@@ -84,6 +86,8 @@ public:
     // FIXME -- could check that the DirichletFermOp block matches this.
     // Pass in Dirichlet FermOp because we really need two dirac operators
     // as double stored gauge fields differ.
+
+    assert(_FermOp.FermionGrid() = _DirichletFermOp.FermionGrid()); // May not be true in future if change communicator scheme
   };
 
   enum Domain { Omega=0, OmegaBar=1 };