Dumux::LinearLocalStiffness< TypeTag, m > Class Template Reference

Base class for linear local assemblers. More...

#include <dumux/porousmediumflow/2p/sequential/diffusion/mimetic/localstiffness.hh>

Inheritance diagram for Dumux::LinearLocalStiffness< TypeTag, m >:

Description

template<class TypeTag, int m>
class Dumux::LinearLocalStiffness< TypeTag, m >

Base class for linear local assemblers.

This class serves as a base class for linear local assemblers. It provides space and access to the local stiffness matrix. The actual assembling is done in a derived class via the virtual assemble method.

Template Parameters

TypeTag	The problem TypeTag
m	number of degrees of freedom per node (system size)

Public Types
using	MBlockType = Dune::FieldMatrix<Scalar, m, m>
using	VBlockType = Dune::FieldVector<Scalar, m>
using	BCBlockType = std::array<BoundaryConditions::Flags, m>
using	BoundaryTypes = GetPropType<TypeTag, Properties::BoundaryTypes>

Public Member Functions
	LinearLocalStiffness ()
virtual	~LinearLocalStiffness ()
virtual void	assemble (const Entity &e, int k=1)=0
	Assembles local stiffness matrix including boundary conditions for given element and order.
virtual void	assemble (const Entity &e, const Dune::BlockVector< VBlockType > &localSolution, int k=1)
	Assembles local stiffness matrix including boundary conditions for given element and order.
virtual void	assembleBoundaryCondition (const Entity &e, int k=1)=0
	assemble only boundary conditions for given element and order
void	print (std::ostream &s, int width, int precision)
	Prints contents of local stiffness matrix.
const MBlockType &	mat (int i, int j) const
	Accesses the local stiffness matrix.
const VBlockType &	rhs (int i) const
	Accesses right hand side.
const BoundaryTypes &	bc (int i) const
	Accesses boundary condition for each dof.
void	setcurrentsize (int s)
	Sets the current size of the local stiffness matrix.
int	currentsize ()
	Gets the current size of the local stiffness matrix.

Protected Attributes
Dune::Matrix< MBlockType >	A
std::vector< VBlockType >	b
std::vector< BoundaryTypes >	bctype

Member Typedef Documentation

BCBlockType

template<class TypeTag, int m>

using Dumux::LinearLocalStiffness< TypeTag, m >::BCBlockType = std::array<BoundaryConditions::Flags, m>

BoundaryTypes

template<class TypeTag, int m>

using Dumux::LocalStiffness< TypeTag, m >::BoundaryTypes = GetPropType<TypeTag, Properties::BoundaryTypes>

inherited

MBlockType

template<class TypeTag, int m>

using Dumux::LinearLocalStiffness< TypeTag, m >::MBlockType = Dune::FieldMatrix<Scalar, m, m>

VBlockType

template<class TypeTag, int m>

using Dumux::LinearLocalStiffness< TypeTag, m >::VBlockType = Dune::FieldVector<Scalar, m>

Constructor & Destructor Documentation

LinearLocalStiffness()

template<class TypeTag, int m>

Dumux::LinearLocalStiffness< TypeTag, m >::LinearLocalStiffness ( )

inline

Initialize local stiffness matrix

~LinearLocalStiffness()

template<class TypeTag, int m>

virtual Dumux::LinearLocalStiffness< TypeTag, m >::~LinearLocalStiffness ( )

inlinevirtual

Member Function Documentation

assemble() [1/2]

template<class TypeTag, int m>

virtual void Dumux::LinearLocalStiffness< TypeTag, m >::assemble ( const Entity & e, const Dune::BlockVector< VBlockType > & localSolution, int k = 1 )

inlinevirtual

Assembles local stiffness matrix including boundary conditions for given element and order.

Since this is a base class for linear assemblers, the local solution will be ignored.

Parameters

e	a codim 0 entity reference
localSolution	The current solution on the entity, which is needed by nonlinear assemblers
k	order of Lagrange basis (default is 1)

Implements Dumux::LocalStiffness< TypeTag, m >.

assemble() [2/2]

template<class TypeTag, int m>

virtual void Dumux::LinearLocalStiffness< TypeTag, m >::assemble ( const Entity & e, int k = 1 )

pure virtual

Assembles local stiffness matrix including boundary conditions for given element and order.

On exit the following things have been done:

• The stiffness matrix for the given entity and polynomial degree has been assembled and is accessible with the mat() method.
• The boundary conditions have been evaluated and are accessible with the bc() method. The boundary conditions are either neumann, process or dirichlet. Neumann indicates that the corresponding node (assuming a nodal basis) is at the Neumann boundary, process indicates that the node is at a process boundary (arising from the parallel decomposition of the mesh). Process boundaries are treated as homogeneous Dirichlet conditions, i.e. the corresponding value in the right hand side is set to 0. Finally, Dirichlet indicates that the node is at the Dirichlet boundary.
• The right hand side has been assembled. It contains either the value of the essential boundary condition or the assembled source term and neumann boundary condition. It is accessible via the rhs() method.

Parameters

e	a codim 0 entity reference
k	order of Lagrange basis (default is 1)

Implements Dumux::LocalStiffness< TypeTag, m >.

assembleBoundaryCondition()

template<class TypeTag, int m>

virtual void Dumux::LocalStiffness< TypeTag, m >::assembleBoundaryCondition ( const Entity & e, int k = 1 )

pure virtualinherited

assemble only boundary conditions for given element and order

On exit the following things have been done:

• The boundary conditions have been evaluated and are accessible with the bc() method. The boundary conditions are either Neumann, process or Dirichlet. Neumann indicates that the corresponding node (assuming a nodal basis) is at the Neumann boundary, process indicates that the node is at a process boundary (arising from the parallel decomposition of the mesh). Process boundaries are treated as homogeneous Dirichlet conditions, i.e. the corresponding value in the right hand side is set to 0. Finally, Dirichlet indicates that the node is at the Dirichlet boundary.
• The right hand side has been assembled as far as boundary conditions are concerned. It contains either the value of the essential boundary condition or the assembled Neumann boundary condition. It is accessible via the rhs() method.

Parameters

e	a codim 0 entity reference
k	order of Lagrange basis (default is 1)

bc()

template<class TypeTag, int m>

const BoundaryTypes & Dumux::LocalStiffness< TypeTag, m >::bc ( int i ) const

inlineinherited

Accesses boundary condition for each dof.

Access boundary condition type for each degree of freedom. Elements are undefined without prior call to the assemble method.

Parameters

i

index i

currentsize()

template<class TypeTag, int m>

int Dumux::LocalStiffness< TypeTag, m >::currentsize ( )

inlineinherited

Gets the current size of the local stiffness matrix.

mat()

template<class TypeTag, int m>

const MBlockType & Dumux::LocalStiffness< TypeTag, m >::mat ( int i, int j ) const

inlineinherited

Accesses the local stiffness matrix.

Access elements of the local stiffness matrix. Elements are undefined without prior call to the assemble method.

Parameters

i	index i
j	index j

print()

template<class TypeTag, int m>

void Dumux::LocalStiffness< TypeTag, m >::print ( std::ostream & s, int width, int precision )

inlineinherited

Prints contents of local stiffness matrix.

Parameters

s	output stream
width	the width
precision	the precision

rhs()

template<class TypeTag, int m>

const VBlockType & Dumux::LocalStiffness< TypeTag, m >::rhs ( int i ) const

inlineinherited

Accesses right hand side.

Access elements of the right hand side vector. Elements are undefined without prior call to the assemble method.

Parameters

i

index i

setcurrentsize()

template<class TypeTag, int m>

void Dumux::LocalStiffness< TypeTag, m >::setcurrentsize ( int s )

inlineinherited

Sets the current size of the local stiffness matrix.

Parameters

s	size which is to be set

Member Data Documentation

A

template<class TypeTag, int m>

Dune::Matrix<MBlockType> Dumux::LocalStiffness< TypeTag, m >::A

protectedinherited

b

template<class TypeTag, int m>

std::vector<VBlockType> Dumux::LocalStiffness< TypeTag, m >::b

protectedinherited

bctype

template<class TypeTag, int m>

std::vector<BoundaryTypes> Dumux::LocalStiffness< TypeTag, m >::bctype

protectedinherited

---

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

• localstiffness.hh