From 8f96c192410a07cba14436c1471b856fb37c3652 Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Sun, 4 Feb 2007 21:21:00 +0000 Subject: Carpet: Call postregrid bin also when only levels have been removed Correct an error in determining whether the grid hierarchy has changed: When only levels are removed, this is now also counted as change. Call postregrid on the newly finest grid in this case. darcs-hash:20070204212111-dae7b-3447e8d8959c437616ed0e0431e13b373aa93f77.gz --- Carpet/Carpet/src/Evolve.cc | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) (limited to 'Carpet/Carpet/src/Evolve.cc') diff --git a/Carpet/Carpet/src/Evolve.cc b/Carpet/Carpet/src/Evolve.cc index e810cc334..46b3ff6fc 100644 --- a/Carpet/Carpet/src/Evolve.cc +++ b/Carpet/Carpet/src/Evolve.cc @@ -224,7 +224,10 @@ namespace Carpet { // Regrid Checkpoint ("Regrid"); + int const oldreflevels = reflevels; bool const did_regrid = Regrid (cctkGH, false); + bool const did_remove_level = reflevels < oldreflevels; + assert (not did_remove_level or did_regrid); if (did_regrid) { BEGIN_META_MODE (cctkGH) { @@ -232,7 +235,7 @@ namespace Carpet { bool const did_recompose = Recompose (cctkGH, rl, true); - if (did_recompose) { + if (did_recompose or (did_remove_level and rl == reflevels - 1)) { BEGIN_MGLEVEL_LOOP (cctkGH) { ENTER_LEVEL_MODE (cctkGH, rl) { do_global_mode = reflevel == reflevels - 1; @@ -285,6 +288,7 @@ namespace Carpet { } LEAVE_LEVEL_MODE; } END_MGLEVEL_LOOP; } // if did_recompose + } // for rl } END_META_MODE; } // if did_regrid -- cgit v1.2.3