Provide wrappers for all boundary conditions provided by thorn

Boundary, so that they can be 'selected for bc' under the new boundary
interface.  The details of default arguments and table keys will be
documented in the thorn guide.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@198 6a38eb6e-646e-4a02-a296-d141613ad6c4

1 files changed, 128 insertions, 8 deletions

diff --git a/src/FlatBoundary.c b/src/FlatBoundary.c
index 63166ef..08657d5 100644
--- a/src/FlatBoundary.c
+++ b/src/FlatBoundary.c
@@ -19,6 +19,9 @@
 #include <string.h>
 
 #include "cctk.h"
+#include "cctk_Faces.h"
+#include "util_Table.h"
+#include "util_ErrorCodes.h"
 #include "cctk_Parameters.h"
 #include "cctk_FortranString.h"
 
@@ -27,12 +30,136 @@
 
 /* the rcs ID and its dummy function to use it */
 static const char *rcsid = "$Header$";
-CCTK_FILEVERSION(CactusBase_Boundary_FlatBoundary_c)
+CCTK_FILEVERSION(CactusBase_Boundary_FlatBoundary_c);
+
+static int ApplyBndFlat (const cGH *GH,
+                         int stencil_dir,
+                         const int *stencil_alldirs,
+                         int dir,
+                         int first_var,
+                         int num_vars);
 
 
 /********************************************************************
  ********************    External Routines   ************************
  ********************************************************************/
+/*@@
+   @routine    BndFlat
+   @date       13 Feb 2003
+   @author     David Rideout
+   @desc
+               Top level function which is registered as handling 
+	       the Flat boundary condition
+   @enddesc
+   @calls      ApplyBndFlat
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      const cGH *
+   @vio        in
+   @endvar
+   @var        num_vars
+   @vdesc      number of variables passed in through var_indices[]
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        var_indices
+   @vdesc      array of variable indicies to which to apply this boundary 
+               condition
+   @vtype      int *
+   @vio        in
+   @endvar
+   @var        faces
+   @vdesc      array of set of faces to which to apply the bc
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        table_handles
+   @vdesc      array of table handles which hold extra arguments
+   @vtype      int
+   @vio        in
+   @endvar
+   @returntype int
+   @returndesc
+               return code of @seeroutine ApplyBndScalar
+   @endreturndesc
+@@*/
+
+int BndFlat(const cGH *GH, int num_vars, int *vars, int *faces, int *tables) 
+{
+  int i, j, k, gdim, max_gdim, err, retval;
+
+  /* variables to pass to ApplyBndFlat */
+  /*int stencil_dir;*/ /* width of stencil in direction dir, not needed */
+  int *stencil_alldirs; /* width of stencil in all directions */
+  int dir; /* direction in which to apply bc */
+
+  retval = 0; stencil_alldirs = NULL;
+
+  /* loop through variables, j at a time */
+  for (i=0; i<num_vars; i+=j) {
+
+    /* find other adjacent vars which are selected for identical bcs */
+    j=1;
+    while (i+j<num_vars && vars[i+j]==vars[i]+j && tables[i+j]==tables[i] && 
+	   faces[i+j]==faces[i])
+    {
+      ++j;
+    }
+
+    /* Check to see if faces specification is valid */
+    if (faces[i] != CCTK_ALL_FACES)
+    {
+      CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING,
+		 "Faces specification %d for Flat boundary conditions on "
+		 "%s is not implemented yet.  "
+		 "Applying radiative bcs to all (external) faces.", faces[i],
+		 CCTK_VarName(vars[i]));
+    }
+    dir = 0; /* apply bc to all faces */
+
+    /* Set up default arguments for ApplyBndFlat */
+    /* Allocate memory and populate stencil width array */
+    gdim = CCTK_GroupDimFromVarI(vars[i]);
+    if (!stencil_alldirs)
+    {
+      stencil_alldirs = (int *) malloc(gdim*sizeof(int));
+      max_gdim = gdim;
+    } else if (gdim > max_gdim)
+    {
+      realloc(stencil_alldirs, gdim*sizeof(int));
+      max_gdim = gdim;
+    }
+    for (k=0; k<gdim; ++k)
+    {
+      stencil_alldirs[k] = 1;
+    }
+
+    /* Look on table for possible non-default arguments
+     * (If any of these table look-ups fail, the value will be unchanged 
+     * from its default value)
+     */
+    /* Stencil width array */
+    err = Util_TableGetIntArray(tables[i], gdim, stencil_alldirs, 
+				"STENCIL WIDTH");
+    if (err == UTIL_ERROR_BAD_HANDLE)
+    {
+      CCTK_VWarn(3, __LINE__, __FILE__, CCTK_THORNSTRING,
+		 "Invalid table handle passed for Flat boundary "
+		 "conditions for %s.  Using all default values.", 
+		 CCTK_VarName(vars[i]));
+    }
+
+    if ((retval = ApplyBndFlat(GH, 0, stencil_alldirs, dir, vars[i], j)) < 0)
+    {
+      CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING,
+		 "ApplyBndFlat() returned %d", retval);
+    }
+  }
+  free(stencil_alldirs);
+
+  return retval;
+}
+
 /* prototypes for external C routines are declared in header Boundary.h
    here only follow the fortran wrapper prototypes */
 void CCTK_FCALL CCTK_FNAME (BndFlatDirGI)
@@ -84,13 +211,6 @@ void CCTK_FCALL CCTK_FNAME (BndFlatVN)
 /********************************************************************
  ********************    Internal Routines   ************************
  ********************************************************************/
-static int ApplyBndFlat (const cGH *GH,
-                         int stencil_dir,
-                         const int *stencil_alldirs,
-                         int dir,
-                         int first_var,
-                         int num_vars);
-
 
 /*@@
    @routine    BndFlatDirGI