NumPy 2.5 eig returns complex numbers even for real matrices—the release note's workaround loses the imaginary part

The NumPy 2.5.0 release notes state that the return values of eig and eigvals will always be complex, and that ...