Useless interpolation removed

Takes too much time on larger grids, and doesn't really change anything
even on relatively small grids.
Multipole expansion replaced by tree code

Now the potential calculation is more accurate, specially at the center
of the galaxy, as well as faster.