5 subroutine init (v1
, v2
, N
)
14 subroutine check (p
, N
)
16 real, parameter :: EPS
= 0.00001
19 diff
= p(i
) - (i
+ 2.0) * (i
- 3.0)
20 if (diff
> EPS
.or
. -diff
> EPS
) call abort
24 subroutine vec_mult (N
)
25 real :: p(N
), v1(N
), v2(N
)
28 !$omp target data map(to: v1, v2) map(from: p)
41 use e_51_1_mod
, only
: vec_mult