From 010f0b88babe2c996ac75280211b5ee564452ef2 Mon Sep 17 00:00:00 2001 From: schnetter Date: Thu, 20 Nov 2003 16:40:34 +0000 Subject: Apply Robin boundary conditions also in the ghost zones and on symmetry boundaties. This makes them consistent with the other boundary conditions. git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/Boundary/trunk@255 6a38eb6e-646e-4a02-a296-d141613ad6c4 --- src/RobinBoundary.c | 13 ++++++++----- 1 file 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); \ } \ -- cgit v1.2.3