MeshReader.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 "geometry/readers/MeshReader.hpp"
 
MeshReader::MeshReader (const CepsString &mesh, const CepsUInt &dim, CepsBool quads) :
    FileReader (mesh),
    m_expectedDim (dim),
    m_quads (quads)
{
  if (m_expectedDim == 0 or m_expectedDim > 3)
  {
    CEPS_ABORT ("Impossible to read meshes with main dimension 0 or >3");
  }
  if (not ceps::fileExists (ceps::getFilename (mesh), ceps::getDir (mesh)))
  {
    CEPS_ABORT ("Cannot find file " << mesh);
  }
 
  m_cellsRead         = 0U;
  m_bdryCellsRead     = 0U;
  m_nbNodes           = 0U;
  m_nbCells           = 0U;
  m_nbBdryCells       = 0U;
  m_nbAttrPerNode     = 0U;
  m_nbAttrPerCell     = 0U;
  m_bdryCellsBegin    = 0U;
  m_cellsBegin        = 0U;
  m_nodesBegin        = 0U;
  m_nbAttrPerBdryCell = 0U;
  m_initialized       = false;
 
  if (m_quads)
  {
    m_nbNodesPerCell     = pow (2, dim);
    m_nbNodesPerBdryCell = pow (2, dim - 1);
  }
  else
  {
    m_nbNodesPerCell     = dim + 1;
    m_nbNodesPerBdryCell = dim;
  }
}
 
MeshReader::~MeshReader ()
{}
 
CepsUInt
MeshReader::getNbCellsRead () const
{
  return m_cellsRead;
}
 
CepsUInt
MeshReader::getNbBdryCellsRead () const
{
  return m_bdryCellsRead;
}
 
CepsUInt
MeshReader::getNbNodes () const
{
  return m_nbNodes;
}
 
CepsUInt
MeshReader::getNbCells () const
{
  return m_nbCells;
}
 
CepsUInt
MeshReader::getNbBdryCells () const
{
  return m_nbBdryCells;
}
 
CepsInt
MeshReader::getNbAttributesPerNode () const
{
  return m_nbAttrPerNode;
}
 
CepsInt
MeshReader::getNbAttributesPerCell () const
{
  return m_nbAttrPerCell;
}
 
CepsInt
MeshReader::getNbAttributesPerBdryCell () const
{
  return m_nbAttrPerBdryCell;
}
 
CepsBool
MeshReader::checkExtension () const
{
  if (ceps::getExtension (m_fileName) != m_extension)
  {
    CEPS_ABORT (
      "Mesh reader created with filename "
      << m_fileName << " but is only compatible with extensions " << m_extension
    );
    return false;
  }
  return true;
}
 
void
MeshReader::checkCanReadNodes (CepsUInt indexStart, CepsUInt indexEnd) const
{
  CEPS_ABORT_IF (not m_initialized,
    "Must call initialize() beforehand"
  );
  CEPS_ABORT_IF (indexStart > indexEnd,
    "index start > index end"
  );
  CEPS_ABORT_IF (indexStart > m_nbNodes or indexEnd > m_nbNodes,
    "index start or index end larger than number of nodes"
  );
}
 
void
MeshReader::checkCanReadCells (CepsUInt indexStart, CepsUInt indexEnd) const
{
  CEPS_ABORT_IF (not m_initialized, 
    "Must call initialize() beforehand"
  );
  CEPS_ABORT_IF (indexStart > indexEnd,
    "index start > index end"
  );
  CEPS_ABORT_IF (indexStart > m_nbCellsTotal or indexEnd > m_nbCellsTotal,
    "Cannot read cells from index #" << indexStart << std::endl
    << "   to index #" << indexEnd << " in file " << m_fileName
  );
}