UPDATE: code cleaned from fictitious 2imp (pam) routines.
Bath variables renamed. Everything is ready for multi-orbital version.
UPDATE: code further simplified to fit the SOLVER+drivers form.
Self-consistency variables removed from global solver memory.
These must be defined and read inside the calling drivers.
Checkout examples in drivers for more informations.
Tested for the PAM on Bethe and 2Dsquare lattice.
Code interfacing with H(k) Hamiltonian from Wannier90 style file.