diff options
author | Erik Schnetter <schnetter@gmail.com> | 2012-06-20 16:23:37 -0700 |
---|---|---|
committer | Barry Wardell <barry.wardell@gmail.com> | 2012-09-11 18:23:26 +0100 |
commit | 10f226c88943f2848b29735b13223819a0e4e9f7 (patch) | |
tree | 59b6afb08ed66f4a2f3bed2cd578822bf1dde7a5 /Carpet/Carpet | |
parent | 592bc3794e9cab179198ddde1f512f74a1a6d3f0 (diff) |
Carpet: Choose epsilon depending on precision
Diffstat (limited to 'Carpet/Carpet')
-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) { |