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>
30 namespace Elasticity {
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)
103 void addElement(
const Dune::FieldMatrix<double,esize,esize>* K,
104 const Dune::FieldVector<double,esize>* S,
134 const std::pair<Direction,NodeValue>& entry);
159 for (MPCMap::iterator it =
mpcs.begin();
160 it !=
mpcs.end();++it)
190 void addDOF(
int row,
int erow,
191 const Dune::FieldMatrix<double,esize,esize>* K,
192 const Dune::FieldVector<double,esize>* S,
211 typedef typename fixMap::iterator
fixIt;
Matrix A
The linear operator.
Definition: asmhandler.hpp:220
int getEquationForDof(int node, int dof)
Get the equation number for a given dof on a given node.
Definition: asmhandler.hpp:73
Helper class with some matrix operations.
MPC * getMPC(int node, int dof)
Look for and return a MPC for a specified node+dof.
Definition: asmhandler_impl.hpp:207
MPCMap mpcs
The set of MPC.
Definition: asmhandler.hpp:199
Representation of multi-point constraint (MPC) equations.
void addMPC(MPC *mpc)
Add a MPC.
Definition: asmhandler_impl.hpp:189
GridType::LeafGridView::template Codim< 0 >::Iterator LeafIterator
An iterator over grid cells.
Definition: asmhandler.hpp:44
~ASMHandler()
Destructor.
Definition: asmhandler.hpp:53
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
GridType::LeafGridView::IndexSet LeafIndexSet
A set of indices.
Definition: asmhandler.hpp:41
AdjacencyPattern adjacencyPattern
Holds the adjacency pattern of the sparse matrix.
Definition: asmhandler.hpp:217
void expandSolution(Vector &result, const Vector &u)
Expand a system vector to a solution vector.
Definition: asmhandler_impl.hpp:147
size_t getEqns() const
Get the number of equations in the system.
Definition: asmhandler.hpp:64
ASMHandler(const GridType &gv_)
The default constructor.
Definition: asmhandler.hpp:48
void determineAdjacencyPattern()
Internal function. Calculate adjacency pattern.
Definition: asmhandler_impl.hpp:377
bool isFixed(int node)
Check if a node is marked as fixed (in any direction)
Definition: asmhandler.hpp:138
void initForAssembly()
This function needs to be called before starting the element assembly.
Definition: asmhandler_impl.hpp:22
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
Vector & getLoadVector()
Obtain a reference to the load vector.
Definition: asmhandler.hpp:87
Dune::FieldVector< double, dim > NodeValue
A vectorial node value.
Definition: asmhandler.hpp:60
Matrix & getOperator()
Obtain a reference to the linear operator.
Definition: asmhandler.hpp:80
Class handling finite element assembly.
Definition: asmhandler.hpp:35
void printOperator() const
Print the current operator.
Definition: asmhandler_impl.hpp:237
void extractValues(Dune::FieldVector< double, comp > &v, const Vector &u, const LeafIterator &it)
Extract values corresponding to cell.
Definition: asmhandler_impl.hpp:114
fixMap::iterator fixIt
Iterator over a fixmap.
Definition: asmhandler.hpp:211
std::pair< Direction, NodeValue > fixEntry
Fixed nodes.
Definition: asmhandler.hpp:205
void preprocess()
Internal function. Generate meqn for registered MPC/fixed nodes.
Definition: asmhandler_impl.hpp:317
Class handling finite element assembly - template implementations.
size_t maxeqn
The number of equations in the system.
Definition: asmhandler.hpp:229
Vector b
The load vector.
Definition: asmhandler.hpp:223
void updateFixedNode(int node, const std::pair< Direction, NodeValue > &entry)
Update/add a fixed node.
Definition: asmhandler_impl.hpp:216
const GridType & gv
A reference to the grid in use.
Definition: asmhandler.hpp:226
void resolveMPCChain(MPC *mpc)
Internal function.
Definition: asmhandler_impl.hpp:252
static const int dim
The dimension of the grid.
Definition: asmhandler.hpp:38
std::vector< int > meqn
Vector of (interleaved) dof<->eqn mapping.
Definition: asmhandler.hpp:202
void printLoadVector() const
Print the current load vector.
Definition: asmhandler_impl.hpp:243
A class for representing a general multi-point constraint equation.
Definition: mpc.hh:58
Logging helper utilities.
void resolveMPCChains()
Resolve chained MPCs.
Definition: asmhandler.hpp:157
std::map< int, fixEntry > fixMap
A mapping from dof to fix value info.
Definition: asmhandler.hpp:208
void nodeAdjacency(const LeafIterator &it, int vertexsize, int row)
Internal function.
Definition: asmhandler_impl.hpp:353
AdjacencyPattern & getAdjacencyPattern()
Access current adjacency pattern.
Definition: asmhandler.hpp:151
fixMap fixedNodes
The map holding information about our fixed nodes.
Definition: asmhandler.hpp:214