SAD_Star: - Working HF solver for l>0 based on the matrix elements from Morten's code. TODO: clean the messy code...