diff options
| author | knarf <knarf@b6729ddc-ac74-4bd1-908c-9dc7244c52a1> | 2010-09-22 14:14:00 +0000 |
|---|---|---|
| committer | knarf <knarf@b6729ddc-ac74-4bd1-908c-9dc7244c52a1> | 2010-09-22 14:14:00 +0000 |
| commit | 63183ce221896a5c6c5db403276405b10a72c06f (patch) | |
| tree | db81d0ad6688d289fd6b0f0191105bf80a895871 | |
| parent | 323f1d8c5ff017626f14872bc157cc5ffdf0b5e2 (diff) | |
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
| -rw-r--r-- | src/Hydro_Analysis.c | 13 |
1 files 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], |