diff options
Diffstat (limited to 'ML_ADM/src/Boundaries.c')
| -rw-r--r-- | ML_ADM/src/Boundaries.c | 1044 |
1 files changed, 1044 insertions, 0 deletions
diff --git a/ML_ADM/src/Boundaries.c b/ML_ADM/src/Boundaries.c new file mode 100644 index 0000000..a79a204 --- /dev/null +++ b/ML_ADM/src/Boundaries.c @@ -0,0 +1,1044 @@ +/* File produced by user eschnett */ +/* Produced with Mathematica Version 6.0 for Mac OS X x86 (32-bit) (April 20, 2007) */ + +/* Mathematica script written by Ian Hinder and Sascha Husa */ + +#include "cctk.h" +#include "cctk_Arguments.h" +#include "cctk_Parameters.h" +#include "cctk_Faces.h" +#include "util_Table.h" +#include "Symmetry.h" + + +/* the boundary treatment is split into 3 steps: */ +/* 1. excision */ +/* 2. symmetries */ +/* 3. "other" boundary conditions, e.g. radiative */ + +/* to simplify scheduling and testing, the 3 steps */ +/* are currently applied in separate functions */ + + +void ML_ADM_CheckBoundaries(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + return; +} + +void ML_ADM_ApplyBoundConds(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS + DECLARE_CCTK_PARAMETERS + + CCTK_INT ierr = 0; + + if (CCTK_EQUALS(curv_bound, "none" ) || + CCTK_EQUALS(curv_bound, "static") || + CCTK_EQUALS(curv_bound, "flat" ) || + CCTK_EQUALS(curv_bound, "zero" ) ) + { + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::curv", curv_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register curv_bound BC for ML_ADM::curv!"); + } + + if (CCTK_EQUALS(lapse_bound, "none" ) || + CCTK_EQUALS(lapse_bound, "static") || + CCTK_EQUALS(lapse_bound, "flat" ) || + CCTK_EQUALS(lapse_bound, "zero" ) ) + { + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::lapse", lapse_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register lapse_bound BC for ML_ADM::lapse!"); + } + + if (CCTK_EQUALS(metric_bound, "none" ) || + CCTK_EQUALS(metric_bound, "static") || + CCTK_EQUALS(metric_bound, "flat" ) || + CCTK_EQUALS(metric_bound, "zero" ) ) + { + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::metric", metric_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register metric_bound BC for ML_ADM::metric!"); + } + + if (CCTK_EQUALS(shift_bound, "none" ) || + CCTK_EQUALS(shift_bound, "static") || + CCTK_EQUALS(shift_bound, "flat" ) || + CCTK_EQUALS(shift_bound, "zero" ) ) + { + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::shift", shift_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register shift_bound BC for ML_ADM::shift!"); + } + + if (CCTK_EQUALS(K11_bound, "none" ) || + CCTK_EQUALS(K11_bound, "static") || + CCTK_EQUALS(K11_bound, "flat" ) || + CCTK_EQUALS(K11_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::K11", K11_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register K11_bound BC for ML_ADM::K11!"); + } + + if (CCTK_EQUALS(K21_bound, "none" ) || + CCTK_EQUALS(K21_bound, "static") || + CCTK_EQUALS(K21_bound, "flat" ) || + CCTK_EQUALS(K21_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::K21", K21_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register K21_bound BC for ML_ADM::K21!"); + } + + if (CCTK_EQUALS(K31_bound, "none" ) || + CCTK_EQUALS(K31_bound, "static") || + CCTK_EQUALS(K31_bound, "flat" ) || + CCTK_EQUALS(K31_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::K31", K31_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register K31_bound BC for ML_ADM::K31!"); + } + + if (CCTK_EQUALS(K22_bound, "none" ) || + CCTK_EQUALS(K22_bound, "static") || + CCTK_EQUALS(K22_bound, "flat" ) || + CCTK_EQUALS(K22_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::K22", K22_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register K22_bound BC for ML_ADM::K22!"); + } + + if (CCTK_EQUALS(K32_bound, "none" ) || + CCTK_EQUALS(K32_bound, "static") || + CCTK_EQUALS(K32_bound, "flat" ) || + CCTK_EQUALS(K32_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::K32", K32_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register K32_bound BC for ML_ADM::K32!"); + } + + if (CCTK_EQUALS(K33_bound, "none" ) || + CCTK_EQUALS(K33_bound, "static") || + CCTK_EQUALS(K33_bound, "flat" ) || + CCTK_EQUALS(K33_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::K33", K33_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register K33_bound BC for ML_ADM::K33!"); + } + + if (CCTK_EQUALS(alpha_bound, "none" ) || + CCTK_EQUALS(alpha_bound, "static") || + CCTK_EQUALS(alpha_bound, "flat" ) || + CCTK_EQUALS(alpha_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::alpha", alpha_bound); + if (ierr < 0) + CCTK_WARN(-1, "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(-1, "Failed to register g11_bound BC for ML_ADM::g11!"); + } + + if (CCTK_EQUALS(g21_bound, "none" ) || + CCTK_EQUALS(g21_bound, "static") || + CCTK_EQUALS(g21_bound, "flat" ) || + CCTK_EQUALS(g21_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g21", g21_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register g21_bound BC for ML_ADM::g21!"); + } + + if (CCTK_EQUALS(g31_bound, "none" ) || + CCTK_EQUALS(g31_bound, "static") || + CCTK_EQUALS(g31_bound, "flat" ) || + CCTK_EQUALS(g31_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g31", g31_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register g31_bound BC for ML_ADM::g31!"); + } + + 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(-1, "Failed to register g22_bound BC for ML_ADM::g22!"); + } + + if (CCTK_EQUALS(g32_bound, "none" ) || + CCTK_EQUALS(g32_bound, "static") || + CCTK_EQUALS(g32_bound, "flat" ) || + CCTK_EQUALS(g32_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::g32", g32_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register g32_bound BC for ML_ADM::g32!"); + } + + 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(-1, "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" ) || + CCTK_EQUALS(beta1_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::beta1", beta1_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register beta1_bound BC for ML_ADM::beta1!"); + } + + if (CCTK_EQUALS(beta2_bound, "none" ) || + CCTK_EQUALS(beta2_bound, "static") || + CCTK_EQUALS(beta2_bound, "flat" ) || + CCTK_EQUALS(beta2_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::beta2", beta2_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register beta2_bound BC for ML_ADM::beta2!"); + } + + if (CCTK_EQUALS(beta3_bound, "none" ) || + CCTK_EQUALS(beta3_bound, "static") || + CCTK_EQUALS(beta3_bound, "flat" ) || + CCTK_EQUALS(beta3_bound, "zero" ) ) + { + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, + "ML_ADM::beta3", beta3_bound); + if (ierr < 0) + CCTK_WARN(-1, "Failed to register beta3_bound BC for ML_ADM::beta3!"); + } + + if (CCTK_EQUALS(curv_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_curv_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_curv_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_curv_bound , curv_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_curv_bound ,curv_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_curv_bound, + "ML_ADM::curv", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::curv!"); + + } + + if (CCTK_EQUALS(lapse_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_lapse_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_lapse_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_lapse_bound , lapse_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_lapse_bound ,lapse_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_lapse_bound, + "ML_ADM::lapse", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::lapse!"); + + } + + if (CCTK_EQUALS(metric_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_metric_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_metric_bound , metric_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_metric_bound ,metric_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_metric_bound, + "ML_ADM::metric", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::metric!"); + + } + + if (CCTK_EQUALS(shift_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_shift_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_shift_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_shift_bound , shift_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_shift_bound ,shift_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_shift_bound, + "ML_ADM::shift", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::shift!"); + + } + + if (CCTK_EQUALS(K11_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_K11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K11_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K11_bound , K11_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_K11_bound ,K11_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K11_bound, + "ML_ADM::K11", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::K11!"); + + } + + if (CCTK_EQUALS(K21_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_K21_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K21_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K21_bound , K21_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_K21_bound ,K21_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K21_bound, + "ML_ADM::K21", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::K21!"); + + } + + if (CCTK_EQUALS(K31_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_K31_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K31_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K31_bound , K31_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_K31_bound ,K31_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K31_bound, + "ML_ADM::K31", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::K31!"); + + } + + if (CCTK_EQUALS(K22_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_K22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K22_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K22_bound , K22_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_K22_bound ,K22_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K22_bound, + "ML_ADM::K22", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::K22!"); + + } + + if (CCTK_EQUALS(K32_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_K32_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K32_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K32_bound , K32_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_K32_bound ,K32_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K32_bound, + "ML_ADM::K32", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::K32!"); + + } + + if (CCTK_EQUALS(K33_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_K33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K33_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K33_bound , K33_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_K33_bound ,K33_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K33_bound, + "ML_ADM::K33", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::K33!"); + + } + + if (CCTK_EQUALS(alpha_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_alpha_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_alpha_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_alpha_bound , alpha_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_alpha_bound ,alpha_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_alpha_bound, + "ML_ADM::alpha", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::alpha!"); + + } + + if (CCTK_EQUALS(g11_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g11_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g11_bound , g11_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g11_bound ,g11_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "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(-1, "Failed to register Radiation BC for ML_ADM::g11!"); + + } + + if (CCTK_EQUALS(g21_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_g21_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g21_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g21_bound , g21_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g21_bound ,g21_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g21_bound, + "ML_ADM::g21", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::g21!"); + + } + + if (CCTK_EQUALS(g31_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_g31_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g31_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g31_bound , g31_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g31_bound ,g31_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g31_bound, + "ML_ADM::g31", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::g31!"); + + } + + if (CCTK_EQUALS(g22_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g22_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g22_bound , g22_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g22_bound ,g22_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "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(-1, "Failed to register Radiation BC for ML_ADM::g22!"); + + } + + if (CCTK_EQUALS(g32_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_g32_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g32_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g32_bound , g32_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g32_bound ,g32_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g32_bound, + "ML_ADM::g32", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::g32!"); + + } + + if (CCTK_EQUALS(g33_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g33_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g33_bound , g33_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_g33_bound ,g33_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "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(-1, "Failed to register Radiation BC for ML_ADM::g33!"); + + } + + if (CCTK_EQUALS(beta1_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_beta1_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta1_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_beta1_bound , beta1_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_beta1_bound ,beta1_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta1_bound, + "ML_ADM::beta1", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::beta1!"); + + } + + if (CCTK_EQUALS(beta2_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_beta2_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta2_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_beta2_bound , beta2_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_beta2_bound ,beta2_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta2_bound, + "ML_ADM::beta2", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::beta2!"); + + } + + if (CCTK_EQUALS(beta3_bound, "radiative")) + { + /* apply radiation boundary condition */ + CCTK_INT handle_beta3_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta3_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_beta3_bound , beta3_bound_limit, "LIMIT") < 0) + CCTK_WARN(-1, "could not set LIMIT value in table!"); + if (Util_TableSetReal(handle_beta3_bound ,beta3_bound_speed, "SPEED") < 0) + CCTK_WARN(-1, "could not set SPEED value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta3_bound, + "ML_ADM::beta3", "Radiation"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Radiation BC for ML_ADM::beta3!"); + + } + + if (CCTK_EQUALS(curv_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_curv_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_curv_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_curv_bound ,curv_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_curv_bound, + "ML_ADM::curv", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Scalar BC for ML_ADM::curv!"); + + } + + if (CCTK_EQUALS(lapse_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_lapse_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_lapse_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_lapse_bound ,lapse_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_lapse_bound, + "ML_ADM::lapse", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Scalar BC for ML_ADM::lapse!"); + + } + + if (CCTK_EQUALS(metric_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_metric_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_metric_bound ,metric_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_metric_bound, + "ML_ADM::metric", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Scalar BC for ML_ADM::metric!"); + + } + + if (CCTK_EQUALS(shift_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_shift_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_shift_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_shift_bound ,shift_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_shift_bound, + "ML_ADM::shift", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Failed to register Scalar BC for ML_ADM::shift!"); + + } + + if (CCTK_EQUALS(K11_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_K11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K11_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K11_bound ,K11_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K11_bound, + "ML_ADM::K11", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::K11!"); + + } + + if (CCTK_EQUALS(K21_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_K21_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K21_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K21_bound ,K21_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K21_bound, + "ML_ADM::K21", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::K21!"); + + } + + if (CCTK_EQUALS(K31_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_K31_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K31_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K31_bound ,K31_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K31_bound, + "ML_ADM::K31", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::K31!"); + + } + + if (CCTK_EQUALS(K22_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_K22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K22_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K22_bound ,K22_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K22_bound, + "ML_ADM::K22", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::K22!"); + + } + + if (CCTK_EQUALS(K32_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_K32_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K32_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K32_bound ,K32_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K32_bound, + "ML_ADM::K32", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::K32!"); + + } + + if (CCTK_EQUALS(K33_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_K33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_K33_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_K33_bound ,K33_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_K33_bound, + "ML_ADM::K33", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::K33!"); + + } + + if (CCTK_EQUALS(alpha_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_alpha_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_alpha_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_alpha_bound ,alpha_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_alpha_bound, + "ML_ADM::alpha", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::alpha!"); + + } + + if (CCTK_EQUALS(g11_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g11_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g11_bound ,g11_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "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(-1, "Error in registering Scalar BC for ML_ADM::g11!"); + + } + + if (CCTK_EQUALS(g21_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_g21_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g21_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g21_bound ,g21_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g21_bound, + "ML_ADM::g21", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::g21!"); + + } + + if (CCTK_EQUALS(g31_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_g31_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g31_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g31_bound ,g31_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g31_bound, + "ML_ADM::g31", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::g31!"); + + } + + if (CCTK_EQUALS(g22_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g22_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g22_bound ,g22_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "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(-1, "Error in registering Scalar BC for ML_ADM::g22!"); + + } + + if (CCTK_EQUALS(g32_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_g32_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g32_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g32_bound ,g32_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g32_bound, + "ML_ADM::g32", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::g32!"); + + } + + if (CCTK_EQUALS(g33_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_g33_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_g33_bound ,g33_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "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(-1, "Error in registering Scalar BC for ML_ADM::g33!"); + + } + + if (CCTK_EQUALS(beta1_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_beta1_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta1_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_beta1_bound ,beta1_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta1_bound, + "ML_ADM::beta1", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::beta1!"); + + } + + if (CCTK_EQUALS(beta2_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_beta2_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta2_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_beta2_bound ,beta2_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta2_bound, + "ML_ADM::beta2", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::beta2!"); + + } + + if (CCTK_EQUALS(beta3_bound, "scalar")) + { + /* apply scalar boundary condition */ + CCTK_INT handle_beta3_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); + if (handle_beta3_bound < 0) CCTK_WARN(-1, "could not create table!"); + if (Util_TableSetReal(handle_beta3_bound ,beta3_bound_scalar, "SCALAR") < 0) + CCTK_WARN(-1, "could not set SCALAR value in table!"); + + ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_beta3_bound, + "ML_ADM::beta3", "scalar"); + + if (ierr < 0) + CCTK_WARN(-1, "Error in registering Scalar BC for ML_ADM::beta3!"); + + } + return; +} + + + +/* template for entries in parameter file: +#$bound$#ML_ADM::curv_bound = "skip" +#$bound$#ML_ADM::curv_bound_speed = 1.0 +#$bound$#ML_ADM::curv_bound_limit = 0.0 +#$bound$#ML_ADM::curv_bound_scalar = 0.0 + +#$bound$#ML_ADM::lapse_bound = "skip" +#$bound$#ML_ADM::lapse_bound_speed = 1.0 +#$bound$#ML_ADM::lapse_bound_limit = 0.0 +#$bound$#ML_ADM::lapse_bound_scalar = 0.0 + +#$bound$#ML_ADM::metric_bound = "skip" +#$bound$#ML_ADM::metric_bound_speed = 1.0 +#$bound$#ML_ADM::metric_bound_limit = 0.0 +#$bound$#ML_ADM::metric_bound_scalar = 0.0 + +#$bound$#ML_ADM::shift_bound = "skip" +#$bound$#ML_ADM::shift_bound_speed = 1.0 +#$bound$#ML_ADM::shift_bound_limit = 0.0 +#$bound$#ML_ADM::shift_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 +#$bound$#ML_ADM::K11_bound_scalar = 0.0 + +#$bound$#ML_ADM::K21_bound = "skip" +#$bound$#ML_ADM::K21_bound_speed = 1.0 +#$bound$#ML_ADM::K21_bound_limit = 0.0 +#$bound$#ML_ADM::K21_bound_scalar = 0.0 + +#$bound$#ML_ADM::K31_bound = "skip" +#$bound$#ML_ADM::K31_bound_speed = 1.0 +#$bound$#ML_ADM::K31_bound_limit = 0.0 +#$bound$#ML_ADM::K31_bound_scalar = 0.0 + +#$bound$#ML_ADM::K22_bound = "skip" +#$bound$#ML_ADM::K22_bound_speed = 1.0 +#$bound$#ML_ADM::K22_bound_limit = 0.0 +#$bound$#ML_ADM::K22_bound_scalar = 0.0 + +#$bound$#ML_ADM::K32_bound = "skip" +#$bound$#ML_ADM::K32_bound_speed = 1.0 +#$bound$#ML_ADM::K32_bound_limit = 0.0 +#$bound$#ML_ADM::K32_bound_scalar = 0.0 + +#$bound$#ML_ADM::K33_bound = "skip" +#$bound$#ML_ADM::K33_bound_speed = 1.0 +#$bound$#ML_ADM::K33_bound_limit = 0.0 +#$bound$#ML_ADM::K33_bound_scalar = 0.0 + +#$bound$#ML_ADM::alpha_bound = "skip" +#$bound$#ML_ADM::alpha_bound_speed = 1.0 +#$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::g21_bound = "skip" +#$bound$#ML_ADM::g21_bound_speed = 1.0 +#$bound$#ML_ADM::g21_bound_limit = 0.0 +#$bound$#ML_ADM::g21_bound_scalar = 0.0 + +#$bound$#ML_ADM::g31_bound = "skip" +#$bound$#ML_ADM::g31_bound_speed = 1.0 +#$bound$#ML_ADM::g31_bound_limit = 0.0 +#$bound$#ML_ADM::g31_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::g32_bound = "skip" +#$bound$#ML_ADM::g32_bound_speed = 1.0 +#$bound$#ML_ADM::g32_bound_limit = 0.0 +#$bound$#ML_ADM::g32_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 +#$bound$#ML_ADM::beta1_bound_scalar = 0.0 + +#$bound$#ML_ADM::beta2_bound = "skip" +#$bound$#ML_ADM::beta2_bound_speed = 1.0 +#$bound$#ML_ADM::beta2_bound_limit = 0.0 +#$bound$#ML_ADM::beta2_bound_scalar = 0.0 + +#$bound$#ML_ADM::beta3_bound = "skip" +#$bound$#ML_ADM::beta3_bound_speed = 1.0 +#$bound$#ML_ADM::beta3_bound_limit = 0.0 +#$bound$#ML_ADM::beta3_bound_scalar = 0.0 + +*/ + |