Hierarchical Eigensolver for Transition Matrices in Spectral Methods