I have implemented a conjugate gradient solver for solving CPHF equation through which I can reproduce SCF dipole polarizabilities results of psi4.
Things to do before this pull request can be merged:
[x] Make the solver a stand alone module and call it in CPHF.py
I have implemented a conjugate gradient solver for solving CPHF equation through which I can reproduce SCF dipole polarizabilities results of psi4. Things to do before this pull request can be merged: