1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
program openmp_24
use omp_lib
integer, parameter :: n = 100000
integer :: i
real, pointer :: b(:)
allocate(b(n))
call omp_set_num_threads(4)
!$omp parallel shared(b) private(i)
!$omp do
do i = 1, n
b(i) = 12.91
end do
!$omp end do
!$omp end parallel
print *, b(1), b(n)
print *, b(1), b(n), b(12), b(12841)
if (abs(b(1) - 12.91) > 1e-8) error stop
if (abs(b(n) - 12.91) > 1e-8) error stop
if (abs(b(12) - 12.91) > 1e-8) error stop
if (abs(b(12841) - 12.91) > 1e-8) error stop
end program
|
