Nested Do Loops

program OMP_Nested_Do

! To run this program: gfortran -fopenmp omp_nested_do.f95; ./a.out

use omp_lib

implicit none

integer ( kind = 4 ) thread_num,num_thread,proc_num,thread_id ! OpenMP Variables

integer:: i,j,Nx,Ny

real*8 sum

! Define Number of Threads.

proc_num = omp_get_num_procs()

thread_num = proc_num

call omp_set_num_threads (thread_num)

write(*,*)

write(*,*) "Number of Processors=", proc_num, "Number of Threads=", thread_num

write(*,*)

Nx = 10

Ny = 5

sum = 0.0d0

!$OMP PARALLEL SHARED(Nx,Ny) PRIVATE(i,j,thread_id)

!$OMP DO

do i = 1,Nx

do j = 1,Ny

thread_id = omp_get_thread_num()

write(*,*) "Thread",thread_id,"says - I am loop index",i,j

enddo

!$OMP END DO

!$OMP END PARALLEL

write(*,*)

!$OMP PARALLEL SHARED(Nx,Ny) PRIVATE(i,j)

!$OMP DO REDUCTION (+:sum)

do i = 1,Nx

do j = 1,Ny

sum = sum + 1.0d0

enddo

!$OMP END PARALLEL

write(*,*) "Total number of points in",Nx,"X",Ny,"grid is", int(sum)

write(*,*)

end program OMP_Nested_Do