CoinFloatEqual.hpp

Go to the documentation of this file.

/* $Id$ */
// 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
*/