FEBase.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 "pde/finiteElements/FEBase.hpp"
#include "pde/finiteElements/FENode.hpp"
 
FEBase::FEBase(GeomCell *cell, CepsBool owned, ReferenceFE *reference) :
    ceps::HoldsGeomObject<GeomCell>(),
    ceps::HoldsNodes<FENode>(),
    ceps::HoldsCells<FEBase>(),
    m_reference(reference)
{
  CEPS_ABORT_IF(ceps::isNullPtr(cell),
    "Cannot create finite element as ptr to geometrical element is null"
  );
  CEPS_ABORT_IF(ceps::isNullPtr(reference),
    "Cannot create finite element as ptr to reference element is null"
  );
  // Check compatibility of reference element and geom cell
  CEPS_ABORT_IF(m_reference->getNumberOfGeomVertices() != cell->getNumberOfNodes()
             or m_reference->getDimension() != cell->getDimension(),
    "Cannot create finite element as the passed geometric cell and reference element\n" <<
    "   do not share the same geometrical properties"
  );
  this->setGeomObject(cell,owned);
}
 
FEBase::~FEBase()
{
}
 
void
FEBase::setGeomObject(GeomCell *cell, CepsBool owned)
{
  CEPS_ABORT_IF(ceps::isNullPtr(cell),
    "Cannot create finite element as ptr to geometrical element is null"
  );
  ceps::HoldsGeomObject<GeomCell>::setGeomObject(cell,owned);
  this->setNodes({});
  this->setCells({});
  return;
}
 
// --------------------------------------------------------------------------------
// Reference Element
// --------------------------------------------------------------------------------
 
ReferenceFE *
FEBase::getReferenceElement() const
{
  return m_reference;
}
 
CepsUInt
FEBase::getNumberOfDofs() const
{
  return getReferenceElement()->getNumberOfDofs();
}
 
const CepsVector<CepsReal3D> &
FEBase::getBasisVertices() const
{
  return getReferenceElement()->getBasisVertices();
}
 
const CepsReal3D &
FEBase::getBasisVertex(CepsUInt i) const
{
  return getReferenceElement()->getBasisVertex(i);
}
 
CepsMathVertex
FEBase::getOrigin() const
{
  return getNodeAt(0)->getProperties()->getCoordinatesForEigen();
}
 
CepsCellType
FEBase::getCellType() const
{
  return getReferenceElement()->getCellType();
}
 
CepsUInt
FEBase::getPolynomialOrder() const
{
  return getReferenceElement()->getPolynomialOrder();
}
 
// --------------------------------------------------------------------------------
// Geometrical
// --------------------------------------------------------------------------------
 
CepsReal
FEBase::getMeasure() const
{
  return getGeomObject()->getMeasure();
}
 
const GeomCell::JacobianMatrixType &
FEBase::getJacobianMatrix() const
{
  return getGeomObject()->getJacobianMatrix();
}
 
CepsReal
FEBase::getJacobianDeterminant() const
{
  return getGeomObject()->getJacobianDeterminant();
}
 
const GeomCell::InverseJacobianMatrixType &
FEBase::getInverseJacobian() const
{
  return getGeomObject()->getInverseJacobianMatrix();
}