diff --git a/src/VCluster/VCluster.cpp b/src/VCluster/VCluster.cpp
index 078a6232c6fd4b0e798b7a3707046d10704ce8f4..baaccc4845fdeb657cf2ac8077f62068a0702022 100644
--- a/src/VCluster/VCluster.cpp
+++ b/src/VCluster/VCluster.cpp
@@ -6,9 +6,13 @@
 #endif
 
 #include "util/print_stack.hpp"
+#include "util/math_util_complex.hpp"
 
 Vcluster * global_v_cluster_private = NULL;
 
+//
+std::vector<int> sieve_spf;
+
 // number of vcluster instances
 size_t n_vcluster = 0;
 bool ofp_initialized = false;
diff --git a/src/VCluster/VCluster.hpp b/src/VCluster/VCluster.hpp
index 7f66ea6840401d319c5834fe10157fa4be0ef630..8b71d56b171fc6b50ebf1f9224484cc9305e9a92 100644
--- a/src/VCluster/VCluster.hpp
+++ b/src/VCluster/VCluster.hpp
@@ -9,8 +9,10 @@
 #define VCLUSTER_HPP
 
 #include <signal.h>
+
 #include "VCluster_base.hpp"
 #include "VCluster_meta_function.hpp"
+#include "util/math_util_complex.hpp"
 
 void bt_sighandler(int sig, siginfo_t * info, void * ctx);
 
@@ -873,6 +875,9 @@ static inline void openfpm_init(int *argc, char ***argv)
 	if (*argc != 0)
 		program_name = std::string(*argv[0]);
 
+	// Initialize math pre-computation tables
+	openfpm::math::init_getFactorization();
+
 	ofp_initialized = true;
 }