Because it is not an index, it really is a mapper.

[[Imported from SVN: r9953]]

[[Imported from SVN: r9936]]

In particular this means having to carry mappers for the local indices around.
This makes the code more verbose, but it also makes it easier to see what is
happening.  While coding this stuff I keep having doubts about whether all
the renumbering works correctly...

[[Imported from SVN: r9925]]

The local copy was to get the thing developed and working.  Now that it is
working we upstream it into dune-grid.

[[Imported from SVN: r9922]]

The IndexSet interface class from dune-grid doesn't have it, and we want
to be like the interface method as much as possible.

[[Imported from SVN: r9921]]

[[Imported from SVN: r9918]]

[[Imported from SVN: r9917]]

[[Imported from SVN: r9909]]

[[Imported from SVN: r9890]]

dune-solvers actually does provide a way to set these matrices directly.
I simply didn't find it until recently.

[[Imported from SVN: r9885]]

[[Imported from SVN: r9881]]

[[Imported from SVN: r9878]]

[[Imported from SVN: r9868]]

[[Imported from SVN: r9854]]

Gauss-Seidel converges equally good/bad as IPOpt, but it is much faster.

[[Imported from SVN: r9839]]

[[Imported from SVN: r9838]]

Previously the energy norm was used, but as we are solving non-convex
problems the Hessians may be indefinite, and then the EnergyNorm is not
a norm.

[[Imported from SVN: r9835]]

[[Imported from SVN: r9817]]

The hasObstacle_ and obstacles_ fields are now single arrays.
The corresponding fields for the entire grid hierarchy is assembled
within dune-solvers now.

[[Imported from SVN: r9816]]

[[Imported from SVN: r9771]]

Previously, the matrix was always assembled using a P1-basis, even if we used
it to measure errors of vectors pertaining to P2 or other spaces.  Nobody noticed
this, because the h1SemiNorm routine would not fail if the vector was larger
than the matrix.   Sigh...

[[Imported from SVN: r9770]]

That makes for slightly simpler code.

[[Imported from SVN: r9768]]

[[Imported from SVN: r9765]]

Even if there is only one processor.  Because even then the local
and global vectors differ in their ordering.

[[Imported from SVN: r9757]]

The Dirichlet nodes are handled directly via the ignoreNodes_ field.
I don't know how the matrix modification code ever got in here.

Removing the code doesn't seem to make any difference.  I find that
a bit surprising: I would expect multigrid and IPOpt to work better
if the matrix is truly symmetric.

[[Imported from SVN: r9742]]

Much easier to use than the previous API.

[[Imported from SVN: r9735]]

[[Imported from SVN: r9734]]

[[Imported from SVN: r9732]]

[[Imported from SVN: r9731]]

[[Imported from SVN: r9729]]

After much back and forth: this is the correct way.  This size cannot
be taken from the transfer operators, because these operators exist
only if there is more than one level.

[[Imported from SVN: r9728]]

[[Imported from SVN: r9727]]

[[Imported from SVN: r9726]]

That works for fe spaces of all orders.

[[Imported from SVN: r9725]]

[[Imported from SVN: r9724]]

[[Imported from SVN: r9723]]

Christian Engwer told me a trick how to get the mpiHelper instance without
known argc and argv.

[[Imported from SVN: r9720]]

[[Imported from SVN: r9710]]

[[Imported from SVN: r9614]]

We have the correct value in the variable 'oldEnergy' anyway.

[[Imported from SVN: r9575]]