Predicates.h

Go to the documentation of this file.

// -*- C++ -*-
//
// This file is part of YODA -- Yet more Objects for Data Analysis
// Copyright (C) 2008-2023 The YODA collaboration (see AUTHORS for details)
//
#ifndef YODA_PREDICATES_H
#define YODA_PREDICATES_H
 
#include <iostream>
 
namespace YODA {
 
 
  struct CmpFloats {
    CmpFloats(double tol=1e-3, double refscale=0.0) : _tol(tol), _refscale(refscale) {}
    bool operator () (const double& a, const double& b) {
      const double div = (_refscale == 0) ? 0.5*(std::abs(a)+std::abs(b)) : _refscale;
      const double dev = (b-a)/div;
      // std::cout << "CmpFloats: " << a << " vs. " << b << " -> dev = " << dev << std::endl;
      return std::abs(dev) < _tol;
    }
    double _tol, _refscale;
  };
 
 
}
 
#endif