/*@@
   @file      ConstantBoundary.c
   @date      Mon Mar 15 15:09:00 1999
   @author    Gerd Lanfermann, Gabrielle Allen
   @desc 
     Constant boundary conditions for a given GF group or single GF
   @enddesc 
 @@*/

/*#define DEBUG_BOUND*/

#include <stdio.h>
#include <assert.h>
#include <stdlib.h>
#include <ctype.h>
#include <stdarg.h>
#include <string.h>

#include "cctk.h"
#include "cctk_Parameters.h"
#include "CactusBase/CartGrid3D/src/Symmetry.h"
#include "cctk_FortranString.h"


 /*@@
   @routine    BndApplyConstant3Di
   @date       Thu Mar  2 11:08:49 2000
   @author     Gerd Lanfermann
   @desc 
   Apply the constant boundary condition for 3d variables
   @enddesc 
   @calls     
   @calledby BndConstantGI,  BndConstantVI,  BndConstantVN,  BndConstantGN,  
   @history 
 
   @endhistory 

@@*/


int BndApplyConstant3Di(cGH *GH,
			int  gdim,
			int *doBC,
			int *lsh,
			int *stencil_size,
			CCTK_REAL var0,
			CCTK_REAL *var)
{
  int i,j,k,sw;

  if (doBC[0] == 1)
  {
    for (k=0;k<lsh[2];k++)
    {
      for (j=0;j<lsh[1];j++)
      {
	for (sw=0;sw<stencil_size[0];sw++)
	{
	  var[CCTK_GFINDEX3D(GH,sw,j,k)]
	    = var0;
	}
      }
    }
  }
  
  if (doBC[1] == 1)
  {
    for (k=0;k<lsh[2];k++)
    {
      for (j=0;j<lsh[1];j++)
      {
	for (sw=0;sw<stencil_size[0];sw++)
	{
	  var[CCTK_GFINDEX3D(GH,lsh[0]-sw-1,j,k)]
	    = var0;
	}
      }
    }
  }

  if (doBC[2] == 1)
  {
    for (k=0;k<lsh[2];k++)
    {
      for (i=0;i<lsh[0];i++)
      {
	for (sw=0;sw<stencil_size[1];sw++)
	{
	  var[CCTK_GFINDEX3D(GH,i,sw,k)]
	    = var0;
	}
      }
    }
  }
  
  if (doBC[3] == 1)
  {
    for (k=0;k<lsh[2];k++)
    {
      for (i=0;i<lsh[0];i++)
      {
	for (sw=0;sw<stencil_size[1];sw++)
	{
	  var[CCTK_GFINDEX3D(GH,i,lsh[1]-sw-1,k)]
	    = var0;
	}
      }
    }
  }

  if (doBC[4] == 1)
  {
    for (j=0;j<lsh[1];j++)
    {
      for (i=0;i<lsh[0];i++)
      {
	for (sw=0;sw<stencil_size[2];sw++)
	{
	  var[CCTK_GFINDEX3D(GH,i,j,sw)]
	    = var0;
	}
      }
    }
  }
  
  if (doBC[5] == 1)
  {
    for (j=0;j<lsh[1];j++)
    {
      for (i=0;i<lsh[0];i++)
      {
	for (sw=0;sw<stencil_size[2];sw++)
	{
	  var[CCTK_GFINDEX3D(GH,i,j,lsh[2]-sw-1)]
	    = var0;
	}
      }
    }
  }
  return(0);
}


 /*@@
   @routine    BndApplyConstant2Di
   @date       Thu Mar  2 11:08:49 2000
   @author     Gerd Lanfermann
   @desc 
   Apply the constant boundary condition for 2d variables
   @enddesc 
   @calls     
   @calledby BndConstantGI,  BndConstantVI,  BndConstantVN,  BndConstantGN,  
   @history 
 
   @endhistory 

@@*/

int BndApplyConstant2Di(cGH *GH,
			int  gdim,
			int *doBC,
			int *lsh,
			int *stencil_size,
			CCTK_REAL var0,
			CCTK_REAL *var)
{
  int i,j,k,sw;
  
  /* lower x */
  if (doBC[0] == 1)
  {
    for (j=0;j<lsh[1];j++)
    {
      for (sw=0;sw<stencil_size[0];sw++)
      {
	var[CCTK_GFINDEX3D(GH,sw,j,k)]
	  = var0;
      }
    }
  }

  /* upper x */
  if (doBC[1] == 1)
  {
    for (j=0;j<lsh[1];j++)
    {
      for (sw=0;sw<stencil_size[0];sw++)
      {
	var[CCTK_GFINDEX3D(GH,lsh[0]-sw-1,j,k)]
	  = var0;
      }
    }
  }

  /* lower y */
  if (doBC[2] == 1)
  {
    for (i=0;i<lsh[0];i++)
    {
      for (sw=0;sw<stencil_size[1];sw++)
      {
	var[CCTK_GFINDEX3D(GH,i,sw,k)]
	  = var0;
      }
    }
  }
  
  /* upper y */
  if (doBC[3] == 1)
  {
    for (i=0;i<lsh[0];i++)
    {
      for (sw=0;sw<stencil_size[1];sw++)
      {
	var[CCTK_GFINDEX3D(GH,i,lsh[1]-sw-1,k)]
	  = var0;
      }
    }
  }
  return(0);
}



/*@@
   @routine    BndApplyConstant1Di
   @date       Thu Mar  2 11:08:49 2000
   @author     Gerd Lanfermann
   @desc 
   Apply the constant boundary condition for 1d variables
   @enddesc 
   @calls     
   @calledby BndConstantGI,  BndConstantVI,  BndConstantVN,  BndConstantGN,  
   @history 
 
   @endhistory 

@@*/

int BndApplyConstant1Di(cGH *GH, 
			int  gdim, 
			int *doBC, 
			int *stencil, 
			int *dstag, 
			CCTK_REAL var0,
			CCTK_REAL *var) 
{ int sw;

  /* lower x */
  if (doBC[0] == 1)
  {
    for (sw=0;sw<stencil[0];sw++)
    {
      var[CCTK_GFINDEX1D(GH,sw)]
	= var0;
    }
  }
  
  /* upper x */
  if (doBC[1] == 1)
  {
    for (sw=0;sw<stencil[0];sw++)
    {
      var[CCTK_GFINDEX1D(GH,GH->cctk_lsh[0]-sw-1)]
	= var0;
    }
  }
  return(0);
}