diff options
Diffstat (limited to 'ML_BSSN/src/RegisterMoL.cc')
| -rw-r--r-- | ML_BSSN/src/RegisterMoL.cc | 20 |
1 files changed, 16 insertions, 4 deletions
diff --git a/ML_BSSN/src/RegisterMoL.cc b/ML_BSSN/src/RegisterMoL.cc index ee001ea..0a76a2f 100644 --- a/ML_BSSN/src/RegisterMoL.cc +++ b/ML_BSSN/src/RegisterMoL.cc @@ -24,7 +24,22 @@ extern "C" void ML_BSSN_RegisterVars(CCTK_ARGUMENTS) ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::Xt1"), CCTK_VarIndex("ML_BSSN::Xt1rhs")); ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::Xt2"), CCTK_VarIndex("ML_BSSN::Xt2rhs")); ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::Xt3"), CCTK_VarIndex("ML_BSSN::Xt3rhs")); - ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::alpha"), CCTK_VarIndex("ML_BSSN::alpharhs")); + + if (CCTK_EQUALS(lapse_evolution_method, "ML_BSSN")) + ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::alpha"), CCTK_VarIndex("ML_BSSN::alpharhs")); + else + ierr += MoLRegisterConstrained(CCTK_VarIndex("ML_BSSN::alpha")); + + if (CCTK_EQUALS(shift_evolution_method, "ML_BSSN")) { + ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::beta1"), CCTK_VarIndex("ML_BSSN::beta1rhs")); + ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::beta2"), CCTK_VarIndex("ML_BSSN::beta2rhs")); + ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::beta3"), CCTK_VarIndex("ML_BSSN::beta3rhs")); + } else { + ierr += MoLRegisterConstrained(CCTK_VarIndex("ML_BSSN::beta1")); + ierr += MoLRegisterConstrained(CCTK_VarIndex("ML_BSSN::beta2")); + ierr += MoLRegisterConstrained(CCTK_VarIndex("ML_BSSN::beta3")); + } + ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::phi"), CCTK_VarIndex("ML_BSSN::phirhs")); ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::gt11"), CCTK_VarIndex("ML_BSSN::gt11rhs")); ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::gt12"), CCTK_VarIndex("ML_BSSN::gt12rhs")); @@ -32,9 +47,6 @@ extern "C" void ML_BSSN_RegisterVars(CCTK_ARGUMENTS) ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::gt22"), CCTK_VarIndex("ML_BSSN::gt22rhs")); ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::gt23"), CCTK_VarIndex("ML_BSSN::gt23rhs")); ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::gt33"), CCTK_VarIndex("ML_BSSN::gt33rhs")); - ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::beta1"), CCTK_VarIndex("ML_BSSN::beta1rhs")); - ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::beta2"), CCTK_VarIndex("ML_BSSN::beta2rhs")); - ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::beta3"), CCTK_VarIndex("ML_BSSN::beta3rhs")); ierr += MoLRegisterEvolved(CCTK_VarIndex("ML_BSSN::trK"), CCTK_VarIndex("ML_BSSN::trKrhs")); /* Register all the evolved Array functions with MoL */ |