Separate the Jacobian and Hessian calls into:

• sparsity computation: compute_jacobian_sparsity(row_indices, column_indices). Will be called once at the beginning.
• evaluation of numerical values: evaluate_jacobian_values(x, jacobian).
• same for Hessian.

cc @amontoison