Always calculate Hpsi

This gets rid of a release compiler warning
Implement Polak-Ribiere update for CG
Add more user friendly error message

When the integralF() results in non zero imaginary part.
Fix a bug, use Venq instead of nenq_pos
Multiply the functional derivative by L^3
Make the FE code work for the H.pseudo.gaussian