program main integer :: n1, n2, i real(8) :: pi, x1, x2 n1 = 1000*1000 n2 = 0 call random_init(.true., .true.) !$omp parallel do reduction(+:n2) do i = 1, n1 call random_number(x1) call random_number(x2) if (x1**2 + x2**2 <= 1.0_8) n2 = n2 + 1 end do pi = 4.0_8*dble(n2)…

Replicable Monte Carlo simulation with OpenMP

PierU January 3, 2023, 1:42pm 5

But it’s probably processor-dependent. Which means that this kind of code is not portable.

Topic		Replies	Views
Random numbers in parallel	22	490	February 18, 2025
Random numbers across systems	13	1542	February 8, 2022
Generate same sequence of random numbers on a program with 1 or N tasks, using `random_number` and MPI	9	836	November 9, 2021
OpenMP stochastic simulation	31	433	September 19, 2024
Initializing RANDOM_NUMBER	9	919	April 2, 2021