Working on the solver, fixed the HF-matrix and setupP()