diff --git a/src/source.cpp b/src/source.cpp
index 3f31541db..2ee8ecf13 100644
--- a/src/source.cpp
+++ b/src/source.cpp
@@ -292,6 +292,7 @@ void initialize_source()
 
   } else {
     // Generation source sites from specified distribution in user input
+    #pragma omp parallel for
     for (int64_t i = 0; i < simulation::work_per_rank; ++i) {
       // initialize random number seed
       int64_t id = simulation::total_gen*settings::n_particles +