#include <deal.II/meshworker/assembler.h>

Public Member Functions
	MGMatrixLocalBlocksToGlobalBlocks (double threshold=1.e-12)

void	initialize (const BlockInfo *block_info, MatrixPtrVector &matrices)

void	initialize (const MGConstrainedDoFs &mg_constrained_dofs)

void	initialize_edge_flux (MatrixPtrVector &up, MatrixPtrVector &down)

void	initialize_interfaces (MatrixPtrVector &interface_in, MatrixPtrVector &interface_out)

template<class DOFINFO >
void	initialize_info (DOFINFO &info, bool face) const

template<class DOFINFO >
void	assemble (const DOFINFO &info)

template<class DOFINFO >
void	assemble (const DOFINFO &info1, const DOFINFO &info2)

Private Member Functions
void	assemble (MatrixType &global, const FullMatrix< number > &local, const unsigned int block_row, const unsigned int block_col, const std::vector< types::global_dof_index > &dof1, const std::vector< types::global_dof_index > &dof2, const unsigned int level1, const unsigned int level2, bool transpose=false)

void	assemble_fluxes (MatrixType &global, const FullMatrix< number > &local, const unsigned int block_row, const unsigned int block_col, const std::vector< types::global_dof_index > &dof1, const std::vector< types::global_dof_index > &dof2, const unsigned int level1, const unsigned int level2)

void	assemble_up (MatrixType &global, const FullMatrix< number > &local, const unsigned int block_row, const unsigned int block_col, const std::vector< types::global_dof_index > &dof1, const std::vector< types::global_dof_index > &dof2, const unsigned int level1, const unsigned int level2)

void	assemble_down (MatrixType &global, const FullMatrix< number > &local, const unsigned int block_row, const unsigned int block_col, const std::vector< types::global_dof_index > &dof1, const std::vector< types::global_dof_index > &dof2, const unsigned int level1, const unsigned int level2)

void	assemble_in (MatrixType &global, const FullMatrix< number > &local, const unsigned int block_row, const unsigned int block_col, const std::vector< types::global_dof_index > &dof1, const std::vector< types::global_dof_index > &dof2, const unsigned int level1, const unsigned int level2)

void	assemble_out (MatrixType &global, const FullMatrix< number > &local, const unsigned int block_row, const unsigned int block_col, const std::vector< types::global_dof_index > &dof1, const std::vector< types::global_dof_index > &dof2, const unsigned int level1, const unsigned int level2)

Private Attributes
MatrixPtrVectorPtr	matrices

MatrixPtrVectorPtr	flux_down

MatrixPtrVectorPtr	flux_up

MatrixPtrVectorPtr	interface_out

MatrixPtrVectorPtr	interface_in

SmartPointer< const BlockInfo, MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number > >	block_info

SmartPointer< const MGConstrainedDoFs, MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number > >	mg_constrained_dofs

const double	threshold

Detailed Description

template<typename MatrixType, typename number = double>
class MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >

A helper class assembling local matrices into global multilevel matrices. This class is the multilevel equivalent of MatrixLocalBlocksToGlobalBlocks and documentation of that class applies here to a large extend.

The global matrices are expected as a vector of pointers to MatrixBlock objects, each containing a MGLevelObject with matrices with a function corresponding to SparseMatrix::add() and information on the block row and column this matrix represents in a block system.

The local matrices are a similar vector of MatrixBlock objects, but containing a FullMatrix.

If local refinement occurs, the Multigrid method needs more matrices, two for continuous elements and another two if numerical fluxes are computed on interfaces. The second set can be added using initialize_edge_flux(). Once added, the contributions in all participating matrices will be assembled from the cell and face matrices automatically.

Author: Guido Kanschat, 2009

Definition at line 302 of file assembler.h.

Constructor & Destructor Documentation

◆ MGMatrixLocalBlocksToGlobalBlocks()

template<typename MatrixType , typename number >

MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::MGMatrixLocalBlocksToGlobalBlocks ( double threshold = 1.e-12 )

inline

Constructor, initializing the threshold, which limits how small numbers may be to be entered into the matrix.

Definition at line 706 of file assembler.h.

Member Function Documentation

◆ initialize() [1/2]

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::initialize	(	const BlockInfo *	block_info,
		MatrixPtrVector &	matrices
	)

inline

Copy the BlockInfo and the matrix pointers into local variables and initialize cell matrix vectors.

Definition at line 715 of file assembler.h.

◆ initialize() [2/2]

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::initialize ( const MGConstrainedDoFs & mg_constrained_dofs )

inline

Initialize the multilevel constraints.

Definition at line 727 of file assembler.h.

◆ initialize_edge_flux()

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::initialize_edge_flux	(	MatrixPtrVector &	up,
		MatrixPtrVector &	down
	)

inline

Multigrid methods on locally refined meshes need additional matrices. For discontinuous Galerkin methods, these are two flux matrices across the refinement edge, which are set by this method.

Definition at line 748 of file assembler.h.

◆ initialize_interfaces()

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::initialize_interfaces	(	MatrixPtrVector &	interface_in,
		MatrixPtrVector &	interface_out
	)

inline

Multigrid methods on locally refined meshes need additional matrices. For discontinuous Galerkin methods, these are two flux matrices across the refinement edge, which are set by this method.

Definition at line 759 of file assembler.h.

◆ initialize_info()

template<typename MatrixType , typename number >

template<class DOFINFO >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::initialize_info	(	DOFINFO &	info,
		bool	face
	)		const

inline

Initialize the local data in the DoFInfo object used later for assembling.

The info object refers to a cell if !face, or else to an interior or boundary face.

Definition at line 737 of file assembler.h.

◆ assemble() [1/3]

template<typename MatrixType , typename number >

template<class DOFINFO >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble ( const DOFINFO & info )

inline

Assemble the local matrices into the global matrices.

Definition at line 1139 of file assembler.h.

◆ assemble() [2/3]

template<typename MatrixType , typename number >

template<class DOFINFO >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble	(	const DOFINFO &	info1,
		const DOFINFO &	info2
	)

inline

Assemble all local matrices into the global matrices.

Definition at line 1174 of file assembler.h.

◆ assemble() [3/3]

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble	(	MatrixType &	global,
		const FullMatrix< number > &	local,
		const unsigned int	block_row,
		const unsigned int	block_col,
		const std::vector< types::global_dof_index > &	dof1,
		const std::vector< types::global_dof_index > &	dof2,
		const unsigned int	level1,
		const unsigned int	level2,
		bool	transpose = `false`
	)

inlineprivate

Assemble a single local matrix into a global one.

Definition at line 770 of file assembler.h.

◆ assemble_fluxes()

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble_fluxes	(	MatrixType &	global,
		const FullMatrix< number > &	local,
		const unsigned int	block_row,
		const unsigned int	block_col,
		const std::vector< types::global_dof_index > &	dof1,
		const std::vector< types::global_dof_index > &	dof2,
		const unsigned int	level1,
		const unsigned int	level2
	)

inlineprivate

Assemble a single local matrix into a global one.

Definition at line 849 of file assembler.h.

◆ assemble_up()

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble_up	(	MatrixType &	global,
		const FullMatrix< number > &	local,
		const unsigned int	block_row,
		const unsigned int	block_col,
		const std::vector< types::global_dof_index > &	dof1,
		const std::vector< types::global_dof_index > &	dof2,
		const unsigned int	level1,
		const unsigned int	level2
	)

inlineprivate

Assemble a single local matrix into a global one.

Definition at line 902 of file assembler.h.

◆ assemble_down()

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble_down	(	MatrixType &	global,
		const FullMatrix< number > &	local,
		const unsigned int	block_row,
		const unsigned int	block_col,
		const std::vector< types::global_dof_index > &	dof1,
		const std::vector< types::global_dof_index > &	dof2,
		const unsigned int	level1,
		const unsigned int	level2
	)

inlineprivate

Assemble a single local matrix into a global one.

Definition at line 955 of file assembler.h.

◆ assemble_in()

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble_in	(	MatrixType &	global,
		const FullMatrix< number > &	local,
		const unsigned int	block_row,
		const unsigned int	block_col,
		const std::vector< types::global_dof_index > &	dof1,
		const std::vector< types::global_dof_index > &	dof2,
		const unsigned int	level1,
		const unsigned int	level2
	)

inlineprivate

Assemble a single local matrix into a global one.

Definition at line 1008 of file assembler.h.

◆ assemble_out()

template<typename MatrixType , typename number >

void MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::assemble_out	(	MatrixType &	global,
		const FullMatrix< number > &	local,
		const unsigned int	block_row,
		const unsigned int	block_col,
		const std::vector< types::global_dof_index > &	dof1,
		const std::vector< types::global_dof_index > &	dof2,
		const unsigned int	level1,
		const unsigned int	level2
	)

inlineprivate

Assemble a single local matrix into a global one.

Definition at line 1073 of file assembler.h.

Member Data Documentation

◆ matrices

template<typename MatrixType , typename number = double>

MatrixPtrVectorPtr MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::matrices

private

The level matrices, stored as a vector of pointers.

Definition at line 441 of file assembler.h.

◆ flux_down

template<typename MatrixType , typename number = double>

MatrixPtrVectorPtr MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::flux_down

private

The flux matrix between the fine and the coarse level at refinement edges.

Definition at line 447 of file assembler.h.

◆ flux_up

template<typename MatrixType , typename number = double>

MatrixPtrVectorPtr MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::flux_up

private

The flux matrix between the coarse and the fine level at refinement edges.

Definition at line 453 of file assembler.h.

◆ interface_out

template<typename MatrixType , typename number = double>

MatrixPtrVectorPtr MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::interface_out

private

The interface matrix between the fine and the coarse level at refinement edges.

Definition at line 459 of file assembler.h.

◆ interface_in

template<typename MatrixType , typename number = double>

MatrixPtrVectorPtr MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::interface_in

private

The interface matrix between the coarse and the fine level at refinement edges.

Definition at line 465 of file assembler.h.

◆ block_info

template<typename MatrixType , typename number = double>

SmartPointer<const BlockInfo, MGMatrixLocalBlocksToGlobalBlocks<MatrixType, number> > MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::block_info

private

A pointer to the object containing the block structure.

Definition at line 470 of file assembler.h.

◆ mg_constrained_dofs

template<typename MatrixType , typename number = double>

SmartPointer<const MGConstrainedDoFs,MGMatrixLocalBlocksToGlobalBlocks<MatrixType, number> > MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::mg_constrained_dofs

private

A pointer to the object containing constraints.

Definition at line 475 of file assembler.h.

◆ threshold

template<typename MatrixType , typename number = double>

const double MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >::threshold

private

The smallest positive number that will be entered into the global matrix. All smaller absolute values will be treated as zero and will not be assembled.

Definition at line 483 of file assembler.h.

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

deal.II/meshworker/assembler.h

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

template<typename MatrixType, typename number = double> class MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >

Constructor & Destructor Documentation

◆ MGMatrixLocalBlocksToGlobalBlocks()

Member Function Documentation

◆ initialize() [1/2]

◆ initialize() [2/2]

◆ initialize_edge_flux()

◆ initialize_interfaces()

◆ initialize_info()

◆ assemble() [1/3]

◆ assemble() [2/3]

◆ assemble() [3/3]

◆ assemble_fluxes()

◆ assemble_up()

◆ assemble_down()

◆ assemble_in()

◆ assemble_out()

Member Data Documentation

◆ matrices

◆ flux_down

◆ flux_up

◆ interface_out

◆ interface_in

◆ block_info

◆ mg_constrained_dofs

◆ threshold

template<typename MatrixType, typename number = double>
class MeshWorker::Assembler::MGMatrixLocalBlocksToGlobalBlocks< MatrixType, number >