diff options
| author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2003-12-03 23:25:24 +0000 |
|---|---|---|
| committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2003-12-03 23:25:24 +0000 |
| commit | cada3378b043d7c8c0e7ee28ecc7087182d311ae (patch) | |
| tree | 9ba537a8768a87bbfc2e4460d3aaf582c628928c /src/driver/mask.cc | |
| parent | 75e72eb8d538dd35632d3228cb4cbd059e3942d2 (diff) | |
At Ian Hawke's request, make the setting of an excision mask specifyable
on a per-horizon basis:
The old parameter set_mask is now replaced by two new parameters:
set_mask_for_all_horizons
set_mask_for_individual_horizon[n]
For most purposes, you can just change set_mask to
set_mask_for_all_horizons in your parameter files. If you
really want to control the setting of the mask on a per-horizon basis,
read section 4.6 in the thorn guide to find out how the new parameters
work.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1222 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/driver/mask.cc')
| -rw-r--r-- | src/driver/mask.cc | 7 |
1 files changed, 5 insertions, 2 deletions
diff --git a/src/driver/mask.cc b/src/driver/mask.cc index 5a36f5f..baa08d4 100644 --- a/src/driver/mask.cc +++ b/src/driver/mask.cc @@ -168,7 +168,7 @@ fp outer_mask_radius(const struct mask_grid_info& mgi, // apparent horizons, to do all this thorn's mask processing. // extern "C" - void AHFinderDirect_do_masks(CCTK_ARGUMENTS) + void AHFinderDirect_maybe_do_masks(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS @@ -177,7 +177,7 @@ const struct verbose_info& verbose_info = state.verbose_info; struct mask_info& mask_info = state.mask_info; // optionally set the mask gridfn based on each horizon's shape -if (mask_info.set_mask) +if (mask_info.set_mask_for_any_horizon) then { // // this setup has to be done each time we're called, and @@ -412,6 +412,9 @@ set_mask_gridfn_to_outside_value(mgi, // for (int hn = 1 ; hn <= N_horizons ; ++hn) { + if (! mask_info.set_mask_for_this_horizon[hn]) + then continue; // *** LOOP CONTROL *** + const struct AH_data& AH_data = *AH_data_array[hn]; if (! AH_data.found_flag) then continue; // *** LOOP CONTROL *** |