aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number > Class Template Reference

Inherits Base< dim, ::LinearAlgebra::distributed::Vector< number > >.

Public Member Functions
	ABlockOperator ()
void	clear () override
void	reinit (const Mapping< dim > &mapping, const DoFHandler< dim > &dof_handler_v, const DoFHandler< dim > &dof_handler_p, const AffineConstraints< number > &constraints_v, const AffineConstraints< number > &constraints_p, std::shared_ptr< MatrixFree< dim, double >> mf_storage, const unsigned int level=numbers::invalid_unsigned_int)
void	set_cell_data (const OperatorCellData< dim, number > &data)
void	compute_diagonal () override
void	set_diagonal (const ::LinearAlgebra::distributed::Vector< number > &diag)

Private Member Functions
void	inner_cell_operation (FEEvaluation< dim, degree_v, degree_v+1, dim, number > &velocity) const
void	cell_operation (FEEvaluation< dim, degree_v, degree_v+1, dim, number > &velocity) const
void	apply_add (::LinearAlgebra::distributed::Vector< number > &dst, const ::LinearAlgebra::distributed::Vector< number > &src) const override
void	local_apply (const ::MatrixFree< dim, number > &data, ::LinearAlgebra::distributed::Vector< number > &dst, const ::LinearAlgebra::distributed::Vector< number > &src, const std::pair< unsigned int, unsigned int > &cell_range) const

Private Attributes
const OperatorCellData< dim, number > *	cell_data

Detailed Description

template<int dim, int degree_v, typename number>
class aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >

Operator for the A block of the Stokes matrix. The same class is used for both active and level mesh operators.

Definition at line 421 of file matrix_free_operators.h.

Constructor & Destructor Documentation

ABlockOperator()

template<int dim, int degree_v, typename number>

aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::ABlockOperator

(

)

Constructor

Member Function Documentation

clear()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::clear ( )

override

Reset the operator.

reinit()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::reinit	(	const Mapping< dim > &	mapping,
		const DoFHandler< dim > &	dof_handler_v,
		const DoFHandler< dim > &	dof_handler_p,
		const AffineConstraints< number > &	constraints_v,
		const AffineConstraints< number > &	constraints_p,
		std::shared_ptr< MatrixFree< dim, double >>	mf_storage,
		const unsigned int	level = numbers::invalid_unsigned_int
	)

Initialize the MatrixFree object given in mf_storage and use that to initialize this operator.

set_cell_data()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::set_cell_data

(

const OperatorCellData< dim, number > &

data

)

Pass in a reference to the problem data.

compute_diagonal()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::compute_diagonal ( )

override

Computes the diagonal of the matrix. Since matrix-free operators have not access to matrix elements, we must apply the matrix-free operator to the unit vectors to recover the diagonal.

set_diagonal()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::set_diagonal

(

const ::LinearAlgebra::distributed::Vector< number > &

diag

)

Manually set the diagonal inside the matrix-free object. This function is needed when using tangential constraints as the function compute_diagonal() cannot handle non-Dirichlet boundary conditions.

inner_cell_operation()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::inner_cell_operation ( FEEvaluation< dim, degree_v, degree_v+1, dim, number > &  velocity ) const

private

Defines the inner-most operator on a single cell batch with the loop over quadrature points.

cell_operation()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::cell_operation ( FEEvaluation< dim, degree_v, degree_v+1, dim, number > &  velocity ) const

private

Defines the operation on a single cell batch including load/store and calls inner_cell_operation().

apply_add()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::apply_add ( ::LinearAlgebra::distributed::Vector< number > &  dst, const ::LinearAlgebra::distributed::Vector< number > &  src  ) const

overrideprivate

Performs the application of the matrix-free operator. This function is called by vmult() functions MatrixFreeOperators::Base.

local_apply()

template<int dim, int degree_v, typename number>

void aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::local_apply ( const ::MatrixFree< dim, number > &  data, ::LinearAlgebra::distributed::Vector< number > &  dst, const ::LinearAlgebra::distributed::Vector< number > &  src, const std::pair< unsigned int, unsigned int > &  cell_range  ) const

private

Defines the application of the cell matrix.

Member Data Documentation

cell_data

template<int dim, int degree_v, typename number>

const OperatorCellData<dim,number>* aspect::MatrixFreeStokesOperators::ABlockOperator< dim, degree_v, number >::cell_data

private

A pointer to the current cell data that contains viscosity and other required parameters per cell.

Definition at line 506 of file matrix_free_operators.h.

---

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

• /home/bob/source/include/aspect/simulator/solver/matrix_free_operators.h