gmDoubleCmp.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/

#ifndef _UTILS_DOUBLE_CMP_
#define _UTILS_DOUBLE_CMP_

#include <math.h>


namespace GmDoubleCmp
{

#define GM_DOUBLECMP_RELTOL (1.0e-8)

#define GM_DOUBLECMP_ABSTOL (1.0E-12)

inline bool equal(double a, double b, double relTol=GM_DOUBLECMP_RELTOL, double absTol=GM_DOUBLECMP_ABSTOL)
{
  // Trata valores muito próximos de zero 
  double diff = fabs(a-b);
  if(diff < absTol) 
    return true; 

  a = fabs(a);
  b = fabs(b);
  return diff <= relTol * (a > b ? a : b);
}

inline bool isZero(double a, double absTol)
{
  return fabs(a) < absTol;
}

inline bool nequal(double a, double b, double relTol=GM_DOUBLECMP_RELTOL, double absTol=GM_DOUBLECMP_ABSTOL)
{
  return !equal(a, b, relTol, absTol);
}  
  
inline bool lt(double a, double b, double relTol=GM_DOUBLECMP_RELTOL, double absTol=GM_DOUBLECMP_ABSTOL)
{
  return (a < b) && !equal(a, b, relTol, absTol);
}

inline bool lte(double a, double b, double relTol=GM_DOUBLECMP_RELTOL, double absTol=GM_DOUBLECMP_ABSTOL)
{
  return (a < b) || equal(a, b, relTol, absTol);
}

inline bool gt(double a, double b, double relTol=GM_DOUBLECMP_RELTOL, double absTol=GM_DOUBLECMP_ABSTOL)
{
  return (a > b) && !equal(a, b, relTol, absTol);
}

inline bool gte(double a, double b, double relTol=GM_DOUBLECMP_RELTOL, double absTol=GM_DOUBLECMP_ABSTOL)
{
  return (a > b) || equal(a, b, relTol, absTol);
}

}

#endif