Stops segfaults when rho iteration is screwed up and generates lots of NaNs
Also trapped some Nan matrices in diagonaliseSvd.