#2286: Add the Fishbone-Moncrief initial data thorn to the ET

Reporter:	Steven R. Brandt
Status:	new
Milestone:	ET_2019_10
Version:
Type:	enhancement
Priority:	major
Component:	EinsteinToolkit thorn

Comment (by Roland Haas):

What you need is something like https://bitbucket.org/einsteintoolkit/einsteininitialdata/src/f2e6accbb09bec42698149cadd069ef470838a0f/TOVSolver/src/tov.c#lines-513 and https://bitbucket.org/einsteintoolkit/einsteininitialdata/src/f2e6accbb09bec42698149cadd069ef470838a0f/TOVSolver/src/tov.c#lines-710 though eg https://bitbucket.org/einsteintoolkit/einsteininitialdata/src/f2e6accbb09bec42698149cadd069ef470838a0f/TO VSolver/src/tov.c#lines-557 is wrong b/c it does not honor the ID flags.

Something like:

do_lapse = CCTK_Equals(initial_lapse, "FishboneMoncriefID");
do_metric = CCTK_Equals(initial_data, "FishboneMoncriefID");

for (...) {
  if(do_lapse) alp[idx] = ...
  if(do_metric) {
    gxx[idx] = ...
  }
}

and if you want to be extra good, you can add something like:

if(!(CCTK_Equals(initial_data_setup_method, "init_sime_levels") ||
  CCTK_Equals(initial_data_setup_method, "init_single_level")) {
  CCTK_VPARAMWARN(
    "initial_data_setup_method muse be either 'foo' or 'bar', but got '%s',
    initial_data_setup_method);
}

to check that the InitiaBase::initial_data_setup_method is set to a setting you support. Admittedly few (none other than Carpet which is the code that actually needs to know) codes do so.

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2286/add-the-fishbone-moncrief-initial-data