just found out that eigenvectors() of hermitian matricies, when you have eigenvalues of >1 multiplicity, the basis construction of those eigenvalues does not ensure that the transpose of the system equals the inverse of the system.  switching to ev(similaritytransform(M), hermitianmatrix=true) instead