Not sure whether useful here but I recently learned about the Galois package, which was useful for my purpose! mhostetter.github.io/galois/latest/ it is quite nicely integrated with numpy.

Thanks for the shout-out! If you liked this one, there are two follow-ups that might be interesting as well :)
Unifying unitary synthesis techniques based on Cartan decompositions: pennylane.ai/qml/demos/tu... and KAK-based compilation of Hamiltonian simulation: pennylane.ai/qml/demos/tu...

If the dynamical Lie algebra of the Hamiltonian is small enough, this gives you a useful time evolution circuit with time-independent(!) depth. Pretty nice 😎 If you like quantum complexity theory more than Lie theory, this is fast forwarding a Hamiltonian in practice 🙂