diff options
Diffstat (limited to 'src/ReductionMin.c')
| -rw-r--r-- | src/ReductionMin.c | 18 |
1 files changed, 16 insertions, 2 deletions
diff --git a/src/ReductionMin.c b/src/ReductionMin.c index d28e24f..ab76dfc 100644 --- a/src/ReductionMin.c +++ b/src/ReductionMin.c @@ -11,6 +11,7 @@ #include <stdlib.h> #include <string.h> +#include <float.h> #include "pugh_reductions.h" @@ -27,6 +28,7 @@ static int ReductionMinVal (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 */], @@ -206,6 +208,7 @@ static int ReductionMinVal (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 */], @@ -322,8 +325,19 @@ static int ReductionMinVal (const cGH *GH, { local_outvals = malloc (total_outvals * sizeof (CCTK_REAL)); - /* outvals[] contains now the local minimum */ - memcpy (local_outvals, outvals, total_outvals * sizeof (CCTK_REAL)); + if (have_local_points) + { + /* outvals[] contains now the local minimum */ + memcpy (local_outvals, outvals, total_outvals * sizeof (CCTK_REAL)); + } + else + { + /* initialize local values to be largest possible fp number */ + for (i = 0; i < total_outvals; i++) + { + local_outvals[i] = DBL_MAX; + } + } if (proc < 0) { CACTUS_MPI_ERROR (MPI_Allreduce (local_outvals, outvals, total_outvals, |