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]]