GeomCell.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 "geometry/base/GeomNode.hpp"
 
class GeomCell :
    public ceps::HoldsGlobalIndex,
    public ceps::HoldsDimension,
    public ceps::HoldsBoundary,
    public ceps::HoldsAttributes,
    public ceps::HoldsProcIds,
    public ceps::HoldsNodes<GeomNode>,
    public ceps::HoldsCells<GeomCell>
{
 
  public:
 
    using JacobianMatrixType        = Eigen::Matrix<CepsReal, 3, Eigen::Dynamic>;
    using InverseJacobianMatrixType = Eigen::Matrix<CepsReal, Eigen::Dynamic, 3>;
  
    GeomCell();
 
    explicit GeomCell(CepsVector<GeomNode*> nodes, const CepsCellGlobalIndex& gID=0);
  
    GeomCell(const GeomCell &) = default;
  
    GeomCell(GeomCell &&) = default;
 
    GeomCell&
    operator=(const GeomCell &) = default;
  
    GeomCell&
    operator=(GeomCell &&) = default;
  
    ~GeomCell() override = default;
  
    // -------------------------------------------------------------------------
    // Nodes manipulation
    // -------------------------------------------------------------------------
 
    GeomNode*
    getNodeWithMinimalIndex() const;
  
    const CepsCellType&
    getCellType() const;
  
    virtual void
    reset();
  
    // --------------------------------------------------------------------------------
    // Jacobian
    // --------------------------------------------------------------------------------
  
    CepsReal
    getJacobianDeterminant() const;
  
    virtual const JacobianMatrixType&
    getJacobianMatrix();
  
    virtual const InverseJacobianMatrixType&
    getInverseJacobianMatrix();
  
    virtual void
    refresh();
 
    void
    setNormal(CepsMathVertex n);
 
    CepsMathVertex 
    getNormal() const;
 
    friend class TestGeomCells;
 
    void
    setNodeSharingCellsIndices(const CepsSet<CepsCellGlobalIndex>& nCellIds);
 
    const CepsSet<CepsCellGlobalIndex>&
    getNodeSharingCellsIndices() const;
 
    CepsReal
    getDiameter();
 
 
  protected:
 
    void
    initialize(CepsUInt dim, CepsCellType cellType, CepsUInt nNodes);
  
    virtual void
    initializeNodes(GeomNode *const *nodes, CepsUInt nNodes);
 
    void
    computeDiameter();
 
  // ===================================================================
  // pure virtual funcs
 
  public:
 
    virtual CepsReal
    getMeasure() const = 0;
 
    virtual CepsUInt
    getNumberOfNodes() const = 0;
 
 
  protected:
 
    virtual void
    computeJacobianMatrix() = 0;
  
    virtual void
    computeInverseJacobianMatrix() = 0;
  
    virtual void
    computeJacobianDeterminant() = 0;
  
    virtual CepsBool
    checkJacobianDeterminant() = 0;
  
 
  protected:
 
    CepsCellType              m_cellType;      
    JacobianMatrixType        m_jacobian;      
    InverseJacobianMatrixType m_invJacobian;   
    CepsReal                  m_detJacobian;   
    CepsBool                  m_invJcomputed;  
    CepsMathVertex            m_normal;        
    CepsReal                  m_diameter;      
    CepsSet<CepsCellGlobalIndex> m_nodeSharingCells; 
};
 
CepsStandardArgs
getStandardArgsFrom(GeomCell* cell);