diff options
Diffstat (limited to 'src/qms.c')
| -rw-r--r-- | src/qms.c | 12 |
1 files changed, 6 insertions, 6 deletions
@@ -1309,7 +1309,7 @@ void qms_mg_solve(CCTK_ARGUMENTS) double *dst = bnd->val + j * bnd->val_stride; for (ptrdiff_t i = -j; i < (ptrdiff_t)cp->solver->local_size[0] + j; i++) { const ptrdiff_t idx = CCTK_GFINDEX3D(cctkGH, ABS(i) + cp->offset_left[0], cp->y_idx, cp->offset_left[1] + cp->solver->local_size[1] - 1 + j); - dst[i] = W_mg[idx]; + dst[i] = W_val[idx]; } } } @@ -1320,7 +1320,7 @@ void qms_mg_solve(CCTK_ARGUMENTS) double *dst = bnd->val + j * bnd->val_stride; for (ptrdiff_t i = -j; i < (ptrdiff_t)cp->solver->local_size[1] + j; i++) { const ptrdiff_t idx = CCTK_GFINDEX3D(cctkGH, cp->offset_left[1] + cp->solver->local_size[0] - 1 + j, cp->y_idx, cp->offset_left[1] + ABS(i)); - dst[i] = W_mg[idx]; + dst[i] = W_val[idx]; } } } @@ -1347,9 +1347,9 @@ void qms_mg_solve(CCTK_ARGUMENTS) start = gettime(); { - double *W_mg_1 = CCTK_VarDataPtr(cctkGH, 1, "QuasiMaximalSlicingMG::W_mg"); - solution_to_grid(cp, cp->solver, W_mg); - memcpy(W_mg_1, W_mg, grid_size * sizeof(*W_mg_1)); + double *W_val_tl1 = CCTK_VarDataPtr(cctkGH, 1, "QuasiMaximalSlicingMG::W_val"); + solution_to_grid(cp, cp->solver, W_val); + memcpy(W_val_tl1, W_val, grid_size * sizeof(*W_val_tl1)); } /* add the solution to the list of past solutions */ @@ -1357,7 +1357,7 @@ void qms_mg_solve(CCTK_ARGUMENTS) memcpy(W_val0, W_val1, sizeof(*W_val0) * grid_size); W_val0_time[reflevel] = W_val1_time[reflevel]; - memcpy(W_val1, W_mg, sizeof(*W_val1) * grid_size); + memcpy(W_val1, W_val, sizeof(*W_val1) * grid_size); W_val1_time[reflevel] = ms->gh->cctk_time; } |