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.
What's the benefit of wrapping non-parallel loops with nogil?