diff options
author | Ian Hinder <ian.hinder@aei.mpg.de> | 2013-09-12 22:09:32 +0200 |
---|---|---|
committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2013-09-12 22:09:32 +0200 |
commit | caadfb1214b142be4512d86994fdfdac9a45206b (patch) | |
tree | b20c527044ad74f0cb953296f456c8c429cb9cb7 /ML_ADM | |
parent | 8c4ac493e00d4cd3978d4504372b3c37b85c64ea (diff) |
Regenerate code (Kranc commit 3237af2566109f099ee78a7effb511a65dc18ae0)
Diffstat (limited to 'ML_ADM')
-rw-r--r-- | ML_ADM/src/Boundaries.cc | 728 | ||||
-rw-r--r-- | ML_ADM/src/RegisterSymmetries.cc | 34 |
2 files changed, 381 insertions, 381 deletions
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,120 +575,6 @@ 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(beta1_bound, "radiative")) { /* select radiation boundary condition */ @@ -632,6 +632,23 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) } + if (CCTK_EQUALS(ML_metric_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!"); + + 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 Scalar BC for ML_ADM::ML_metric!"); + + } + if (CCTK_EQUALS(ML_curv_bound, "scalar")) { /* select scalar boundary condition */ @@ -666,23 +683,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(ML_metric_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!"); - - 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 Scalar BC for ML_ADM::ML_metric!"); - - } - if (CCTK_EQUALS(ML_shift_bound, "scalar")) { /* select scalar boundary condition */ @@ -700,6 +700,108 @@ 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(K11_bound, "scalar")) { /* select scalar boundary condition */ @@ -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; |