CoinFloatEqual.hpp

Go to the documentation of this file.

/* $Id: CoinFloatEqual.hpp 2083 2019-01-06 19:38:09Z unxusr $ */
// Copyright (C) 2000, International Business Machines
// Corporation and others.  All Rights Reserved.
// This code is licensed under the terms of the Eclipse Public License (EPL).

#ifndef CoinFloatEqual_H
#define CoinFloatEqual_H

#include <algorithm>
#include <cmath>

#include "CoinFinite.hpp"

class CoinAbsFltEq {
public:

  inline bool operator()(const double f1, const double f2) const

  {
    if (CoinIsnan(f1) || CoinIsnan(f2))
      return false;
    if (f1 == f2)
      return true;
    return (fabs(f1 - f2) < epsilon_);
  }


  CoinAbsFltEq()
    : epsilon_(1.e-10)
  {
  }


  CoinAbsFltEq(const double epsilon)
    : epsilon_(epsilon)
  {
  }


  virtual ~CoinAbsFltEq() {}


  CoinAbsFltEq(const CoinAbsFltEq &src)
    : epsilon_(src.epsilon_)
  {
  }


  CoinAbsFltEq &operator=(const CoinAbsFltEq &rhs)

  {
    if (this != &rhs)
      epsilon_ = rhs.epsilon_;
    return (*this);
  }


private:


  double epsilon_;

};

class CoinRelFltEq {
public:

  inline bool operator()(const double f1, const double f2) const

  {
    if (CoinIsnan(f1) || CoinIsnan(f2))
      return false;
    if (f1 == f2)
      return true;
    if (!CoinFinite(f1) || !CoinFinite(f2))
      return false;

    double tol = (fabs(f1) > fabs(f2)) ? fabs(f1) : fabs(f2);

    return (fabs(f1 - f2) <= epsilon_ * (1 + tol));
  }


#ifndef COIN_FLOAT

  CoinRelFltEq()
    : epsilon_(1.e-10)
  {
  }
#else

  CoinRelFltEq()
    : epsilon_(1.e-6) {}; // as float
#endif


  CoinRelFltEq(const double epsilon)
    : epsilon_(epsilon)
  {
  }


  virtual ~CoinRelFltEq() {}


  CoinRelFltEq(const CoinRelFltEq &src)
    : epsilon_(src.epsilon_)
  {
  }


  CoinRelFltEq &operator=(const CoinRelFltEq &rhs)

  {
    if (this != &rhs)
      epsilon_ = rhs.epsilon_;
    return (*this);
  }


private:


  double epsilon_;

};

#endif

/* vi: softtabstop=2 shiftwidth=2 expandtab tabstop=2
*/