2022-11-04 09:35:42 -07:00

24 lines
387 B
Fortran

! @@name: reduction.7
! @@type: F-free
! @@operation: compile
! @@expect: success
! @@version: pre_omp_3.0
program array_red
integer,parameter :: n=100
integer :: j
real :: a(n), b(n,n)
call init(n,b)
a(:) = 0.0e0
!$omp parallel do reduction(+:a)
do j = 1, n
a(:) = a(:) + b(:,j)
end do
print*, " a(1) a(n): ", a(1), a(n)
end program