It is so uncompileable that I am not able to fix it quickly.  Since I don't
actually need it right now I let it wait for another day.

This initial commit simply brings a cleaned-up version of the harmonicmaps program.
The first step will be to make it compute curve-shortening flow on the sphere.

Here is a new program that does nothing but compute discretization errors:
that is L2 and H1 norms of differences between two functions.  The reference
function can be either discrete or given in closed form (the latter is not
enabled yet).

So far I have only used this program to measure discretization errors of
Cosserat shell models, and I did *not* get the error behavior I expected.
So either the error behavior of such shells is weird, or this new measurement
program is still buggy.

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

[[Imported from SVN: r10095]]

While we are there, set all the flags in a loop, so we don't have
to keep copying boilerplate code as the number of programs to build
increases.

[[Imported from SVN: r10074]]

[[Imported from SVN: r10061]]