/*@@
   @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"

#include "Symmetry.h"
#include "SpaceMask.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    SetupSpaceMaskRegistry
   @date       November 2002
   @author     Denis Pollney
   @desc 
               Allocate and initialise the global spacemask registry.
   @enddesc 
@@*/
void SetupSpaceMaskRegistry (void)
{
  SpaceMask_Registry* spacemask_registry;

  spacemask_registry = (SpaceMask_Registry*) malloc
    (sizeof(SpaceMask_Registry));

  spacemask_registry->ntypes=0;
  spacemask_registry->type_list = NULL;

  return;
}

/*@@
   @routine    MaskSym
   @date       October 2002
   @author     Denis Pollney
   @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::space_mask");

  return;
}

/*@@
   @routine    MaskSym_emask
   @date       Fri 3 May 2002
   @author     Gabrielle Allen
   @desc 
   Scheduled routine to set symmetries for mask
   @enddesc 
   @calls     
   @calledby   
   @history 
 
   @endhistory 
 
@@*/
void MaskSym_emask(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;
}

/*@@
   @routine    MaskZero
   @date       
   @author     Denis Pollney
   @desc 
               Initialise the mask to zero.
   @enddesc 
   @calls     
   @calledby   
   @history 
 
   @endhistory 
 
@@*/
void MaskZero(CCTK_ARGUMENTS)
{
  DECLARE_CCTK_ARGUMENTS

  int i;

  for(i = 0; i < cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
    space_mask[i] = 0;

  return;
}

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