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
39
40
41
42
43
44
45
46
47
48
49
50
|
#include <cctk.h>
#include <cctk_Arguments.h>
#include <cctk_Parameters.h>
#include <carpet.hh>
#include <loopcontrol.h>
namespace Carpet {
using namespace std;
void
CarpetUnusedMask (CCTK_ARGUMENTS)
{
DECLARE_CCTK_ARGUMENTS;
DECLARE_CCTK_PARAMETERS;
dh const & dd = *vdd.AT(map);
ibbox const & ext = dd.light_boxes.AT(mglevel).AT(reflevel).AT(component).exterior;
ibset const & unused_region = dd.local_boxes.AT(mglevel).AT(reflevel).AT(local_component).unused_region;
assert(dim == 3);
// Zero out
LOOP_OVER_BSET (cctkGH, ext, box, imin, imax) {
#pragma omp parallel
CCTK_LOOP3(unused_mask_zero, i, j, k,
imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
CCTK_INT i3D = CCTK_GFINDEX3D(cctkGH, i, j, k);
carpet_unusedpoints_mask[i3D] = 0;
} CCTK_ENDLOOP3(unused_mask_zero);
} END_LOOP_OVER_BSET;
// Set it where unused
LOOP_OVER_BSET (cctkGH, unused_region, box, imin, imax) {
#pragma omp parallel
CCTK_LOOP3(unused_mask_set, i, j, k,
imin[0],imin[1],imin[2], imax[0],imax[1],imax[2],
cctk_ash[0],cctk_ash[1],cctk_ash[2])
{
CCTK_INT i3D = CCTK_GFINDEX3D(cctkGH, i, j, k);
carpet_unusedpoints_mask[i3D] = 1;
} CCTK_ENDLOOP3(unused_mask_set);
} END_LOOP_OVER_BSET;
} // CarpetUnusedMask
}
|