feq.h

#ifndef UTIL_FEQ_H
#define UTIL_FEQ_H
 
/*
* Util Package - C++ Utilities for Scientific Computation
*
* Copyright 2010 - 2017, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
#include <cmath>
 
namespace Util
{
 
   inline bool feq(double x, double y, double eps = 1.0E-10)
   {
      const double c = 1.0E-5;
      double diff = fabs(x - y) * c / (fabs(x) + fabs(y) + c);
      return (diff < eps);
   }
 
}
#endif