Lanczos determination of ground (and excited) states updated with better handling of
special cases. if #required eigenvalues/vectors exceeds the dimension of the sector
or so does the block number then the hamiltonian is usually small (<10) and is solved
by direct diagonalization.
(P)ARPACK does not handle this case and produces an error so no eigenvalues/vectors
are given on output. This is because of the condition:
nev <= ncv < n
which is not respected.
This solution allows immediate extension to finite-temperature case as soon as I
find an ideal criteria to select the low-lying states.
some updates imported from devel_lanczos_multi-band branch, mainly about renaming of
some variables (the impurity GF), move of chi arrays and other small things.
LANCZOS modules are defined. PLAIN_LANCZOS has complex extension.
FULLED routines are moved to included files to make room to LANCED routines.
updated ED_BATH with new print format for ground-state and hamiltonian
updated drivers fulled_lda1b
Some adjustments in the fulled_lda1b code.
Removed double timer print for Norb*Nspin=1 case
BUG: updated the calculation of the spin susceptibility.
The dynamical response function is defined as the trace of the commutator of the operators, so that the susceptibility
corresponds to the retarded component and NOT to the time-ordered one, as calculated up to now.
Added calculation of the Matsubara component, to check Kubo theory from extrapolation.
The real-frequency Chi_spin is printed to a new file Chi_w.ed so not to overwrite existing file and allow comparisons.
Merge branch 'devel'