Hello, [...] I would suggest to remove the phrase stating that the 2-by-2 diagonal blocks in H (corresponding to complex conjugate pairs of eigenvalues) are returned in standard form. That seems not to be true in my experiments. Best regards, Vasile