As discussed in #296, defining the following operators would allow some backends to optimize the joint computation: - `value_derivative_and_second_derivative` - `value_gradient_and_hessian` Main issue: with what kind of scenario do we test those?