#include #include "rng-omp.h" int main() { int const n_threads = 4; omp_set_num_threads( n_threads ); RNG rng( n_threads ); int const DIM_ENSEMBLE = 10; #pragma omp parallel for for(int i = 0; i < DIM_ENSEMBLE; ++i) { #pragma omp critical std::cout << "hello from thread " << omp_get_thread_num() << " with random number " << rng() << std::endl; } }