Generalized boundary condition routines for applying to

arbitrary CCTK data types (except CCTK_COMPLEX).

Added/completed grdoc.


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

1 files changed, 799 insertions, 771 deletions

diff --git a/src/CopyBoundary.c b/src/CopyBoundary.c
index 117b742..37a59ad 100644
--- a/src/CopyBoundary.c
+++ b/src/CopyBoundary.c
@@ -2,956 +2,984 @@
    @file      CopyBoundary.c
    @date      Mon Mar 15 15:09:00 1999
    @author    Gerd Lanfermann, Gabrielle Allen
-   @desc 
-     Wrappers for calling routines for applying copying-boudnary
-     conditions
-   @enddesc 
+   @desc
+              Routines for applying copying-boundary conditions
+   @enddesc
+   @history
+   @hdate     Sun 25 Feb 2001
+   @hauthor   Thomas Radke
+   @hdesc     BC routines generalized for applying to arbitrary CCTK data types
+   @endhistory
+   @version   $Id$
  @@*/
 
-#include <stdio.h>
-#include <assert.h>
 #include <stdlib.h>
-#include <ctype.h>
-#include <stdarg.h>
 #include <string.h>
 
 #include "cctk.h"
 #include "cctk_FortranString.h"
-#include "cctk_Parameters.h"
 
-#include "BoundarySymmetries.h"
+#include "Symmetry.h"
 #include "Boundary.h"
 
-/* Internal routine prototypes */
+/* the rcs ID and its dummy function to use it */
+static char *rcsid = "$Header$";
+CCTK_FILEVERSION(CactusBase_Boundary_CopyBoundary_c)
 
-static int BndApplyCopy3Di(cGH *GH,
-			   int doBC,
-			   int *lssh,
-			   int *stencil_size,
-			   CCTK_REAL *var,
-			   CCTK_REAL *varp);
-
-static int BndApplyCopy2Di(cGH *GH,
-			   int doBC,
-			   int *lssh,
-			   int *stencil_size,
-			   CCTK_REAL *var,
-			   CCTK_REAL *varp);
-
-static int BndApplyCopy1Di(cGH *GH,
-			   int doBC,
-			   int *lssh,
-			   int *stencil_size,
-			   CCTK_REAL *var,
-			   CCTK_REAL *varp);
-
-static int ApplyBndCopy(cGH *GH, 
-			int stencil,
-			int *stencil_array, 
-			int first_var,
-			int second_var,
-			int num_vars,
-			int dir);
-
-/* Local variables */
 
+/********************************************************************
+ ********************    External Routines   ************************
+ ********************************************************************/
+/* prototypes for external C routines are declared in header Boundary.h
+   here only follow the fortran wrapper prototypes */
+void CCTK_FCALL CCTK_FNAME (BndCopyDirVI)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            const int *vi_to,
+                            const int *vi_from);
+void CCTK_FCALL CCTK_FNAME (BndCopyVI)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            const int *vi_to,
+                            const int *vi_from);
+void CCTK_FCALL CCTK_FNAME (BndCopyDirGI)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            const int *gi_to,
+                            const int *gi_from);
+void CCTK_FCALL CCTK_FNAME (BndCopyGI)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            const int *gi_to,
+                            const int *gi_from);
+void CCTK_FCALL CCTK_FNAME (BndCopyDirVN)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            TWO_FORTSTRINGS_ARGS);
+void CCTK_FCALL CCTK_FNAME (BndCopyVN)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            TWO_FORTSTRINGS_ARGS);
+void CCTK_FCALL CCTK_FNAME (BndCopyDirGN)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            TWO_FORTSTRINGS_ARGS);
+void CCTK_FCALL CCTK_FNAME (BndCopyGN)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            TWO_FORTSTRINGS_ARGS);
 
 
 /********************************************************************
- ********************    External Routines   ************************
+ ********************    Internal Routines   ************************
  ********************************************************************/
+static int ApplyBndCopy (cGH *GH,
+                         int stencil_dir,
+                         int stencil_alldirs[],
+                         int dir,
+                         int first_var_to,
+                         int first_var_from,
+                         int num_vars);
+
 
 /*@@
    @routine    BndCopyDirVI
-   @date       Sat Mar 20
+   @date       Sat Jan 20 2001
    @author     Gabrielle Allen
-   @desc 
-   Apply copy boundary routines by var index in given direction
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundary routines by var index in given direction
+   @enddesc
+   @calls      ApplyBndCopy
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil_size
+   @vdesc      stencil size in this direction
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        dir
+   @vdesc      direction to copy
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        vi_to
+   @vdesc      index of variable to copy boundaries to
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        vi_from
+   @vdesc      index of variable to copy boundaries from
+   @vtype      int
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine ApplyBndCopy <BR>
+               -1 if invalid variable indices are given
+   @endreturndesc
 @@*/
-
-int BndCopyDirVI(cGH *GH, 
-		 int stencil_size,
-		 int dir,
-		 int vi_to, 
-		 int vi_from) 
+int BndCopyDirVI (cGH *GH,
+                  int stencil_size,
+                  int dir,
+                  int vi_to,
+                  int vi_from)
 {
-  return ApplyBndCopy(GH,stencil_size,NULL,vi_to,vi_from,1,dir);
+  int retval, num_vars;
+
+
+  num_vars = CCTK_NumVars ();
+  if (vi_to >= 0 && vi_to < num_vars && vi_from >= 0 && vi_from < num_vars)
+  {
+    retval = ApplyBndCopy (GH, stencil_size, NULL, dir, vi_to, vi_from, 1);
+  }
+  else
+  {
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid variable indices %d and/or %d in BndCopyDirVI",
+                vi_to, vi_from);
+    retval = -1;
+  }
+
+  return (retval);
 }
 
-void CCTK_FCALL CCTK_FNAME(BndCopyDirVI)
-     (int *ierr, 
-      cGH *GH, 
-      int *stencil_size, 
-      int *dir,
-      int *vi_to, 
-      int *vi_from)
+void CCTK_FCALL CCTK_FNAME (BndCopyDirVI)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            const int *vi_to,
+                            const int *vi_from)
 {
-  *ierr = BndCopyDirVI(GH, *stencil_size, *dir, *vi_to, *vi_from);
+  *ierr = BndCopyDirVI (GH, *stencil_size, *dir, *vi_to, *vi_from);
 }
 
 
-
 /*@@
    @routine    BndCopyVI
    @date       Thu Mar  2 11:02:10 2000
    @author     Gerd Lanfermann
-   @desc 
-   Apply copy boundary routines by var index
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundary routines by var index
+   @enddesc
+   @calls      ApplyBndCopy
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil
+   @vdesc      stencil width
+   @vtype      int [ dimension of variable ]
+   @vio        in
+   @endvar
+   @var        vi_to
+   @vdesc      index of variable to copy boundaries to
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        vi_from
+   @vdesc      index of variable to copy boundaries from
+   @vtype      int
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine ApplyBndCopy <BR>
+               -1 if invalid variable indices are given
+   @endreturndesc
 @@*/
+int BndCopyVI (cGH *GH,
+               int stencil[],
+               int vi_to,
+               int vi_from)
+{
+  int retval, num_vars;
+
+
+  num_vars = CCTK_NumVars ();
+  if (vi_to >= 0 && vi_to < num_vars && vi_from >= 0 && vi_from < num_vars)
+  {
+    retval = ApplyBndCopy (GH, -1, stencil, 0, vi_to, vi_from, 1);
+  }
+  else
+  {
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid variable indices %d and/or %d in BndCopyVI",
+                vi_to, vi_from);
+    retval = -1;
+  }
+
+  return (retval);
+}
 
-int BndCopyVI(cGH *GH, 
-	      int *stencil, 
-	      int vi_to, 
-	      int vi_from) 
+void CCTK_FCALL CCTK_FNAME (BndCopyVI)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            const int *vi_to,
+                            const int *vi_from)
 {
-  return ApplyBndCopy(GH,-1,stencil,vi_to,vi_from,1,0);
+  *ierr = BndCopyVI (GH, stencil, *vi_to, *vi_from);
 }
 
+
 /* ====================================================== */
 
  /*@@
    @routine    BndCopyDirGI
    @date       Sat Jan 20 2001
    @author     Gabrielle Allen
-   @desc 
-   Apply copy boundaries by group index in given direction
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundaries by group index in given direction
+   @enddesc
+   @calls      ApplyBndCopy
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil_size
+   @vdesc      stencil size in this direction
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        dir
+   @vdesc      direction to copy boundaries
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        gi_to
+   @vdesc      index of group to copy boundaries to
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        gi_from
+   @vdesc      index of group to copy boundaries from
+   @vtype      int
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine ApplyBndCopy <BR>
+               -1 if invalid group indices are given
+   @endreturndesc
 @@*/
-
-int BndCopyDirGI(cGH *GH, 
-		 int stencil_size,
-		 int dir,
-		 int gi_to, 
-		 int gi_from) 
+int BndCopyDirGI (cGH *GH,
+                  int stencil_size,
+                  int dir,
+                  int gi_to,
+                  int gi_from)
 {
-  int first_vi_to,first_vi_from;
-  int numvars;
+  int first_vi_to, first_vi_from, retval;
 
-  first_vi_to  = CCTK_FirstVarIndexI(gi_to);
-  first_vi_from= CCTK_FirstVarIndexI(gi_from);
-  numvars      = CCTK_NumVarsInGroupI(gi_to);
 
-  return ApplyBndCopy(GH,stencil_size,NULL,first_vi_to,first_vi_from,numvars,dir);
+  first_vi_to   = CCTK_FirstVarIndexI (gi_to);
+  first_vi_from = CCTK_FirstVarIndexI (gi_from);
+  if (first_vi_to >= 0 && first_vi_from >= 0)
+  {
+    retval = ApplyBndCopy (GH, stencil_size, NULL, dir, first_vi_to,
+                           first_vi_from, CCTK_NumVarsInGroupI (gi_to));
+  }
+  else
+  {
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid group indices %d and/or %d in BndCopyDirGI",
+                gi_to, gi_from);
+    retval = -1;
+  }
+
+  return (retval);
 }
 
-void CCTK_FCALL CCTK_FNAME(BndCopyDirGI)
-     (int *ierr, 
-      cGH *GH, 
-      int *stencil_size,
-      int *dir,
-      int *gi_to, 
-      int *gi_from)
+void CCTK_FCALL CCTK_FNAME (BndCopyDirGI)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            const int *gi_to,
+                            const int *gi_from)
 {
-  *ierr = BndCopyDirGI(GH, *stencil_size, *dir, *gi_to, *gi_from);
-  return;
+  *ierr = BndCopyDirGI (GH, *stencil_size, *dir, *gi_to, *gi_from);
 }
 
+
  /*@@
    @routine    BndCopyGI
    @date       Thu Mar  2 11:07:11 2000
    @author     Gerd Lanfermann
-   @desc 
-   Apply copy boundaries by group index
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundaries by group index
+   @enddesc
+   @calls      ApplyBndCopy
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil
+   @vdesc      stencil width
+   @vtype      int [ dimension of group ]
+   @vio        in
+   @endvar
+   @var        gi_to
+   @vdesc      index of group to copy boundaries to
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        gi_from
+   @vdesc      index of group to copy boundaries from
+   @vtype      int
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine ApplyBndCopy <BR>
+               -1 if invalid group indices are given
+   @endreturndesc
 @@*/
-
-int BndCopyGI(cGH *GH, 
-	      int *stencil, 
-	      int gi_to, 
-	      int gi_from) 
+int BndCopyGI (cGH *GH,
+               int stencil[],
+               int gi_to,
+               int gi_from)
 {
-  int first_vi_to,first_vi_from;
-  int numvars;
+  int first_vi_to, first_vi_from, retval;
+
 
-  first_vi_to  = CCTK_FirstVarIndexI(gi_to);
-  first_vi_from= CCTK_FirstVarIndexI(gi_from);
-  numvars      = CCTK_NumVarsInGroupI(gi_to);
+  first_vi_to   = CCTK_FirstVarIndexI (gi_to);
+  first_vi_from = CCTK_FirstVarIndexI (gi_from);
+  if (first_vi_to >= 0 && first_vi_from >= 0)
+  {
+    retval = ApplyBndCopy (GH, -1, stencil, 0, first_vi_to, first_vi_from,
+                           CCTK_NumVarsInGroupI (gi_to));
+  }
+  else
+  {
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid group indices %d and/or %d in BndCopyGI",
+                gi_to, gi_from);
+    retval = -1;
+  }
 
-  return ApplyBndCopy(GH,-1,stencil,first_vi_to,first_vi_from,numvars,0);
+  return (retval);
 }
 
-void CCTK_FCALL CCTK_FNAME(BndCopyGI)
-     (int *ierr, 
-      cGH *GH, 
-      int *stencil, 
-      int *gi_to, 
-      int *gi_from)
+void CCTK_FCALL CCTK_FNAME (BndCopyGI)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            const int *gi_to,
+                            const int *gi_from)
 {
-  *ierr = BndCopyGI(GH, stencil, *gi_to, *gi_from);
+  *ierr = BndCopyGI (GH, stencil, *gi_to, *gi_from);
 }
 
+
 /* ======================================================= */
 
 /*@@
    @routine    BndCopyDirGN
    @date       Sat Jan 20 2001
    @author     Gabrielle Allen
-   @desc 
-   Apply copy boundary routines by group name in given direction
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundary routines by group name in given direction
+   @enddesc
+   @calls      BndCopyDirGI
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil_size
+   @vdesc      stencil size in this direction
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        dir
+   @vdesc      direction to copy boundaries
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        gn_to
+   @vdesc      name of group to copy boundaries to
+   @vtype      char []
+   @vio        in
+   @endvar
+   @var        gn_from
+   @vdesc      name of group to copy boundaries from
+   @vtype      char []
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine BndCopyDirGI <BR>
+               -1 if invalid group names are given
+   @endreturndesc
 @@*/
-
-int BndCopyDirGN(cGH *GH, 
-		 int stencil_size,
-		 int dir,
-		 const char *gn_to, 
-		 const char *gn_from)
+int BndCopyDirGN (cGH *GH,
+                  int stencil_size,
+                  int dir,
+                  char gn_to[],
+                  char gn_from[])
 {
-  int retval;
-  int gi_to, gi_from;
-
-  gi_to   = CCTK_GroupIndex(gn_to);
-  gi_from = CCTK_GroupIndex(gn_from);
-  
-  if ((gi_to>-1)&&(gi_from>-1))
-    retval = BndCopyDirGI(GH, stencil_size, dir, gi_to, gi_from);
-  else 
+  int gi_to, gi_from, num_groups, retval;
+
+
+  gi_to      = CCTK_GroupIndex (gn_to);
+  gi_from    = CCTK_GroupIndex (gn_from);
+  num_groups = CCTK_NumGroups ();
+
+  if (gi_to >= 0 && gi_to < num_groups && gi_from >= 0 && gi_from < num_groups)
+  {
+    retval = BndCopyDirGI (GH, stencil_size, dir, gi_to, gi_from);
+  }
+  else
   {
-    CCTK_VWarn(2,__LINE__,__FILE__,"Boundary",
-	       "BndCopyDirGN: Cannot find group indices for %s, %s ",
-	       gn_to,gn_from);
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid group names '%s' and/or '%s' in BndCopyDirGN",
+                gn_to, gn_from);
     retval = -1;
   }
 
-  return retval;
-
+  return (retval);
 }
 
-void CCTK_FCALL CCTK_FNAME(BndCopyDirGN)
-     (int *ierr, 
-      cGH *GH, 
-      int *stencil_size,
-      int *dir,
-      TWO_FORTSTRINGS_ARGS)
+void CCTK_FCALL CCTK_FNAME (BndCopyDirGN)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            TWO_FORTSTRINGS_ARGS)
 {
-  TWO_FORTSTRINGS_CREATE(gn_to, gn_from)
-  *ierr = BndCopyDirGN(GH, *stencil_size, *dir, gn_to, gn_from);
-  free(gn_to);
-  free(gn_from);
-  return;
+  TWO_FORTSTRINGS_CREATE (gn_to, gn_from)
+  *ierr = BndCopyDirGN (GH, *stencil_size, *dir, gn_to, gn_from);
+  free (gn_to);
+  free (gn_from);
 }
 
+
 /*@@
    @routine    BndCopyGN
    @date       Thu Mar  2 11:02:10 2000
    @author     Gerd Lanfermann
-   @desc 
-   Apply copy boundary routines by group name
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundary routines by group name
+   @enddesc
+   @calls      BndCopyGI
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil
+   @vdesc      stencil width
+   @vtype      int [ dimension of group ]
+   @vio        in
+   @endvar
+   @var        gn_to
+   @vdesc      name of group to copy boundaries to
+   @vtype      char []
+   @vio        in
+   @endvar
+   @var        gn_from
+   @vdesc      name of group to copy boundaries from
+   @vtype      char []
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine BndCopyGI <BR>
+               -1 if invalid group names are given
+   @endreturndesc
 @@*/
-
-int BndCopyGN(cGH *GH, 
-	      int *stencil, 
-	      const char *gn_to, 
-	      const char *gn_from)
+int BndCopyGN (cGH *GH,
+               int stencil[],
+               char gn_to[],
+               char gn_from[])
 {
-  int retval;
-  int gi_to, gi_from;
-
-  gi_to   = CCTK_GroupIndex(gn_to);
-  gi_from = CCTK_GroupIndex(gn_from);
-  
-  if ((gi_to>-1)&&(gi_from>-1))
-    retval = BndCopyGI(GH, stencil, gi_to, gi_from);
-  else 
+  int gi_to, gi_from, num_groups, retval;
+
+
+  gi_to      = CCTK_GroupIndex (gn_to);
+  gi_from    = CCTK_GroupIndex (gn_from);
+  num_groups = CCTK_NumGroups ();
+
+  if (gi_to >= 0 && gi_to < num_groups && gi_from >= 0 && gi_from < num_groups)
+  {
+    retval = BndCopyGI (GH, stencil, gi_to, gi_from);
+  }
+  else
   {
-    CCTK_VWarn(2,__LINE__,__FILE__,"Boundary",
-	       "BndCopyGN: Cannot find group indices for %s, %s ",
-	       gn_to,gn_from);
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid group names '%s' and/or '%s' in BndCopyGN",
+                gn_to, gn_from);
     retval = -1;
   }
 
-  return retval;
-
+  return (retval);
 }
 
-void CCTK_FCALL CCTK_FNAME(BndCopyGN)
-     (int *ierr, 
-      cGH *GH, 
-      int *stencil, 
-      TWO_FORTSTRINGS_ARGS)
+void CCTK_FCALL CCTK_FNAME (BndCopyGN)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            TWO_FORTSTRINGS_ARGS)
 {
-  TWO_FORTSTRINGS_CREATE(gn_to, gn_from)
-  *ierr = BndCopyGN(GH, stencil, gn_to, gn_from);
-  free(gn_to);
-  free(gn_from);
-  return;
+  TWO_FORTSTRINGS_CREATE (gn_to, gn_from)
+  *ierr = BndCopyGN (GH, stencil, gn_to, gn_from);
+  free (gn_to);
+  free (gn_from);
 }
 
-/* ================================================================ */
+
+/* ======================================================= */
 
 /*@@
    @routine    BndCopyDirVN
    @date       Sat Jan 20 2001
    @author     Gabrielle Allen
-   @desc 
-   Apply copy boundary routines by variable name in given direction
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundary routines by group name in given direction
+   @enddesc
+   @calls      BndCopyDirVI
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil_size
+   @vdesc      stencil size in this direction
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        dir
+   @vdesc      direction to copy boundaries
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        vn_to
+   @vdesc      name of variable to copy boundaries to
+   @vtype      char []
+   @vio        in
+   @endvar
+   @var        vn_from
+   @vdesc      name of variable to copy boundaries from
+   @vtype      char []
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine BndCopyDirVI <BR>
+               -1 if invalid variable names are given
+   @endreturndesc
 @@*/
-
-int BndCopyDirVN(cGH *GH, 
-		 int stencil_size,
-		 int dir,
-		 const char *vn_to, 
-		 const char *vn_from)
+int BndCopyDirVN (cGH *GH,
+                  int stencil_size,
+                  int dir,
+                  char vn_to[],
+                  char vn_from[])
 {
-  int vi_to, vi_from;
-  int retval;
-
-  vi_to   = CCTK_VarIndex(vn_to);
-  vi_from = CCTK_VarIndex(vn_from);
-  
-  if ((vi_to>-1)&&(vi_from>-1))
-    retval = BndCopyDirVI(GH, stencil_size, dir, vi_to, vi_from);
-  else 
+  int vi_to, vi_from, num_vars, retval;
+
+
+  vi_to    = CCTK_VarIndex (vn_to);
+  vi_from  = CCTK_VarIndex (vn_from);
+  num_vars = CCTK_NumVars ();
+
+  if (vi_to >= 0 && vi_to < num_vars && vi_from >= 0 && vi_from < num_vars)
+  {
+    retval = BndCopyDirVI (GH, stencil_size, dir, vi_to, vi_from);
+  }
+  else
   {
-    CCTK_VWarn(2,__LINE__,__FILE__,"Boundary",
-	       "BndCopyDirVN: Cannot find variable indices for %s, %s",
-	       vn_to,vn_from);
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid variable names '%s' and/or '%s' in BndCopyDirVN",
+                vn_to, vn_from);
     retval = -1;
   }
 
-  return retval;
+  return (retval);
 }
 
-void CCTK_FCALL CCTK_FNAME(BndCopyDirVN)
-     (int *ierr, 
-      cGH *GH, 
-      int *stencil_size, 
-      int *dir,
-      TWO_FORTSTRINGS_ARGS)
+void CCTK_FCALL CCTK_FNAME (BndCopyDirVN)
+                           (int *ierr,
+                            cGH *GH,
+                            const int *stencil_size,
+                            const int *dir,
+                            TWO_FORTSTRINGS_ARGS)
 {
-  TWO_FORTSTRINGS_CREATE(vn_to, vn_from)
-  *ierr = BndCopyDirVN(GH, *stencil_size, *dir, vn_to, vn_from);
-  free(vn_to);
-  free(vn_from);
-  return;
+  TWO_FORTSTRINGS_CREATE (vn_to, vn_from)
+  *ierr = BndCopyDirVN (GH, *stencil_size, *dir, vn_to, vn_from);
+  free (vn_to);
+  free (vn_from);
 }
 
+
 /*@@
    @routine    BndCopyVN
    @date       Thu Mar  2 11:02:10 2000
    @author     Gerd Lanfermann
-   @desc 
-   Apply copy boundary routines by variable name
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundary routines by variable name
+   @enddesc
+   @calls      BndCopyVI
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil
+   @vdesc      stencil width
+   @vtype      int [ dimension of variable ]
+   @vio        in
+   @endvar
+   @var        vn_to
+   @vdesc      name of variable to copy boundaries to
+   @vtype      char []
+   @vio        in
+   @endvar
+   @var        vn_from
+   @vdesc      name of variable to copy boundaries from
+   @vtype      char []
+   @vio        in
+   @endvar
+
+   @returntype int
+   @returndesc
+               return code of @seeroutine BndCopyVI <BR>
+               -1 if invalid variable names are given
+   @endreturndesc
 @@*/
-
-int BndCopyVN(cGH *GH, 
-	      int *stencil, 
-	      const char *vn_to, 
-	      const char *vn_from)
+int BndCopyVN (cGH *GH,
+               int stencil[],
+               char vn_to[],
+               char vn_from[])
 {
-  int vi_to, vi_from;
-  int retval;
-
-  vi_to   = CCTK_VarIndex(vn_to);
-  vi_from = CCTK_VarIndex(vn_from);
-  
-  if ((vi_to>-1)&&(vi_from>-1))
-    retval = BndCopyVI(GH, stencil, vi_to, vi_from);
-  else 
+  int vi_to, vi_from, num_vars, retval;
+
+
+  vi_to    = CCTK_VarIndex (vn_to);
+  vi_from  = CCTK_VarIndex (vn_from);
+  num_vars = CCTK_NumVars ();
+
+  if (vi_to >= 0 && vi_to < num_vars && vi_from >= 0 && vi_from < num_vars)
+  {
+    retval = BndCopyVI (GH, stencil, vi_to, vi_from);
+  }
+  else
   {
-    CCTK_VWarn(2,__LINE__,__FILE__,"Boundary",
-	       "BndCopyVN: Cannot find variable indices for %s, %s",
-	       vn_to,vn_from);
+    CCTK_VWarn (2, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Invalid variable names '%s' and/or '%s' in BndCopyVN",
+                vn_to, vn_from);
     retval = -1;
   }
 
-  return retval;
+  return (retval);
 }
 
-void CCTK_FCALL CCTK_FNAME(BndCopyVN)
-     (int *ierr, 
-      cGH *GH, 
-      int *stencil, 
-      TWO_FORTSTRINGS_ARGS)
+void CCTK_FCALL CCTK_FNAME (BndCopyVN)
+                           (int *ierr,
+                            cGH *GH,
+                            int stencil[],
+                            TWO_FORTSTRINGS_ARGS)
 {
-  TWO_FORTSTRINGS_CREATE(vn_to, vn_from)
-  *ierr = BndCopyVN(GH, stencil, vn_to, vn_from);
-  free(vn_to);
-  free(vn_from);
-  return;
+  TWO_FORTSTRINGS_CREATE (vn_to, vn_from)
+  *ierr = BndCopyVN (GH, stencil, vn_to, vn_from);
+  free (vn_to);
+  free (vn_from);
 }
 
 
-
 /********************************************************************
  *********************     Local Routines   *************************
  ********************************************************************/
 
+/* maximum dimension we can deal with */
+#define MAXDIM  3
+
+/* macro to compute the linear index of a 3D point */
+#define INDEX_3D(lsh, i, j, k)      ((i) + (lsh)[0]*((j) + (lsh)[1]*(k)))
+
+/*@@
+   @routine    COPY_BOUNDARY
+   @date       Sat 20 Jan 2001
+   @author     Thomas Radke
+   @desc
+               Macro to apply copy boundary conditions to a variable
+               Currently it is limited up to 3D variables only.
+   @enddesc
+
+   @var        doBC
+   @vdesc      flag telling whether to apply boundary conditions or not
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        iend, jend, kend
+   @vdesc      upper ranges for the loopers
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        ii, jj, kk
+   @vdesc      indices of the current grid point
+   @vtype      int
+   @vio        in
+   @endvar
+@@*/
+#define COPY_BOUNDARY(doBC,                                                   \
+                      iend, jend, kend,                                       \
+                      ii, jj, kk)                                             \
+{                                                                             \
+  if (doBC)                                                                   \
+  {                                                                           \
+    for (k = 0; k < kend; k++)                                                \
+    {                                                                         \
+      for (j = 0; j < jend; j++)                                              \
+      {                                                                       \
+        for (i = 0; i < iend; i++)                                            \
+        {                                                                     \
+          int _index;                                                         \
+                                                                              \
+                                                                              \
+          _index = INDEX_3D (lsh, ii, jj, kk) * vtypesize;                    \
+          memcpy ((char *) GH->data[var_to][timelvl_to] + _index,             \
+                  (char *) GH->data[var_from][timelvl_from] + _index,         \
+                  vtypesize);                                                 \
+        }                                                                     \
+      }                                                                       \
+    }                                                                         \
+  }                                                                           \
+}
+
+
 /*@@
    @routine    ApplyBndCopy
    @date       Thu Mar  2 11:02:10 2000
    @author     Gerd Lanfermann
-   @desc 
-   Apply copy boundary routines by group index
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
+   @desc
+               Apply copy boundary conditions to a group of grid functions
+               given by their indices
+               This routine is called by the various BndCopyXXX wrappers.
+
+               Although it is currently limited to handle 1D, 2D, or 3D
+               variables only it can easily be extended for higher dimensions
+               by adapting the appropriate macros.
+   @enddesc
+
+   @var        GH
+   @vdesc      Pointer to CCTK grid hierarchy
+   @vtype      cGH *
+   @vio        in
+   @endvar
+   @var        stencil_dir
+   @vdesc      stencil width in direction dir
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        stencil_alldirs
+   @vdesc      stencil widths for all directions
+   @vtype      int [ dimension of variable(s) ]
+   @vio        in
+   @endvar
+   @var        dir
+   @vdesc      direction to copy boundaries (0 for copying all directions)
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        first_var_to
+   @vdesc      index of first variable to copy boundaries to
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        first_var_from
+   @vdesc      index of first variable to copy boundaries from
+   @vtype      int
+   @vio        in
+   @endvar
+   @var        num_vars
+   @vdesc      number of variables
+   @vtype      int
+   @vio        in
+   @endvar
+
+   @calls      CCTK_GroupIndexFromVarI
+               CCTK_GroupDimI
+               CCTK_VarTypeI
+               CCTK_GroupStaggerDirArrayGI
+               CCTK_NumTimeLevelsFromVarI
+               COPY_BOUNDARY
+   @history
+   @hdate      Sat 20 Jan 2001
+   @hauthor    Thomas Radke
+   @hdesc      Merged separate routines for 1D, 2D, and 3D
+               into a single generic routine
+   @endhistory
+
+   @returntype int
+   @returndesc
+                0 for success
+               -1 if dimension is not supported
+               -2 if direction parameter is invalid
+               -3 if stencil width array parameter is NULL
+   @endreturndesc
 @@*/
-
-static int ApplyBndCopy(cGH *GH, 
-			int stencil,
-			int *stencil_array,
-			int first_var,
-			int second_var,
-			int num_vars,
-			int dir)
+static int ApplyBndCopy (cGH *GH,
+                         int stencil_dir,
+                         int stencil_alldirs[],
+                         int dir,
+                         int first_var_to,
+                         int first_var_from,
+                         int num_vars)
 {
-  int symmetry_handle;  /* handle for the optional symmetry structure */
-  int vi, vi2, gi, dim;
-  int idim;
-  int retval=0;
-  int dirloop1,dirloop2;
-  int doBC,count;
-  int berr,ierr;
-  int *dstag;
-  int *lssh;
-  int *sw;
-  int timelevel;        /* timelevel that condition applied on */
+  int i, j, k;
+  int timelvl_to, timelvl_from;
+  int gindex, gdim;
+  int var_to, var_from, vtypesize;
+  int doBC[2*MAXDIM], dstag[MAXDIM], lsh[MAXDIM], lssh[MAXDIM], stencil[MAXDIM];
   SymmetryGHex *sGHex;
 
-#ifdef DEBUG_BOUNDARY
-  printf("Input arguments for ApplyBndCopy:\n");
-  printf("GH = %x\n",GH);
-  printf("stencil = %d\n",stencil);
-  printf("stencil_array = %x\n",stencil_array);
-  printf("first_var = %d\n",first_var);
-  printf("second_var = %d\n",second_var);
-  printf("num_vars = %d\n",num_vars); 
-  printf("dir = %d\n",dir);
-  printf("-----------------------------------\n");
-#endif
-
-  /* Get group dimension */
-  dim = CCTK_GroupDimFromVarI(first_var);
-
-  if (dim<1 || dim>3)
-  {
-    CCTK_VWarn(1,__LINE__,__FILE__,CCTK_THORNSTRING,
-	       "ApplyBndCopy: Invalid grid function dimension %d",
-	       dim);
-    retval=-1;
-  }
-  else
-  {
-
-    /* See if we have a symmetry array */
-    symmetry_handle = CCTK_GHExtensionHandle("Symmetry");
-    if (symmetry_handle < 0)
-    {
-      sGHex = NULL;
-    }
-    else
-    {
-      sGHex  = (SymmetryGHex*)GH->extensions[symmetry_handle];
-    }
-
-    /* Set up stencil width array if needed */
-    if (stencil_array && dir == 0)
-    {
-      sw = stencil_array;
-    }
-    else if (dir==0)
-    {
-      CCTK_WARN(0,"ApplyBndCopy: Direction 0 invalid for directional BC");
-    }
-    else if (dir>dim||-dir>dim)
-    {
-      CCTK_VWarn(0,__LINE__,__FILE__,CCTK_THORNSTRING,
-		 "ApplyBndCopy: Direction %d greater than dimension %d",
-		 dir,dim);
-    }
-    else
-    {
-      int i;
-      /* Should be single direction */
-      sw  = (int *)malloc(dim*sizeof(int));
-      if (sw)
-      {
-	for (i=0;i<dim;i++)
-	{
-	  sw[i] = 0;
-	}
-	if (dir > 0)
-	{
-	  sw[dir-1] = stencil;
-	}
-	else
-	{
-	  sw[-dir-1] = stencil;
-	}
-      }
-      else
-      {
-	CCTK_WARN(0,"ApplyBndCopy: Malloc sw failed");
-      }
-    }
-
-    /* allocate arrays */
-    dstag  = (int *)malloc(dim*sizeof(int));
-    lssh   = (int *)malloc(dim*sizeof(int));
-
-    /* get the directional staggering of the group */
-    gi   = CCTK_GroupIndexFromVarI(first_var);
-    ierr = CCTK_GroupStaggerDirArrayGI(dstag, dim, gi);
-
-    /* get the current timelevel */
-    timelevel = CCTK_NumTimeLevelsFromVarI (first_var) - 1;
-
-    /* Treat all boundaries or just one? */
-    if (dir>0)
-    {
-      dirloop1 = 2*(dir-1)+1;
-      dirloop2 = dirloop1+1;
-    }
-    else if (dir<0)
-    {
-      dirloop1 = 2*(-dir-1);
-      dirloop2 = dirloop1+1;
-    }
-    else
-    {
-      dirloop1 = 0;
-      dirloop2 = 2*dim;
-    }
-  
-    for (vi=first_var; vi<first_var+num_vars; vi++) 
-    {
-      vi2 = (vi-first_var+second_var);
-
-      /* Apply condition if:
-	 + boundary is not a symmetry boundary 
-           (no symmetry or unset(=unsed))
-	 + boundary is a physical boundary
-	 + have enough grid points             */
-
-
-      for (idim=0;idim<dim;idim++)
-      {
-	lssh[idim] = (int)GH->cctk_lssh[CCTK_LSSH_IDX(dstag[idim],idim)];
-      }
-
-      for (count=dirloop1;count<dirloop2;count++)
-      {
-	if (sGHex)
-	{
-	  doBC = (
-		  ( (sGHex->GFSym[vi][count]==GFSYM_NOSYM) ||
-		    (sGHex->GFSym[vi][count]==GFSYM_UNSET) ) &&
-		  GH->cctk_lsh[count/2]>1 && GH->cctk_bbox[count]) ?
-	    count : -1;
-	}
-	else
-	{
-	  doBC  = (GH->cctk_lsh[count/2]>1 && GH->cctk_bbox[count]) ?
-	    count : -1 ;
-	}
-
-
-	if (doBC >=0)
-	{
-	  switch (dim)
-	  {
-	  case 1: 
-	    berr = BndApplyCopy1Di(GH, 
-				     doBC, 
-				     lssh,
-				     sw,
-				     GH->data[vi][timelevel],
-				     GH->data[vi2][timelevel]); 
-	    break;
-	  case 2: 
-	    berr = BndApplyCopy2Di(GH, 
-				     doBC, 
-				     lssh,
-				     sw,
-				     GH->data[vi][timelevel],
-				     GH->data[vi2][timelevel]); 
-	    break;
-	  case 3: 
-	    berr = BndApplyCopy3Di(GH, 
-				     doBC,
-				     lssh,
-				     sw,
-				     GH->data[vi][timelevel],
-				     GH->data[vi2][timelevel]); 
-	    break;     
-	  default : 
-	    berr = -1; 
-	    CCTK_VWarn(1,__LINE__,__FILE__,"Boundary", 
-		       "No copy boundary for dim>3: grid variable '%s'",
-		       CCTK_VarName(vi));
-	  }
-	}
-      }
-
-      /*      retval = (berr>-1) ? 0 : -1;*/
-    }
-    
-    free(dstag);
-    free(lssh);
-
-    if (!stencil_array)
-    {
-      free(sw);
-    }
-  }
-
-  return retval;
-}
 
+  /* get the group index of the target variable */
+  gindex = CCTK_GroupIndexFromVarI (first_var_to);
 
+  /* get the number of dimensions and the size of the variable's type */
+  gdim      = CCTK_GroupDimI (gindex);
+  vtypesize = CCTK_VarTypeSize (CCTK_VarTypeI (first_var_to));
 
- /*@@
-   @routine    BndApplyCopy3Di
-   @date       Thu Mar  2 11:04:52 2000
-   @author     Gerd Lanfermann
-   @desc 
-      Apply the copy boundary routine, internal routines
-      for the 3d case
-   @enddesc 
-   @calls     
-   @calledby
-   @history 
- 
-   @endhistory 
-
-@@*/
-
-static int BndApplyCopy3Di(cGH *GH,
-			   int doBC,
-			   int *lssh,
-			   int *stencil_size,
-			   CCTK_REAL *var,
-			   CCTK_REAL *varp)
-{
-  int i,j,k,sw,lin;  
-
-  /* lower x */
-  if (doBC == 0)
+  /* make sure we can deal with this number of dimensions */
+  if (gdim > MAXDIM)
   {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (lower x)\n");
-#endif
-    for (k=0;k<lssh[2];k++)
-    {
-      for (j=0;j<lssh[1];j++)
-      {
-	for (sw=0;sw<stencil_size[0];sw++)
-	{
-	  lin = CCTK_GFINDEX3D(GH,sw,j,k);
-	  var[lin] = varp[lin];
-	}
-      }
-    }
+    CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "Variable dimension of %d not supported", gdim);
+    return (-1);
   }
-  
-  /* upper x */
-  if (doBC == 1)
+
+  /* check the direction parameter */
+  if (abs (dir) > gdim)
   {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (upper x)\n");
-#endif
-    for (k=0;k<lssh[2];k++)
-    {
-      for (j=0;j<lssh[1];j++)
-      {
-	for (sw=0;sw<stencil_size[0];sw++)
-	{
-	  lin = CCTK_GFINDEX3D(GH,lssh[0]-sw-1,j,k);
-	  var[lin] = varp[lin];
-	}
-      }
-    }
+    CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING,
+                "ApplyBndCopy: direction %d greater than dimension %d",
+                dir, gdim);
+    return (-2);
   }
 
-  /* lower y */
-  if (doBC == 2)
+  /* set up stencil width array */
+  if (dir)
   {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (lower y)\n");
-#endif
-    for (k=0;k<lssh[2];k++)
-    {
-      for (i=0;i<lssh[0];i++)
-      {
-	for (sw=0;sw<stencil_size[1];sw++)
-	{
-	  lin = CCTK_GFINDEX3D(GH,i,sw,k);
-	  var[lin] = varp[lin];
-	}
-      }
-    }
+    stencil[abs (dir) - 1] = stencil_dir;
   }
-
-  /* upper y */
-  if (doBC == 3)
+  else if (stencil_alldirs)
   {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (upper y)\n");
-#endif
-    for (k=0;k<lssh[2];k++)
-    {
-      for (i=0;i<lssh[0];i++)
-      {
-	for (sw=0;sw<stencil_size[1];sw++)
-	{
-	  lin = CCTK_GFINDEX3D(GH,i,lssh[1]-sw-1,k);
-	  var[lin] = varp[lin];
-	}
-      }
-    }
+    memcpy (stencil, stencil_alldirs, gdim * sizeof (int));
   }
-
-  /* lower z */
-  if (doBC == 4)
+  else
   {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (lower z)\n");
-#endif
-    for (j=0;j<lssh[1];j++)
-    {
-      for (i=0;i<lssh[0];i++)
-      {
-	for (sw=0;sw<stencil_size[2];sw++)
-	{
-	  lin = CCTK_GFINDEX3D(GH,i,j,sw);
-	  var[lin] = varp[lin];
-	}
-      }
-    }
+    CCTK_WARN (1, "ApplyBndCopy: NULL pointer passed for stencil width array");
+    return (-3);
   }
 
-  /* upper z */  
-  if (doBC == 5)
+  /* initialize arrays for variables with less dimensions than MAXDIM
+     so that we can use the INDEX_3D macro later on */
+  for (i = gdim; i < MAXDIM; i++)
   {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (upper z)\n");
-#endif
-    for (j=0;j<lssh[1];j++)
-    {
-      for (i=0;i<lssh[0];i++)
-      {
-	for (sw=0;sw<stencil_size[2];sw++)
-	{
-	  lin = CCTK_GFINDEX3D(GH,i,j,lssh[2]-sw-1);
-	  var[lin] = varp[lin];
-	}
-      }
-    }
+    lssh[i] = 1;
+    lsh[i]  = 0;
   }
 
-  return(0);
-}
+  /* get the directional staggering of the group */
+  CCTK_GroupStaggerDirArrayGI (dstag, gdim, gindex);
 
+  /* get the current timelevel */
+  timelvl_to   = CCTK_NumTimeLevelsFromVarI (first_var_to) - 1;
+  timelvl_from = CCTK_NumTimeLevelsFromVarI (first_var_from) - 1;
 
+  /* see if we have a symmetry array */
+  sGHex = (SymmetryGHex *) CCTK_GHExtension (GH, "Symmetry");
 
- /*@@
-   @routine    BndApplyCopy2Di
-   @date       Thu Mar  2 11:04:52 2000
-   @author     Gerd Lanfermann
-   @desc 
-      Apply the copy boundary routine, internal routines
-      for the 2d case
-   @enddesc 
-   @calls     
-   @calledby BndCopyGI,  BndCopyVN,  BndCopyGN,  BndCopyVI,  
-   @history 
- 
-   @endhistory 
-
-@@*/
-
-static int BndApplyCopy2Di(cGH *GH,
-			   int doBC,
-			   int *lssh,
-			   int *stencil_size,
-			   CCTK_REAL *var,
-			   CCTK_REAL *varp)
-
-{
-  int i,j,sw,lin;  
-
-  /* lower x */
-  if (doBC == 0)
-  {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (lower x)\n");
-#endif
-    for (j=0;j<lssh[1];j++)
-    {
-      for (sw=0;sw<stencil_size[0];sw++)
-      {
-	lin=CCTK_GFINDEX2D(GH,sw,j);
-	var[lin] = varp[lin];
-      }
-    }
-  }
-   
-  /* upper x */
-  if (doBC == 1)
+  /* now loop over all variables */
+  for (var_to = first_var_to, var_from = first_var_from;
+       var_to < first_var_to + num_vars;
+       var_to++, var_from++)
   {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (upper x)\n");
-#endif
-    for (j=0;j<lssh[1];j++)
+    /* Apply condition if:
+       + boundary is not a symmetry boundary (no symmetry or unset(=unsed))
+       + boundary is a physical boundary
+       + have enough grid points
+    */
+    memset (doBC, 1, sizeof (doBC));
+    if (sGHex)
     {
-      for (sw=0;sw<stencil_size[0];sw++)
+      for (i = 0; i < 2 * gdim; i++)
       {
-	lin=CCTK_GFINDEX2D(GH,lssh[0]-sw-1,j);
-	var[lin] = varp[lin];
+        doBC[i] = sGHex->GFSym[var_to][i] == GFSYM_NOSYM ||
+                  sGHex->GFSym[var_to][i] == GFSYM_UNSET;
       }
     }
-  }
-
-  /* lower y */
-  if (doBC == 2)
-  {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (lower y)\n");
-#endif
-    for (i=0;i<lssh[0];i++)
+    for (i = 0; i < gdim; i++)
     {
-      for (sw=0;sw<stencil_size[1];sw++)
+      lsh[i]       = GH->cctk_lsh[i];
+      lssh[i]      = GH->cctk_lssh[CCTK_LSSH_IDX (dstag[i], i)];
+      doBC[i*2]   &= GH->cctk_lsh[i] > 1 && GH->cctk_bbox[i*2];
+      doBC[i*2+1] &= GH->cctk_lsh[i] > 1 && GH->cctk_bbox[i*2+1];
+      if (dir != 0)
       {
-	lin = CCTK_GFINDEX2D(GH,i,sw);
-	var[lin] = varp[lin];
+        doBC[i*2]   &= (dir < 0 && (i + 1 == abs (dir)));
+        doBC[i*2+1] &= (dir > 0 && (i + 1 == abs (dir)));
       }
     }
-  }
 
-  /* upper y */
-  if (doBC == 3)
-  {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (upper y)\n");
-#endif
-    for (i=0;i<lssh[0];i++)
+    /* now copy the boundaries face by face */
+    if (gdim > 0)
     {
-      for (sw=0;sw<stencil_size[1];sw++)
-      {
-	lin = CCTK_GFINDEX2D(GH,i,lssh[1]-sw-1);
-	var[lin] = varp[lin];
-      }
+      /* lower x */
+      COPY_BOUNDARY (doBC[0], stencil[0], lssh[1], lssh[2], i, j, k);
+      /* upper x */
+      COPY_BOUNDARY (doBC[1], stencil[0], lssh[1], lssh[2], lssh[0]-i-1, j, k);
     }
-  }
-
-  return(0);
-}
-
-
-
- /*@@
-   @routine    BndApplyCopy1Di
-   @date       Thu Mar  2 11:04:52 2000
-   @author     Gerd Lanfermann
-   @desc 
-      Apply the copy boundary routine, internal routines
-      for the 1d case
-   @enddesc 
-   @calls     
-   @calledby BndCopyGI,  BndCopyVN,  BndCopyGN,  BndCopyVI,  
-   @history 
- 
-   @endhistory 
-
-@@*/
-
-static int BndApplyCopy1Di(cGH *GH,
-			   int doBC,
-			   int *lssh,
-			   int *stencil_size,
-			   CCTK_REAL *var,
-			   CCTK_REAL *varp)
-
-{
-  int sw,lin;  
-
-
-  /* lower x */
-  if (doBC == 0)
-  {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (lower x)\n");
-#endif
-    for (sw=0;sw<stencil_size[0];sw++)
+    if (gdim > 1)
     {
-      lin=CCTK_GFINDEX1D(GH,sw);
-      var[lin] = varp[lin];
+      /* lower y */
+      COPY_BOUNDARY (doBC[2], lssh[0], stencil[1], lssh[2], i, j, k);
+      /* upper y */
+      COPY_BOUNDARY (doBC[3], lssh[0], stencil[1], lssh[2], i, lssh[1]-j-1, k);
     }
-  }
-   
-  /* upper x */
-  if (doBC == 1)
-  {
-#ifdef DEBUG_BOUNDARY
-    printf("Applying copy boundary (upper x)\n");
-#endif
-    for (sw=0;sw<stencil_size[0];sw++)
+    if (gdim > 2)
     {
-      lin=CCTK_GFINDEX1D(GH,lssh[0]-sw-1);
-      var[lin] = varp[lin];
+      /* lower z */
+      COPY_BOUNDARY (doBC[4], lssh[0], lssh[1], stencil[2], i, j, k);
+      /* upper z */
+      COPY_BOUNDARY (doBC[5], lssh[0], lssh[1], stencil[2], i, j, lssh[2]-k-1);
     }
   }
 
   return(0);
 }
-
-
-