Drop commented out obsolete code.

1 files changed, 0 insertions, 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;