Computes the Galerkin LSE for the Hodge Laplacian of the whitney two form. More...

#include </home/nico/bildung/SemVI/thesis/lehrfempp/projects/hldo_sphere/operators/whitney_two_hodge_laplacian.h>

Public Member Functions
	WhitneyTwoHodgeLaplace ()
	Constructor initializes basic mesh (Octaeder with radius 1.0) initializes zerovalued function f. More...

void	Compute ()
	Computes the Galerkin LSE. More...

void	SetMesh (std::shared_ptr< const lf::mesh::Mesh > mesh_p)
	Sets the mesh and creates dof_handler. More...

void	SetLoadFunction (std::function< SCALAR(const Eigen::Matrix< double, 3, 1 > &)> &f)
	Sets the load function. More...

Eigen::Matrix< SCALAR, Eigen::Dynamic, 1 >	GetLoadVector ()
	returns the Loadvector More...

lf::assemble::COOMatrix< SCALAR >	GetGalerkinMatrix ()
	returns the Galerkin Matrix More...

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

std::function< SCALAR(const Eigen::Matrix< double, 3, 1 > &)>	f_

lf::assemble::COOMatrix< SCALAR >	coo_matrix_

Eigen::Matrix< SCALAR, Eigen::Dynamic, 1 >	phi_

Detailed Description

template<typename SCALAR>
class projects::hldo_sphere::operators::WhitneyTwoHodgeLaplace< SCALAR >

Computes the Galerkin LSE for the Hodge Laplacian of the whitney two form.

\[ -\Delta_2 = -\text{div}_{\Gamma} \circ \mathbf{grad}_{\Gamma} \\ -\Delta_2 u = f \]

Basis functions used are the rotated Whitney 1-forms and the cellwise constant functions

Details regarding the mathematical derivations can be found in the thesis Hodge-Laplacians and Dirac Operators on the Surface of the 3-Sphere section 4.3.3.

Note: Only triangular meshes are supported

Definition at line 54 of file whitney_two_hodge_laplacian.h.

Constructor & Destructor Documentation

◆ WhitneyTwoHodgeLaplace()

template<typename SCALAR >

projects::hldo_sphere::operators::WhitneyTwoHodgeLaplace< SCALAR >::WhitneyTwoHodgeLaplace ( )

inline

Constructor initializes basic mesh (Octaeder with radius 1.0) initializes zerovalued function f.

Definition at line 62 of file whitney_two_hodge_laplacian.h.

Member Function Documentation

◆ Compute()

template<typename SCALAR >

void projects::hldo_sphere::operators::WhitneyTwoHodgeLaplace< SCALAR >::Compute ( )

inline

Computes the Galerkin LSE.

The Galerkin matrix will be accessable with GetGalerkinMatrix() The load vector will be accessable with GetLoadVector()

Definition at line 88 of file whitney_two_hodge_laplacian.h.

References lf::assemble::AssembleVectorLocally(), lf::base::RefEl::kSegment(), lf::base::RefEl::kTria(), lf::assemble::DofHandler::NumDofs(), lf::assemble::COOMatrix< SCALAR >::setZero(), and lf::assemble::COOMatrix< SCALAR >::triplets().

Referenced by projects::hldo_sphere::operators::HodgeLaplaciansSourceProblems< SCALAR >::Compute().

◆ GetGalerkinMatrix()

template<typename SCALAR >

lf::assemble::COOMatrix< SCALAR > projects::hldo_sphere::operators::WhitneyTwoHodgeLaplace< SCALAR >::GetGalerkinMatrix ( )

inline

returns the Galerkin Matrix

This is the Matrix of the LSE containing the negative laplacian

Note: The Galerkin matrix must be computed with Compute() before calling this funciton

Definition at line 234 of file whitney_two_hodge_laplacian.h.

Referenced by projects::hldo_sphere::operators::HodgeLaplaciansSourceProblems< SCALAR >::Compute().

◆ GetLoadVector()

template<typename SCALAR >

Eigen::Matrix< SCALAR, Eigen::Dynamic, 1 > projects::hldo_sphere::operators::WhitneyTwoHodgeLaplace< SCALAR >::GetLoadVector ( )

inline

returns the Loadvector

This is the righthandside of the LSE

Note: The loadvector must be computed with Compute() before calling this function

Definition at line 223 of file whitney_two_hodge_laplacian.h.

Referenced by projects::hldo_sphere::operators::HodgeLaplaciansSourceProblems< SCALAR >::Compute().

◆ SetLoadFunction()

template<typename SCALAR >

void projects::hldo_sphere::operators::WhitneyTwoHodgeLaplace< SCALAR >::SetLoadFunction ( std::function< SCALAR(const Eigen::Matrix< double, 3, 1 > &)> & f )

inline

Sets the load function.

Parameters