From 2bd1e6a9820538a186c685a894d49c9a8bb85de5 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Mon, 5 Nov 2012 11:10:50 +0100 Subject: Zero the last element of gyz as well. --- src/trumpet.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src') diff --git a/src/trumpet.c b/src/trumpet.c index e2f13ed..1e36b2f 100644 --- a/src/trumpet.c +++ b/src/trumpet.c @@ -187,7 +187,7 @@ void trumpet_data(CCTK_ARGUMENTS) gsl_spline_init(spline, r_iso, r_areal, size); acc = gsl_interp_accel_alloc(); - memset(gyz, 0, sizeof(*gxy)*CCTK_GFINDEX3D(cctkGH, cctk_lsh[0]-1, cctk_lsh[1]-1, cctk_lsh[2]-1)); + memset(gyz, 0, sizeof(*gxy)*CCTK_GFINDEX3D(cctkGH, cctk_lsh[0]-1, cctk_lsh[1]-1, cctk_lsh[2]-1) + 1); #pragma omp parallel for for (int k = 0; k < cctk_lsh[2]; k++) -- cgit v1.2.3