1! { dg-do run }
2
3module e_51_5_mod
4contains
5  subroutine init (v1, v2, N)
6    integer :: i, N
7    real :: v1(N), v2(N)
8    do i = 1, N
9      v1(i) = i + 2.0
10      v2(i) = i - 3.0
11    end do
12  end subroutine
13
14  subroutine check (p, N)
15    integer :: i, N
16    real, parameter :: EPS = 0.00001
17    real :: diff, p(N)
18    do i = 1, N
19      diff = p(i) - (i + 2.0) * (i - 3.0)
20      if (diff > EPS .or. -diff > EPS) call abort
21    end do
22  end subroutine
23
24  subroutine foo (p, v1, v2, N)
25    real, dimension(:) :: p, v1, v2
26    integer :: N
27    call init (v1, v2, N)
28    !$omp target data map(to: v1, v2, N) map(from: p)
29      call vec_mult (p, v1, v2, N)
30    !$omp end target data
31    call check (p, N)
32  end subroutine
33
34  subroutine vec_mult (p, v1, v2, N)
35    real, dimension(:) :: p, v1, v2
36    integer :: i, N
37    !$omp target map(to: v1, v2, N) map(from: p)
38      !$omp parallel do
39      do i = 1, N
40        p(i) = v1(i) * v2(i)
41      end do
42    !$omp end target
43  end subroutine
44end module
45
46program e_51_5
47  use e_51_5_mod, only : foo
48  integer, parameter :: N = 1024
49  real, allocatable, dimension(:) :: p, v1, v2
50  allocate(p(N), v1(N), v2(N))
51  call foo (p, v1, v2, N)
52  deallocate (p, v1, v2)
53end program
54