From c528242a98091617cf14018b695c313c578c7a1a Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Wed, 2 May 2018 14:04:48 +0200 Subject: Add an option to use the maximal lapse from libteukolskydata --- src/teukolsky.c | 7 +++++++ 1 file changed, 7 insertions(+) (limited to 'src') 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