Extends CROperatorBase by a generic methods to assemble global stiffness matrix from local stiffness matrices.
More...
template<class TypeTag>
class Dumux::CROperatorAssemblerTwoP< TypeTag >
Extends CROperatorBase by a generic methods to assemble global stiffness matrix from local stiffness matrices.
A class for mapping a CR function to a CR function This class sets up a compressed row storage matrix with connectivity for CR elements. This class does not fill any entries into the matrix.
The template parameter TypeTag describes what kind of Assembler we are. There two choices:
LevelTag We assemble on a grid level.
LeafTag We assemble on the leaf entities of the grid
- Template Parameters
-
TypeTag | The problem Type Tag |
Assembles global stiffness matrix.
This method takes an object that can compute local stiffness matrices and assembles the global linear system Au=f.
- Parameters
-
loc | the local assembler providing element stiffness and boundary conditions for all elements |
u | solution, contains initial values on input, Dirichlet values are set. The type of boundary condition for a node is inferred from the values returned by the local assembler. A node is of Neumann type if all elements referring to that node report a Neumann boundary condition, it is set to Dirichlet if a least one element reports a process or Dirichlet boundary condition. The difference between process and Dirichlet is that process always denotes a homogeneous Dirichlet value. |
f | right hand side is filled by this method |
Note that the rows corresponding to nodes at the Dirichlet boundary are filled with trivial equations of the form
\[1\cdot u_i = f_i \]
where \(u_i\) and \(f_i\) are both set to the Dirichlet value at the \(i\)th node.