diff options
author | Erik Schnetter <schnetter@gmail.com> | 2013-01-16 14:41:41 -0500 |
---|---|---|
committer | Erik Schnetter <schnetter@gmail.com> | 2013-01-16 14:41:41 -0500 |
commit | 0441c40fe9ad49d9e25d360b0da7f69d26ca82e6 (patch) | |
tree | 8b384ef0efedc8250b551fdb93069e7cc929f52b /Carpet | |
parent | ba9d75ac3f15bf6285f763be96531dc807ca7a0a (diff) |
CarpetReduce: Clean up code defining reduction weight masks
Diffstat (limited to 'Carpet')
-rw-r--r-- | Carpet/CarpetReduce/src/mask_coords.c | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/Carpet/CarpetReduce/src/mask_coords.c b/Carpet/CarpetReduce/src/mask_coords.c index 121bbe7f4..0a07bc854 100644 --- a/Carpet/CarpetReduce/src/mask_coords.c +++ b/Carpet/CarpetReduce/src/mask_coords.c @@ -164,26 +164,26 @@ CoordBase_SetupMask (CCTK_ARGUMENTS) } /* Loop over the points next to boundary */ - if (verbose) { - CCTK_VInfo (CCTK_THORNSTRING, - "Setting non-staggered boundary points in direction %d face %d to weight 1/2 on level %d", d, f, reflevel); - } unsigned const bits = BMSK(cctk_dim); unsigned bmask = 0; for (unsigned b=0; b<bits; ++b) { - if (BGET(b,d) == f) { + if (BGET(b,d) == !f) { bmask = BSET(bmask, b); } } assert (BCNT(bmask) == bits/2); + if (verbose) { + CCTK_VInfo (CCTK_THORNSTRING, + "Setting non-staggered boundary points on level %d in direction %d face %d to bmask 0x%x", reflevel, d, f, bmask); + } #pragma omp parallel CCTK_LOOP3(CoordBase_SetupMask_boundary2, i,j,k, bmin[0],bmin[1],bmin[2], bmax[0],bmax[1],bmax[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_ash[0],cctk_ash[1],cctk_ash[2]) { int const ind = CCTK_GFINDEX3D (cctkGH, i, j, k); - iweight[ind] &= ~bmask; + iweight[ind] &= bmask; } CCTK_ENDLOOP3(CoordBase_SetupMask_boundary2); } /* if the domain is not empty */ |