diff options
| author | Anton Khirnov <anton@khirnov.net> | 2019-03-15 14:16:24 +0100 |
|---|---|---|
| committer | Anton Khirnov <anton@khirnov.net> | 2019-03-15 14:39:57 +0100 |
| commit | b50edf6ff11001edc99c42991f7c50988e551648 (patch) | |
| tree | 77330cf3e6f7755cbfc5c583e8486bbb78690788 | |
| parent | 76985513479fb58020a5c87a7acd25c37715ad73 (diff) | |
Fill the corner ghost points for fine-level boundaries.
| -rw-r--r-- | src/maximal_slicing_axi_mg.c | 9 |
1 files changed, 4 insertions, 5 deletions
diff --git a/src/maximal_slicing_axi_mg.c b/src/maximal_slicing_axi_mg.c index b229551..e9f14fd 100644 --- a/src/maximal_slicing_axi_mg.c +++ b/src/maximal_slicing_axi_mg.c @@ -574,14 +574,13 @@ void msa_mg_eval(CCTK_ARGUMENTS) start = gettime(); for (ptrdiff_t idx_z = 0; idx_z < cp->offset_right[1]; idx_z++) - for (ptrdiff_t idx_x = 0; idx_x < cp->boundary_val_stride[0]; idx_x++) { - const ptrdiff_t idx = CCTK_GFINDEX3D(cctkGH, idx_x + cp->offset_left[0], cp->y_idx, cctkGH->cctk_lsh[2] - cp->offset_right[1] + idx_z); + for (ptrdiff_t idx_x = 0; idx_x < cctkGH->cctk_lsh[0]; idx_x++) { + const ptrdiff_t idx = CCTK_GFINDEX3D(cctkGH, idx_x, cp->y_idx, cctkGH->cctk_lsh[2] - cp->offset_right[1] + idx_z); lapse_mg_eval[idx] = fact0 * lapse_prev0[idx] + fact1 * lapse_prev1[idx]; } - for (int idx_z = 0; idx_z < cp->boundary_val_stride[1]; idx_z++) + for (int idx_z = 0; idx_z < cctkGH->cctk_lsh[2]; idx_z++) for (int idx_x = 0; idx_x < cp->offset_right[0]; idx_x++) { - const ptrdiff_t idx = CCTK_GFINDEX3D(cctkGH, cctkGH->cctk_lsh[0] - cp->offset_right[0] + idx_x, cp->y_idx, - cp->offset_left[1] + idx_z); + const ptrdiff_t idx = CCTK_GFINDEX3D(cctkGH, cctkGH->cctk_lsh[0] - cp->offset_right[0] + idx_x, cp->y_idx, idx_z); lapse_mg_eval[idx] = fact0 * lapse_prev0[idx] + fact1 * lapse_prev1[idx]; } |