diff options
| author | schnetter <schnetter@6a38eb6e-646e-4a02-a296-d141613ad6c4> | 2003-11-20 16:40:34 +0000 |
|---|---|---|
| committer | schnetter <schnetter@6a38eb6e-646e-4a02-a296-d141613ad6c4> | 2003-11-20 16:40:34 +0000 |
| commit | 010f0b88babe2c996ac75280211b5ee564452ef2 (patch) | |
| tree | 0e2baf854d44132bdc9fcc2e547e157ffbfff22a | |
| parent | 50462e7aa66393795aadeba16fb0b58871951e23 (diff) | |
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
| -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); \ } \ |