adds option `generate_jacobian_apply`

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.