Precision.hpp

Go to the documentation of this file.

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   This file is part of CEPS.
 
   CEPS is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.
 
   CEPS is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with CEPS (see file LICENSE at root of project).
   If not, see <https://www.gnu.org/licenses/>.
 
 
   Copyright 2019-2024 Inria, Universite de Bordeaux
 
   Authors, in alphabetical order:
 
     Pierre-Elliott BECUE, Florian CARO, Yves COUDIERE(*), Andjela DAVIDOVIC, 
     Charlie DOUANLA-LONTSI, Marc FUENTES, Mehdi JUHOOR, Michael LEGUEBE(*), 
     Pauline MIGERDITICHAN, Valentin PANNETIER(*), Nejib ZEMZEMI.
     * : currently active authors
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#pragma once
 
#include <cmath>
#include <limits>
 
#include "common/core/CepsTypes.hpp"
 
namespace ceps
{
 
  inline CepsBool
  equals(CepsReal a, CepsReal b, CepsReal error = 1.0)
  {
    return (std::fabs(a-b)<=std::abs(std::min(a,b))*std::numeric_limits<CepsReal>::epsilon()*error);
  }
  
  inline CepsBool
  approxEquals(CepsReal a, CepsReal b, CepsReal epsilon)
  {
    return (std::fabs(a-b)<epsilon);
  }
  
  inline CepsBool
  approxEquals(CepsReal *a, CepsReal *b, CepsUInt n, CepsReal epsilon)
  {
    CepsBool res = true;
    for (CepsUInt i=0U; i<n; i++)
      res = res and (std::fabs(a[i]-b[i])>=epsilon);
    return res;
  }
  
  inline CepsBool
  greaterThan(CepsReal a, CepsReal b, CepsReal epsilon = FLOATING_POINT_EPSILON)
  {
    return ((a-b) > ((std::fabs(a)<std::fabs(b)) ? std::fabs(b) : std::fabs(a))*epsilon);
  }
  
  inline CepsBool
  lesserThan(CepsReal a, CepsReal b, CepsReal epsilon = FLOATING_POINT_EPSILON)
  {
    return ((b-a) >((std::fabs(a) < std::fabs(b)) ? std::fabs(b) : std::fabs(a))*epsilon);
  }
 
}