This options enables to generate jacobian apply methods regardless
of the `matrix_free` option. If `matrix_free` is set to true,
`generate_jacobian_apply` is also set to true. This option replaces
`matrix_free` in all contexts except in the selection of the solver.

This is currently only used for level 4 instrumentation,
as these start and stop markers were messing heavily with
dependencies in the sum factorization code. This way, the necessary but
ugly work of adding instrumentation is separated much more
nicely.

It is not necessary anymore

The control operator doesn't accumulate in the residual. Instead the
localoperator has a member vector where the value of the derivative is
accumulated.

In that case you need two GridFunctions, which can be
bound to the two cells.

After consulting Steffen