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