/*@@
   @file      MaskInit.c
   @date      Fri Apr 26 17:20:13 2002
   @author    Miguel Alcubierre
   @desc 
   Initialise the mask (I just copied LapseInits.c)
   @enddesc 
   @version $Header$
 @@*/

#include "cctk.h"

#include "cctk_Parameters.h"
#include "cctk_Arguments.h"

static const char *rcsid = "$Header$";

CCTK_FILEVERSION(CactusEinstein_Einstein_MaskInit_c)


/********************************************************************
 *********************     Local Data Types   ***********************
 ********************************************************************/

/********************************************************************
 ********************* Local Routine Prototypes *********************
 ********************************************************************/

/********************************************************************
 ***************** Scheduled Routine Prototypes *********************
 ********************************************************************/

void MaskOne(CCTK_ARGUMENTS);

/********************************************************************
 ********************* Other Routine Prototypes *********************
 ********************************************************************/

/********************************************************************
 *********************     Local Data   *****************************
 ********************************************************************/

/********************************************************************
 *********************     External Routines   **********************
 ********************************************************************/

 /*@@
   @routine    MaskSym
   @date       Fri 3 May 2002
   @author     Gabrielle Allen
   @desc 
   Scheduled routine to set symmetries for mask
   @enddesc 
   @calls     
   @calledby   
   @history 
 
   @endhistory 
 
@@*/
void MaskSym(CCTK_ARGUMENTS)
{
  DECLARE_CCTK_ARGUMENTS

  int one;
  int sym[3];

  one = 1;
  
  sym[0] =  one;
  sym[1] =  one;
  sym[2] =  one;

  SetCartSymVN(cctkGH, sym,"spacemask::emask");

  return;
}

 /*@@
   @routine    MaskOne
   @date       
   @author     Miguel Alcubierre
   @desc 
   Scheduled routine to initialise the mask to one.
   @enddesc 
   @calls     
   @calledby   
   @history 
 
   @endhistory 
 
@@*/
void MaskOne(CCTK_ARGUMENTS)
{
  int i;
  DECLARE_CCTK_ARGUMENTS

  for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
  {
    emask[i] = 1.0;
  }

  return;
}

/********************************************************************
 *********************     Local Routines   *************************
 ********************************************************************/