Added a new C function to find the position and width of all cells in an Octree
Tidy up SPH importers and split construct_octree into two main functions, which allows the octree to be defined using a different set of particles compared to the one used for the density discretization.