Lomont.h

Go to the documentation of this file.

#ifndef GAUDIKERNEL_LOMONT_H 
#define GAUDIKERNEL_LOMONT_H 1
// ============================================================================
// Include file
// ============================================================================
// STD & STL 
// ============================================================================
#include <functional>
// ============================================================================
// GaudiKernel
// ============================================================================
#include "GaudiKernel/Kernel.h"
// ============================================================================
namespace Gaudi
{
  // ==========================================================================
  namespace Math 
  {
    // ========================================================================
    template <class TYPE> class Lomont ;
    // ========================================================================
    GAUDI_API 
    bool lomont_compare_float
    ( const float          af      , 
      const float          bf      , 
      const unsigned short maxULPs ) ;
    // ========================================================================
    GAUDI_API 
    bool lomont_compare_double
    ( const double         af      , 
      const double         bf      , 
      const unsigned int   maxULPs ) ;
    // ========================================================================    
    template <>
    class Lomont<float> : public std::binary_function<float,float,bool>
    {
      // ======================================================================
    public:
      // ======================================================================
      Lomont ( const unsigned short ulps ) : m_ulps ( ulps ) {}
      // ======================================================================
    public:
      // ======================================================================
      inline bool operator () ( const float a , const float b ) const 
      { return lomont_compare_float ( a , b , m_ulps ) ; }
      // ======================================================================
    private: 
      // ======================================================================
      Lomont() ;                         // the default constructor is disabled
      // ======================================================================
    private: 
      // ======================================================================
      unsigned short m_ulps ;         // the precision in "units in last place"
      // ======================================================================      
    };
    // ========================================================================
    template <>
    class Lomont<double> : public std::binary_function<double,double,bool>
    {
      // ======================================================================
    public:
      // ======================================================================
      Lomont ( const unsigned int ulps ) : m_ulps ( ulps ) {}
      // ======================================================================
    public:
      // ======================================================================
      inline bool operator () ( const double a , const double b ) const 
      { return lomont_compare_double ( a , b , m_ulps ) ; }
      // ======================================================================
    private: 
      // ======================================================================
      Lomont () ;                        // the default constructor is disabled
      // ======================================================================
    private: 
      // ======================================================================
      unsigned int m_ulps ;           // the precision in "units in last place"
      // ======================================================================      
    };
    // ========================================================================
    GAUDI_API 
    float next_float ( const float af , const short ulps ) ;
    // ========================================================================
    GAUDI_API 
    double next_double ( const double  af , const short ulps ) ;
    // ========================================================================
  } //                                             end of namespace Gaudi::Math
  // ==========================================================================
} //                                                     end of namespace Gaudi
// ============================================================================
// The END 
// ============================================================================
#endif // GAUDIKERNEL_LOMONT_H
// ============================================================================