make sure the fd hessian is symmetric; hardwire Dexp(0) for more speed; use symmetry when computing the differential of the geodesic interpolation

[[Imported from SVN: r2016]]