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

It has been removed from dune-functions in the last release
(at least).

I want the interpolation rule to appear in the log files.
It is to easy to use the wrong one by accident.

Otherwise certain matrix-copying methods will not work anymore.

Makes the code a bit shorter and easier to read.

Only cmake is supported from now on.

Make them a parameter instead.  This allows to switch more easily between
geodesic finite elements and projected finite elements.

Make them a parameter instead.  This allows to switch more easily between
geodesic finite elements and projected finite elements.

Otherwise, the ParaView 'glyph' filter won't recognize the direction fields.

Such measurements are now centralized in the program compute-disc-error.cc.

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.

It is still no function in the dune-functions sense, but at least the
basis parameter needs to be a dune-functions basis now, rather than
a deprecated dune-fufem basis.

For those, it is necessary to have separate measurement for the displacement and orientation components.

Also, check whether the input files have the correct length

This creates a field with topological charge 2.