Specific preparation for in-place operators?

[gdalle]

We could introduce a prepare_jacobian! in addition to prepare_jacobian.
Advantages: we signal that mutation is allowed, which allows us to generate fewer basis vectors (and possibly exploit a sparsity pattern).