File: openmp_32.f90

package info (click to toggle)
lfortran 0.45.0-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 46,332 kB
  • sloc: cpp: 137,068; f90: 51,260; python: 6,444; ansic: 4,277; yacc: 2,285; fortran: 806; sh: 524; makefile: 30; javascript: 15
file content (20 lines) | stat: -rw-r--r-- 408 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
program openmp_32
    use omp_lib
    implicit none
    real :: phi
    integer :: j, i

    phi = 10124.142

    call omp_set_num_threads(4)
    do i = 1, 10
        !$omp parallel do private(j) reduction(+:phi)
        do j = 1, 100
            phi = phi + 1.0
        end do
        !$omp end parallel do
    end do

    print *, phi
    if (abs(phi - 11124.1416) > 1e-8) error stop
end program openmp_32