Add an option to use the maximal lapse from libteukolskydata

1 files changed, 7 insertions, 0 deletions

diff --git a/src/teukolsky.c b/src/teukolsky.c
index 60cef99..234fedb 100644
--- a/src/teukolsky.c
+++ b/src/teukolsky.c
@@ -144,6 +144,13 @@ void teukolsky_data(CCTK_ARGUMENTS)
         }
     }
 
+    if (!strcmp(initial_lapse, "teukolsky_max")) {
+        ret = td_eval_lapse(td, cctk_lsh[2] * cctk_lsh[0], r_val, theta_val,
+                            (const unsigned int [2]){ 0, 0}, alp);
+        if (ret)
+            CCTK_WARN(0, "Error evaluating the lapse\n");
+    }
+
     free(r_val);
     free(theta_val);
     free(psi_val);