Rewrite using a new method.

Based on 'Adaptive mesh refinement for coupled elliptic-hyperbolic
problems' by Pretorius and Choptuik, 2006.

1 files changed, 25 insertions, 2 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 246aa63..0b992b9 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -1,10 +1,26 @@
 # Schedule definitions for thorn MaximalSlicingAxiMG
 #
 if (CCTK_Equals(lapse_evolution_method, "maximal_axi_mg")) {
+    SCHEDULE msa_mg_inithist in CCTK_INITIAL {
+        LANG: C
+    } ""
 
-    SCHEDULE maximal_slicing_axi_mg IN ML_BSSN_evolCalcGroup BEFORE ML_BSSN_RHS {
+    SCHEDULE msa_mg_init IN CCTK_BASEGRID {
         LANG: C
-    } "Maximal slicing in axisymmetry"
+    } ""
+
+    SCHEDULE msa_mg_eval IN MoL_CalcRHS BEFORE ML_BSSN_evolCalcGroup {
+        LANG: C
+    } ""
+
+    SCHEDULE msa_mg_sync IN CCTK_POSTSTEP BEFORE msa_mg_solve {
+        SYNC: lapse_mg
+        LANG: C
+    } ""
+
+    SCHEDULE msa_mg_solve IN CCTK_POSTSTEP {
+        LANG: C
+    } ""
 
     #SCHEDULE maximal_slicing_axi IN MoL_PostStepModify {
     #    LANG: C
@@ -21,4 +37,11 @@ if (CCTK_Equals(lapse_evolution_method, "maximal_axi_mg")) {
     SCHEDULE maximal_slicing_axi_mg_register_mol IN MoL_Register {
         LANG: C
     } ""
+
+    STORAGE: lapse_mg[2]
+    STORAGE: lapse_mg_eval
+    STORAGE: lapse_prev0
+    STORAGE: lapse_prev1
+    STORAGE: lapse_prev0_time
+    STORAGE: lapse_prev1_time
 }