method to compute the actual Hesse matrix. For some strange reason this slows down convergence, rather than speeding it up. to be investigated [[Imported from SVN: r7173]]
