instead of a singular value decomposition.  This makes computing
a Hesse matrix roughly 20% faster.

[[Imported from SVN: r7232]]