From e6a0b26fdbbf2d7669de937c67c8da48ce07cc56 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Thu, 23 May 2019 14:57:53 +0200 Subject: Drop commented out obsolete code. --- src/maximal_slicing_axi_mg.c | 66 -------------------------------------------- 1 file changed, 66 deletions(-) diff --git a/src/maximal_slicing_axi_mg.c b/src/maximal_slicing_axi_mg.c index cae7197..69832cf 100644 --- a/src/maximal_slicing_axi_mg.c +++ b/src/maximal_slicing_axi_mg.c @@ -875,72 +875,6 @@ void msa_mg_sync(CCTK_ARGUMENTS) return; } -#if 0 -static void maximal_slicing_axi_mg_old(CCTK_ARGUMENTS) -{ - CoordPatch *cp; - - DECLARE_CCTK_ARGUMENTS; - DECLARE_CCTK_PARAMETERS; - - double *coeffs = NULL; - int i, ret; - - fprintf(stderr, "ms mg solve: level %d time %g\n", ctz(ms->gh->cctk_levfac[0]), ms->gh->cctk_time); - - cp = get_coord_patch(ms); - - fill_eq_coeffs(ms->gh, cp); - - CCTK_TimerStart("MaximalSlicingAxiMG_Solve"); - ret = mg2d_solve(cp->solver); - if (ret < 0) - CCTK_WARN(0, "Error solving the maximal slicing equation"); - CCTK_TimerStop("MaximalSlicingAxiMG_Solve"); - - - CCTK_TimerStart("MaximalSlicingAxiMG_Copy"); -//#pragma omp parallel for - for (int j = 0; j < cp->solver->domain_size; j++) - for (int i = 0; i < cp->solver->domain_size; i++) { - const int idx_dst = CCTK_GFINDEX3D(ms->gh, i + cp->offset_left[0], cp->y_idx, j + cp->offset_left[1]); - const int idx_src = j * cp->solver->u_stride + i; - alpha[idx_dst] = 1.0 + cp->solver->u[idx_src]; - } - if (cp->level > 0) { - /* fill in the interpolated outer buffer points */ - for (int idx_z = 0; idx_z < cp->offset_right[1]; idx_z++) - for (int idx_x = 0; idx_x < cp->boundary_val_stride[0]; idx_x++) { - const ptrdiff_t idx_src = idx_z * cp->boundary_val_stride[0] + idx_x; - const ptrdiff_t idx_dst = CCTK_GFINDEX3D(ms->gh, cp->offset_left[0] + idx_x, cp->y_idx, ms->gh->cctk_lsh[2] - cp->offset_right[1] + idx_z); - alpha[idx_dst] = cp->boundary_val[0][idx_src]; - } - for (int idx_z = 0; idx_z < cp->boundary_val_stride[1]; idx_z++) - for (int idx_x = 0; idx_x < cp->offset_right[0]; idx_x++) { - const ptrdiff_t idx_src = idx_x * cp->boundary_val_stride[1] + idx_z; - const ptrdiff_t idx_dst = CCTK_GFINDEX3D(ms->gh, ms->gh->cctk_lsh[0] - cp->offset_right[0] + idx_x, cp->y_idx, - cp->offset_left[1] + idx_z); - alpha[idx_dst] = cp->boundary_val[1][idx_src]; - } - } - - /* fill in the axial symmetry ghostpoints by mirroring */ - for (int idx_z = cp->offset_left[1]; idx_z < ms->gh->cctk_lsh[2]; idx_z++) - for (int idx_x = 0; idx_x < cp->offset_left[0]; idx_x++) { - const ptrdiff_t idx_dst = CCTK_GFINDEX3D(ms->gh, idx_x, cp->y_idx, idx_z); - const ptrdiff_t idx_src = CCTK_GFINDEX3D(ms->gh, 2 * cp->offset_left[0] - idx_x, cp->y_idx, idx_z); - alpha[idx_dst] = alpha[idx_src]; - } - for (int idx_z = 0; idx_z < cp->offset_left[1]; idx_z++) - for (int idx_x = 0; idx_x < ms->gh->cctk_lsh[0]; idx_x++) { - const ptrdiff_t idx_dst = CCTK_GFINDEX3D(ms->gh, idx_x, cp->y_idx, idx_z); - const ptrdiff_t idx_src = CCTK_GFINDEX3D(ms->gh, idx_x, cp->y_idx, 2 * cp->offset_left[1] - idx_z); - alpha[idx_dst] = alpha[idx_src]; - } - CCTK_TimerStop("MaximalSlicingAxiMG_Copy"); -} -#endif - void msa_mg_init(CCTK_ARGUMENTS) { DECLARE_CCTK_ARGUMENTS; -- cgit v1.2.3