Only pass *sequential* (in var index) selected group members to ApplyBnd<BC>.

Fixes failing testsuite bh_shift_rad in BSSN_MoL.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@226 6a38eb6e-646e-4a02-a296-d141613ad6c4

1 files changed, 8 insertions, 4 deletions

diff --git a/src/RadiationBoundary.c b/src/RadiationBoundary.c
index 2033658..4ce0f34 100644
--- a/src/RadiationBoundary.c
+++ b/src/RadiationBoundary.c
@@ -193,10 +193,14 @@ int BndRadiative(const cGH *GH, int num_vars, int *vars, int *faces,
        can do is find variables of the same group which are selected
        for identical bcs.  If all GFs had the same staggering then we
        could groups many GFs together. */
-    gi = CCTK_GroupIndexFromVarI (vars[i]);    
-    while (i+j<num_vars && CCTK_GroupIndexFromVarI(vars[i+j])==gi && 
-	   tables[i+j]==tables[i] && faces[i+j]==faces[i] && 
-	   widths[i+j]==widths[i])
+    gi = CCTK_GroupIndexFromVarI (vars[i]);
+#ifdef DEBUG
+    printf("starting increment computation with group %d:\n", gi);
+    printf("this group had %d members\n", CCTK_NumVarsInGroupI(gi));
+#endif
+    while (i+j<num_vars && vars[i+j]==vars[i]+j && 
+	   CCTK_GroupIndexFromVarI(vars[i+j])==gi && tables[i+j]==tables[i] 
+	   && faces[i+j]==faces[i] && widths[i+j]==widths[i])
     {
       ++j;
     }