diff options
| -rw-r--r-- | src/patch/ghost_zone.cc | 12 |
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!) |