diff options
-rw-r--r-- | src/volume_integral.c | 13 |
1 files changed, 10 insertions, 3 deletions
diff --git a/src/volume_integral.c b/src/volume_integral.c index a4d3748..b6fa223 100644 --- a/src/volume_integral.c +++ b/src/volume_integral.c @@ -89,15 +89,22 @@ void ADMMass_Volume(CCTK_ARGUMENTS) ghost = *(const int *)ghost_ptr; } - for(i=1; i<cctk_lsh[0]-1; i++) - for(j=1; j<cctk_lsh[1]-1; j++) - for(k=1; k<cctk_lsh[2]-1; k++) + 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_VolumeMass_GF[ijk] = 0.0; ADMMass_VolumeMass_pot_x[ijk] = 0.0; ADMMass_VolumeMass_pot_y[ijk] = 0.0; ADMMass_VolumeMass_pot_z[ijk] = 0.0; + } + + for(i=1; i<cctk_lsh[0]-1; i++) + for(j=1; j<cctk_lsh[1]-1; j++) + for(k=1; k<cctk_lsh[2]-1; k++) + { + ijk = CCTK_GFINDEX3D(cctkGH, i, j, k); #include "CactusEinstein/ADMMacros/src/macro/UPPERMET_guts.h" u[0][0] = UPPERMET_UXX; |