That's already done in adolcnamespaceinjections.hh, in dune-fufem.

Capturing any polymorphic object by value is wrong,
and gcc 8 started to warn about this in the case
of exceptions.

The latter is deprecated as of dune-functions 2.6.

Otherwise some dune-common stuff will not work properly,
for example copying of FieldMatrices that use adouble
as their number type.

In particular, assemblers now hand LocalView objects around, rather than pairs of
elements and LocalFiniteElements.

[[Imported from SVN: r9997]]

[[Imported from SVN: r9913]]

[[Imported from SVN: r9912]]

[[Imported from SVN: r9877]]

That should be no problem now that we use a 50-decimal digits type
for the FD approximation.

[[Imported from SVN: r9876]]

[[Imported from SVN: r9875]]

It offers a 'double' interface, but we don't really need that.

[[Imported from SVN: r9874]]

[[Imported from SVN: r9873]]

[[Imported from SVN: r9872]]

[[Imported from SVN: r9870]]

This code compares scalar and vector AD, and a finite difference
approximation.

[[Imported from SVN: r9869]]

[[Imported from SVN: r9681]]

[[Imported from SVN: r9680]]

[[Imported from SVN: r9614]]

This still fails -- there must still be a bug in the Hessian computation

[[Imported from SVN: r9427]]

[[Imported from SVN: r9420]]

[[Imported from SVN: r9419]]

Hence many input configurations are tested now

[[Imported from SVN: r9418]]

They are easier to debug

[[Imported from SVN: r9417]]

[[Imported from SVN: r9413]]

[[Imported from SVN: r9411]]

[[Imported from SVN: r9410]]

[[Imported from SVN: r9409]]

[[Imported from SVN: r9408]]