diff options
Diffstat (limited to 'src/ReductionAvg.c')
-rw-r--r-- | src/ReductionAvg.c | 14 |
1 files changed, 12 insertions, 2 deletions
diff --git a/src/ReductionAvg.c b/src/ReductionAvg.c index 792202b..a2bc35d 100644 --- a/src/ReductionAvg.c +++ b/src/ReductionAvg.c @@ -27,6 +27,7 @@ static int ReductionAvg (const cGH *GH, int iterator[/* dim */], const int points_per_dim[/* dim */], int num_points, + int have_local_points, int num_inarrays, const int intypes[/* num_inarrays */], const void *const inarrays[/* num_inarrays */], @@ -205,6 +206,7 @@ static int ReductionAvg (const cGH *GH, int iterator[/* dim */], const int points_per_dim[/* dim */], int num_points, + int have_local_points, int num_inarrays, const int intypes[/* num_inarrays */], const void *const inarrays[/* num_inarrays */], @@ -321,8 +323,16 @@ static int ReductionAvg (const cGH *GH, { local_outvals = malloc (total_outvals * sizeof (CCTK_REAL)); - /* outvals[] contains now the local average */ - memcpy (local_outvals, outvals, total_outvals * sizeof (CCTK_REAL)); + if (have_local_points) + { + /* outvals[] contains now the local average */ + memcpy (local_outvals, outvals, total_outvals * sizeof (CCTK_REAL)); + } + else + { + /* initialize local values to zero */ + memset (local_outvals, 0, total_outvals * sizeof (CCTK_REAL)); + } if (proc < 0) { CACTUS_MPI_ERROR (MPI_Allreduce (local_outvals, outvals, total_outvals, |