ReferenceFE.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/math/CepsPolynomial.hpp"
#include "pde/finiteElements/FEType.hpp"
 
class ReferenceFE
{
 
 
  public:
    ReferenceFE() = delete;
  
    ReferenceFE(CepsUInt dimension, CepsFeLagrangeType type);
  
    ReferenceFE(const ReferenceFE& that) = default;
  
    ReferenceFE &
    operator=(const ReferenceFE& that) = default;
  
    ~ReferenceFE();
  
    CepsCellType
    getCellType() const;
  
    CepsUInt
    getPolynomialOrder() const;
  
    CepsUInt
    getNumberOfDofs() const;
  
    CepsMathDynamic1D
    computeBasisFunctions(const CepsReal3D& point);
  
    CepsMathDynamic2D
    computeBasisFunctionDerivatives(const CepsReal3D& point);
  
    const CepsVector<CepsReal3D>&
    getBasisVertices() const;
  
    const CepsReal3D&
    getBasisVertex(CepsUInt i) const;
  
    CepsInt
    getBasisVertexGeomCellIndex(CepsUInt i) const;
  
    GeomNode*
    getNodeOfOwnerForBasisVertex(CepsUInt i, GeomCell* cell) const;
  
    CepsReal3D
    getTransformedBasisVertex(CepsUInt i, GeomCell* cell) const;
  
    CepsUInt
    getNumberOfGeomVertices() const;
  
    CepsUInt
    getDimension() const;
  
    CepsUInt
    getNumberOfBasisVertices() const;
  
    CepsUInt
    getNumberOfBasisFunctions() const;
  
    CepsFeLagrangeType
    getFEType() const;
  
    CepsReal
    computeBasisFunction(CepsUInt iPhi, const CepsReal3D& point);
  
    CepsReal3D
    computeBasisFunctionDerivative(CepsUInt iPhi, const CepsReal3D& point);
 
 
  protected:
  
    void
    buildBasisFunction();
 
 
  protected :
 
    CepsFeLagrangeType         m_type;            
    CepsUInt                   m_dim;             
    CepsUInt                   m_order;           
    CepsUInt                   m_nbDofs;          
    CepsVector<Polynomial<3U>> m_basisFunctions;  
    CepsVector<CepsReal3D>     m_vertices;        
    CepsVector<CepsInt>        m_geomCellIndex;
    
    CepsUInt                   m_nGeomVertices;   
 
    CepsVector<CepsSet<CepsUInt>> m_ownerDeductionRule;
};