From b50edf6ff11001edc99c42991f7c50988e551648 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Fri, 15 Mar 2019 14:16:24 +0100 Subject: Fill the corner ghost points for fine-level boundaries. --- src/maximal_slicing_axi_mg.c | 9 ++++----- 1 file 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]; } -- cgit v1.2.3