Use symmetric eigensolver for diagonalizing rate matrix

With my BLAS, and a 500x500 matrix, the hermetian eigensolver is about 5x to 10x faster than the nonsymmetric one.
Maybe just add a comment here, that's all.