Use HVP linear map instead of Hessian

[gdalle]

Once DifferentiationInterface is integrated (#51), I would recommend that you use DI.hvp instead of DI.hessian to get the lazy, basis-free Hessian-vector product