CepsTypeTools.cpp

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
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#include "common/core/CepsTypeTools.hpp"
#include "common/core/CepsIoMacros.hpp"
#include "common/core/CepsParallelTools.hpp"
 
CepsBool
operator<(const CepsHash3& h31, const CepsHash3& h32)
{
  return h31.h1<h32.h1 or (h31.h1==h32.h1 and 
        (h31.h2<h32.h2 or (h31.h2==h32.h2 and 
         h31.h3<h32.h3)));
}
 
std::ostream&
operator<<(std::ostream& os, const CepsHash3& h)
{
  os << "(" << h.h1 << ", " << h.h2 << ", " << h.h3 << ")";
  return os;
}
 
CepsHash3
ceps::hash::getHash3(const CepsReal3D& xyz)
{
  CepsHash3 res;
  res.h1 = CepsSHash(std::round(xyz[0]/CEPS_HASH3_PRECISION));
  res.h2 = CepsSHash(std::round(xyz[1]/CEPS_HASH3_PRECISION));
  res.h3 = CepsSHash(std::round(xyz[2]/CEPS_HASH3_PRECISION));
  return res;
}
 
void
ceps::checkNanOrInf(CepsReal value, CepsString message)
{
  CEPS_ABORT_IF(std::isnan(value) or std::isinf(value),
    "Detected NaN or Infty: " << message
  );
}