From 63183ce221896a5c6c5db403276405b10a72c06f Mon Sep 17 00:00:00 2001 From: knarf Date: Wed, 22 Sep 2010 14:14:00 +0000 Subject: complain when multiple maxima have been found git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Hydro_Analysis/trunk@112 b6729ddc-ac74-4bd1-908c-9dc7244c52a1 --- src/Hydro_Analysis.c | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) diff --git a/src/Hydro_Analysis.c b/src/Hydro_Analysis.c index d02b64c..565043d 100644 --- a/src/Hydro_Analysis.c +++ b/src/Hydro_Analysis.c @@ -127,10 +127,15 @@ void Hydro_Analysis_LocationSearch(CCTK_ARGUMENTS) if (CCTK_ReduceLocArrayToArray1D(cctkGH, -1, handle_sum, &local_rho_max_loc, level_rho_max_loc, 4, CCTK_VARIABLE_REAL)) CCTK_WARN(0, "Error while reducing local_rho_max_loc"); - if (level_rho_max_loc[3]) { - Hydro_Analysis_rho_max_loc[0] = level_rho_max_loc[0]; - Hydro_Analysis_rho_max_loc[1] = level_rho_max_loc[1]; - Hydro_Analysis_rho_max_loc[2] = level_rho_max_loc[2]; + if (fabs(level_rho_max_loc[3]) < 1.e-14) { + if (fabs(level_rho_max_loc[3]-1.) < 1.e-14) + { + Hydro_Analysis_rho_max_loc[0] = level_rho_max_loc[0]; + Hydro_Analysis_rho_max_loc[1] = level_rho_max_loc[1]; + Hydro_Analysis_rho_max_loc[2] = level_rho_max_loc[2]; + } else + CCTK_WARN(1, "Found more than one identical maximum, not setting anything."); + } /* CCTK_VInfo(CCTK_THORNSTRING, "New location: %g,%g,%g", Hydro_Analysis_rho_max_loc[0], -- cgit v1.2.3