compute_jacobian