From 7c2737faa45cc68ea9059be525a14d2613573909 Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Sun, 23 Jun 2019 15:02:43 +0200
Subject: Add an option to extrapolate on the finest level as well.

---
 param.ccl                    | 6 ++++++
 src/maximal_slicing_axi_mg.c | 2 +-
 2 files changed, 7 insertions(+), 1 deletion(-)

diff --git a/param.ccl b/param.ccl
index 55b52a6..03bf6b6 100644
--- a/param.ccl
+++ b/param.ccl
@@ -7,6 +7,12 @@ EXTENDS KEYWORD lapse_evolution_method
   "maximal_axi_mg" :: "Maximal slicing for an axisymmetric spacetime"
 }
 
+RESTRICTED:
+CCTK_INT fine_solve "Whether to solve on each sub-step of the finest level evolution"
+{
+    0:1 :: ""
+} 1
+
 RESTRICTED:
 CCTK_INT fd_stencil "finite differencing stencil"
 {
diff --git a/src/maximal_slicing_axi_mg.c b/src/maximal_slicing_axi_mg.c
index 11770d5..9d32120 100644
--- a/src/maximal_slicing_axi_mg.c
+++ b/src/maximal_slicing_axi_mg.c
@@ -960,7 +960,7 @@ void msa_mg_eval(CCTK_ARGUMENTS)
         fact1 = (t - lapse_prev0_time[reflevel]) / time_interp_step;
     }
 
-    if (reflevel < ms->nb_levels - 1) {
+    if (!fine_solve || reflevel < ms->nb_levels - 1) {
         /* on coarse levels use extrapolated past solutions */
         int64_t extrap_start = gettime();
 
-- 
cgit v1.2.3