GeometryReader.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 "common/CepsCommon.hpp"
#include "geometry/base/GeomSimplex.hpp"
#include "geometry/readers/CoupledNodesReader.hpp"
#include "geometry/readers/GmshMeshReader.hpp"
#include "geometry/readers/MeditMeshReader.hpp"
#include "geometry/readers/TetgenMeshReader.hpp"
#include "geometry/readers/VtkMeshReader.hpp"
 
namespace MeshReaderType
{
  enum Type
  {
    VOLUMIC  = 0,  
    SURFACIC = 1,  
    CABLE    = 2   
  };
}  // namespace MeshReaderType
 
// --------------------------------------------------------
// The next class can be later completed with meshes that have quads
// Add corresponding vectors of volumicQuadMeshes, surfacicMeshes, etc...
// volumicTriMeshes, ...
// and rewrite computeCellMemory with an option for quads
// Normally, the mesh readers are ready for quads.
// --------------------------------------------------------
 
class GeometryReader
{
public:
  GeometryReader (
    const CepsVector<CepsString> &volumicMeshes,
    const CepsVector<CepsString> &surfacicMeshes,
    const CepsVector<CepsString> &cableMeshes,
    const CepsString             &coupledNodesFile = ""
  );
 
  ~GeometryReader ();
 
  // Cell/node count
 
  void
  getCellOffsetsPerMeshType (CepsUInt cellOffsets[4]);
 
  void
  getCellOffsetPerMesh (CepsVector<CepsUInt> &offsets);
 
  void
  getNodeOffsetsPerMeshType (CepsUInt offsets[4]);
 
  void
  getNodeOffsetPerMesh (CepsVector<CepsUInt> &offsets);
 
  CepsUInt
  getNbNodes ();
 
  CepsUInt
  getNbCells ();
 
  CepsUInt
  getNbBdryCells ();
 
  CepsUInt
  readNodes (
    CepsUInt                   indexStart,
    CepsUInt                   nbNodesToRead,
    CepsVector<CepsReal>      &coords,
    CepsVector<CepsAttribute> &attr,
    CepsVector<CepsUInt>      &attrOffset,
    CepsUInt                  &attrSize
  );
 
  CepsUInt
  readCells (
    CepsUInt                         indexStart,
    CepsUInt                         nbCellsToRead,
    CepsVector<CepsNodeGlobalIndex> &nodeIndices,
    CepsVector<CepsUInt>            &nodeIndicesOffset,
    CepsVector<CepsAttribute>       &attr,
    CepsVector<CepsUInt>            &attrOffset,
    CepsVector<CepsChar>            &isBoundary,
    CepsUInt                        &eindSize,
    CepsUInt                        &attrSize
  );
 
  const CepsMultiMap<CepsNodeGlobalIndex, CepsNodeGlobalIndex> &
  getCoupledNodes ();
 
private:
  MeshReader *
  getMeshReaderFor (const CepsString &meshName, CepsInt index, CepsUInt dim);
 
  void
  initializeDataStructures ();
 
  void
  readCoupledNodes (const CepsString &fileName);
 
  CepsUInt
  computeNumberOfNodeIndices (CepsUInt nbCells[3], CepsUInt nbBoundary[3]);
 
protected:
  CepsUInt m_nbMeshes;  
 
  // Mesh readers
  MeshReader           **m_readers;  
  CepsVector<CepsString> m_meshFileNames;  
  CepsUInt              *m_readerType;     
 
  // Coupling meshes
  CepsBool m_hasCoupledNodes;  
  CepsMultiMap<CepsNodeGlobalIndex, CepsNodeGlobalIndex> m_coupledNodes;  
 
  // Mesh data
  CepsUInt *m_meshNbNodes;      
  CepsUInt *m_nodeOffsets;      
  CepsUInt *m_meshNbCells;      
  CepsUInt *m_meshNbBdryCells;  
  CepsUInt *m_cellOffsets;      
 
  // All meshes data
  CepsUInt m_totalNbNodes;      
  CepsUInt m_totalNbCells;      
  CepsUInt m_totalNbBdryCells;  
};