block_diagonalize fails to recognize unsimplified sympy Hamiltonians as diagonal

block_diagonalize raises ValueError: The off-diagonal elements of the unperturbed Hamiltonian must be zero when provided with a sympy Hamiltonian and subspace_eigenvectors. Seems to be because of unsimplified sympy expressions for the off-diagonal elements.

issue_offdiag.py

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information