12#ifndef ASMHANDLER_HPP_
13#define ASMHANDLER_HPP_
15#include <opm/common/utility/platform_dependent/disable_warnings.h>
17#include <dune/geometry/referenceelements.hh>
18#include <dune/common/fmatrix.hh>
19#include <dune/istl/bvector.hh>
20#include <dune/common/fvector.hh>
22#include <opm/common/utility/platform_dependent/reenable_warnings.h>
34 template<
class Gr
idType>
38 static const int dim = GridType::dimension;
44 typedef typename GridType::LeafGridView::template Codim<0>::Iterator
LeafIterator;
55 for (MPCMap::iterator it=
mpcs.begin(); it !=
mpcs.end();++it)
109 void addElement(
const Dune::FieldMatrix<double,esize,esize>* K,
110 const Dune::FieldVector<double,esize>* S,
140 const std::pair<Direction,NodeValue>& entry);
165 for (MPCMap::iterator it =
mpcs.begin();
166 it !=
mpcs.end();++it)
196 void addDOF(
int row,
int erow,
197 const Dune::FieldMatrix<double,esize,esize>* K,
198 const Dune::FieldVector<double,esize>* S,
217 typedef typename fixMap::iterator
fixIt;
Class handling finite element assembly - template implementations.
Class handling finite element assembly.
Definition asmhandler.hpp:35
static const int dim
The dimension of the grid.
Definition asmhandler.hpp:38
std::pair< Direction, NodeValue > fixEntry
Fixed nodes.
Definition asmhandler.hpp:211
bool isFixed(int node)
Check if a node is marked as fixed (in any direction)
Definition asmhandler.hpp:144
void initForAssembly()
This function needs to be called before starting the element assembly.
Definition asmhandler_impl.hpp:22
MPCMap mpcs
The set of MPC.
Definition asmhandler.hpp:205
void resolveMPCChain(MPC *mpc)
Internal function.
Definition asmhandler_impl.hpp:250
~ASMHandler()
Destructor.
Definition asmhandler.hpp:53
const GridType & gv
A reference to the grid in use.
Definition asmhandler.hpp:232
fixMap fixedNodes
The map holding information about our fixed nodes.
Definition asmhandler.hpp:220
ASMHandler & operator=(const ASMHandler &)=delete
No copying of this class.
Dune::FieldVector< double, dim > NodeValue
A vectorial node value.
Definition asmhandler.hpp:66
size_t maxeqn
The number of equations in the system.
Definition asmhandler.hpp:235
AdjacencyPattern adjacencyPattern
Holds the adjacency pattern of the sparse matrix.
Definition asmhandler.hpp:223
std::vector< int > meqn
Vector of (interleaved) dof<->eqn mapping.
Definition asmhandler.hpp:208
fixMap::iterator fixIt
Iterator over a fixmap.
Definition asmhandler.hpp:217
Vector & getLoadVector()
Obtain a reference to the load vector.
Definition asmhandler.hpp:93
GridType::LeafGridView::IndexSet LeafIndexSet
A set of indices.
Definition asmhandler.hpp:41
int getEquationForDof(int node, int dof)
Get the equation number for a given dof on a given node.
Definition asmhandler.hpp:79
ASMHandler(const ASMHandler &)=delete
No copying of this class.
void addElement(const Dune::FieldMatrix< double, esize, esize > *K, const Dune::FieldVector< double, esize > *S, const LeafIterator &cell, Vector *b=NULL)
Add an element matrix/vector to the system.
Definition asmhandler_impl.hpp:85
void resolveMPCChains()
Resolve chained MPCs.
Definition asmhandler.hpp:163
Vector b
The load vector.
Definition asmhandler.hpp:229
ASMHandler(const GridType &gv_)
The default constructor.
Definition asmhandler.hpp:48
std::map< int, fixEntry > fixMap
A mapping from dof to fix value info.
Definition asmhandler.hpp:214
Matrix A
The linear operator.
Definition asmhandler.hpp:226
void updateFixedNode(int node, const std::pair< Direction, NodeValue > &entry)
Update/add a fixed node.
Definition asmhandler_impl.hpp:214
void extractValues(Dune::FieldVector< double, comp > &v, const Vector &u, const LeafIterator &it)
Extract values corresponding to cell.
Definition asmhandler_impl.hpp:114
void expandSolution(Vector &result, const Vector &u)
Expand a system vector to a solution vector.
Definition asmhandler_impl.hpp:145
void addMPC(MPC *mpc)
Add a MPC.
Definition asmhandler_impl.hpp:187
AdjacencyPattern & getAdjacencyPattern()
Access current adjacency pattern.
Definition asmhandler.hpp:157
void preprocess()
Internal function. Generate meqn for registered MPC/fixed nodes.
Definition asmhandler_impl.hpp:315
void nodeAdjacency(const LeafIterator &it, int vertexsize, int row)
Internal function.
Definition asmhandler_impl.hpp:351
GridType::LeafGridView::template Codim< 0 >::Iterator LeafIterator
An iterator over grid cells.
Definition asmhandler.hpp:44
void printLoadVector() const
Print the current load vector.
Definition asmhandler_impl.hpp:241
Matrix & getOperator()
Obtain a reference to the linear operator.
Definition asmhandler.hpp:86
MPC * getMPC(int node, int dof)
Look for and return a MPC for a specified node+dof.
Definition asmhandler_impl.hpp:205
size_t getEqns() const
Get the number of equations in the system.
Definition asmhandler.hpp:70
void determineAdjacencyPattern()
Internal function. Calculate adjacency pattern.
Definition asmhandler_impl.hpp:375
void addDOF(int row, int erow, const Dune::FieldMatrix< double, esize, esize > *K, const Dune::FieldVector< double, esize > *S, const LeafIndexSet &set, const LeafIterator &cell, Vector *b, double scale=1.f)
Internal function.
Definition asmhandler_impl.hpp:52
void printOperator() const
Print the current operator.
Definition asmhandler_impl.hpp:235
A class for representing a general multi-point constraint equation.
Definition mpc.hh:59
Logging helper utilities.
Helper class with some matrix operations.
Dune::BCRSMatrix< Dune::FieldMatrix< double, 1, 1 > > Matrix
A sparse matrix holding our operator.
Definition matrixops.hpp:27
Dune::BlockVector< Dune::FieldVector< double, 1 > > Vector
A vector holding our RHS.
Definition matrixops.hpp:33
std::vector< std::set< int > > AdjacencyPattern
For storing matrix adjacency/sparsity patterns.
Definition matrixops.hpp:30
Representation of multi-point constraint (MPC) equations.
std::map< int, MPC * > MPCMap
A mapping from dof to MPCs.
Definition mpc.hh:156
Inverting small matrices.
Definition ImplicitAssembly.hpp:43