fix setup_scalar_gridfn_in_ghost_zone()

to properly respect  want_*  flags


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@132 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 11 insertions, 1 deletions

diff --git a/src/patch/ghost_zone.cc b/src/patch/ghost_zone.cc
index 30911b3..c0dbca4 100644
--- a/src/patch/ghost_zone.cc
+++ b/src/patch/ghost_zone.cc
@@ -131,6 +131,15 @@ void symmetry_ghost_zone::setup_scalar_gridfn_in_ghost_zone
 	{
 	for (int ipar = min_ipar(iperp) ; ipar <= max_ipar(iperp) ; ++ipar)
 	{
+	// are we supposed to do this point?
+	if (    (want_min_par_corner && ipar_is_in_min_ipar_corner(ipar))
+	     || (want_non_corner     && ipar_is_in_non_corner(ipar))
+	     || (want_max_par_corner && ipar_is_in_max_ipar_corner(ipar))    )
+	   then {
+		// we want to do this point ==> nop here
+		}
+	   else continue;				// *** LOOP CONTROL ***
+
 	// where is the symmetry-image of this point?
 	const int sym_iperp = symmetry_iperp_of_iperp(iperp);
 	const int sym_ipar  = symmetry_ipar_of_ipar  (ipar );
@@ -140,7 +149,8 @@ void symmetry_ghost_zone::setup_scalar_gridfn_in_ghost_zone
 			       .isigma_of_iperp_ipar(sym_iperp,sym_ipar);
 
 	// get the gridfn there
-	const fp sym_gridfn = symmetry_patch().gridfn(gfn, sym_irho,sym_isigma);
+	const fp sym_gridfn = symmetry_patch()
+			      .gridfn(gfn, sym_irho,sym_isigma);
 
 	// store the gridfn here
 	// (here's where we use the fact that it's a scalar!)