diff options
-rw-r--r-- | src/RobinBoundary.c | 13 |
1 files changed, 8 insertions, 5 deletions
diff --git a/src/RobinBoundary.c b/src/RobinBoundary.c index f90252c..3dbff30 100644 --- a/src/RobinBoundary.c +++ b/src/RobinBoundary.c @@ -671,11 +671,6 @@ void CCTK_FCALL CCTK_FNAME (BndRobinVN) /* outermost loop over all z points */ \ for (k = 0; k < GH->cctk_lsh[2]; k++) \ { \ - if ((k == 0 && ! doBC[4]) || (k == GH->cctk_lsh[2]-1 && ! doBC[5])) \ - { \ - continue; \ - } \ - \ dz = 0; \ if (k == 0 && doBC[4]) \ { \ @@ -700,9 +695,13 @@ void CCTK_FCALL CCTK_FNAME (BndRobinVN) } \ \ /* lower x */ \ + dx = 0; \ if (doBC[0]) \ { \ dx = +1; \ + } \ + if (dx || dy || dz) \ + { \ SET_LINEAR_INDICES (0); \ ROBIN_BOUNDARY_TYPED_3D (cctk_type); \ } \ @@ -719,9 +718,13 @@ void CCTK_FCALL CCTK_FNAME (BndRobinVN) } \ \ /* upper x */ \ + dx = 0; \ if (doBC[1]) \ { \ dx = -1; \ + } \ + if (dx || dy || dz) \ + { \ SET_LINEAR_INDICES (GH->cctk_lsh[0]-1); \ ROBIN_BOUNDARY_TYPED_3D (cctk_type); \ } \ |