Fix the computation of the Jacobian

[amontoison]

The function jac_coord_residual! returns the wrong result.
We need to use AD on the function residual! if we want accurate Jacobians for now.