Follow-up from "WIP: automatic Peierls phase calculation"
The following discussion from !219 (merged) should be addressed:
-
@anton-akhmerov started a discussion: With not a lot of extra effort we could also support non-Abelian gauge fields. To do that we would need:
- To operate with exponentiated values (so to always compute
scipy.linalg.expm(1j * flux)
if flux-is matrix-valued, andnp.exp(1j * flux)
for scalars). - To use multiplication instead of addition in phase computation
- To operate with exponentiated values (so to always compute