diff options
Diffstat (limited to 'src/maximal_slicing_axi_mg.c')
-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]; } |