Carpet/CarpetEvolutionMask/src/evolution_mask_init.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

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

#include "util_Table.h"


void
EvolutionMaskBase_InitEvolutionMask (CCTK_ARGUMENTS)
{
  DECLARE_CCTK_ARGUMENTS;
  DECLARE_CCTK_PARAMETERS;
  
  int i, j, k;
  
  
  /* Initialise the mask to 1 everywhere */
  if (verbose) {
    CCTK_INFO ("Initialising the mask to 1");
  }
  
  for (k=0; k<cctk_lsh[2]; ++k) {
    for (j=0; j<cctk_lsh[1]; ++j) {
      for (i=0; i<cctk_lsh[0]; ++i) {
        
        int const ind = CCTK_GFINDEX3D (cctkGH, i, j, k);
        evolution_mask[ind] = 1.0;
        if (provide_buffer_mask) {
          buffer_mask[ind] = 0;
        }
        
      }
    }
  }

}