Judge unitary-ness relative to scaling.
We have some flaky tests when checking if eigenvectors are all orthonormal, since our error increases with the scaling (or largest eigenvalue) of the original matrix. Modify the test to increase the error tolerance.
Edited by Antonio Sánchez