gcc/libgomp/testsuite/libgomp.fortran/reverse-offload-5a.f90
Tobias Burnus ea4b23d9c8 libgomp: Handle OpenMP's reverse offloads
This commit enabled reverse offload for nvptx such that gomp_target_rev
actually gets called.  And it fills the latter function to do all of
the following: finding the host function to the device func ptr and
copying the arguments to the host, processing the mapping/firstprivate,
calling the host function, copying back the data and freeing as needed.

The data handling is made easier by assuming that all host variables
either existed before (and are in the mapping) or that those are
devices variables not yet available on the host. Thus, the reverse
mapping can do without refcounts etc. Note that the spec disallows
inside a target region device-affecting constructs other than target
plus ancestor device-modifier and it also limits the clauses permitted
on this construct.

For the function addresses, an additional splay tree is used; for
the lookup of mapped variables, the existing splay-tree is used.
Unfortunately, its data structure requires a full walk of the tree;
Additionally, the just mapped variables are recorded in a separate
data structure an extra lookup. While the lookup is slow, assuming
that only few variables get mapped in each reverse offload construct
and that reverse offload is the exception and not performance critical,
this seems to be acceptable.

libgomp/ChangeLog:

	* libgomp.h (struct target_mem_desc): Predeclare; move
	below after 'reverse_splay_tree_node' and add rev_array
	member.
	(struct reverse_splay_tree_key_s, reverse_splay_compare): New.
	(reverse_splay_tree_node, reverse_splay_tree,
	reverse_splay_tree_key): New typedef.
	(struct gomp_device_descr): Add mem_map_rev member.
	* oacc-host.c (host_dispatch): NULL init .mem_map_rev.
	* plugin/plugin-nvptx.c (GOMP_OFFLOAD_get_num_devices): Claim
	support for GOMP_REQUIRES_REVERSE_OFFLOAD.
	* splay-tree.h (splay_tree_callback_stop): New typedef; like
	splay_tree_callback but returning int not void.
	(splay_tree_foreach_lazy): Define; like splay_tree_foreach but
	taking splay_tree_callback_stop as argument.
	* splay-tree.c (splay_tree_foreach_internal_lazy,
	splay_tree_foreach_lazy): New; but early exit if callback returns
	nonzero.
	* target.c: Instatiate splay_tree_c with splay_tree_prefix 'reverse'.
	(gomp_map_lookup_rev): New.
	(gomp_load_image_to_device): Handle reverse-offload function
	lookup table.
	(gomp_unload_image_from_device): Free devicep->mem_map_rev.
	(struct gomp_splay_tree_rev_lookup_data, gomp_splay_tree_rev_lookup,
	gomp_map_rev_lookup, struct cpy_data, gomp_map_cdata_lookup_int,
	gomp_map_cdata_lookup): New auxiliary structs and functions for
	gomp_target_rev.
	(gomp_target_rev): Implement reverse offloading and its mapping.
	(gomp_target_init): Init current_device.mem_map_rev.root.
	* testsuite/libgomp.fortran/reverse-offload-2.f90: New test.
	* testsuite/libgomp.fortran/reverse-offload-3.f90: New test.
	* testsuite/libgomp.fortran/reverse-offload-4.f90: New test.
	* testsuite/libgomp.fortran/reverse-offload-5.f90: New test.
	* testsuite/libgomp.fortran/reverse-offload-5a.f90: New test without
	mapping of on-device allocated variables.
2022-12-10 13:42:08 +01:00

98 lines
2.4 KiB
Fortran

! { dg-additional-options "-foffload-options=nvptx-none=-misa=sm_35" { target { offload_target_nvptx } } }
! Copying on-device allocated memory fails with cuMemcpyDtoHAsync error: invalid argument
! Hence, reverse-offload-5.f90 has been copied to *this* file,
! which uses on-host allocated vars - but only on the device side stack variables
implicit none
!$omp requires reverse_offload
integer, allocatable :: A(:), A2(:), s1, s2
integer :: i,D(4)
logical :: shared_mem
shared_mem = .false.
a = [1,2,3,4]
a2 = [8,7,6,5]
s1 = 532
s2 = 55
!$omp target map(to: shared_mem)
shared_mem = .true.
!$omp end target
!$omp target map(to: A, A2, s1, s2)
block
! stack variables:
integer :: ai(4), ai2(4), si1, si2
a = a * 2
a2 = a2 * 3
s1 = s1 * 4
s2 = s2 * 5
ai = [23,35,86,43]
ai2 = [8,4,7,1]
si1 = 64
si2 = 765
!$omp target device (ancestor:1) &
!$omp& map(to: A, s1, ai, si1) map(always, to: a2, s2) &
!$omp& map(tofrom: ai2, si2)
if (shared_mem) then
if (any (a /= 2 * [1,2,3,4])) stop 1
if (s1 /= 4 * 532) stop 2
else
if (any (a /= [1,2,3,4])) stop 3
if (s1 /= 532) stop 4
endif
if (any (a2 /= 3 * [8,7,6,5])) stop 5
if (s2 /= 5 * 55) stop 6
if (any (ai /= [23,35,86,43])) stop 7
if (any (ai2 /= [8,4,7,1])) stop 8
if (si1 /= 64) stop 9
if (si2 /= 765) stop 10
a = a*3
a2 = a2*7
s1 = s1*11
s2 = s2*5
ai = ai*13
ai2 = ai2*21
si1 = si1*27
si2 = si2*31
!$omp end target
if (shared_mem) then
if (any (a /= 3 * 2 * [1,2,3,4])) stop 11
if (any (a2 /= 7 * 3 * [8,7,6,5])) stop 12
if (s1 /= 11 * 4 * 532) stop 13
if (s2 /= 5 * 5 * 55) stop 14
if (any (ai /= 13 * [23,35,86,43])) stop 15
if (si1 /= 27 * 64) stop 16
else
if (any (a /= 2 * [1,2,3,4])) stop 17
if (any (a2 /= 3 * [8,7,6,5])) stop 18
if (s1 /= 4 * 532) stop 19
if (s2 /= 5 * 55) stop 20
if (any (ai /= [23,35,86,43])) stop 22
if (si1 /= 64) stop 23
endif
if (any (ai2 /= 21 * [8,4,7,1])) stop 24
if (si2 /= 31 * 765) stop 25
end block
if (shared_mem) then
if (any (a /= 3 * 2 * [1,2,3,4])) stop 30
if (any (a2 /= 7 * 3 * [8,7,6,5])) stop 31
if (s1 /= 11 * 4 * 532) stop 32
if (s2 /= 5 * 5 * 55) stop 33
else
if (any (a /= 3 * [1,2,3,4])) stop 34
if (any (a2 /= 3 * 7 * [8,7,6,5])) stop 35
if (s1 /= 11 * 532) stop 36
if (s2 /= 5 * 5 * 55) stop 37
endif
deallocate (a, a2, s1, s2)
end