Used for calculating vibronic eigenvalues and eigenvectors in a Jahn-Teller active electronic state.  Can be used with and without spin-orbit coupling.