From c1fb1d45924159fa8fe96cfbfdcd6cdc5e6d0bc3 Mon Sep 17 00:00:00 2001 From: Christian Reisswig Date: Tue, 7 Jun 2011 11:22:57 -0700 Subject: Erik's patch to fix a reduction operator problem. Without this fix, reductions may depend on the number of processors used. --- Carpet/CarpetReduce/src/mask_coords.c | 12 ++++++++---- 1 file changed, 8 insertions(+), 4 deletions(-) (limited to 'Carpet/CarpetReduce') diff --git a/Carpet/CarpetReduce/src/mask_coords.c b/Carpet/CarpetReduce/src/mask_coords.c index 814f1ecfb..c324671d1 100644 --- a/Carpet/CarpetReduce/src/mask_coords.c +++ b/Carpet/CarpetReduce/src/mask_coords.c @@ -109,10 +109,12 @@ CoordBase_SetupMask (CCTK_ARGUMENTS) } if (f==0) { /* lower face */ - bmax[d] = imin[d] + bnd_points[2*d+f]; + //bmax[d] = imin[d] + bnd_points[2*d+f]; + bmax[d] = imin[d] + cctk_nghostzones[d]; } else { /* upper face */ - bmin[d] = imax[d] - bnd_points[2*d+f]; + //bmin[d] = imax[d] - bnd_points[2*d+f]; + bmin[d] = imax[d] - cctk_nghostzones[d]; } /* Loop over the boundary */ @@ -201,10 +203,12 @@ CoordBase_SetupMask (CCTK_ARGUMENTS) } if (f==0) { /* lower face */ - bmax[d] = imin[d] + bnd_points[2*d+f]; + //bmax[d] = imin[d] + bnd_points[2*d+f]; + bmax[d] = imin[d] + cctk_nghostzones[d]; } else { /* upper face */ - bmin[d] = imax[d] - bnd_points[2*d+f]; + //bmin[d] = imax[d] - bnd_points[2*d+f]; + bmin[d] = imax[d] - cctk_nghostzones[d]; } /* Loop over the boundary */ -- cgit v1.2.3