update:
found a bug in the stopping condition for iloop > nloop. fixed in all drivers.

update.
- evaluates spin-spin and density-density correlations <n_a,n_b>
- introduced numbering in the columns.info file from greater readibility.
bugs fixed:
- now evaluates spin-spin <sz_a,sz_b> rather than moment-moment <m_a,m_b> correlation
function.
- found a bug in the printing routine for the irreducible bath chi^2 fit:
always printed G0 field rather than G0 or \Delta depending on the value of
cg_scheme.
- new bath initialization: generate a smart bath in the range -Hband:Hband (input of
ed_init_solver, default=2) and optionally read bath from file so that
restarting from different bath levels is possible. bath center is at 1.d-4.
Code entirely written.
First debug successfull and tested against Massimo Capone lanczos code (MClanc).
Results at first DMFT loop for U=0 and U=1 give very same spectrum.
Green's functions are found in very good agreement and even better
fullfillment of the symmetry (i.e. Re_G or Re_Sigma identically zero).

Initial results for fully converged DMFT calculation obtained for U=2.
Order parameter is found to be on top of MClanc code.

To be fixed:
- Code to be cleaned.
- Delta functions routines in ED_BATH to be replaced with a unique subroutine
for real or matsubara frequencies.
- not interacting Grenn's function (impG0 and impF0) of the Anderson problem looks wrong.
Inversion to be check out.

Test against F. Petocchi results for bethe lattice in order.
Code cleaned. Ready to be merged with the +devel branch.
Delta functions routines have been fixed: there are now 4routines:
2 to get normal component Delta (matsubara,real),
2 to get anomalous component FDelta (matsubara,real).

bug fixed: a bug in the calculation of delta_imp observables was found
and fixed in ED_OBSERVABLES.

Benchmarks against M.Capone and F.Petocchi data are very satisfactory.
update:
code little revisioned, some bug fixed.
- order parameter +phi was evaluated as twice is value
- an sign error in the definition of the anomalous hybridization functio on the real-axis
(fdelta_bath_real) caused a wrong evaluation of the self-enrergy (impSAreal) and thus
of the spectral density at the self-consistency.
- bug in the evaluation of the G0,F0 has been fixed (wrong definition of the determinant).

Test against IPT real-axis code show very good results at zero temperature.
Keep testing at finite temperature.
Merge branch 'devel_sc' into merge_devel_and_develSC

Conflicts:
	.gitignore
	ED_AUX_FUNX.f90
	ED_CHI2FIT.f90
	ED_HAMILTONIAN.f90
	ED_OBSERVABLES.f90
	drivers/ed_tddpam_bethe.f90
Merge branch 'devel_sc_gm' of github.com:aamaricci/Exact-Diagonalization-DMFT into devel_sc_gm

Conflicts:
	ED_BATH.f90
	ED_CHI2FIT.f90
	ED_GREENS_FUNCTIONS.f90
	ED_OBSERVABLES.f90
	drivers/ed_ahm_bethe.f90
	drivers/ed_tddpam_bethe.f90
Pull and merge from devel_sc_gm in GITHUB.
Some minor modifications to compile the code
Code slightly changed and tested against Dualla (Giacomo Mazza Version of the code)
with identical results.

We noticed indeed a small change with respect to a previous version of SciFor library.
(not necessarily an error).

We use this version to test the devel branch.