diff options
author | knarf <knarf@54511f98-0e4f-0410-826e-eb8b393f5a1e> | 2011-08-23 18:58:36 +0000 |
---|---|---|
committer | knarf <knarf@54511f98-0e4f-0410-826e-eb8b393f5a1e> | 2011-08-23 18:58:36 +0000 |
commit | e943a66a52228669b89c7403ee8c1eb1b56900a0 (patch) | |
tree | ee4ed0a6e0ef926da1077e56f85c10242c0f4c10 | |
parent | 72b795598c4af9577fea8c9cc980d3a518b6178f (diff) |
initialize ADMMass_SurfaceMass_GF completely
git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/ADMMass/trunk@41 54511f98-0e4f-0410-826e-eb8b393f5a1e
-rw-r--r-- | src/surface_integral.c | 8 |
1 files changed, 7 insertions, 1 deletions
diff --git a/src/surface_integral.c b/src/surface_integral.c index fa14602..6886981 100644 --- a/src/surface_integral.c +++ b/src/surface_integral.c @@ -190,12 +190,18 @@ void ADMMass_Surface(CCTK_ARGUMENTS) z_max_k = find_closest(cctkGH, cctk_lsh, cctk_delta_space, ghost, z, *ADMMass_box_z_max, 2); + for(i=0; i<cctk_lsh[0]; i++) + for(j=0; j<cctk_lsh[1]; j++) + for(k=0; k<cctk_lsh[2]; k++) + { + ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); + ADMMass_SurfaceMass_GF[ijk] = 0.0; + } for(i=ghost; i<cctk_lsh[0]-ghost; i++) for(j=ghost; j<cctk_lsh[1]-ghost; j++) for(k=ghost; k<cctk_lsh[2]-ghost; k++) { ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); - ADMMass_SurfaceMass_GF[ijk] = 0.0; /* delimit the cube on whose surface we want to integrate */ if ((i >= x_min_i) && |