#include <disc_constraint_fvcr.h>

Inheritance diagram for ug::DiscConstraintFVCR< TGridFunction >:

Public Types
typedef TAlgebra	algebra_type
	Algebra type. More...

typedef Attachment< dimMat >	AMathDimMat

typedef Attachment< std::vector< MathVector< dim > > >	ANumberArray

typedef PeriodicAttachmentAccessor< side_type, AMathDimMat >	aSideDimMat

typedef PeriodicAttachmentAccessor< side_type, ANumber >	aSideNumber

typedef PeriodicAttachmentAccessor< side_type, ANumberArray >	aSideNumberArray

typedef PeriodicAttachmentAccessor< side_type, ASizetArray >	aSideSizetArray

typedef Attachment< std::vector< DoFIndex > >	ASizetArray

typedef MathMatrix< dim, dim >	dimMat

typedef TDomain	domain_type
	Type of Domain. More...

typedef TGridFunction::template dim_traits< dim >::grid_base_object	elem_type
	element type More...

typedef TGridFunction::template dim_traits< dim >::const_iterator	ElemIterator
	element iterator More...

typedef face_type_traits< dim >::face_type0	face_type0

typedef face_type_traits< dim >::face_type1	face_type1

typedef domain_traits< TDomain::dim >::grid_base_object	grid_base_object
	Type of geometric base object. More...

typedef domain_type::grid_type	grid_type
	grid type More...

typedef std::pair< MathVector< dim >, MathVector< dim > >	MathVector_Pair

typedef algebra_type::matrix_type	matrix_type
	Type of algebra matrix. More...

typedef domain_type::position_accessor_type	position_accessor_type
	position accessor More...

typedef elem_type::side	side_type
	side type More...

typedef TGridFunction::template traits< side_type >::const_iterator	SideIterator
	side iterator More...

typedef TGridFunction::algebra_type	TAlgebra

typedef TGridFunction::domain_type	TDomain

typedef algebra_type::vector_type	vector_type
	Type of algebra vector. More...

typedef std::vector< std::pair< DoFIndex, MathVector< dim > > >	vIndexPosPair

typedef std::vector< std::vector< std::pair< DoFIndex, MathVector< dim > > > >	vvIndexPosPair

Public Types inherited from ug::IDomainConstraint< TGridFunction::domain_type, TGridFunction::algebra_type >
typedef TAlgebra	algebra_type

typedef TDomain	domain_type

typedef algebra_type::matrix_type	matrix_type

typedef algebra_type::vector_type	vector_type

Public Types inherited from ug::IConstraint< class >
typedef TAlgebra	algebra_type

typedef algebra_type::matrix_type	matrix_type

typedef algebra_type::vector_type	vector_type

Public Member Functions
virtual void	add_defect (vector_type &d, const vector_type &u, ConstSmartPtr< DoFDistribution > dd, const number time=0.0, const number s_a=1.0)

virtual void	add_pressure_defect (vector_type &d, const vector_type &u, ConstSmartPtr< DoFDistribution > dd, const number time=0.0, const number s_a=1.0)

virtual void	adjust_jacobian (matrix_type &J, const vector_type &u, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0, ConstSmartPtr< VectorTimeSeries< vector_type > > vSol=NULL, const number s_a0=1.0)

void	compute_grad_shapes ()
	compute gradient shapes for velocity on rotated elements More...

	DiscConstraintFVCR (SmartPtr< TGridFunction > u)
	constructor More...

	DiscConstraintFVCR (SmartPtr< TGridFunction > u, bool bLinUpConvDefect, bool bLinUpConvJacobian, bool bLinPressureDefect, bool bLinPressureJacobian, bool bAdaptive)

	DiscConstraintFVCR (SmartPtr< TGridFunction > u, bool bLinUpConvDefect, bool bLinUpConvJacobian, bool bLinPressureDefect, bool bLinPressureJacobian, bool bAdaptive, bool bLimiter)

	DiscConstraintFVCR (SmartPtr< TGridFunction > u, bool bLinUpConvDefect, bool bLinUpConvJacobian, bool bLinPressureDefect, bool bLinPressureJacobian, bool bAdaptive, bool bLimiter, const char *subsets)

	DiscConstraintFVCR (SmartPtr< TGridFunction > u, bool bLinUpConvDefect, bool bLinUpConvJacobian, bool bLinPressureDefect, bool bLinPressureJacobian, bool bAdaptive, const char *subsets)

void	init (SmartPtr< TGridFunction > u, bool bLinUpConvDefect, bool bLinUpConvJacobian, bool bLinPressureDefect, bool bLinPressureJacobian, bool bAdaptive, bool bLimiter)

void	set_limiter (bool bLimiter)

void	set_zero_grad_bnd (const char *subsets)
	set boundaries, in associated elements there is no linear pressure and no linear velocity upwind More...

bool	zeroGradBndElem (typename grid_type::template traits< side_type >::secure_container sides)

	~DiscConstraintFVCR ()
	destructor More...


virtual void	adjust_defect (vector_type &d, const vector_type &u, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0, ConstSmartPtr< VectorTimeSeries< vector_type > > vSol=SPNULL, const std::vector< number > vScaleMass=NULL, const std::vector< number > vScaleStiff=NULL)


virtual void	adjust_linear (matrix_type &mat, vector_type &rhs, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0)


virtual void	adjust_rhs (vector_type &rhs, const vector_type &u, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0)


virtual void	adjust_solution (vector_type &u, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0)

virtual int	type () const
	returns the type of the constraints More...

Public Member Functions inherited from ug::IDomainConstraint< TGridFunction::domain_type, TGridFunction::algebra_type >
SmartPtr< ApproximationSpace< TDomain > >	approximation_space ()

ConstSmartPtr< ApproximationSpace< TDomain > >	approximation_space () const

virtual SmartPtr< IErrEstData< TDomain > >	err_est_data ()

bool	err_est_enabled () const

	IDomainConstraint ()

virtual void	set_approximation_space (SmartPtr< ApproximationSpace< TDomain > > approxSpace)

void	set_ass_tuner (ConstSmartPtr< AssemblingTuner< TAlgebra > > spAssemblingTuner=NULL)

void	set_error_estimator (SmartPtr< IErrEstData< TDomain > > ee)

Public Member Functions inherited from ug::IConstraint< class >
virtual void	adjust_correction (vector_type &c, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0)

virtual void	adjust_error (const vector_type &u, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0, ConstSmartPtr< VectorTimeSeries< vector_type > > vSol=SPNULL, const std::vector< number > vScaleMass=NULL, const std::vector< number > vScaleStiff=NULL)

virtual void	adjust_linear_residual (vector_type &d, const vector_type &u, ConstSmartPtr< DoFDistribution > dd, int type, number time=0.0)

virtual void	adjust_prolongation (matrix_type &P, ConstSmartPtr< DoFDistribution > ddFine, ConstSmartPtr< DoFDistribution > ddCoarse, int type, number time=0.0)

virtual void	adjust_prolongation (vector_type &uFine, GridLevel fineLvl, const vector_type &uCoarse, GridLevel coarseLvl, int type)

virtual void	adjust_restriction (matrix_type &R, ConstSmartPtr< DoFDistribution > ddCoarse, ConstSmartPtr< DoFDistribution > ddFine, int type, number time=0.0)

virtual void	adjust_restriction (vector_type &uCoarse, GridLevel coarseLvl, const vector_type &uFine, GridLevel fineLvl, int type)

virtual void	modify_solution (SmartPtr< VectorTimeSeries< vector_type > > vSolMod, ConstSmartPtr< VectorTimeSeries< vector_type > > vSol, ConstSmartPtr< DoFDistribution > dd, int type)

virtual void	modify_solution (vector_type &uMod, const vector_type &u, ConstSmartPtr< DoFDistribution > dd, int type)

virtual	~IConstraint ()

Public Attributes
aSideDimMat	acGrad

aSideNumberArray	acGradSh

aSideSizetArray	acGradShInd

aSideNumber	acVol

AMathDimMat	aGrad

ANumberArray	aGradSh

ASizetArray	aGradShInd

ANumber	aVol

Static Public Attributes
static const size_t	_P_ = dim

static const int	blockSize = algebra_type::blockSize
	blockSize of used algebra More...

static const int	dim = TDomain::dim
	world Dimension More...

static const size_t	maxShapeSize = 2*DimCRFVGeometry<dim>::maxNumSCV-1

Private Attributes
bool	m_bAdaptive

bool	m_bLinPressureDefect

bool	m_bLinPressureJacobian

bool	m_bLinUpConvDefect

bool	m_bLinUpConvJacobian

grid_type *	m_grid

ISubsetHandler *	m_ish

bool	m_limiter

SmartPtr< TGridFunction >	m_u

SubsetGroup	m_zeroGradSg

Additional Inherited Members
Protected Member Functions inherited from ug::IDomainConstraint< TGridFunction::domain_type, TGridFunction::algebra_type >
ConstSmartPtr< DoFDistribution >	dd (const GridLevel &gl) const

Protected Attributes inherited from ug::IDomainConstraint< TGridFunction::domain_type, TGridFunction::algebra_type >
SmartPtr< ApproximationSpace< TDomain > >	m_spApproxSpace

ConstSmartPtr< AssemblingTuner< TAlgebra > >	m_spAssTuner

SmartPtr< IErrEstData< TDomain > >	m_spErrEstData

Member Typedef Documentation

◆ algebra_type

template<typename TGridFunction >

typedef TAlgebra ug::DiscConstraintFVCR< TGridFunction >::algebra_type

Algebra type.

◆ AMathDimMat

template<typename TGridFunction >

typedef Attachment<dimMat> ug::DiscConstraintFVCR< TGridFunction >::AMathDimMat

◆ ANumberArray

template<typename TGridFunction >

typedef Attachment<std::vector< MathVector<dim> > > ug::DiscConstraintFVCR< TGridFunction >::ANumberArray

◆ aSideDimMat

template<typename TGridFunction >

typedef PeriodicAttachmentAccessor<side_type,AMathDimMat > ug::DiscConstraintFVCR< TGridFunction >::aSideDimMat

◆ aSideNumber

template<typename TGridFunction >

typedef PeriodicAttachmentAccessor<side_type,ANumber > ug::DiscConstraintFVCR< TGridFunction >::aSideNumber

◆ aSideNumberArray

template<typename TGridFunction >

typedef PeriodicAttachmentAccessor<side_type,ANumberArray> ug::DiscConstraintFVCR< TGridFunction >::aSideNumberArray

◆ aSideSizetArray

template<typename TGridFunction >

typedef PeriodicAttachmentAccessor<side_type,ASizetArray> ug::DiscConstraintFVCR< TGridFunction >::aSideSizetArray

◆ ASizetArray

template<typename TGridFunction >

typedef Attachment<std::vector< DoFIndex > > ug::DiscConstraintFVCR< TGridFunction >::ASizetArray

◆ dimMat

template<typename TGridFunction >

typedef MathMatrix<dim,dim> ug::DiscConstraintFVCR< TGridFunction >::dimMat

◆ domain_type

template<typename TGridFunction >

typedef TDomain ug::DiscConstraintFVCR< TGridFunction >::domain_type

Type of Domain.

◆ elem_type

template<typename TGridFunction >

typedef TGridFunction::template dim_traits<dim>::grid_base_object ug::DiscConstraintFVCR< TGridFunction >::elem_type

element type

◆ ElemIterator

template<typename TGridFunction >

typedef TGridFunction::template dim_traits<dim>::const_iterator ug::DiscConstraintFVCR< TGridFunction >::ElemIterator

element iterator

◆ face_type0

template<typename TGridFunction >

typedef face_type_traits<dim>::face_type0 ug::DiscConstraintFVCR< TGridFunction >::face_type0

◆ face_type1

template<typename TGridFunction >

typedef face_type_traits<dim>::face_type1 ug::DiscConstraintFVCR< TGridFunction >::face_type1

◆ grid_base_object

template<typename TGridFunction >

typedef domain_traits<TDomain::dim>::grid_base_object ug::DiscConstraintFVCR< TGridFunction >::grid_base_object

Type of geometric base object.

◆ grid_type

template<typename TGridFunction >

typedef domain_type::grid_type ug::DiscConstraintFVCR< TGridFunction >::grid_type

grid type

◆ MathVector_Pair

template<typename TGridFunction >

typedef std::pair<MathVector<dim>, MathVector<dim> > ug::DiscConstraintFVCR< TGridFunction >::MathVector_Pair

◆ matrix_type

template<typename TGridFunction >

typedef algebra_type::matrix_type ug::DiscConstraintFVCR< TGridFunction >::matrix_type

Type of algebra matrix.

◆ position_accessor_type

template<typename TGridFunction >

typedef domain_type::position_accessor_type ug::DiscConstraintFVCR< TGridFunction >::position_accessor_type

position accessor

◆ side_type

template<typename TGridFunction >

typedef elem_type::side ug::DiscConstraintFVCR< TGridFunction >::side_type

side type

◆ SideIterator

template<typename TGridFunction >

typedef TGridFunction::template traits<side_type>::const_iterator ug::DiscConstraintFVCR< TGridFunction >::SideIterator

side iterator

◆ TAlgebra

template<typename TGridFunction >

typedef TGridFunction::algebra_type ug::DiscConstraintFVCR< TGridFunction >::TAlgebra

◆ TDomain

template<typename TGridFunction >

typedef TGridFunction::domain_type ug::DiscConstraintFVCR< TGridFunction >::TDomain

◆ vector_type

template<typename TGridFunction >

typedef algebra_type::vector_type ug::DiscConstraintFVCR< TGridFunction >::vector_type

Type of algebra vector.

◆ vIndexPosPair

template<typename TGridFunction >

typedef std::vector<std::pair<DoFIndex, MathVector<dim> > > ug::DiscConstraintFVCR< TGridFunction >::vIndexPosPair

◆ vvIndexPosPair

template<typename TGridFunction >

typedef std::vector<std::vector<std::pair<DoFIndex, MathVector<dim> > > > ug::DiscConstraintFVCR< TGridFunction >::vvIndexPosPair

Constructor & Destructor Documentation

◆ DiscConstraintFVCR() [1/5]

template<typename TGridFunction >

ug::DiscConstraintFVCR< TGridFunction >::DiscConstraintFVCR ( SmartPtr< TGridFunction > u )

inline

constructor

References ug::DiscConstraintFVCR< TGridFunction >::init().

◆ DiscConstraintFVCR() [2/5]

template<typename TGridFunction >

ug::DiscConstraintFVCR< TGridFunction >::DiscConstraintFVCR	(	SmartPtr< TGridFunction >	u,
		bool	bLinUpConvDefect,
		bool	bLinUpConvJacobian,
		bool	bLinPressureDefect,
		bool	bLinPressureJacobian,
		bool	bAdaptive
	)

inline

References ug::DiscConstraintFVCR< TGridFunction >::init().

◆ DiscConstraintFVCR() [3/5]

template<typename TGridFunction >

ug::DiscConstraintFVCR< TGridFunction >::DiscConstraintFVCR	(	SmartPtr< TGridFunction >	u,
		bool	bLinUpConvDefect,
		bool	bLinUpConvJacobian,
		bool	bLinPressureDefect,
		bool	bLinPressureJacobian,
		bool	bAdaptive,
		bool	bLimiter
	)

inline

References ug::DiscConstraintFVCR< TGridFunction >::init().

◆ DiscConstraintFVCR() [4/5]

template<typename TGridFunction >

ug::DiscConstraintFVCR< TGridFunction >::DiscConstraintFVCR	(	SmartPtr< TGridFunction >	u,
		bool	bLinUpConvDefect,
		bool	bLinUpConvJacobian,
		bool	bLinPressureDefect,
		bool	bLinPressureJacobian,
		bool	bAdaptive,
		const char *	subsets
	)

inline

References ug::DiscConstraintFVCR< TGridFunction >::init(), and ug::DiscConstraintFVCR< TGridFunction >::set_zero_grad_bnd().

◆ DiscConstraintFVCR() [5/5]

template<typename TGridFunction >

ug::DiscConstraintFVCR< TGridFunction >::DiscConstraintFVCR	(	SmartPtr< TGridFunction >	u,
		bool	bLinUpConvDefect,
		bool	bLinUpConvJacobian,
		bool	bLinPressureDefect,
		bool	bLinPressureJacobian,
		bool	bAdaptive,
		bool	bLimiter,
		const char *	subsets
	)

inline

References ug::DiscConstraintFVCR< TGridFunction >::init(), and ug::DiscConstraintFVCR< TGridFunction >::set_zero_grad_bnd().

◆ ~DiscConstraintFVCR()

template<typename TGridFunction >

ug::DiscConstraintFVCR< TGridFunction >::~DiscConstraintFVCR ( )

inline

destructor

Member Function Documentation

◆ add_defect()

template<typename TGridFunction >

virtual void ug::DiscConstraintFVCR< TGridFunction >::add_defect	(	vector_type &	d,
		const vector_type &	u,
		ConstSmartPtr< DoFDistribution >	dd,
		const number	time = `0.0`,
		const number	s_a = `1.0`
	)

inlinevirtual

Referenced by ug::DiscConstraintFVCR< TGridFunction >::adjust_defect().

◆ add_pressure_defect()

template<typename TGridFunction >

virtual void ug::DiscConstraintFVCR< TGridFunction >::add_pressure_defect	(	vector_type &	d,
		const vector_type &	u,
		ConstSmartPtr< DoFDistribution >	dd,
		const number	time = `0.0`,
		const number	s_a = `1.0`
	)

inlinevirtual

Referenced by ug::DiscConstraintFVCR< TGridFunction >::adjust_defect().

◆ adjust_defect()

template<typename TGridFunction >

virtual void ug::DiscConstraintFVCR< TGridFunction >::adjust_defect	(	vector_type &	d,
		const vector_type &	u,
		ConstSmartPtr< DoFDistribution >	dd,
		int	type,
		number	time = `0.0`,
		ConstSmartPtr< VectorTimeSeries< vector_type > >	vSol = `SPNULL`,
		const std::vector< number > *	vScaleMass = `NULL`,
		const std::vector< number > *	vScaleStiff = `NULL`
	)

inlinevirtual

adapts defect to enforce constraints

Implements ug::IConstraint< class >.

References ug::DiscConstraintFVCR< TGridFunction >::add_defect(), ug::DiscConstraintFVCR< TGridFunction >::add_pressure_defect(), ug::IDomainConstraint< TGridFunction::domain_type, TGridFunction::algebra_type >::dd(), ug::DiscConstraintFVCR< TGridFunction >::m_bLinUpConvDefect, and SPNULL.

◆ adjust_jacobian()

template<typename TGridFunction >

virtual void ug::DiscConstraintFVCR< TGridFunction >::adjust_jacobian	(	matrix_type &	J,
		const vector_type &	u,
		ConstSmartPtr< DoFDistribution >	dd,
		int	type,
		number	time = `0.0`,
		ConstSmartPtr< VectorTimeSeries< vector_type > >	vSol = `NULL`,
		const number	s_a0 = `1.0`
	)

inlinevirtual

handle convection

handle pressure

Implements ug::IConstraint< class >.

◆ adjust_linear()

template<typename TGridFunction >

virtual void ug::DiscConstraintFVCR< TGridFunction >::adjust_linear	(	matrix_type &	mat,
		vector_type &	rhs,
		ConstSmartPtr< DoFDistribution >	dd,
		int	type,
		number	time = `0.0`
	)

inlinevirtual

adapts matrix and rhs (linear case) to enforce constraints

Implements ug::IConstraint< class >.

◆ adjust_rhs()

template<typename TGridFunction >

virtual void ug::DiscConstraintFVCR< TGridFunction >::adjust_rhs	(	vector_type &	rhs,
		const vector_type &	u,
		ConstSmartPtr< DoFDistribution >	dd,
		int	type,
		number	time = `0.0`
	)

inlinevirtual

adapts a rhs to enforce constraints

Implements ug::IConstraint< class >.

◆ adjust_solution()

template<typename TGridFunction >

virtual void ug::DiscConstraintFVCR< TGridFunction >::adjust_solution	(	vector_type &	u,
		ConstSmartPtr< DoFDistribution >	dd,
		int	type,
		number	time = `0.0`
	)

inlinevirtual

sets the constraints in a solution vector

Implements ug::IConstraint< class >.

◆ compute_grad_shapes()

template<typename TGridFunction >

void ug::DiscConstraintFVCR< TGridFunction >::compute_grad_shapes ( )

inline

compute gradient shapes for velocity on rotated elements

References ug::DiscConstraintFVCR< TGridFunction >::acGradSh, ug::DiscConstraintFVCR< TGridFunction >::acGradShInd, ug::DiscConstraintFVCR< TGridFunction >::acVol, ug::CollectCornerCoordinates(), ug::DiscConstraintFVCR< TGridFunction >::dim, SmartPtr< TGridFunction >::get(), ug::DimCRFVGeometry< TDim, TWorldDim >::SCV::global_grad(), ug::PeriodicBoundaryManager::is_slave(), ug::DiscConstraintFVCR< TGridFunction >::m_grid, ug::DiscConstraintFVCR< TGridFunction >::m_u, ug::DimCRFVGeometry< TDim, TWorldDim >::num_scv(), num_subsets(), s, ug::DimCRFVGeometry< TDim, TWorldDim >::scv(), ug::SetAttachmentValues(), UG_ASSERT, ug::DimCRFVGeometry< TDim, TWorldDim >::update(), and ug::DimCRFVGeometry< TDim, TWorldDim >::SCV::volume().

Referenced by ug::DiscConstraintFVCR< TGridFunction >::adjust_jacobian(), and ug::DiscConstraintFVCR< TGridFunction >::init().

◆ init()

template<typename TGridFunction >

void ug::DiscConstraintFVCR< TGridFunction >::init	(	SmartPtr< TGridFunction >	u,
		bool	bLinUpConvDefect,
		bool	bLinUpConvJacobian,
		bool	bLinPressureDefect,
		bool	bLinPressureJacobian,
		bool	bAdaptive,
		bool	bLimiter
	)