by keeping in memory psi,hpsi,spsi and making linear-algebra operations
to keep them updated rather just applying the routine.

adding the calculation of nhpsi in cegterg (optionally output), and some similar printout in paro_gamma

call from cb_paro_gamma_main modified accordingly

rm specific solver_param module (just containing DP and stdout) in favuor of  util_param in UtilXlib

adding a separate directory for DENSE matrix diagonalization which is used by other solvers as prelimiary or finale stage

addition of references for paro_gamma
update to run_all_examples script

change to pcg.f90 in ParO to force the residual and search vectors to be orthogonal to the occupied manyfold.
Not totally sure of the S|psi0Xpsi0| or |psi0Xpsi0|S order