Comment 4 for bug 1043704

OK, so I can see currently that the PointSource::apply() method calls the Mesh::intersected_cell() method which in turn will call the appropriate method in the appropriate IntersectionOperator class which is dependent on the primitive type of the mesh.

So your suggestion is to move the logic you submitted in the patch to the Mesh::intersected_cell(const Point& point) method so that it will be used consistently across DOLFIN.