diff --git a/include/openmc/timer.h b/include/openmc/timer.h
index 9149c4ef7..b2d5f0005 100644
--- a/include/openmc/timer.h
+++ b/include/openmc/timer.h
@@ -21,6 +21,7 @@ extern Timer time_finalize;
 extern Timer time_inactive;
 extern Timer time_initialize;
 extern Timer time_read_xs;
+extern Timer time_sample_source;
 extern Timer time_tallies;
 extern Timer time_total;
 extern Timer time_transport;
diff --git a/src/output.cpp b/src/output.cpp
index c8b91c4dc..35cb9c00e 100644
--- a/src/output.cpp
+++ b/src/output.cpp
@@ -474,6 +474,8 @@ void print_runtime()
     show_time("Time synchronizing fission bank", time_bank.elapsed(), 1);
     show_time("Sampling source sites", time_bank_sample.elapsed(), 2);
     show_time("SEND/RECV source sites", time_bank_sendrecv.elapsed(), 2);
+  } else {
+    show_time("Time sampling source", time_sample_source.elapsed(), 1);
   }
   show_time("Time accumulating tallies", time_tallies.elapsed(), 1);
   show_time("Total time for finalization", time_finalize.elapsed());
diff --git a/src/simulation.cpp b/src/simulation.cpp
index eaed2badf..39c140a0d 100644
--- a/src/simulation.cpp
+++ b/src/simulation.cpp
@@ -462,7 +462,9 @@ void finalize_generation()
 
   } else if (settings::run_mode == RUN_MODE_FIXEDSOURCE) {
     // For fixed-source mode, we need to sample the external source
+    simulation::time_sample_source.start();
     fill_source_bank_fixedsource();
+    simulation::time_sample_source.stop();
   }
 }
 
diff --git a/src/source.cpp b/src/source.cpp
index d2e13a173..6314eb23f 100644
--- a/src/source.cpp
+++ b/src/source.cpp
@@ -328,6 +328,7 @@ void free_memory_source()
 void fill_source_bank_fixedsource()
 {
   if (settings::path_source.empty()) {
+    #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 + overall_generation()) *
diff --git a/src/state_point.cpp b/src/state_point.cpp
index 4c159c167..1e72049bf 100644
--- a/src/state_point.cpp
+++ b/src/state_point.cpp
@@ -283,6 +283,8 @@ openmc_statepoint_write(const char* filename, bool* write_source)
       write_dataset(runtime_group, "synchronizing fission bank", time_bank.elapsed());
       write_dataset(runtime_group, "sampling source sites", time_bank_sample.elapsed());
       write_dataset(runtime_group, "SEND-RECV source sites", time_bank_sendrecv.elapsed());
+    } else {
+      write_dataset(runtime_group, "sampling source sites", time_sample_source.elapsed());
     }
     write_dataset(runtime_group, "accumulating tallies", time_tallies.elapsed());
     write_dataset(runtime_group, "total", time_total.elapsed());
diff --git a/src/timer.cpp b/src/timer.cpp
index 2303c6090..616140aa9 100644
--- a/src/timer.cpp
+++ b/src/timer.cpp
@@ -16,6 +16,7 @@ Timer time_finalize;
 Timer time_inactive;
 Timer time_initialize;
 Timer time_read_xs;
+Timer time_sample_source;
 Timer time_tallies;
 Timer time_total;
 Timer time_transport;
@@ -68,6 +69,7 @@ void reset_timers()
   simulation::time_inactive.reset();
   simulation::time_initialize.reset();
   simulation::time_read_xs.reset();
+  simulation::time_sample_source.reset();
   simulation::time_tallies.reset();
   simulation::time_total.reset();
   simulation::time_transport.reset();