diff options
Diffstat (limited to 'ML_BSSN_O2/src/Boundaries.cc')
-rw-r--r-- | ML_BSSN_O2/src/Boundaries.cc | 104 |
1 files changed, 0 insertions, 104 deletions
diff --git a/ML_BSSN_O2/src/Boundaries.cc b/ML_BSSN_O2/src/Boundaries.cc index 2a6334b..e711921 100644 --- a/ML_BSSN_O2/src/Boundaries.cc +++ b/ML_BSSN_O2/src/Boundaries.cc @@ -120,17 +120,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) CCTK_WARN(0, "Failed to register ML_shift_bound BC for ML_BSSN_O2::ML_shift!"); } - if (CCTK_EQUALS(ML_Theta_bound, "none" ) || - CCTK_EQUALS(ML_Theta_bound, "static") || - CCTK_EQUALS(ML_Theta_bound, "flat" ) || - CCTK_EQUALS(ML_Theta_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_BSSN_O2::ML_Theta", ML_Theta_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register ML_Theta_bound BC for ML_BSSN_O2::ML_Theta!"); - } - if (CCTK_EQUALS(ML_trace_curv_bound, "none" ) || CCTK_EQUALS(ML_trace_curv_bound, "static") || CCTK_EQUALS(ML_trace_curv_bound, "flat" ) || @@ -406,17 +395,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) CCTK_WARN(0, "Failed to register beta3_bound BC for ML_BSSN_O2::beta3!"); } - if (CCTK_EQUALS(Theta_bound, "none" ) || - CCTK_EQUALS(Theta_bound, "static") || - CCTK_EQUALS(Theta_bound, "flat" ) || - CCTK_EQUALS(Theta_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_BSSN_O2::Theta", Theta_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Theta_bound BC for ML_BSSN_O2::Theta!"); - } - if (CCTK_EQUALS(trK_bound, "none" ) || CCTK_EQUALS(trK_bound, "static") || CCTK_EQUALS(trK_bound, "flat" ) || @@ -580,25 +558,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(ML_Theta_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_ML_Theta_bound = -1; - if (handle_ML_Theta_bound < 0) handle_ML_Theta_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_ML_Theta_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_ML_Theta_bound , ML_Theta_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_ML_Theta_bound ,ML_Theta_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_Theta_bound, - "ML_BSSN_O2::ML_Theta", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_BSSN_O2::ML_Theta!"); - - } - if (CCTK_EQUALS(ML_trace_curv_bound, "radiative")) { /* select radiation boundary condition */ @@ -1074,25 +1033,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(Theta_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Theta_bound = -1; - if (handle_Theta_bound < 0) handle_Theta_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Theta_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Theta_bound , Theta_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Theta_bound ,Theta_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Theta_bound, - "ML_BSSN_O2::Theta", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_BSSN_O2::Theta!"); - - } - if (CCTK_EQUALS(trK_bound, "radiative")) { /* select radiation boundary condition */ @@ -1248,23 +1188,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(ML_Theta_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_ML_Theta_bound = -1; - if (handle_ML_Theta_bound < 0) handle_ML_Theta_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_ML_Theta_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_ML_Theta_bound ,ML_Theta_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_Theta_bound, - "ML_BSSN_O2::ML_Theta", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for ML_BSSN_O2::ML_Theta!"); - - } - if (CCTK_EQUALS(ML_trace_curv_bound, "scalar")) { /* select scalar boundary condition */ @@ -1690,23 +1613,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) } - if (CCTK_EQUALS(Theta_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Theta_bound = -1; - if (handle_Theta_bound < 0) handle_Theta_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Theta_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Theta_bound ,Theta_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Theta_bound, - "ML_BSSN_O2::Theta", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_BSSN_O2::Theta!"); - - } - if (CCTK_EQUALS(trK_bound, "scalar")) { /* select scalar boundary condition */ @@ -1769,11 +1675,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) #$bound$#ML_BSSN_O2::ML_shift_bound_limit = 0.0 #$bound$#ML_BSSN_O2::ML_shift_bound_scalar = 0.0 -#$bound$#ML_BSSN_O2::ML_Theta_bound = "skip" -#$bound$#ML_BSSN_O2::ML_Theta_bound_speed = 1.0 -#$bound$#ML_BSSN_O2::ML_Theta_bound_limit = 0.0 -#$bound$#ML_BSSN_O2::ML_Theta_bound_scalar = 0.0 - #$bound$#ML_BSSN_O2::ML_trace_curv_bound = "skip" #$bound$#ML_BSSN_O2::ML_trace_curv_bound_speed = 1.0 #$bound$#ML_BSSN_O2::ML_trace_curv_bound_limit = 0.0 @@ -1899,11 +1800,6 @@ extern "C" void ML_BSSN_O2_SelectBoundConds(CCTK_ARGUMENTS) #$bound$#ML_BSSN_O2::beta3_bound_limit = 0.0 #$bound$#ML_BSSN_O2::beta3_bound_scalar = 0.0 -#$bound$#ML_BSSN_O2::Theta_bound = "skip" -#$bound$#ML_BSSN_O2::Theta_bound_speed = 1.0 -#$bound$#ML_BSSN_O2::Theta_bound_limit = 0.0 -#$bound$#ML_BSSN_O2::Theta_bound_scalar = 0.0 - #$bound$#ML_BSSN_O2::trK_bound = "skip" #$bound$#ML_BSSN_O2::trK_bound_speed = 1.0 #$bound$#ML_BSSN_O2::trK_bound_limit = 0.0 |