! { dg-do run }

module target4
contains
  subroutine foo (a,m,n)
    integer :: m,n,i,j
    double precision :: a(m, n), t
    !$omp target data map(a) map(to: m, n)
    do i=1,n
      t = 0.0d0
      !$omp target map(t)
        !$omp parallel do reduction(+:t)
          do j=1,m
            t = t + a(j,i) * a(j,i)
          end do
      !$omp end target
      t = 2.0d0 * t
      !$omp target
        !$omp parallel do
          do j=1,m
            a(j,i) = a(j,i) * t
          end do
      !$omp end target
    end do
    !$omp end target data
  end subroutine foo
end module target4
  use target4, only : foo
  integer :: i, j
  double precision :: a(8, 9), res(8, 9)
  do i = 1, 8
    do j = 1, 9
      a(i, j) = i + j
    end do
  end do
  call foo (a, 8, 9)
  res = reshape ((/ 1136.0d0, 1704.0d0, 2272.0d0, 2840.0d0, 3408.0d0, 3976.0d0, &
&   4544.0d0, 5112.0d0, 2280.0d0, 3040.0d0, 3800.0d0, 4560.0d0, 5320.0d0, 6080.0d0, &
&   6840.0d0, 7600.0d0, 3936.0d0, 4920.0d0, 5904.0d0, 6888.0d0, 7872.0d0, 8856.0d0, &
&   9840.0d0, 10824.0d0, 6200.0d0, 7440.0d0, 8680.0d0, 9920.0d0, 11160.0d0, 12400.0d0, &
&   13640.0d0, 14880.0d0, 9168.0d0, 10696.0d0, 12224.0d0, 13752.0d0, 15280.0d0, 16808.0d0, &
&   18336.0d0, 19864.0d0, 12936.0d0, 14784.0d0, 16632.0d0, 18480.0d0, 20328.0d0, 22176.0d0, &
&   24024.0d0, 25872.0d0, 17600.0d0, 19800.0d0, 22000.0d0, 24200.0d0, 26400.0d0, 28600.0d0, &
&   30800.0d0, 33000.0d0, 23256.0d0, 25840.0d0, 28424.0d0, 31008.0d0, 33592.0d0, 36176.0d0, &
&   38760.0d0, 41344.0d0, 30000.0d0, 33000.0d0, 36000.0d0, 39000.0d0, 42000.0d0, 45000.0d0, &
&   48000.0d0, 51000.0d0 /), (/ 8, 9 /))
  if (any (a /= res)) stop 1
end