Skip to content
Snippets Groups Projects
Commit 608f0ce1 authored by Jonathan Youett's avatar Jonathan Youett
Browse files

Improve documentation

parent 91217e95
No related branches found
No related tags found
No related merge requests found
...@@ -22,10 +22,6 @@ restrict(const std::vector<BoxConstraint<field_type, blocksize> >&f, ...@@ -22,10 +22,6 @@ restrict(const std::vector<BoxConstraint<field_type, blocksize> >&f,
for (size_t i=0; i<t.size(); i++) for (size_t i=0; i<t.size(); i++)
t[i].clear(); t[i].clear();
// The Mandel restriction computes the maximum $m$ of all nodal base coefficients
// in the support of a given coarse grid base function \f$ \phi_c \f$. If the
// coefficient of \f$ \phi_c \f$ is smaller than \f$ m \f$ , it is set to \f$ m \f$.
// We use the multigrid transfer matrix to find all relevant fine grid dofs // We use the multigrid transfer matrix to find all relevant fine grid dofs
typedef typename Dune::BCRSMatrix<Dune::FieldMatrix<field_type,1,1> >::row_type RowType; typedef typename Dune::BCRSMatrix<Dune::FieldMatrix<field_type,1,1> >::row_type RowType;
typedef typename RowType::ConstIterator ColumnIterator; typedef typename RowType::ConstIterator ColumnIterator;
......
...@@ -8,6 +8,14 @@ ...@@ -8,6 +8,14 @@
#include <dune/solvers/transferoperators/multigridtransfer.hh> #include <dune/solvers/transferoperators/multigridtransfer.hh>
#include <dune/solvers/common/boxconstraint.hh> #include <dune/solvers/common/boxconstraint.hh>
/** \brief Mandel restriction to compute coarse grid obstacles, needed within montone multigrid methods.
*
* The Mandel restriction computes the maximum $m$ of all nodal base coefficients
* in the support of a given coarse grid base function \f$ \phi_c \f$. If the
* coefficient of \f$ \phi_c \f$ is smaller than \f$ m \f$ , it is set to \f$ m \f$.
*
* \tparam VectorType The block vector type of the algebraic problem.
*/
template <class VectorType> template <class VectorType>
class MandelObstacleRestrictor : public ObstacleRestrictor<VectorType> class MandelObstacleRestrictor : public ObstacleRestrictor<VectorType>
{ {
...@@ -17,6 +25,7 @@ class MandelObstacleRestrictor : public ObstacleRestrictor<VectorType> ...@@ -17,6 +25,7 @@ class MandelObstacleRestrictor : public ObstacleRestrictor<VectorType>
public: public:
//! Compute coarse grid obstacles using the mandel restriction
virtual void restrict(const std::vector<BoxConstraint<field_type,blocksize> >& f, virtual void restrict(const std::vector<BoxConstraint<field_type,blocksize> >& f,
std::vector<BoxConstraint<field_type,blocksize> >& t, std::vector<BoxConstraint<field_type,blocksize> >& t,
const Dune::BitSetVector<blocksize>& fHasObstacle, const Dune::BitSetVector<blocksize>& fHasObstacle,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment