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).

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.

By using the same mechanism as in the shell script:  piping everything through
the 'tee' program.  This may not be the most elegant way from a Python standpoint,
but it seems to work.

Since now my Python testsuite seems to do everything the shell version did,
this commit also removes the shell version.

Having it in bash was starting to be a pain.  This particular testsuite uses
a bit of arithmetic, which was difficult to do in bash, but trivial in Python.

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 still tested for the dune-fufem wrapper of a first-order basis, while nowadays
the RiemannianTRSolver class is instantiated with dune-functions bases.

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.

This is needed to be able to multiply FieldMatrices with Quaternions.

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