diff --git a/numerical_methods/qr_eigen_values.cpp b/numerical_methods/qr_eigen_values.cpp
index 95ace03a7..4ced8591a 100644
--- a/numerical_methods/qr_eigen_values.cpp
+++ b/numerical_methods/qr_eigen_values.cpp
@@ -28,7 +28,6 @@ using qr_algorithm::operator<<;
 void create_matrix(std::valarray<std::valarray<double>> *A) {
     int i, j, tmp, lim2 = LIMS >> 1;
     int N = A->size();
-    std::srand(std::time(nullptr));
 
 #ifdef _OPENMP
 #pragma omp for
@@ -242,6 +241,9 @@ int main(int argc, char **argv) {
         return -1;
     }
 
+    // initialize random number generator
+    std::srand(std::time(nullptr));
+
     int i, rows = mat_size, columns = mat_size;
 
     std::valarray<std::valarray<double>> A(rows);