Use the polarization identity from Absil, Mahony, Sepulchre (page 107) to compute an fd approximation of the Riemannian Hessian

[[Imported from SVN: r6431]]