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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
subroutine initialize_array_kernel_1(n, a)
integer, intent(in) :: n
real, intent(out) :: a(n)
integer :: i
!$omp parallel shared(a) private(i)
!$omp do
do i = 1, n
a(i) = 12.91
end do
!$omp end do
!$omp end parallel
print *, a(1), a(n), a(12), a(12841)
if (abs(a(1) - 12.91) > 1e-8) error stop
if (abs(a(n) - 12.91) > 1e-8) error stop
if (abs(a(12) - 12.91) > 1e-8) error stop
if (abs(a(12841) - 12.91) > 1e-8) error stop
end subroutine
subroutine initialize_array(n, b)
use omp_lib
implicit none
integer, intent(in) :: n
real, dimension(n), intent(out) :: b
integer :: i
call initialize_array_kernel_1(n, b)
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 subroutine
program openmp_13
use omp_lib
integer, parameter :: n = 100000
real, dimension(n) :: c
print *, c(2)
call omp_set_num_threads(4)
call initialize_array(n, c)
print *, c(1), c(n), c(12), c(12841)
if (abs(c(1) - 12.91) > 1e-8) error stop
if (abs(c(n) - 12.91) > 1e-8) error stop
if (abs(c(12) - 12.91) > 1e-8) error stop
if (abs(c(12841) - 12.91) > 1e-8) error stop
end program
|
