We can now completely control the approximation order by a const int,
which is much nicer.

This data is used by the compute-disc-error program to measure discretization errors.

Because these two classes are for the most part identical.

Currently, this merged version involves some trickery, including some template
specializations and multiple inheritances.  I hope to get rid of that, eventually.

This is an important step towards making the 'mixed' Cosserat code more like
the regular one.  The mid-term goal is to merge the two implementations, because
there is way too much code duplication here.

... with values expected to be given in the parameter file.

Scaling the trust region is the key to success, because in general the midsurface
deformation and the rotations will have values in different ranges.

Unfortunately, this bases needs to be switched on and off using preprocessor
commands.  Hopefully I'll find a way to make that prettier in the future.

It was only in the dune-gfe module for historic reasons.

The inline expression in the parameter file that we used previously is
too short and limited.  In particular, we want to be able to use the
same expression for error measurements.  However, for that we need
a function with a derivative.

[[Imported from SVN: r10124]]

[[Imported from SVN: r10123]]

[[Imported from SVN: r10122]]

[[Imported from SVN: r10120]]

[[Imported from SVN: r10119]]

[[Imported from SVN: r10113]]

That really was a hack: the build system should do it.

[[Imported from SVN: r10111]]

[[Imported from SVN: r10110]]

[[Imported from SVN: r10109]]

[[Imported from SVN: r10108]]

[[Imported from SVN: r10107]]

[[Imported from SVN: r10106]]

[[Imported from SVN: r10104]]

[[Imported from SVN: r10103]]

[[Imported from SVN: r10102]]

[[Imported from SVN: r10098]]

[[Imported from SVN: r10096]]

[[Imported from SVN: r10095]]

[[Imported from SVN: r10094]]

[[Imported from SVN: r10093]]