WHITESAPCES for previous commit: Add the option to use the Riemannian Proximal Newton solver also for problems with dim = dimwold

Use cholmod if the dune-solvers version >= 2.8, use umfpack for older versions

Now all combinations of dim = dimworld or dim != dimworld and MIXED_SPACE = 0 or MIXED_SPACE = 1 at least compile.

This has to be dimworld and not dim.

The StressFreeStateGridFunction is the actual parametrization from the grid to the
nonplanar Cosserat shell in stress-free state.
The geometries of the StressFreeStateGridFunction will then be used for all calculations
instead of the linear geometries of the piecewise linear grid.

In the 'mixed' case the volume term was simply forgotten, and not
actually handed to the assembler.

When dune-curvedgeometry is not found, the normal Derivative is set to zero.

Virtual Function is deprecated after dune 2.7.

Simon Praetorius' dune-vtk module has much better support for reading
and writing than my old hacks in dune-gfe.  Let's start adopting
dune-vtk for VTK I/O.

In this initial step, dune-vtk becomes a dependency of dune-gfe,
but only an optional one.  This may change.

Also adjust CosseratEnergyStiffness and MixedGFEAssembler for the use with a CompositeBasis instead of a PowerBasis

Up until dune-common 2.7, the Dune namespace contained shared_ptr:
First as an actual implementation, later simply as a namespace
import.  Dune-common post-2.7 has removed that import,
and users of shared_ptr really have to explicitly specify
the std:: namespace prefix.

Previously, they all inherited from LocalGeodesicFEStiffness,
which was weird, because none of them implemented first or
second derivatives as you would expect from something called
'stiffness'.

These two files are very similar, and merging them will reduce
maintenance effort in the future.  As of this patch, there is
now quite a number of preprocessor directives in the new merged
cosserat-continuum.cc.  That is not pretty, but the new state
is still much easier to maintain than two completely separate
files.  What is more, I hope to be able to reduce some of that
preprocessor trickery by smarter use of the dune-functions
basis trees.

This used to work only with dune-fufem bases, but the method
constructBoundaryDofs recently learned to also handle the bases
of dune-functions.

This used to be a yes/no, parameter, now we require 'cube'
rather that 'yes' or 'true'.

The code it contains is adolcnamespaceinjections.hh.

Capturing any polymorphic object by value is wrong,
and gcc 8 started to warn about this in the case
of exceptions.

The latter is deprecated as of dune-functions 2.6.

This is needed to build with a recent dune-common and dune-fufem

gcc warns about some of those nowadays.

Otherwise some dune-common stuff will not work properly,
for example copying of FieldMatrices that use adouble
as their number type.

The method indexSet has been removed from dune-functions
a long time ago.

This is planned to be used to start a fine-grid computation on the coarse-grid
solution, or to start a high-order computation from a low-order one.  This is
helpful, in particular, for measurements of the discretization error, which
take a lot of time.