diff options
Diffstat (limited to 'ML_WaveToy/src/Boundaries.c')
| -rw-r--r-- | ML_WaveToy/src/Boundaries.c | 249 |
1 files changed, 0 insertions, 249 deletions
diff --git a/ML_WaveToy/src/Boundaries.c b/ML_WaveToy/src/Boundaries.c deleted file mode 100644 index 995ed75..0000000 --- a/ML_WaveToy/src/Boundaries.c +++ /dev/null @@ -1,249 +0,0 @@ -/* File produced by Kranc */ - -#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_WaveToy_CheckBoundaries(CCTK_ARGUMENTS) -{ - DECLARE_CCTK_ARGUMENTS; - DECLARE_CCTK_PARAMETERS; - - return; -} - -void ML_WaveToy_SelectBoundConds(CCTK_ARGUMENTS) -{ - DECLARE_CCTK_ARGUMENTS; - DECLARE_CCTK_PARAMETERS; - - CCTK_INT ierr = 0; - - if (CCTK_EQUALS(WT_rho_bound, "none" ) || - CCTK_EQUALS(WT_rho_bound, "static") || - CCTK_EQUALS(WT_rho_bound, "flat" ) || - CCTK_EQUALS(WT_rho_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_WaveToy::WT_rho", WT_rho_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register WT_rho_bound BC for ML_WaveToy::WT_rho!"); - } - - if (CCTK_EQUALS(WT_u_bound, "none" ) || - CCTK_EQUALS(WT_u_bound, "static") || - CCTK_EQUALS(WT_u_bound, "flat" ) || - CCTK_EQUALS(WT_u_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_WaveToy::WT_u", WT_u_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register WT_u_bound BC for ML_WaveToy::WT_u!"); - } - - if (CCTK_EQUALS(rho_bound, "none" ) || - CCTK_EQUALS(rho_bound, "static") || - CCTK_EQUALS(rho_bound, "flat" ) || - CCTK_EQUALS(rho_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_WaveToy::rho", rho_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register rho_bound BC for ML_WaveToy::rho!"); - } - - if (CCTK_EQUALS(u_bound, "none" ) || - CCTK_EQUALS(u_bound, "static") || - CCTK_EQUALS(u_bound, "flat" ) || - CCTK_EQUALS(u_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "ML_WaveToy::u", u_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register u_bound BC for ML_WaveToy::u!"); - } - - if (CCTK_EQUALS(WT_rho_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_WT_rho_bound = -1; - if (handle_WT_rho_bound < 0) handle_WT_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_WT_rho_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_WT_rho_bound , WT_rho_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_WT_rho_bound ,WT_rho_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_rho_bound, - "ML_WaveToy::WT_rho", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::WT_rho!"); - - } - - if (CCTK_EQUALS(WT_u_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_WT_u_bound = -1; - if (handle_WT_u_bound < 0) handle_WT_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_WT_u_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_WT_u_bound , WT_u_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_WT_u_bound ,WT_u_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_u_bound, - "ML_WaveToy::WT_u", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::WT_u!"); - - } - - if (CCTK_EQUALS(rho_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_rho_bound = -1; - if (handle_rho_bound < 0) handle_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_rho_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_rho_bound , rho_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_rho_bound ,rho_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_rho_bound, - "ML_WaveToy::rho", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::rho!"); - - } - - if (CCTK_EQUALS(u_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_u_bound = -1; - if (handle_u_bound < 0) handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_u_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_u_bound , u_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_u_bound ,u_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_u_bound, - "ML_WaveToy::u", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for ML_WaveToy::u!"); - - } - - if (CCTK_EQUALS(WT_rho_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_WT_rho_bound = -1; - if (handle_WT_rho_bound < 0) handle_WT_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_WT_rho_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_WT_rho_bound ,WT_rho_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_rho_bound, - "ML_WaveToy::WT_rho", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for ML_WaveToy::WT_rho!"); - - } - - if (CCTK_EQUALS(WT_u_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_WT_u_bound = -1; - if (handle_WT_u_bound < 0) handle_WT_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_WT_u_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_WT_u_bound ,WT_u_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_WT_u_bound, - "ML_WaveToy::WT_u", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for ML_WaveToy::WT_u!"); - - } - - if (CCTK_EQUALS(rho_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_rho_bound = -1; - if (handle_rho_bound < 0) handle_rho_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_rho_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_rho_bound ,rho_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_rho_bound, - "ML_WaveToy::rho", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_WaveToy::rho!"); - - } - - if (CCTK_EQUALS(u_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_u_bound = -1; - if (handle_u_bound < 0) handle_u_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_u_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_u_bound ,u_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_u_bound, - "ML_WaveToy::u", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for ML_WaveToy::u!"); - - } - return; -} - - - -/* template for entries in parameter file: -#$bound$#ML_WaveToy::WT_rho_bound = "skip" -#$bound$#ML_WaveToy::WT_rho_bound_speed = 1.0 -#$bound$#ML_WaveToy::WT_rho_bound_limit = 0.0 -#$bound$#ML_WaveToy::WT_rho_bound_scalar = 0.0 - -#$bound$#ML_WaveToy::WT_u_bound = "skip" -#$bound$#ML_WaveToy::WT_u_bound_speed = 1.0 -#$bound$#ML_WaveToy::WT_u_bound_limit = 0.0 -#$bound$#ML_WaveToy::WT_u_bound_scalar = 0.0 - -#$bound$#ML_WaveToy::rho_bound = "skip" -#$bound$#ML_WaveToy::rho_bound_speed = 1.0 -#$bound$#ML_WaveToy::rho_bound_limit = 0.0 -#$bound$#ML_WaveToy::rho_bound_scalar = 0.0 - -#$bound$#ML_WaveToy::u_bound = "skip" -#$bound$#ML_WaveToy::u_bound_speed = 1.0 -#$bound$#ML_WaveToy::u_bound_limit = 0.0 -#$bound$#ML_WaveToy::u_bound_scalar = 0.0 - -*/ - |