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.