Changeset 499

Timestamp:

Jan 18, 2006, 5:26:09 PM (17 years ago)

Author:

Peter

Message:

set tolerance to round_error

File:

• trunk/test/averager_test.cc (modified) (3 diffs)

trunk/test/averager_test.cc

@@ -8 +8 @@
 
 #include <fstream>
+#include <limits>
 #include <iostream>
 
@@ -71 +72 @@
       
 
-  if ( a.variance() != a.std()*a.std() ){
+  if ( fabs(a.variance() - a.std()*a.std())>
+       std::numeric_limits<double>().round_error() ){
     ok=false;
     *error << "error: std squared should be variance" << std::endl;
-    *error << "std: " << a.std() << std::endl;
+    *error << "std2: " << a.std()*a.std() << std::endl;
     *error << "variance: " << a.variance() << std::endl;
+    *error << "difference is: " << a.std()*a.std()-a.variance() << std::endl;
   }
   
@@ -97 +100 @@
   a.reset();
   a.add(x);
-  const double tol=0.0001;
+  const double tol=std::numeric_limits<double>().round_error();
   if (!equal(a,aw,tol,error)){
     *error << "error: AveragerWeighted with unitary weights should "