Reference documentation for deal.II version 9.6.0
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
Loading...
Searching...
No Matches
MeshWorker::Assembler Namespace Reference

Classes

class  CellsAndFaces
 
class  Functional
 
class  GnuplotPatch
 
class  MatrixLocalBlocksToGlobalBlocks
 
class  MatrixSimple
 
class  MGMatrixLocalBlocksToGlobalBlocks
 
class  MGMatrixSimple
 
class  ResidualLocalBlocksToGlobalBlocks
 
class  ResidualSimple
 
class  SystemSimple
 

Detailed Description

The namespace containing objects that can be used to assemble data computed on cells and faces into global objects. This can reach from collecting the total error estimate from cell and face contributions to assembling matrices and multilevel matrices.

Data models

The class chosen from this namespace determines which data model is used. For the local as well as the global objects, we have the choice between two models:

The comprehensive data model

This is the structure set up by the FESystem class. Globally, this means, data is assembled into one residual vector and into one matrix. These objects may be block vectors and block matrices, but the process of assembling ignores this fact.

Similarly, there is only a single cell vector and cell matrix, respectively, which is indexed by all degrees of freedom of the FESystem. When building the cell matrix, it is necessary to distinguish between the different components of the system and select the right operator for each pair.

The blocked data model

Here, all the blocks are treated separately (in spite of using FESystem for its convenience in other places). For instance, no block matrix is assembled, but a list of blocks, which can be combined later by BlockMatrixArray. Locally, this means, that each matrix block of a system is generated separately and assembled into the corresponding global block.

This approach is advantageous, if the number of matrices for each block position in the global system is different. For instance, block preconditioners for the Oseen problem require 3 pressure matrices, but only one divergence and one advection-diffusion operator for velocities.

Additionally, this approach enables the construction of a system of equations from building blocks for each equation and coupling operator.

Nevertheless, since a separate FEValues object must be created for each base element, it is not quite clear a priori, which data model is more efficient.