From f931e591a2f1c74cee121289776902be5fcca6e8 Mon Sep 17 00:00:00 2001
From: Erik Schnetter <schnetter@gmail.com>
Date: Wed, 20 Jun 2012 16:23:37 -0700
Subject: Carpet: Choose epsilon depending on precision

---
 Carpet/Carpet/src/Evolve.cc | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

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) {
-- 
cgit v1.2.3