From caadfb1214b142be4512d86994fdfdac9a45206b Mon Sep 17 00:00:00 2001 From: Ian Hinder Date: Thu, 12 Sep 2013 22:09:32 +0200 Subject: Regenerate code (Kranc commit 3237af2566109f099ee78a7effb511a65dc18ae0) --- ML_ADM/src/Boundaries.cc | 736 +++++++++++++++++++-------------------- ML_ADM/src/RegisterSymmetries.cc | 34 +- 2 files changed, 385 insertions(+), 385 deletions(-) (limited to 'ML_ADM') diff --git a/ML_ADM/src/Boundaries.cc b/ML_ADM/src/Boundaries.cc index f150f48..bc3f0e5 100644 --- a/ML_ADM/src/Boundaries.cc +++ b/ML_ADM/src/Boundaries.cc @@ -32,6 +32,17 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0; + if (CCTK_EQUALS(ML_metric_bound, "none" ) || + CCTK_EQUALS(ML_metric_bound, "static") || + CCTK_EQUALS(ML_metric_bound, "flat" ) || + CCTK_EQUALS(ML_metric_bound, "zero" ) ) + { + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::ML_metric", ML_metric_bound); + if (ierr < 0) + CCTK_WARN(0, "Failed to register ML_metric_bound BC for ML_ADM::ML_metric!"); + } + if (CCTK_EQUALS(ML_curv_bound, "none" ) || CCTK_EQUALS(ML_curv_bound, "static") || CCTK_EQUALS(ML_curv_bound, "flat" ) || @@ -54,17 +65,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) CCTK_WARN(0, "Failed to register ML_lapse_bound BC for ML_ADM::ML_lapse!"); } - if (CCTK_EQUALS(ML_metric_bound, "none" ) || - CCTK_EQUALS(ML_metric_bound, "static") || - CCTK_EQUALS(ML_metric_bound, "flat" ) || - CCTK_EQUALS(ML_metric_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_ADM::ML_metric", ML_metric_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register ML_metric_bound BC for ML_ADM::ML_metric!"); - } - if (CCTK_EQUALS(ML_shift_bound, "none" ) || CCTK_EQUALS(ML_shift_bound, "static") || CCTK_EQUALS(ML_shift_bound, "flat" ) || @@ -76,6 +76,72 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) CCTK_WARN(0, "Failed to register ML_shift_bound BC for ML_ADM::ML_shift!"); } + if (CCTK_EQUALS(g11_bound, "none" ) || + CCTK_EQUALS(g11_bound, "static") || + CCTK_EQUALS(g11_bound, "flat" ) || + CCTK_EQUALS(g11_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g11", g11_bound); + if (ierr < 0) + CCTK_WARN(0, "Failed to register g11_bound BC for ML_ADM::g11!"); + } + + if (CCTK_EQUALS(g12_bound, "none" ) || + CCTK_EQUALS(g12_bound, "static") || + CCTK_EQUALS(g12_bound, "flat" ) || + CCTK_EQUALS(g12_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g12", g12_bound); + if (ierr < 0) + CCTK_WARN(0, "Failed to register g12_bound BC for ML_ADM::g12!"); + } + + if (CCTK_EQUALS(g13_bound, "none" ) || + CCTK_EQUALS(g13_bound, "static") || + CCTK_EQUALS(g13_bound, "flat" ) || + CCTK_EQUALS(g13_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g13", g13_bound); + if (ierr < 0) + CCTK_WARN(0, "Failed to register g13_bound BC for ML_ADM::g13!"); + } + + if (CCTK_EQUALS(g22_bound, "none" ) || + CCTK_EQUALS(g22_bound, "static") || + CCTK_EQUALS(g22_bound, "flat" ) || + CCTK_EQUALS(g22_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g22", g22_bound); + if (ierr < 0) + CCTK_WARN(0, "Failed to register g22_bound BC for ML_ADM::g22!"); + } + + if (CCTK_EQUALS(g23_bound, "none" ) || + CCTK_EQUALS(g23_bound, "static") || + CCTK_EQUALS(g23_bound, "flat" ) || + CCTK_EQUALS(g23_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g23", g23_bound); + if (ierr < 0) + CCTK_WARN(0, "Failed to register g23_bound BC for ML_ADM::g23!"); + } + + if (CCTK_EQUALS(g33_bound, "none" ) || + CCTK_EQUALS(g33_bound, "static") || + CCTK_EQUALS(g33_bound, "flat" ) || + CCTK_EQUALS(g33_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g33", g33_bound); + if (ierr < 0) + CCTK_WARN(0, "Failed to register g33_bound BC for ML_ADM::g33!"); + } + if (CCTK_EQUALS(K11_bound, "none" ) || CCTK_EQUALS(K11_bound, "static") || CCTK_EQUALS(K11_bound, "flat" ) || @@ -153,72 +219,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) CCTK_WARN(0, "Failed to register alpha_bound BC for ML_ADM::alpha!"); } - if (CCTK_EQUALS(g11_bound, "none" ) || - CCTK_EQUALS(g11_bound, "static") || - CCTK_EQUALS(g11_bound, "flat" ) || - CCTK_EQUALS(g11_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_ADM::g11", g11_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register g11_bound BC for ML_ADM::g11!"); - } - - if (CCTK_EQUALS(g12_bound, "none" ) || - CCTK_EQUALS(g12_bound, "static") || - CCTK_EQUALS(g12_bound, "flat" ) || - CCTK_EQUALS(g12_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_ADM::g12", g12_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register g12_bound BC for ML_ADM::g12!"); - } - - if (CCTK_EQUALS(g13_bound, "none" ) || - CCTK_EQUALS(g13_bound, "static") || - CCTK_EQUALS(g13_bound, "flat" ) || - CCTK_EQUALS(g13_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_ADM::g13", g13_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register g13_bound BC for ML_ADM::g13!"); - } - - if (CCTK_EQUALS(g22_bound, "none" ) || - CCTK_EQUALS(g22_bound, "static") || - CCTK_EQUALS(g22_bound, "flat" ) || - CCTK_EQUALS(g22_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_ADM::g22", g22_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register g22_bound BC for ML_ADM::g22!"); - } - - if (CCTK_EQUALS(g23_bound, "none" ) || - CCTK_EQUALS(g23_bound, "static") || - CCTK_EQUALS(g23_bound, "flat" ) || - CCTK_EQUALS(g23_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_ADM::g23", g23_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register g23_bound BC for ML_ADM::g23!"); - } - - if (CCTK_EQUALS(g33_bound, "none" ) || - CCTK_EQUALS(g33_bound, "static") || - CCTK_EQUALS(g33_bound, "flat" ) || - CCTK_EQUALS(g33_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_ADM::g33", g33_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register g33_bound BC for ML_ADM::g33!"); - } - if (CCTK_EQUALS(beta1_bound, "none" ) || CCTK_EQUALS(beta1_bound, "static") || CCTK_EQUALS(beta1_bound, "flat" ) || @@ -252,6 +252,25 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) CCTK_WARN(0, "Failed to register beta3_bound BC for ML_ADM::beta3!"); } + if (CCTK_EQUALS(ML_metric_bound, "radiative")) + { + /* select radiation boundary condition */ + static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_ML_metric_bound , ML_metric_bound_limit, "LIMIT") < 0) + CCTK_WARN(0, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_speed, "SPEED") < 0) + CCTK_WARN(0, "could not set SPEED value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound, + "ML_ADM::ML_metric", "Radiation"); + + if (ierr < 0) + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::ML_metric!"); + + } + if (CCTK_EQUALS(ML_curv_bound, "radiative")) { /* select radiation boundary condition */ @@ -290,25 +309,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(ML_metric_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_ML_metric_bound , ML_metric_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound, - "ML_ADM::ML_metric", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::ML_metric!"); - - } - if (CCTK_EQUALS(ML_shift_bound, "radiative")) { /* select radiation boundary condition */ @@ -328,6 +328,120 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) } + if (CCTK_EQUALS(g11_bound, "radiative")) + { + /* select radiation boundary condition */ + static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g11_bound , g11_bound_limit, "LIMIT") < 0) + CCTK_WARN(0, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g11_bound ,g11_bound_speed, "SPEED") < 0) + CCTK_WARN(0, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound, + "ML_ADM::g11", "Radiation"); + + if (ierr < 0) + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g11!"); + + } + + if (CCTK_EQUALS(g12_bound, "radiative")) + { + /* select radiation boundary condition */ + static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g12_bound , g12_bound_limit, "LIMIT") < 0) + CCTK_WARN(0, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g12_bound ,g12_bound_speed, "SPEED") < 0) + CCTK_WARN(0, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound, + "ML_ADM::g12", "Radiation"); + + if (ierr < 0) + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g12!"); + + } + + if (CCTK_EQUALS(g13_bound, "radiative")) + { + /* select radiation boundary condition */ + static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g13_bound , g13_bound_limit, "LIMIT") < 0) + CCTK_WARN(0, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g13_bound ,g13_bound_speed, "SPEED") < 0) + CCTK_WARN(0, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound, + "ML_ADM::g13", "Radiation"); + + if (ierr < 0) + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g13!"); + + } + + if (CCTK_EQUALS(g22_bound, "radiative")) + { + /* select radiation boundary condition */ + static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g22_bound , g22_bound_limit, "LIMIT") < 0) + CCTK_WARN(0, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g22_bound ,g22_bound_speed, "SPEED") < 0) + CCTK_WARN(0, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound, + "ML_ADM::g22", "Radiation"); + + if (ierr < 0) + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g22!"); + + } + + if (CCTK_EQUALS(g23_bound, "radiative")) + { + /* select radiation boundary condition */ + static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g23_bound , g23_bound_limit, "LIMIT") < 0) + CCTK_WARN(0, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g23_bound ,g23_bound_speed, "SPEED") < 0) + CCTK_WARN(0, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound, + "ML_ADM::g23", "Radiation"); + + if (ierr < 0) + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g23!"); + + } + + if (CCTK_EQUALS(g33_bound, "radiative")) + { + /* select radiation boundary condition */ + static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g33_bound , g33_bound_limit, "LIMIT") < 0) + CCTK_WARN(0, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g33_bound ,g33_bound_speed, "SPEED") < 0) + CCTK_WARN(0, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound, + "ML_ADM::g33", "Radiation"); + + if (ierr < 0) + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g33!"); + + } + if (CCTK_EQUALS(K11_bound, "radiative")) { /* select radiation boundary condition */ @@ -461,242 +575,230 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(g11_bound, "radiative")) + if (CCTK_EQUALS(beta1_bound, "radiative")) { /* select radiation boundary condition */ - static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g11_bound , g11_bound_limit, "LIMIT") < 0) + static CCTK_INT handle_beta1_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_beta1_bound < 0) handle_beta1_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta1_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_beta1_bound , beta1_bound_limit, "LIMIT") < 0) CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_g11_bound ,g11_bound_speed, "SPEED") < 0) + if (Util_TableSetReal(handle_beta1_bound ,beta1_bound_speed, "SPEED") < 0) CCTK_WARN(0, "could not set SPEED value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound, - "ML_ADM::g11", "Radiation"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta1_bound, + "ML_ADM::beta1", "Radiation"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g11!"); + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::beta1!"); } - if (CCTK_EQUALS(g12_bound, "radiative")) + if (CCTK_EQUALS(beta2_bound, "radiative")) { /* select radiation boundary condition */ - static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g12_bound , g12_bound_limit, "LIMIT") < 0) + static CCTK_INT handle_beta2_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_beta2_bound < 0) handle_beta2_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta2_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_beta2_bound , beta2_bound_limit, "LIMIT") < 0) CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_g12_bound ,g12_bound_speed, "SPEED") < 0) + if (Util_TableSetReal(handle_beta2_bound ,beta2_bound_speed, "SPEED") < 0) CCTK_WARN(0, "could not set SPEED value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound, - "ML_ADM::g12", "Radiation"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta2_bound, + "ML_ADM::beta2", "Radiation"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g12!"); + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::beta2!"); } - if (CCTK_EQUALS(g13_bound, "radiative")) + if (CCTK_EQUALS(beta3_bound, "radiative")) { /* select radiation boundary condition */ - static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g13_bound , g13_bound_limit, "LIMIT") < 0) + static CCTK_INT handle_beta3_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_beta3_bound < 0) handle_beta3_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta3_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_beta3_bound , beta3_bound_limit, "LIMIT") < 0) CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_g13_bound ,g13_bound_speed, "SPEED") < 0) + if (Util_TableSetReal(handle_beta3_bound ,beta3_bound_speed, "SPEED") < 0) CCTK_WARN(0, "could not set SPEED value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound, - "ML_ADM::g13", "Radiation"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta3_bound, + "ML_ADM::beta3", "Radiation"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g13!"); + CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::beta3!"); } - if (CCTK_EQUALS(g22_bound, "radiative")) + if (CCTK_EQUALS(ML_metric_bound, "scalar")) { - /* select radiation boundary condition */ - static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g22_bound , g22_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_g22_bound ,g22_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); + /* select scalar boundary condition */ + static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound, - "ML_ADM::g22", "Radiation"); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound, + "ML_ADM::ML_metric", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g22!"); + CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_metric!"); } - if (CCTK_EQUALS(g23_bound, "radiative")) + if (CCTK_EQUALS(ML_curv_bound, "scalar")) { - /* select radiation boundary condition */ - static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g23_bound , g23_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_g23_bound ,g23_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); + /* select scalar boundary condition */ + static CCTK_INT handle_ML_curv_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_ML_curv_bound < 0) handle_ML_curv_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_ML_curv_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_ML_curv_bound ,ML_curv_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound, - "ML_ADM::g23", "Radiation"); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_curv_bound, + "ML_ADM::ML_curv", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g23!"); + CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_curv!"); } - if (CCTK_EQUALS(g33_bound, "radiative")) + if (CCTK_EQUALS(ML_lapse_bound, "scalar")) { - /* select radiation boundary condition */ - static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g33_bound , g33_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_g33_bound ,g33_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); + /* select scalar boundary condition */ + static CCTK_INT handle_ML_lapse_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_ML_lapse_bound < 0) handle_ML_lapse_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_ML_lapse_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_ML_lapse_bound ,ML_lapse_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound, - "ML_ADM::g33", "Radiation"); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_lapse_bound, + "ML_ADM::ML_lapse", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g33!"); + CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_lapse!"); } - if (CCTK_EQUALS(beta1_bound, "radiative")) + if (CCTK_EQUALS(ML_shift_bound, "scalar")) { - /* select radiation boundary condition */ - static CCTK_INT handle_beta1_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_beta1_bound < 0) handle_beta1_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_beta1_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_beta1_bound , beta1_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_beta1_bound ,beta1_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); + /* select scalar boundary condition */ + static CCTK_INT handle_ML_shift_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_ML_shift_bound < 0) handle_ML_shift_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_ML_shift_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_ML_shift_bound ,ML_shift_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta1_bound, - "ML_ADM::beta1", "Radiation"); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_shift_bound, + "ML_ADM::ML_shift", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::beta1!"); + CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_shift!"); } - if (CCTK_EQUALS(beta2_bound, "radiative")) + if (CCTK_EQUALS(g11_bound, "scalar")) { - /* select radiation boundary condition */ - static CCTK_INT handle_beta2_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_beta2_bound < 0) handle_beta2_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_beta2_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_beta2_bound , beta2_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_beta2_bound ,beta2_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); + /* select scalar boundary condition */ + static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g11_bound ,g11_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta2_bound, - "ML_ADM::beta2", "Radiation"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound, + "ML_ADM::g11", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::beta2!"); + CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g11!"); } - if (CCTK_EQUALS(beta3_bound, "radiative")) + if (CCTK_EQUALS(g12_bound, "scalar")) { - /* select radiation boundary condition */ - static CCTK_INT handle_beta3_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_beta3_bound < 0) handle_beta3_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_beta3_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_beta3_bound , beta3_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_beta3_bound ,beta3_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); + /* select scalar boundary condition */ + static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g12_bound ,g12_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta3_bound, - "ML_ADM::beta3", "Radiation"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound, + "ML_ADM::g12", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::beta3!"); + CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g12!"); } - if (CCTK_EQUALS(ML_curv_bound, "scalar")) + if (CCTK_EQUALS(g13_bound, "scalar")) { /* select scalar boundary condition */ - static CCTK_INT handle_ML_curv_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_ML_curv_bound < 0) handle_ML_curv_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_ML_curv_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_ML_curv_bound ,ML_curv_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); + static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g13_bound ,g13_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_curv_bound, - "ML_ADM::ML_curv", "scalar"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound, + "ML_ADM::g13", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_curv!"); + CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g13!"); } - if (CCTK_EQUALS(ML_lapse_bound, "scalar")) + if (CCTK_EQUALS(g22_bound, "scalar")) { /* select scalar boundary condition */ - static CCTK_INT handle_ML_lapse_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_ML_lapse_bound < 0) handle_ML_lapse_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_ML_lapse_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_ML_lapse_bound ,ML_lapse_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); + static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g22_bound ,g22_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_lapse_bound, - "ML_ADM::ML_lapse", "scalar"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound, + "ML_ADM::g22", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_lapse!"); + CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g22!"); } - if (CCTK_EQUALS(ML_metric_bound, "scalar")) + if (CCTK_EQUALS(g23_bound, "scalar")) { /* select scalar boundary condition */ - static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); + static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g23_bound ,g23_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound, - "ML_ADM::ML_metric", "scalar"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound, + "ML_ADM::g23", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_metric!"); + CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g23!"); } - if (CCTK_EQUALS(ML_shift_bound, "scalar")) + if (CCTK_EQUALS(g33_bound, "scalar")) { /* select scalar boundary condition */ - static CCTK_INT handle_ML_shift_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_ML_shift_bound < 0) handle_ML_shift_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_ML_shift_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_ML_shift_bound ,ML_shift_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); + static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1; + if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!"); + if (Util_TableSetReal(handle_g33_bound ,g33_bound_scalar, "SCALAR") < 0) + CCTK_WARN(0, "could not set SCALAR value in table!"); - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_shift_bound, - "ML_ADM::ML_shift", "scalar"); + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound, + "ML_ADM::g33", "scalar"); if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_shift!"); + CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g33!"); } @@ -819,108 +921,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(g11_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g11_bound ,g11_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound, - "ML_ADM::g11", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g11!"); - - } - - if (CCTK_EQUALS(g12_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g12_bound ,g12_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound, - "ML_ADM::g12", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g12!"); - - } - - if (CCTK_EQUALS(g13_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g13_bound ,g13_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound, - "ML_ADM::g13", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g13!"); - - } - - if (CCTK_EQUALS(g22_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g22_bound ,g22_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound, - "ML_ADM::g22", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g22!"); - - } - - if (CCTK_EQUALS(g23_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g23_bound ,g23_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound, - "ML_ADM::g23", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g23!"); - - } - - if (CCTK_EQUALS(g33_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1; - if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_g33_bound ,g33_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound, - "ML_ADM::g33", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g33!"); - - } - if (CCTK_EQUALS(beta1_bound, "scalar")) { /* select scalar boundary condition */ @@ -977,6 +977,11 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) /* template for entries in parameter file: +#$bound$#ML_ADM::ML_metric_bound = "skip" +#$bound$#ML_ADM::ML_metric_bound_speed = 1.0 +#$bound$#ML_ADM::ML_metric_bound_limit = 0.0 +#$bound$#ML_ADM::ML_metric_bound_scalar = 0.0 + #$bound$#ML_ADM::ML_curv_bound = "skip" #$bound$#ML_ADM::ML_curv_bound_speed = 1.0 #$bound$#ML_ADM::ML_curv_bound_limit = 0.0 @@ -987,16 +992,41 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) #$bound$#ML_ADM::ML_lapse_bound_limit = 0.0 #$bound$#ML_ADM::ML_lapse_bound_scalar = 0.0 -#$bound$#ML_ADM::ML_metric_bound = "skip" -#$bound$#ML_ADM::ML_metric_bound_speed = 1.0 -#$bound$#ML_ADM::ML_metric_bound_limit = 0.0 -#$bound$#ML_ADM::ML_metric_bound_scalar = 0.0 - #$bound$#ML_ADM::ML_shift_bound = "skip" #$bound$#ML_ADM::ML_shift_bound_speed = 1.0 #$bound$#ML_ADM::ML_shift_bound_limit = 0.0 #$bound$#ML_ADM::ML_shift_bound_scalar = 0.0 +#$bound$#ML_ADM::g11_bound = "skip" +#$bound$#ML_ADM::g11_bound_speed = 1.0 +#$bound$#ML_ADM::g11_bound_limit = 0.0 +#$bound$#ML_ADM::g11_bound_scalar = 0.0 + +#$bound$#ML_ADM::g12_bound = "skip" +#$bound$#ML_ADM::g12_bound_speed = 1.0 +#$bound$#ML_ADM::g12_bound_limit = 0.0 +#$bound$#ML_ADM::g12_bound_scalar = 0.0 + +#$bound$#ML_ADM::g13_bound = "skip" +#$bound$#ML_ADM::g13_bound_speed = 1.0 +#$bound$#ML_ADM::g13_bound_limit = 0.0 +#$bound$#ML_ADM::g13_bound_scalar = 0.0 + +#$bound$#ML_ADM::g22_bound = "skip" +#$bound$#ML_ADM::g22_bound_speed = 1.0 +#$bound$#ML_ADM::g22_bound_limit = 0.0 +#$bound$#ML_ADM::g22_bound_scalar = 0.0 + +#$bound$#ML_ADM::g23_bound = "skip" +#$bound$#ML_ADM::g23_bound_speed = 1.0 +#$bound$#ML_ADM::g23_bound_limit = 0.0 +#$bound$#ML_ADM::g23_bound_scalar = 0.0 + +#$bound$#ML_ADM::g33_bound = "skip" +#$bound$#ML_ADM::g33_bound_speed = 1.0 +#$bound$#ML_ADM::g33_bound_limit = 0.0 +#$bound$#ML_ADM::g33_bound_scalar = 0.0 + #$bound$#ML_ADM::K11_bound = "skip" #$bound$#ML_ADM::K11_bound_speed = 1.0 #$bound$#ML_ADM::K11_bound_limit = 0.0 @@ -1032,36 +1062,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) #$bound$#ML_ADM::alpha_bound_limit = 0.0 #$bound$#ML_ADM::alpha_bound_scalar = 0.0 -#$bound$#ML_ADM::g11_bound = "skip" -#$bound$#ML_ADM::g11_bound_speed = 1.0 -#$bound$#ML_ADM::g11_bound_limit = 0.0 -#$bound$#ML_ADM::g11_bound_scalar = 0.0 - -#$bound$#ML_ADM::g12_bound = "skip" -#$bound$#ML_ADM::g12_bound_speed = 1.0 -#$bound$#ML_ADM::g12_bound_limit = 0.0 -#$bound$#ML_ADM::g12_bound_scalar = 0.0 - -#$bound$#ML_ADM::g13_bound = "skip" -#$bound$#ML_ADM::g13_bound_speed = 1.0 -#$bound$#ML_ADM::g13_bound_limit = 0.0 -#$bound$#ML_ADM::g13_bound_scalar = 0.0 - -#$bound$#ML_ADM::g22_bound = "skip" -#$bound$#ML_ADM::g22_bound_speed = 1.0 -#$bound$#ML_ADM::g22_bound_limit = 0.0 -#$bound$#ML_ADM::g22_bound_scalar = 0.0 - -#$bound$#ML_ADM::g23_bound = "skip" -#$bound$#ML_ADM::g23_bound_speed = 1.0 -#$bound$#ML_ADM::g23_bound_limit = 0.0 -#$bound$#ML_ADM::g23_bound_scalar = 0.0 - -#$bound$#ML_ADM::g33_bound = "skip" -#$bound$#ML_ADM::g33_bound_speed = 1.0 -#$bound$#ML_ADM::g33_bound_limit = 0.0 -#$bound$#ML_ADM::g33_bound_scalar = 0.0 - #$bound$#ML_ADM::beta1_bound = "skip" #$bound$#ML_ADM::beta1_bound_speed = 1.0 #$bound$#ML_ADM::beta1_bound_limit = 0.0 diff --git a/ML_ADM/src/RegisterSymmetries.cc b/ML_ADM/src/RegisterSymmetries.cc index 0d05f97..0f972fb 100644 --- a/ML_ADM/src/RegisterSymmetries.cc +++ b/ML_ADM/src/RegisterSymmetries.cc @@ -19,67 +19,67 @@ extern "C" void ML_ADM_RegisterSymmetries(CCTK_ARGUMENTS) sym[0] = 1; sym[1] = 1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::K11"); + SetCartSymVN(cctkGH, sym, "ML_ADM::g11"); sym[0] = -1; sym[1] = -1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::K12"); + SetCartSymVN(cctkGH, sym, "ML_ADM::g12"); sym[0] = -1; sym[1] = 1; sym[2] = -1; - SetCartSymVN(cctkGH, sym, "ML_ADM::K13"); + SetCartSymVN(cctkGH, sym, "ML_ADM::g13"); sym[0] = 1; sym[1] = 1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::K22"); + SetCartSymVN(cctkGH, sym, "ML_ADM::g22"); sym[0] = 1; sym[1] = -1; sym[2] = -1; - SetCartSymVN(cctkGH, sym, "ML_ADM::K23"); - - sym[0] = 1; - sym[1] = 1; - sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::K33"); + SetCartSymVN(cctkGH, sym, "ML_ADM::g23"); sym[0] = 1; sym[1] = 1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::alpha"); + SetCartSymVN(cctkGH, sym, "ML_ADM::g33"); sym[0] = 1; sym[1] = 1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::g11"); + SetCartSymVN(cctkGH, sym, "ML_ADM::K11"); sym[0] = -1; sym[1] = -1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::g12"); + SetCartSymVN(cctkGH, sym, "ML_ADM::K12"); sym[0] = -1; sym[1] = 1; sym[2] = -1; - SetCartSymVN(cctkGH, sym, "ML_ADM::g13"); + SetCartSymVN(cctkGH, sym, "ML_ADM::K13"); sym[0] = 1; sym[1] = 1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::g22"); + SetCartSymVN(cctkGH, sym, "ML_ADM::K22"); sym[0] = 1; sym[1] = -1; sym[2] = -1; - SetCartSymVN(cctkGH, sym, "ML_ADM::g23"); + SetCartSymVN(cctkGH, sym, "ML_ADM::K23"); sym[0] = 1; sym[1] = 1; sym[2] = 1; - SetCartSymVN(cctkGH, sym, "ML_ADM::g33"); + SetCartSymVN(cctkGH, sym, "ML_ADM::K33"); + + sym[0] = 1; + sym[1] = 1; + sym[2] = 1; + SetCartSymVN(cctkGH, sym, "ML_ADM::alpha"); sym[0] = -1; sym[1] = 1; -- cgit v1.2.3