From 959291633ec89448f4123c5ecbf16db876f29429 Mon Sep 17 00:00:00 2001 From: hinder Date: Mon, 13 Sep 2010 11:18:30 +0000 Subject: Update automatically generated code git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@61 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843 --- src/Boundaries.c | 208 ---------------------------------------------------- src/psis_calc_2nd.c | 15 ++++ src/psis_calc_4th.c | 15 ++++ 3 files changed, 30 insertions(+), 208 deletions(-) (limited to 'src') diff --git a/src/Boundaries.c b/src/Boundaries.c index b6e6cd0..ef1a009 100644 --- a/src/Boundaries.c +++ b/src/Boundaries.c @@ -31,219 +31,11 @@ void WeylScal4_SelectBoundConds(CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS; CCTK_INT ierr = 0; - - if (CCTK_EQUALS(Psi4i_group_bound, "none" ) || - CCTK_EQUALS(Psi4i_group_bound, "static") || - CCTK_EQUALS(Psi4i_group_bound, "flat" ) || - CCTK_EQUALS(Psi4i_group_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4i_group", Psi4i_group_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4i_group_bound BC for WeylScal4::Psi4i_group!"); - } - - if (CCTK_EQUALS(Psi4r_group_bound, "none" ) || - CCTK_EQUALS(Psi4r_group_bound, "static") || - CCTK_EQUALS(Psi4r_group_bound, "flat" ) || - CCTK_EQUALS(Psi4r_group_bound, "zero" ) ) - { - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4r_group", Psi4r_group_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4r_group_bound BC for WeylScal4::Psi4r_group!"); - } - - if (CCTK_EQUALS(Psi4i_bound, "none" ) || - CCTK_EQUALS(Psi4i_bound, "static") || - CCTK_EQUALS(Psi4i_bound, "flat" ) || - CCTK_EQUALS(Psi4i_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4i", Psi4i_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4i_bound BC for WeylScal4::Psi4i!"); - } - - if (CCTK_EQUALS(Psi4r_bound, "none" ) || - CCTK_EQUALS(Psi4r_bound, "static") || - CCTK_EQUALS(Psi4r_bound, "flat" ) || - CCTK_EQUALS(Psi4r_bound, "zero" ) ) - { - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1, - "WeylScal4::Psi4r", Psi4r_bound); - if (ierr < 0) - CCTK_WARN(0, "Failed to register Psi4r_bound BC for WeylScal4::Psi4r!"); - } - - if (CCTK_EQUALS(Psi4i_group_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4i_group_bound = -1; - if (handle_Psi4i_group_bound < 0) handle_Psi4i_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_group_bound , Psi4i_group_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4i_group_bound ,Psi4i_group_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_group_bound, - "WeylScal4::Psi4i_group", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4i_group!"); - - } - - if (CCTK_EQUALS(Psi4r_group_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4r_group_bound = -1; - if (handle_Psi4r_group_bound < 0) handle_Psi4r_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_group_bound , Psi4r_group_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4r_group_bound ,Psi4r_group_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_group_bound, - "WeylScal4::Psi4r_group", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4r_group!"); - - } - - if (CCTK_EQUALS(Psi4i_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4i_bound = -1; - if (handle_Psi4i_bound < 0) handle_Psi4i_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_bound , Psi4i_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4i_bound ,Psi4i_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_bound, - "WeylScal4::Psi4i", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4i!"); - - } - - if (CCTK_EQUALS(Psi4r_bound, "radiative")) - { - /* select radiation boundary condition */ - static CCTK_INT handle_Psi4r_bound = -1; - if (handle_Psi4r_bound < 0) handle_Psi4r_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_bound , Psi4r_bound_limit, "LIMIT") < 0) - CCTK_WARN(0, "could not set LIMIT value in table!"); - if (Util_TableSetReal(handle_Psi4r_bound ,Psi4r_bound_speed, "SPEED") < 0) - CCTK_WARN(0, "could not set SPEED value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_bound, - "WeylScal4::Psi4r", "Radiation"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Radiation BC for WeylScal4::Psi4r!"); - - } - - if (CCTK_EQUALS(Psi4i_group_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4i_group_bound = -1; - if (handle_Psi4i_group_bound < 0) handle_Psi4i_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_group_bound ,Psi4i_group_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_group_bound, - "WeylScal4::Psi4i_group", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for WeylScal4::Psi4i_group!"); - - } - - if (CCTK_EQUALS(Psi4r_group_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4r_group_bound = -1; - if (handle_Psi4r_group_bound < 0) handle_Psi4r_group_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_group_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_group_bound ,Psi4r_group_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_group_bound, - "WeylScal4::Psi4r_group", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Failed to register Scalar BC for WeylScal4::Psi4r_group!"); - - } - - if (CCTK_EQUALS(Psi4i_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4i_bound = -1; - if (handle_Psi4i_bound < 0) handle_Psi4i_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4i_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4i_bound ,Psi4i_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4i_bound, - "WeylScal4::Psi4i", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for WeylScal4::Psi4i!"); - - } - - if (CCTK_EQUALS(Psi4r_bound, "scalar")) - { - /* select scalar boundary condition */ - static CCTK_INT handle_Psi4r_bound = -1; - if (handle_Psi4r_bound < 0) handle_Psi4r_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE); - if (handle_Psi4r_bound < 0) CCTK_WARN(0, "could not create table!"); - if (Util_TableSetReal(handle_Psi4r_bound ,Psi4r_bound_scalar, "SCALAR") < 0) - CCTK_WARN(0, "could not set SCALAR value in table!"); - - ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_Psi4r_bound, - "WeylScal4::Psi4r", "scalar"); - - if (ierr < 0) - CCTK_WARN(0, "Error in registering Scalar BC for WeylScal4::Psi4r!"); - - } return; } /* template for entries in parameter file: -#$bound$#WeylScal4::Psi4i_group_bound = "skip" -#$bound$#WeylScal4::Psi4i_group_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4i_group_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4i_group_bound_scalar = 0.0 - -#$bound$#WeylScal4::Psi4r_group_bound = "skip" -#$bound$#WeylScal4::Psi4r_group_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4r_group_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4r_group_bound_scalar = 0.0 - -#$bound$#WeylScal4::Psi4i_bound = "skip" -#$bound$#WeylScal4::Psi4i_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4i_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4i_bound_scalar = 0.0 - -#$bound$#WeylScal4::Psi4r_bound = "skip" -#$bound$#WeylScal4::Psi4r_bound_speed = 1.0 -#$bound$#WeylScal4::Psi4r_bound_limit = 0.0 -#$bound$#WeylScal4::Psi4r_bound_scalar = 0.0 - */ diff --git a/src/psis_calc_2nd.c b/src/psis_calc_2nd.c index 193c552..ecdc0e0 100644 --- a/src/psis_calc_2nd.c +++ b/src/psis_calc_2nd.c @@ -19,6 +19,21 @@ #define CUB(x) ((x) * (x) * (x)) #define QAD(x) ((x) * (x) * (x) * (x)) +void psis_calc_2nd_SelectBCs(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + CCTK_INT ierr = 0; + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4i_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4i_group."); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4r_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4r_group."); + return; +} + void psis_calc_2nd_Body(cGH const * restrict const cctkGH, int const dir, int const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], int const min[3], int const max[3], int const n_subblock_gfs, CCTK_REAL * restrict const subblock_gfs[]) { DECLARE_CCTK_ARGUMENTS; diff --git a/src/psis_calc_4th.c b/src/psis_calc_4th.c index 32fd7c2..ca493ac 100644 --- a/src/psis_calc_4th.c +++ b/src/psis_calc_4th.c @@ -19,6 +19,21 @@ #define CUB(x) ((x) * (x) * (x)) #define QAD(x) ((x) * (x) * (x) * (x)) +void psis_calc_4th_SelectBCs(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_ARGUMENTS; + DECLARE_CCTK_PARAMETERS; + + CCTK_INT ierr = 0; + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4i_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4i_group."); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "WeylScal4::Psi4r_group","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for WeylScal4::Psi4r_group."); + return; +} + void psis_calc_4th_Body(cGH const * restrict const cctkGH, int const dir, int const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], int const min[3], int const max[3], int const n_subblock_gfs, CCTK_REAL * restrict const subblock_gfs[]) { DECLARE_CCTK_ARGUMENTS; -- cgit v1.2.3