Finite Element Space that supports arbitrary, local degrees. More...

#include <lf/fe/fe.h>

Inheritance diagram for lf::fe::HierarchicScalarFESpace< SCALAR >:

Public Types
using	Scalar = SCALAR

Public Types inherited from lf::fe::ScalarFESpace< SCALAR >
using	Scalar = SCALAR

Public Member Functions
	HierarchicScalarFESpace ()=delete

	HierarchicScalarFESpace (const HierarchicScalarFESpace &)=delete

	HierarchicScalarFESpace (HierarchicScalarFESpace &&) noexcept=default

HierarchicScalarFESpace &	operator= (const HierarchicScalarFESpace &)=delete

HierarchicScalarFESpace &	operator= (HierarchicScalarFESpace &&) noexcept=default

	HierarchicScalarFESpace (const std::shared_ptr< const lf::mesh::Mesh > &mesh_p, unsigned degree)
	Construct a new Hierarchic FESpace with uniform polynomial degree. More...

template<class F , class = std::enable_if_t< std::is_invocable_v<F, const mesh::Entity &>>>
	HierarchicScalarFESpace (const std::shared_ptr< const lf::mesh::Mesh > &mesh_p, F &&degree_functor)
	Construct a new Hierarchic FESpace with (possibly) varying polynomial degrees. More...

std::shared_ptr< const lf::mesh::Mesh >	Mesh () const override
	access to underlying mesh More...

const lf::assemble::DofHandler &	LocGlobMap () const override
	access to associated local-to-global map More...

ScalarReferenceFiniteElement< SCALAR > const *	ShapeFunctionLayout (const lf::mesh::Entity &entity) const override
	access to shape function layout for cells More...

size_type	NumRefShapeFunctions (const lf::mesh::Entity &entity) const override
	number of interior shape functions associated to entities of various types More...

	~HierarchicScalarFESpace () override=default

Public Member Functions inherited from lf::fe::ScalarFESpace< SCALAR >
virtual std::shared_ptr< const lf::mesh::Mesh >	Mesh () const =0
	acess to underlying mesh More...

virtual const lf::assemble::DofHandler &	LocGlobMap () const =0
	access to associated local-to-global map More...

virtual ScalarReferenceFiniteElement< SCALAR > const *	ShapeFunctionLayout (const lf::mesh::Entity &entity) const =0
	access to shape function layout for mesh entities More...

virtual size_type	NumRefShapeFunctions (const lf::mesh::Entity &entity) const =0
	number of interior shape functions associated to a particular mesh entity. More...

virtual	~ScalarFESpace ()=default
	No special destructor. More...

Private Member Functions
template<class F >
void	Init (F &&degree_functor)

Private Attributes
std::shared_ptr< const lf::mesh::Mesh >	mesh_p_

quad::QuadRuleCache	qr_cache_

lf::mesh::utils::AllCodimMeshDataSet< std::variant< std::monostate, FePoint< SCALAR >, FeHierarchicSegment< SCALAR >, FeHierarchicTria< SCALAR >, FeHierarchicQuad< SCALAR > > >	ref_el_

lf::assemble::DynamicFEDofHandler	dofh_

Additional Inherited Members
Protected Member Functions inherited from lf::fe::ScalarFESpace< SCALAR >
	ScalarFESpace ()=default
	default constructor, needed by std::vector More...

	ScalarFESpace (const ScalarFESpace &)=default

	ScalarFESpace (ScalarFESpace &&) noexcept=default

ScalarFESpace &	operator= (const ScalarFESpace &)=default

ScalarFESpace &	operator= (ScalarFESpace &&) noexcept=default

Detailed Description

template<typename SCALAR>
class lf::fe::HierarchicScalarFESpace< SCALAR >

Finite Element Space that supports arbitrary, local degrees.

Template Parameters

SCALAR underlying scalar type, usually either double or complex<double>

This FE Space contains hierarchic Finite Elements, meaning that the function spaces of lower polynomial degree are contained in the higher order function spaces.

The polynomial degree can vary from entity to entity, i.e. local \(p\)-refinement is supported.

Example usage

The following code snippet computes the solution of the BVP

\begin{align} - \Delta u &= 1 && \text{on }\Omega := [0,1]^2 \\ u &= 0 && \text{on }\partial \Omega \end{align}

  // Get a mesh
  auto mesh_p = lf::mesh::test_utils::GenerateHybrid2DTestMesh();
  // Create HierarchicalFESpace with uniform degree 5.
  const unsigned degree = 5;
  const auto fe_space =
      std::make_shared<lf::fe::HierarchicScalarFESpace<double>>(mesh_p, degree);
 
  // define diffusion coefficient
  const lf::mesh::utils::MeshFunctionConstant mf_alpha(1);
  // define rhs load
  auto mf_load = lf::mesh::utils::MeshFunctionConstant(1.);
 
  // Assemble the system matrix and right hand side
  Eigen::VectorXd rhs = Eigen::VectorXd::Zero(fe_space->LocGlobMap().NumDofs());
  const assemble::DofHandler &dofh = fe_space->LocGlobMap();
  assemble::COOMatrix<double> A_COO(dofh.NumDofs(), dofh.NumDofs());
 
  DiffusionElementMatrixProvider element_matrix_provider(fe_space, mf_alpha);
  AssembleMatrixLocally(0, dofh, dofh, element_matrix_provider, A_COO);
  ScalarLoadElementVectorProvider element_vector_provider(fe_space, mf_load);
  AssembleVectorLocally(0, dofh, element_vector_provider, rhs);
 
  // Enforce zero dirichlet boundary conditions
 
  const auto boundary = lf::mesh::utils::flagEntitiesOnBoundary(mesh_p);
  const auto selector = [&](unsigned int idx) -> std::pair<bool, double> {
    const auto &entity = dofh.Entity(idx);
    return {entity.Codim() > 0 && boundary(entity), 0};
  };
  FixFlaggedSolutionComponents(selector, A_COO, rhs);
 
  // Solve the LSE using the cholesky decomposition
  Eigen::SparseMatrix<double> A = A_COO.makeSparse();
  Eigen::SimplicialLDLT<Eigen::SparseMatrix<double>> solver(A);
  const Eigen::VectorXd solution = solver.solve(rhs);
 
  // visualize the solution:
  io::VtkWriter vtk(mesh_p, "solution.vtk", 0, 5);
  vtk.WritePointData("solution", MeshFunctionFE(fe_space, solution));

Definition at line 43 of file hierarchic_scalar_fe_space.h.

Member Typedef Documentation

◆ Scalar

template<typename SCALAR >

using lf::fe::HierarchicScalarFESpace< SCALAR >::Scalar = SCALAR

Definition at line 45 of file hierarchic_scalar_fe_space.h.

Constructor & Destructor Documentation

◆ HierarchicScalarFESpace() [1/5]

template<typename SCALAR >

lf::fe::HierarchicScalarFESpace< SCALAR >::HierarchicScalarFESpace ( )

delete

◆ HierarchicScalarFESpace() [2/5]

template<typename SCALAR >

lf::fe::HierarchicScalarFESpace< SCALAR >::HierarchicScalarFESpace ( const HierarchicScalarFESpace< SCALAR > & )

delete

◆ HierarchicScalarFESpace() [3/5]

template<typename SCALAR >

lf::fe::HierarchicScalarFESpace< SCALAR >::HierarchicScalarFESpace ( HierarchicScalarFESpace< SCALAR > && )

defaultnoexcept

◆ HierarchicScalarFESpace() [4/5]

template<typename SCALAR >

lf::fe::HierarchicScalarFESpace< SCALAR >::HierarchicScalarFESpace	(	const std::shared_ptr< const lf::mesh::Mesh > &	mesh_p,
		unsigned	degree
	)

inlineexplicit

Construct a new Hierarchic FESpace with uniform polynomial degree.

Parameters

mesh_p	A shared pointer to the underlying mesh (immutable)
degree	The uniform polynomial degree.

Definition at line 59 of file hierarchic_scalar_fe_space.h.

References lf::base::RefEl::kPoint().

◆ HierarchicScalarFESpace() [5/5]

template<typename SCALAR >

template<class F , class = std::enable_if_t< std::is_invocable_v<F, const mesh::Entity &>>>

lf::fe::HierarchicScalarFESpace< SCALAR >::HierarchicScalarFESpace	(	const std::shared_ptr< const lf::mesh::Mesh > &	mesh_p,
		F &&	degree_functor
	)

inlineexplicit

Construct a new Hierarchic FESpace with (possibly) varying polynomial degrees.

Template Parameters

F	type of the `degree_functor`

Parameters

mesh_p	A shared pointer to the underlying mesh (immutable)
degree_functor	A function object that assigns a polynomial degree to every entity in the mesh. See below for more info.

Degree Functor

The degree_functor must overload the call operator as follows:

unsigned degree_functor(const lf::mesh::Entity& e)

lf::mesh::Entity

Interface class representing a topological entity in a cellular complex

Definition: entity.h:39

and should return the polynomial degree for the respective entity. The degree functor will be called for all entities (i.e. edges, triangles, quadrilaterals) in the mesh except for points.

Definition at line 88 of file hierarchic_scalar_fe_space.h.

◆ ~HierarchicScalarFESpace()

template<typename SCALAR >

lf::fe::HierarchicScalarFESpace< SCALAR >::~HierarchicScalarFESpace ( )

overridedefault

Member Function Documentation

◆ Init()

template<typename SCALAR >

template<class F >

void lf::fe::HierarchicScalarFESpace< SCALAR >::Init ( F && degree_functor )

inlineprivate

Definition at line 173 of file hierarchic_scalar_fe_space.h.

References lf::base::RefEl::kQuad(), and lf::base::RefEl::kTria().

◆ LocGlobMap()

template<typename SCALAR >

const lf::assemble::DofHandler & lf::fe::HierarchicScalarFESpace< SCALAR >::LocGlobMap ( ) const

inlineoverridevirtual

access to associated local-to-global map

Returns: a reference to the lf::assemble::DofHandler object (immutable)

Implements lf::fe::ScalarFESpace< SCALAR >.

Definition at line 124 of file hierarchic_scalar_fe_space.h.

◆ Mesh()

template<typename SCALAR >

std::shared_ptr< const lf::mesh::Mesh > lf::fe::HierarchicScalarFESpace< SCALAR >::Mesh ( ) const

inlineoverridevirtual

access to underlying mesh

Returns: a shared pointer to the mesh

Implements lf::fe::ScalarFESpace< SCALAR >.

Definition at line 116 of file hierarchic_scalar_fe_space.h.

◆ NumRefShapeFunctions()

template<typename SCALAR >

size_type lf::fe::HierarchicScalarFESpace< SCALAR >::NumRefShapeFunctions ( const lf::mesh::Entity & entity ) const

inlineoverridevirtual

number of interior shape functions associated to entities of various types

Implements lf::fe::ScalarFESpace< SCALAR >.

Definition at line 153 of file hierarchic_scalar_fe_space.h.

◆ operator=() [1/2]

template<typename SCALAR >

HierarchicScalarFESpace & lf::fe::HierarchicScalarFESpace< SCALAR >::operator= ( const HierarchicScalarFESpace< SCALAR > & )

delete

◆ operator=() [2/2]

template<typename SCALAR >

HierarchicScalarFESpace & lf::fe::HierarchicScalarFESpace< SCALAR >::operator= ( HierarchicScalarFESpace< SCALAR > && )

defaultnoexcept

◆ ShapeFunctionLayout()

template<typename SCALAR >

ScalarReferenceFiniteElement< SCALAR > const * lf::fe::HierarchicScalarFESpace< SCALAR >::ShapeFunctionLayout ( const lf::mesh::Entity & entity ) const

inlineoverridevirtual

access to shape function layout for cells

access to shape function layout for mesh entities

Parameters

entity The entity to get the reference element for

Warning: NULL pointers may be returned by this method in case a finite element specification was not given for a particular topological type of entity.