From ebca0ef99f27521d07f419ac6695fc0a2bda22d2 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Fri, 15 Feb 2019 18:23:59 +0100 Subject: Adapt to mg2d boundary specification changes. Rename FIXDIFF to REFLECT. Use the new FALLOFF BC on the physical boundary. --- src/maximal_slicing_axi_mg.c | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/maximal_slicing_axi_mg.c b/src/maximal_slicing_axi_mg.c index c7d6d18..44fe1ba 100644 --- a/src/maximal_slicing_axi_mg.c +++ b/src/maximal_slicing_axi_mg.c @@ -237,10 +237,10 @@ static CoordPatch *get_coord_patch(MSMGContext *ms, int level) cp->solver->fd_stencil = ms->fd_stencil; - cp->solver->boundaries[MG2D_BOUNDARY_0L]->type = MG2D_BC_TYPE_FIXDIFF; - cp->solver->boundaries[MG2D_BOUNDARY_1L]->type = MG2D_BC_TYPE_FIXDIFF; - cp->solver->boundaries[MG2D_BOUNDARY_0U]->type = MG2D_BC_TYPE_FIXVAL; - cp->solver->boundaries[MG2D_BOUNDARY_1U]->type = MG2D_BC_TYPE_FIXVAL; + cp->solver->boundaries[MG2D_BOUNDARY_0L]->type = MG2D_BC_TYPE_REFLECT; + cp->solver->boundaries[MG2D_BOUNDARY_1L]->type = MG2D_BC_TYPE_REFLECT; + cp->solver->boundaries[MG2D_BOUNDARY_0U]->type = level ? MG2D_BC_TYPE_FIXVAL : MG2D_BC_TYPE_FALLOFF; + cp->solver->boundaries[MG2D_BOUNDARY_1U]->type = level ? MG2D_BC_TYPE_FIXVAL : MG2D_BC_TYPE_FALLOFF; cp->solver->maxiter = ms->maxiter; cp->solver->tol = ms->tol_residual; -- cgit v1.2.3