Inheritance diagram for Esfem::Impl::MCF_op:

Collaboration diagram for Esfem::Impl::MCF_op:

Public Types
using	Scalar_fef = Esfem::Grid::Scal_FEfun::Dune_FEfun
	$f\colon \R^3 \to \R$

using	Vector_fef = Esfem::Grid::Vec_FEfun::Dune_FEfun
	$f\colon \R^3 \to \R^3$

template<typename T >
using	Local_function = typename T::LocalFunctionType
	Finite element function restricted to a triangle.

template<typename T >
using	Domain = typename Local_function< T >::DomainType

template<typename T >
using	Range = typename Local_function< T >::RangeType

using	RangeField = typename Local_function< Vector_fef >::RangeFieldType
	Simply $\R$ .

using	Geometry = Vector_fef::DiscreteFunctionSpaceType::IteratorType::Entity::Geometry
	Geometry means access to the finite element.

using	Grid_part = Vector_fef::GridPartType
	Auxiliary template argument.

using	Quadrature = Dune::Fem::CachingQuadrature< Grid_part, 0 >

template<typename T >
using	Jacobian_range = typename Local_function< T >::JacobianRangeType
	Jacobian of the local finite element function. More...

Public Member Functions
	MCF_op (const Io::Parameter &, const Grid::Grid_and_time &, const Scalar_fef &u_input)

	MCF_op (const MCF_op &)=delete

MCF_op &	operator= (const MCF_op &)=delete

void	operator() (const Vector_fef &rhs, Vector_fef &lhs) const override
	Applying the mean curvature operator. The cg solver uses this. More...

void	brusselator_rhs (const Vector_fef &rhs, Vector_fef &lhs) const
	Generates rhs for the linear system. More...

Static Public Attributes
Variable containing dimensions
static constexpr int	dim_vec_domain = Domain<Vector_fef>::dimension

static constexpr int	dim_vec_range = Range<Vector_fef>::dimension

static constexpr int	dim_scalar_domain = Domain<Scalar_fef>::dimension

static constexpr int	dim_scalar_range = Range<Scalar_fef>::dimension

Private Member Functions
void	mcf_lhs_matrixFree_assembly (const Geometry &, const Quadrature &, const Local_function< Vector_fef > &, Local_function< Vector_fef > &) const
	Used by operator() More...

void	mcf_rhs_matrixFree_assembly (const Geometry &, const Quadrature &, const Local_function< Vector_fef > &, const Local_function< Scalar_fef > &u_loc, Local_function< Vector_fef > &) const
	Used by rhs() More...

Range< Vector_fef >	surface_normal (const Geometry &) const
	Calculates the outwards pointing normal vector.

Private Attributes
const double	alpha
	Velocity Laplace regularization parameter by Lubich.

const double	delta
	Tumor growth parameter.

const double	epsilon
	Mean curvature regularization parameter by Elliott.

const double	eps
	Generic tolerance.

const Dune::Fem::TimeProviderBase &	tp
	Time step provider.

const Scalar_fef &	u
	Concentration of the growth promoting chemical substance.

Detailed Description

Member Typedef Documentation

template<typename T >

Jacobian of the local finite element function.

Generates rhs for the linear system.

The new value of the finite element function will be

$\begin{equation*} (M_3^n + \alpha A_3^n) \nodalValue{X}^n + \tau \delta M_3^n(\nodalValue{u}^n, \nodalValue{\surfaceNormal}). \end{equation*}$

void MCF_op::mcf_lhs_matrixFree_assembly	(	const Geometry &	g,
		const Quadrature &	q,
		const Local_function< Vector_fef > &	cf,
		Local_function< Vector_fef > &	f
	)		const

private

Used by operator()

The matrix vector formulation reads as

$\begin{equation*} \parentheses[\big]{M + (\alpha + \varepsilon \tau) A} \nodalValue{X} \end{equation*}$

void MCF_op::mcf_rhs_matrixFree_assembly	(	const Geometry &	g,
		const Quadrature &	q,
		const Local_function< Vector_fef > &	cf,
		const Local_function< Scalar_fef > &	u_loc,
		Local_function< Vector_fef > &	f
	)		const