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.
changed to `build_ratemat`