CoinFloatEqual.hpp

Go to the documentation of this file.

/* $Id: CoinFloatEqual.hpp 1416 2011-04-17 09:57:29Z stefan $ */
// 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