1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// SPDX-FileCopyrightText: Copyright Contributors to the Kokkos project
#include <Kokkos_Macros.hpp>
#ifdef KOKKOS_ENABLE_EXPERIMENTAL_CXX20_MODULES
import kokkos.core;
#else
#include <Kokkos_Core.hpp>
#endif
namespace Test {
#ifndef KOKKOS_ENABLE_SYCL
__global__ void offset_streams(int* p) {
int idx = blockIdx.x * blockDim.x + threadIdx.x;
if (idx < 100) {
p[idx] += idx;
}
}
#endif
template <typename MemorySpace>
struct FunctorRange {
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> a;
FunctorRange(
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>>
a_)
: a(a_) {}
KOKKOS_INLINE_FUNCTION
void operator()(const int i) const { a(i) += 1; }
};
template <typename MemorySpace>
struct FunctorRangeReduce {
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> a;
FunctorRangeReduce(
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>>
a_)
: a(a_) {}
KOKKOS_INLINE_FUNCTION
void operator()(const int i, int& lsum) const { lsum += a(i); }
};
template <typename MemorySpace>
struct FunctorMDRange {
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> a;
FunctorMDRange(
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>>
a_)
: a(a_) {}
KOKKOS_INLINE_FUNCTION
void operator()(const int i, const int j) const { a(i * 10 + j) += 1; }
};
template <typename MemorySpace>
struct FunctorMDRangeReduce {
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> a;
FunctorMDRangeReduce(
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>>
a_)
: a(a_) {}
KOKKOS_INLINE_FUNCTION
void operator()(const int i, const int j, int& lsum) const {
lsum += a(i * 10 + j);
}
};
template <typename MemorySpace, typename ExecutionSpace>
struct FunctorTeam {
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> a;
FunctorTeam(
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>>
a_)
: a(a_) {}
KOKKOS_INLINE_FUNCTION
void operator()(
typename Kokkos::TeamPolicy<ExecutionSpace>::member_type const& team)
const {
int i = team.league_rank();
Kokkos::parallel_for(Kokkos::TeamThreadRange(team, 10),
[&](const int j) { a(i * 10 + j) += 1; });
}
};
template <typename MemorySpace, typename ExecutionSpace>
struct FunctorTeamReduce {
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>> a;
FunctorTeamReduce(
Kokkos::View<int*, MemorySpace, Kokkos::MemoryTraits<Kokkos::Unmanaged>>
a_)
: a(a_) {}
KOKKOS_INLINE_FUNCTION
void operator()(
typename Kokkos::TeamPolicy<ExecutionSpace>::member_type const& team,
int& lsum) const {
int i = team.league_rank();
int team_sum;
Kokkos::parallel_reduce(
Kokkos::TeamThreadRange(team, 10),
[&](const int j, int& tsum) { tsum += a(i * 10 + j); }, team_sum);
Kokkos::single(Kokkos::PerTeam(team), [&]() { lsum += team_sum; });
}
};
} // namespace Test
|
