Speed up facility sampling (select_facility())

[Hussein-Mahfouz]

select_facility is slow. Script 3.3_assign_facility_all.py takes up the large majority of the pipeline runtime. The select_facility function should be improved. Parallelisation should probably be the first approach to try