We may want to override options from the parameter file.

[[Imported from SVN: r9946]]

[[Imported from SVN: r9945]]

[[Imported from SVN: r9944]]

[[Imported from SVN: r9943]]

To really make sure I am doing the infamous L-shape example
like the people before me did I need point loads.  Being a
mathematician I am no particular fan of point loads, so they
may go out again eventually.

[[Imported from SVN: r9942]]

[[Imported from SVN: r9941]]

This is the interpolation method used by Ingo Münch and Wolfgang Müller for
Cosserat materials.  My implementation is only able to do all that I needed
for the illustrations in the GAMM Rundbrief article.  I'm sure there'll be
more fixes needed if you actually want to do finite elements with this.

[[Imported from SVN: r9940]]

For example, I'll use this for finite-strain elasticity problems.

Some of this should be in dune-fufem and/or dune-pdelab, but for the time
being it's easier to have it here.

[[Imported from SVN: r9939]]

[[Imported from SVN: r9938]]

This form allows to set the shear correction factor kappa.  It is still
not ruled out that we need to twiddle with this factor in order to
weaken the material a little bit.

[[Imported from SVN: r9937]]

[[Imported from SVN: r9936]]

[[Imported from SVN: r9935]]

[[Imported from SVN: r9934]]

[[Imported from SVN: r9932]]

[[Imported from SVN: r9931]]

[[Imported from SVN: r9930]]

Gauss-Seidel is simply to slow, the Newton convergence suffers too much.

[[Imported from SVN: r9929]]

[[Imported from SVN: r9928]]

[[Imported from SVN: r9927]]

That code was buggy, and we only didn't see that because the wrong sizes
were always overwritten by the correct ones before being used.
(That's so as long as you call a reducing method before 'scatter').

Removing the code leaves the situation as bad as before, but it
allows us to get rid of the nOwnedLocalEntity() method from the global
mapper.

[[Imported from SVN: r9926]]

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 global indexing infrastructure does not currently support P3 spaces.

[[Imported from SVN: r9924]]

[[Imported from SVN: r9923]]

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: r9920]]

That's what this method is called in the dune-grid IndexSet
interface.

[[Imported from SVN: r9919]]

[[Imported from SVN: r9918]]

[[Imported from SVN: r9917]]

Yes it is to be moved to dune-grid, but for the transition it is
easier to use the one from this module.

[[Imported from SVN: r9916]]

[[Imported from SVN: r9915]]

[[Imported from SVN: r9914]]

[[Imported from SVN: r9913]]

[[Imported from SVN: r9912]]

That code was written back in the days when we were only using
first order spaces.  Now, with other spaces, we have to properly
distinguish between grid vertices and Lagrange nodes.

[[Imported from SVN: r9911]]

[[Imported from SVN: r9910]]

[[Imported from SVN: r9909]]

Use MPIHelper instead.

[[Imported from SVN: r9908]]

[[Imported from SVN: r9907]]

[[Imported from SVN: r9906]]