diff options
-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], |