FLOATassert.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TESTTOOLS_FLOATASSERT_H
#define TESTTOOLS_FLOATASSERT_H
 
#include <iostream>
#include <cassert>
#include <cfloat>
#include <cmath>
 
 
namespace Athena_test {
  inline
  bool floatEQ(float lhs, float rhs) {
    return fabs(lhs-rhs)<=FLT_EPSILON;
  }
  inline
  bool floatNEQ(float lhs, float rhs) {
    return fabs(lhs-rhs)>FLT_EPSILON;
  }
 
  inline
  bool isEqual (double x1, double x2, double thresh = 1e-6)
  {
    double den = std::abs(x1+x2);
    if (den < thresh) return true;
    double diff = std::abs (x1-x2) / den;
    if (diff < thresh)
      return true;
    std::cout << "Match failure: " << x1 << " " << x2 << " (" << diff << ")\n";
    return false;
  }
 
}
 
#define FLOAT_NEQassert( LHS, RHS ) assert(Athena_test::floatNEQ(LHS, RHS));    
#define FLOAT_EQassert( LHS, RHS ) assert(Athena_test::floatEQ(LHS, RHS));  
 
#endif