Support MPI parallelism with BVH ray intersection

- [x] The `intersect` method of BVH must return global results by gathering all local results. Support also some special cases/optimization : only sequential results, closest or any intersection   
- [x] Implement a first MPI version,  a vector of ray distributed over all process. NOTE: this solution will be improved, here we apply an all_gather of all rays