`lazy_jacobian` is less generic than `pullback_function`

[gdalle]

The result J of lazy_jacobian can supposedly be used as v' * J for vector-Jacobian products, but the tranpose operator doesn't work if v is, for instance, a 3D array