Refactoring and cleaning the code.
Added a class for calualtion of the Coulom matrix. The Coulomb matrix
elements are now stored in a matrix instead of calculated on the fly.