H-div Projection for discontinuous discrete functions. The projection is described in detail in: More...

#include <dune/fem/operator/projection/hdivprojection.hh>

Inheritance diagram for Dune::Fem::HdivProjection< DiscreteFunctionType >:

Public Types
typedef DiscreteFunctionType	DestinationType
	type of argument and destination

typedef DiscreteFunctionSpaceType	SpaceType
	convenience typedef for space type

typedef BaseType::RangeFunctionType	RangeFunctionType

typedef BaseType::DomainFunctionType	DomainFunctionType

typedef Dune::FieldTraits< RangeFieldType >::real_type	RealType

typedef BaseType::JacobianOperatorType	JacobianOperatorType

typedef RangeFunctionType::DiscreteFunctionSpaceType	RangeSpaceType

typedef DomainFunctionType::DiscreteFunctionSpaceType	DomainSpaceType

Public Member Functions
	HdivProjection (const DiscreteFunctionSpaceType &space)
	constructor taking space

	HdivProjection (const HdivProjection &org)

const DiscreteFunctionSpaceType &	space () const
	return reference to space

void	setTime (double)
	set time for operators

double	timeStepEstimate () const
	estimate maximum time step

virtual void	operator() (const DiscreteFunctionType &arg, DiscreteFunctionType &dest) const
	application operator projection arg to H-div space

virtual int	size () const
	return size of discrete function space, i.e. number of unknowns

virtual void	initializeTimeStepSize (const DestinationType &U0) const
	call operator once to calculate initial time step size

virtual bool	hasLimiter () const
	return true if limit method is implemented

virtual void	limit (const DestinationType &arg, DestinationType &dest) const
	limiter application operator

virtual void	applyLimiter (DestinationType &U) const
	limiter application operator

virtual const DestinationType *	destination () const
	return reference to pass's local memory

virtual void	jacobian (const DomainFunctionType &u, JacobianOperatorType &jOp) const

virtual void	jacobian (const DomainFunctionType &u, JacobianOperatorType &jOp) const =0
	obtain linearization

virtual void	operator() (const DomainFunctionType &u, RangeFunctionType &w) const=0
	application operator

virtual void	finalize ()
	finalization of operator

virtual bool	nonlinear () const

Static Public Member Functions
static double	normalJump (const DiscreteFunctionType &discFunc, const int polyOrder=-1)
	return sum of jumps of discrete function normal to intersection

template<class AdaptationType >
static void	estimator (const DiscreteFunctionType &velo, AdaptationType &adaptation)

Protected Attributes
std::unique_ptr< DestinationType >	uTmp_

Detailed Description

template<class DiscreteFunctionType>
class Dune::Fem::HdivProjection< DiscreteFunctionType >

H-div Projection for discontinuous discrete functions. The projection is described in detail in:

P. Bastian and B. Riviere. Superconvergence and H(div)-projection for discontinuous Galerkin methods. Int. J. Numer. Meth. Fluids., 42:1043-1057, 2003.

(see homepage of Peter Bastian: http://hal.iwr.uni-heidelberg.de/~peter/Papers/BDMpaper.pdf )

Note:

This projection only works for polynomial order 1 and the following spaces:

SimplexGrids + DiscontinuousGalerkinSpace
CubeGrids + LegendreDiscontinuousGalerkinSpace

Member Typedef Documentation

◆ DestinationType

typedef DiscreteFunctionType Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::DestinationType

inherited

type of argument and destination

◆ DomainFunctionType

template<class DomainFunction , class RangeFunction = DomainFunction, class LinearOperator = AutomaticDifferenceLinearOperator< DomainFunction, RangeFunction >>

typedef BaseType::DomainFunctionType Dune::Fem::AutomaticDifferenceOperator< DomainFunction, RangeFunction, LinearOperator >::DomainFunctionType

inherited

◆ DomainSpaceType

template<class DomainFunction , class RangeFunction = DomainFunction, class LinearOperator = AutomaticDifferenceLinearOperator< DomainFunction, RangeFunction >>

typedef DomainFunctionType::DiscreteFunctionSpaceType Dune::Fem::AutomaticDifferenceOperator< DomainFunction, RangeFunction, LinearOperator >::DomainSpaceType

inherited

◆ JacobianOperatorType

template<class DomainFunction , class RangeFunction = DomainFunction, class LinearOperator = AutomaticDifferenceLinearOperator< DomainFunction, RangeFunction >>

typedef BaseType::JacobianOperatorType Dune::Fem::AutomaticDifferenceOperator< DomainFunction, RangeFunction, LinearOperator >::JacobianOperatorType

inherited

◆ RangeFunctionType

template<class DomainFunction , class RangeFunction = DomainFunction, class LinearOperator = AutomaticDifferenceLinearOperator< DomainFunction, RangeFunction >>

typedef BaseType::RangeFunctionType Dune::Fem::AutomaticDifferenceOperator< DomainFunction, RangeFunction, LinearOperator >::RangeFunctionType

inherited

◆ RangeSpaceType

template<class DomainFunction , class RangeFunction = DomainFunction, class LinearOperator = AutomaticDifferenceLinearOperator< DomainFunction, RangeFunction >>

typedef RangeFunctionType::DiscreteFunctionSpaceType Dune::Fem::AutomaticDifferenceOperator< DomainFunction, RangeFunction, LinearOperator >::RangeSpaceType

inherited

◆ RealType

template<class DomainFunction , class RangeFunction = DomainFunction, class LinearOperator = AutomaticDifferenceLinearOperator< DomainFunction, RangeFunction >>

typedef Dune::FieldTraits<RangeFieldType>::real_type Dune::Fem::AutomaticDifferenceOperator< DomainFunction, RangeFunction, LinearOperator >::RealType

inherited

◆ SpaceType

typedef DiscreteFunctionSpaceType Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::SpaceType

inherited

convenience typedef for space type

Constructor & Destructor Documentation

◆ HdivProjection() [1/2]

template<class DiscreteFunctionType >

Dune::Fem::HdivProjection< DiscreteFunctionType >::HdivProjection ( const DiscreteFunctionSpaceType & space )

inline

constructor taking space

◆ HdivProjection() [2/2]

template<class DiscreteFunctionType >

Dune::Fem::HdivProjection< DiscreteFunctionType >::HdivProjection ( const HdivProjection< DiscreteFunctionType > & org )

inline

Member Function Documentation

◆ applyLimiter()

virtual void Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::applyLimiter ( DestinationType & U ) const

inlinevirtualinherited

limiter application operator

Parameters

[in,out] U argument and destination to apply Limiter(u), needs internal copying

Note: : Default implementation is to do nothing (hasLimiter == false)

◆ destination()

virtual const DestinationType * Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::destination ( ) const

inlinevirtualinherited

return reference to pass's local memory

◆ estimator()

template<class DiscreteFunctionType >

template<class AdaptationType >

static void Dune::Fem::HdivProjection< DiscreteFunctionType >::estimator	(	const DiscreteFunctionType &	velo,
		AdaptationType &	adaptation
	)

inlinestatic

◆ finalize()

virtual void Dune::Fem::Operator< JacobianOperator::DomainFunctionType , JacobianOperator::RangeFunctionType >::finalize ( )

inlinevirtualinherited

finalization of operator

Note: The default implementation is empty.

◆ hasLimiter()

virtual bool Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::hasLimiter ( ) const

inlinevirtualinherited

return true if limit method is implemented

Returns: true if limit is implemented

◆ initializeTimeStepSize()

void Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::initializeTimeStepSize ( const DestinationType & U0 ) const

inlinevirtualinherited

call operator once to calculate initial time step size

Parameters

U0	initial data to compute initial time step size

◆ jacobian() [1/2]

template<class DomainFunction , class RangeFunction = DomainFunction, class LinearOperator = AutomaticDifferenceLinearOperator< DomainFunction, RangeFunction >>

virtual void Dune::Fem::AutomaticDifferenceOperator< DomainFunction, RangeFunction, LinearOperator >::jacobian	(	const DomainFunctionType &	u,
		JacobianOperatorType &	jOp
	)		const

inlinevirtualinherited

◆ jacobian() [2/2]

template<class JacobianOperator >

virtual void Dune::Fem::DifferentiableOperator< JacobianOperator >::jacobian	(	const DomainFunctionType &	u,
		JacobianOperatorType &	jOp
	)		const

pure virtualinherited

obtain linearization

Parameters

[in]	u	argument discrete function
[out]	jOp	destination Jacobian operator

Note: This method has to be implemented by all derived classes.

Implemented in Dune::Fem::DifferentiableGalerkinOperator< Integrands, JacobianOperator >, Dune::Fem::MassLumpingDifferentiableOperator< Integrands, MassIntegrands, JacobianOperator >, and Dune::Fem::MOLDifferentiableGalerkinOperator< Integrands, JacobianOperator >.

◆ limit()

virtual void Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::limit	(	const DestinationType &	arg,
		DestinationType &	dest
	)		const

inlinevirtualinherited

limiter application operator

Parameters

arg	argument, u
dest	destination, Limiter(u)

Note: : Default implementation is to copy arg into dest.

◆ nonlinear()

virtual bool Dune::Fem::Operator< JacobianOperator::DomainFunctionType , JacobianOperator::RangeFunctionType >::nonlinear ( ) const

inlinevirtualinherited

Return true if the Operator is nonlinear and false otherwise (default is true).

Reimplemented in DirichletWrapperOperator< Operator, Constraints >, Dune::Fem::DifferentiableGalerkinOperator< Integrands, JacobianOperator >, and Dune::Fem::DifferentiableGalerkinOperator< ModelIntegrands, LinearOperator >.

◆ normalJump()

template<class DiscreteFunctionType >

static double Dune::Fem::HdivProjection< DiscreteFunctionType >::normalJump	(	const DiscreteFunctionType &	discFunc,
		const int	polyOrder = `-1`
	)

inlinestatic

return sum of jumps of discrete function normal to intersection

◆ operator()() [1/2]

template<class DiscreteFunctionType >

virtual void Dune::Fem::HdivProjection< DiscreteFunctionType >::operator()	(	const DiscreteFunctionType &	arg,
		DiscreteFunctionType &	dest
	)		const

inlinevirtual

application operator projection arg to H-div space

◆ operator()() [2/2]

virtual void Dune::Fem::Operator< JacobianOperator::DomainFunctionType , JacobianOperator::RangeFunctionType >::operator()	(	const DomainFunctionType &	u,
		RangeFunctionType &	w
	)		const

pure virtualinherited

application operator

Parameters

[in]	u	argument discrete function
[out]	w	destination discrete function

Note: This method has to be implemented by all derived classes.

◆ setTime()

template<class DiscreteFunctionType >

void Dune::Fem::HdivProjection< DiscreteFunctionType >::setTime ( double time )

inlinevirtual

set time for operators

Parameters

time	current time of evaluation

Reimplemented from Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType >.

◆ size()

virtual int Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::size ( ) const

inlinevirtualinherited

return size of discrete function space, i.e. number of unknowns

◆ space()

template<class DiscreteFunctionType >

const DiscreteFunctionSpaceType & Dune::Fem::HdivProjection< DiscreteFunctionType >::space ( ) const

inlinevirtual

return reference to space

Implements Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType >.

◆ timeStepEstimate()

template<class DiscreteFunctionType >

double Dune::Fem::HdivProjection< DiscreteFunctionType >::timeStepEstimate ( ) const

inlinevirtual

estimate maximum time step

For an explicit time discretization, the time step has to be limited. An estimate for the maximum time step of an explicit Euler scheme is returned by this function. Maximum time steps for higher order Runge Kutta schemes can be derived from this value.

Reimplemented from Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType >.

Member Data Documentation

◆ uTmp_

std::unique_ptr< DestinationType > Dune::Fem::SpaceOperatorInterface< DiscreteFunctionType , Fem::AutomaticDifferenceOperator< DiscreteFunctionType > >::uTmp_

mutableprotectedinherited

The documentation for this class was generated from the following file:

hdivprojection.hh

Public Types

Public Member Functions

Static Public Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

◆ DestinationType

◆ DomainFunctionType

◆ DomainSpaceType

◆ JacobianOperatorType

◆ RangeFunctionType

◆ RangeSpaceType

◆ RealType

◆ SpaceType

Constructor & Destructor Documentation

◆ HdivProjection() [1/2]

◆ HdivProjection() [2/2]

Member Function Documentation

◆ applyLimiter()

◆ destination()

◆ estimator()

◆ finalize()

◆ hasLimiter()

◆ initializeTimeStepSize()

◆ jacobian() [1/2]

◆ jacobian() [2/2]

◆ limit()

◆ nonlinear()

◆ normalJump()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ setTime()

◆ size()

◆ space()

◆ timeStepEstimate()

Member Data Documentation

◆ uTmp_