diff options
| author | Erik Schnetter <schnetter@gmail.com> | 2012-06-20 16:23:37 -0700 |
|---|---|---|
| committer | Erik Schnetter <schnetter@gmail.com> | 2012-06-20 16:23:37 -0700 |
| commit | f931e591a2f1c74cee121289776902be5fcca6e8 (patch) | |
| tree | b9a268a8a8a9766e10a0ab80ee9e7f0898e80521 | |
| parent | 55fdb0854e7af97103d04f5e8a23b15f2d15f715 (diff) | |
Carpet: Choose epsilon depending on precision
| -rw-r--r-- | Carpet/Carpet/src/Evolve.cc | 5 |
1 files changed, 4 insertions, 1 deletions
diff --git a/Carpet/Carpet/src/Evolve.cc b/Carpet/Carpet/src/Evolve.cc index 5fd146991..19e25fc9b 100644 --- a/Carpet/Carpet/src/Evolve.cc +++ b/Carpet/Carpet/src/Evolve.cc @@ -1,8 +1,10 @@ #include <algorithm> #include <cassert> +#include <cmath> #include <cstdio> #include <cstdlib> #include <iostream> +#include <limits> #include <map> #include <string> #include <sstream> @@ -104,7 +106,8 @@ namespace Carpet { { Timer timer("CheckLevelTimes"); timer.start(); - CCTK_REAL const eps = 1.0e-12; + CCTK_REAL const eps = + pow(numeric_limits<CCTK_REAL>::epsilon(), CCTK_REAL(0.75)); assert (abs (cctkGH->cctk_time - global_time) <= eps * global_time); for (int ml=0; ml<mglevels; ++ml) { for (int rl=0; rl<reflevels; ++rl) { |