diff options
Diffstat (limited to 'src/ScalarBoundary.c')
-rw-r--r-- | src/ScalarBoundary.c | 27 |
1 files changed, 11 insertions, 16 deletions
diff --git a/src/ScalarBoundary.c b/src/ScalarBoundary.c index e978344..eb84636 100644 --- a/src/ScalarBoundary.c +++ b/src/ScalarBoundary.c @@ -954,29 +954,29 @@ void CCTK_FCALL CCTK_FNAME (BndScalarVN) if (gdim > 0) \ { \ /* lower x */ \ - SCALAR_BOUNDARY_TYPED (doBC[0], widths[0], lssh[1], lssh[2], \ + SCALAR_BOUNDARY_TYPED (doBC[0], widths[0], lsh[1], lsh[2], \ i, j, k, left_cctk_type, right_cctk_type); \ /* upper x */ \ - SCALAR_BOUNDARY_TYPED (doBC[1], widths[1], lssh[1], lssh[2], \ - lssh[0]-i-1, j, k, left_cctk_type, right_cctk_type);\ + SCALAR_BOUNDARY_TYPED (doBC[1], widths[1], lsh[1], lsh[2], \ + lsh[0]-i-1, j, k, left_cctk_type, right_cctk_type);\ } \ if (gdim > 1) \ { \ /* lower y */ \ - SCALAR_BOUNDARY_TYPED (doBC[2], lssh[0], widths[2], lssh[2], \ + SCALAR_BOUNDARY_TYPED (doBC[2], lsh[0], widths[2], lsh[2], \ i, j, k, left_cctk_type, right_cctk_type); \ /* upper y */ \ - SCALAR_BOUNDARY_TYPED (doBC[3], lssh[0], widths[3], lssh[2], \ - i, lssh[1]-j-1, k, left_cctk_type, right_cctk_type);\ + SCALAR_BOUNDARY_TYPED (doBC[3], lsh[0], widths[3], lsh[2], \ + i, lsh[1]-j-1, k, left_cctk_type, right_cctk_type);\ } \ if (gdim > 2) \ { \ /* lower z */ \ - SCALAR_BOUNDARY_TYPED (doBC[4], lssh[0], lssh[1], widths[4], \ + SCALAR_BOUNDARY_TYPED (doBC[4], lsh[0], lsh[1], widths[4], \ i, j, k, left_cctk_type, right_cctk_type); \ /* upper z */ \ - SCALAR_BOUNDARY_TYPED (doBC[5], lssh[0], lssh[1], widths[5], \ - i, j, lssh[2]-k-1, left_cctk_type, right_cctk_type);\ + SCALAR_BOUNDARY_TYPED (doBC[5], lsh[0], lsh[1], widths[5], \ + i, j, lsh[2]-k-1, left_cctk_type, right_cctk_type);\ } \ } @@ -1062,7 +1062,7 @@ static int ApplyBndScalar (const cGH *GH, int i, j, k; int gindex, gdim; int var, timelvl; - int doBC[2*MAXDIM], dstag[MAXDIM], lsh[MAXDIM], lssh[MAXDIM]; + int doBC[2*MAXDIM], lsh[MAXDIM]; CCTK_INT widths[2*MAXDIM]; CCTK_INT symtable; CCTK_INT symbnd[2*MAXDIM]; @@ -1111,13 +1111,9 @@ static int ApplyBndScalar (const cGH *GH, so that we can use the INDEX_3D macro later on */ for (i = gdim; i < MAXDIM; i++) { - lsh[i] = 0; - lssh[i] = 1; + lsh[i] = 0; } - /* get the directional staggering of the group */ - CCTK_GroupStaggerDirArrayGI (dstag, gdim, gindex); - /* get the current timelevel */ timelvl = 0; @@ -1146,7 +1142,6 @@ static int ApplyBndScalar (const cGH *GH, for (i = 0; i < gdim; i++) { lsh[i] = GH->cctk_lsh[i]; - lssh[i] = GH->cctk_lssh[CCTK_LSSH_IDX (dstag[i], i)]; doBC[i*2] &= GH->cctk_lsh[i] > widths[i*2] && GH->cctk_bbox[i*2]; doBC[i*2+1] &= GH->cctk_lsh[i] > widths[i*2+1] && GH->cctk_bbox[i*2+1]; if (dir != 0) |