The entire code has been extended to multi-orbital case in the 'irreducible' representation of the Bath.
This requires that NO correlated orbital-orbital 'hybridizations' are present, which would requires
mixed gf calculations and full bath representation. [this is actually under development in the _nomix branch].
The code was tested for single and two-band case in the Bethe Lattice (easiest) for both FULLED and LANCED channels.
Results are perfect, though some care is needed when solving multi-band with Lanczos: converged is harder to get!
In general I noticed that a large fake-Temperature is needed to converge (that is the code find it hard to resolve
finer energy scales).

I coded the full Kanamori type Hamiltonian with:
-Uloc = same orbital, opposite spins interaction term
-Ust  = different orbital, opposite spins interaction term
-Und  = (Ust-J) different orbital, parallel spins interaction term (Hund's effect)
-Jh   = Spin-Exchange and Pair-Hopping terms, turned off with Jhflag=F
NB: I DID NOT TEST ALL THE INTERACTIONS  TERM, THIS REQUIRE ADDITIONAL NECESSARY WORK.

-->  FURTHER TETS REQUIRED BEFORE TO MERGE WITH STABLE-DEVEL  <--
I would like to perform some additional tests for stability and reliability before
merging with the stable development branch and then to stable branch master.
This is a paranoid-security procedure, to ensure that this codes produces the very same results of
the _devel branch code (one-band).