... with values expected to be given in the parameter file.

Scaling the trust region is the key to success, because in general the midsurface
deformation and the rotations will have values in different ranges.

That really was a hack: the build system should do it.

[[Imported from SVN: r10111]]

[[Imported from SVN: r10103]]

[[Imported from SVN: r10098]]

Seems to be a bug in IPOpt, or in our build system

[[Imported from SVN: r10086]]

This is much cleaner than handing over the Grid only, and then trying
to guess the basis using the preprocessor.  After all, the
RiemannianTrustRegionSolver calls the assembler, so it may know the
function space basis.

[[Imported from SVN: r10035]]

[[Imported from SVN: r9989]]

[[Imported from SVN: r9988]]

[[Imported from SVN: r9983]]

We may want to override options from the parameter file.

[[Imported from SVN: r9946]]

[[Imported from SVN: r9945]]

That code was written back in the days when we were only using
first order spaces.  Now, with other spaces, we have to properly
distinguish between grid vertices and Lagrange nodes.

[[Imported from SVN: r9911]]

[[Imported from SVN: r9910]]

[[Imported from SVN: r9904]]

[[Imported from SVN: r9888]]

[[Imported from SVN: r9879]]

[[Imported from SVN: r9862]]

For some reason, when numpy is used, I can't use ctrl-c to stop
the program.  But numpy is not really used anyway, so let's get
rid of it.

[[Imported from SVN: r9861]]

[[Imported from SVN: r9829]]

That makes it easier to keep them consistent.  Also, the code is shorter.

[[Imported from SVN: r9822]]

[[Imported from SVN: r9821]]

[[Imported from SVN: r9820]]

This used to be 'int', because dunepython didn't support the python-to-C
conversion for bool.

[[Imported from SVN: r9819]]

In particular, I can use math.cos in the initial values.

[[Imported from SVN: r9808]]

[[Imported from SVN: r9803]]

[[Imported from SVN: r9802]]

Newer versions of ADOL-C do not define 'overwrite' anymore, hence there
is no risk of macro clashing.

[[Imported from SVN: r9798]]

[[Imported from SVN: r9797]]

[[Imported from SVN: r9796]]

That way, we can now control what vertices are Dirichlet and Neumann without having
to recompile the code.

[[Imported from SVN: r9794]]

This will eventually allow to not hard-wired the Dirichlet vertices
into the C++ code, and change them without recompiling.  This also
means we can remove some preprocessor switches.

Currently, though, the Python code is still hard-wired, so we are
only half-way there.

Note that this patch may break several test problems.  I did not test
all of them.

[[Imported from SVN: r9793]]

[[Imported from SVN: r9789]]

[[Imported from SVN: r9787]]

[[Imported from SVN: r9780]]

[[Imported from SVN: r9775]]