Update automatically generated code

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@61 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843

3 files changed, 30 insertions, 208 deletions

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;