first debug completed.
Hubbard model FULLED is working. Crosscheck with devel branch same code, gives
perfect comparison, but the egs. This depends most probably on the fit,given
that the hamiltonian.restart files are very similar but not identical, whereas
all the local and impurity GF are identical within many digits.
I need to debug a little the CHI**2 fit routines.
Go on with more tests.
testing the PAM case with pd-impurity.
First results on the laptop look ok, using a not-so-weird-in-the-end
self-consistency in which all the components of the Delta function but the
pp one are set to zero. This is clear from Bethe lattice self-consistency of
the tdd_PAM (todo): the other components have the form Delta_dd,pd = (at_pp)**2/a**2t_ppG_dd,pd
that in the case of PAM (a=0) disappear.