bugfix: we need the derivative of the geodesic interpolation in the gradient computation (chain rule), not the derivative of the velocity interpolation

[[Imported from SVN: r4066]]