From 3e448435d397f78742153920fe9df0465b5c51db Mon Sep 17 00:00:00 2001
From: Peter Boyle <paboyle@ph.ed.ac.uk>
Date: Thu, 23 Nov 2023 18:23:29 -0500
Subject: [PATCH] Restrict to interior

---
 .../multigrid/GeneralCoarsenedMatrixMultiRHS.h     | 14 ++++++++------
 1 file changed, 8 insertions(+), 6 deletions(-)

diff --git a/Grid/algorithms/multigrid/GeneralCoarsenedMatrixMultiRHS.h b/Grid/algorithms/multigrid/GeneralCoarsenedMatrixMultiRHS.h
index fed84a33..fcd79d5f 100644
--- a/Grid/algorithms/multigrid/GeneralCoarsenedMatrixMultiRHS.h
+++ b/Grid/algorithms/multigrid/GeneralCoarsenedMatrixMultiRHS.h
@@ -149,13 +149,15 @@ public:
 	  int32_t point= bp/nbasis;
 	  int32_t b    = bp%nbasis;
 	  auto SE  = Stencil_v.GetEntry(point,ss);
-	  int32_t snbr= SE->_offset;
-	  auto nbr = coalescedReadGeneralPermute(in_v[snbr],SE->_permute,Nd);
-	  auto res = Aview_p[point][ss](0,b)*nbr(0);
-	  for(int bb=1;bb<nbasis;bb++) {
-	    res = res + Aview_p[point][ss](bb,b)*nbr(bb);
+	  if ( SE->_permute == 0 ) { 
+	    int32_t snbr= SE->_offset;
+	    auto nbr = coalescedReadGeneralPermute(in_v[snbr],SE->_permute,Nd);
+	    auto res = Aview_p[point][ss](0,b)*nbr(0);
+	    for(int bb=1;bb<nbasis;bb++) {
+	      res = res + Aview_p[point][ss](bb,b)*nbr(bb);
+	    }
+	    coalescedWrite(Vview_p[point][ss](b),res);
 	  }
-	  coalescedWrite(Vview_p[point][ss](b),res);
       });
       accelerator_for(sb, osites*nbasis, Nsimd, {
 	  int ss = sb/nbasis;