[[Imported from SVN: r9526]]

[[Imported from SVN: r9525]]

and move from GenericReferenceElement to ReferenceElement

[[Imported from SVN: r9524]]

The Weingarten map is block-diagonal for a product space.

[[Imported from SVN: r9517]]

It used to be a compile-term error, but now that we have ADOL-C support
the missing gradient is not a problem anymore.

[[Imported from SVN: r9516]]

[[Imported from SVN: r9510]]

[[Imported from SVN: r9509]]

[[Imported from SVN: r9508]]

This patch brings it back.

[[Imported from SVN: r9477]]

[[Imported from SVN: r9476]]

[[Imported from SVN: r9467]]

[[Imported from SVN: r9451]]

[[Imported from SVN: r9450]]

[[Imported from SVN: r9449]]

[[Imported from SVN: r9448]]

I can't make up my mind whether I want 'ctype' or 'field_type'

[[Imported from SVN: r9447]]

ADOL-C doesn't play nice with OpenMP.  Hence we have to disable
the latter for the time being.

[[Imported from SVN: r9446]]

[[Imported from SVN: r9445]]

When using the classes with ADOL-C, T becomes an 'adouble',
which is a non-trivial type.  However, the convexityField
is constexpr (and with a reason -- otherwise we could write
the initialization together with the declaration), and
constexpr and non-trivial types don't work together.

[[Imported from SVN: r9444]]

[[Imported from SVN: r9443]]

[[Imported from SVN: r9442]]

[[Imported from SVN: r9441]]

And set the inner loop solver verbosity to QUIET

[[Imported from SVN: r9435]]

This number is nowadays infered from the number of transfer operators.

[[Imported from SVN: r9432]]

Since ADOL-C is not multi-threaded we will not be using OpenMP
in the foreseeable future.  And on a single-threaded code the
OpenMP timer and the Dune timer return the same values anyway.

[[Imported from SVN: r9431]]

[[Imported from SVN: r9430]]

Offer setInitialIterate instead

[[Imported from SVN: r9429]]

[[Imported from SVN: r9428]]

Because ADOL-C doesn't like calls to sqrt(0), obviously.
The use of sqrt(x) for small x can be avoided completely
using modified series expansions, because the overall
function exp is differentiable even at x=0.

[[Imported from SVN: r9426]]

Makes for a more precise Hessian matrix when using automatic differentiation.

[[Imported from SVN: r9425]]

[[Imported from SVN: r9424]]

[[Imported from SVN: r9422]]

[[Imported from SVN: r9421]]

[[Imported from SVN: r9415]]

Because it may contain a projection onto the manifold, which
is part of the energy evaluation if we view the energy as
being defined on all of R^n (as ADOL-C does).

[[Imported from SVN: r9414]]

Not actually working yet, though.

[[Imported from SVN: r9412]]

[[Imported from SVN: r9405]]

[[Imported from SVN: r9404]]

[[Imported from SVN: r9403]]

[[Imported from SVN: r9402]]