//~ TODO
//-
/**
 *
 [-] Use arena for matrices instead of malloc maybe
 [DONE] Factor l-dependent matrix
 [DONE] Print eigenfunctions for plotting
 [DONE] Solve for more shells
 [-] Make some good orbital/atom abstraction
 [-] Then try to make hartree fock work maybe.
 [-] Need to make a Complex double (Z64) version if I want to go relativistic
 [-] Smarter way to multiply bsplines and do gauss integration?
 *
 */