So far: 2d only

This will be used in the paper with Hanne, Philipp and Markus
on projection-based FE.

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.

In particular, this patch removes tests for the gradient
of the Cosserat energy and related things.  The
CosseratEnergyLocalStiffness class stopped computing
these things back in 2014 (in 63e42ece )
when it was decided that algorithmic differentiation was
the way to go.

With this patch, the test builds again, but it doesn't run yet.
I am afraid that fixing that may require a larger investment
of time...

- use more shared_ptr
- New MCMGMapper interface from dune-grid 2.6
- basis.indexSet() has been removed from dune-functions

Makes it faster, but only by 1 or 2 percent.

The new value is what you need for linear target spaces.
I am not convinced that this is sufficient, but it matches
what we do for simplices.

The whole quadrature business will need some theoretical
investigations eventually.

Makes the code a bit shorter and easier to read.

Only cmake is supported from now on.

A lot of raw pointers have been replaced by std::shared_ptr.

A lot of data members are inaccessible now, and getters have
to be used.

A small patch is still needed for dune-solvers: The iterationStep_
member is still needed to be public, because the new getter only
grants access to the object, not the pointer to it.

It seemed strange to me to demonstrate a second-order tangent field on a
first-order base function.  People may think that this is a necessary restriction.

Also, this patch includes a bugfix:  The basis of the tangent space was hard-wired
to (0,1,0) and (0,0,1).  But this only works if we are at the tangent space at (1,0,0).