From c528242a98091617cf14018b695c313c578c7a1a Mon Sep 17 00:00:00 2001
From: Anton Khirnov <anton@khirnov.net>
Date: Wed, 2 May 2018 14:04:48 +0200
Subject: Add an option to use the maximal lapse from libteukolskydata

---
 param.ccl       | 5 +++++
 src/teukolsky.c | 7 +++++++
 2 files changed, 12 insertions(+)

diff --git a/param.ccl b/param.ccl
index 7c3b340..cd81760 100644
--- a/param.ccl
+++ b/param.ccl
@@ -5,6 +5,11 @@ EXTENDS KEYWORD initial_data
   "teukolskydata" :: "Teukolsky wave initial data"
 }
 
+EXTENDS KEYWORD initial_lapse
+{
+  "teukolsky_max" :: "Teukolsky maximal slicing lapse"
+}
+
 RESTRICTED:
 CCTK_REAL amplitude "Wave amplitude A."
 {
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);
-- 
cgit v1.2.3