Bug 854: fix numerical issue in SelfAdjointEigenSolver::computeDirect for 3x3 matrices. The tolerance to detect stable cross products was too optimistic.
Add respective unit tests.
Do not apply the preconditioner before starting the iterations as this might destroy a very good initial guess.