#line "./omp.f90" 1
      program main
      implicit none
      integer :: nkgmax, nr, i, j, k
      real(kind=8), allocatable, dimension(:) :: cont_wave
      real(kind=8), allocatable, dimension(:,:) :: vmat2D
      real(kind=8) :: tmp
      integer :: time1, time2, dt, count_rate, count_max
      real(kind=8) :: secs_acc

      nkgmax = 10000
      nr = 10000

      allocate(cont_wave(1:nkgmax*nr))


      cont_wave(:) = 0.d0
      do i = 1, nkgmax
        do j = 1, nr
            cont_wave((i-1)*nr + j) = dble(i - j) / dble(i + j)
        end do
      end do


      allocate(vmat2D(1:nkgmax, 1:nkgmax))


      call system_clock(count_max=count_max, count_rate=count_rate)
      call system_clock(time1)

      print*, 'Before OpenMP Offload'



    !$omp target teams distribute parallel do num_teams(220) map(to: con&
!$omp&t_wave) map(from: vmat2D)
      do i = 1, nkgmax

        do j = 1, nkgmax
            if (j > i) cycle
            tmp = 0.d0

            do k = 1, nr
                tmp = tmp + cont_wave((k-1)*nkgmax + i) * cont_wave((k-1&
     &)*nkgmax + j)
#line 44
            end do
            vmat2D(i, j) = tmp
            if (i /= j) vmat2D(j, i) = vmat2D(i, j)
        end do
      end do
    !$omp end target teams distribute parallel do



      call system_clock(time2)
      dt = time2 - time1
      secs_acc = real(dt) / real(count_rate)
      print*, 'Time in secs with OpenMP Offload: ', secs_acc

      print*, 'Min value in vmat2D: ', minval(vmat2D(1:nkgmax, 1:nkgmax)&
     &)
#line 59
      print*, 'Max value in vmat2D: ', maxval(vmat2D(1:nkgmax, 1:nkgmax)&
     &)
#line 60
      print*, 'Mean value in vmat2D: ', sum(vmat2D(1:nkgmax, 1:nkgmax)) &
     &/ dble(nkgmax * nkgmax)
      end program main