UPDATE: code updated with new version of PARSE_INPUT and default values of CG fit adapted
to improved versino of OPTIMIZE_MINIMIZE.
Best CG fit is obtained using i) analytic evaluation of the gradient (i.e. for Delta fit mode), and
ii) keeping the tolerance around 1.d-5.

BUG [SOLVED]: a micro-bug hidden in the arpack interface +lanczos_arpack.f90

setting info=0 and resid=0.d0 forced the toss of a random vector INSIDE the
arpack routine, using the LAPACK interface. This made the user lose contact
with the starting, random vector.
This caused serious problems in the RDMFT application and in parallel calculations,
where the random vector was most proabably slightly different from call to call, giving rise
to slightly different values in the eigen-vectors, ultimately leading to DIFFERENT impurity
Green's functions with the same input bath (!!)

FIXED: we generate an identical input vector resid from fortran random number generator (KISS)
using seed=1234567 (the good old hirsch-fye qmc). This ensures that all processed, threads or sites
have identical inputy vectors, notwithstanding their situation.

small update: BHZ and BHZ-AFM updated drivers.