use which=LA/SA instead of LM/SM in for hermitian sparse matrices (added temporary debug printouts). This seems to be required to get the expected sorting of real eigenvalues of hermitian matrices