diff options
| author | Anton Khirnov <anton@khirnov.net> | 2019-02-15 18:23:59 +0100 |
|---|---|---|
| committer | Anton Khirnov <anton@khirnov.net> | 2019-02-15 18:23:59 +0100 |
| commit | ebca0ef99f27521d07f419ac6695fc0a2bda22d2 (patch) | |
| tree | 3468c6ba927fc15d90b00c39e06f78a9a4525f31 | |
| parent | 0a532d47e46a0d7780b7b1bda8b57e762dda88ef (diff) | |
Adapt to mg2d boundary specification changes.
Rename FIXDIFF to REFLECT.
Use the new FALLOFF BC on the physical boundary.
| -rw-r--r-- | src/maximal_slicing_axi_mg.c | 8 |
1 files 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; |