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.
Is it just me or is it confusing to be calling `buildK()` with "S" mode but dumping into a variable K?