1 | // -*- C++ -*-
|
---|
2 | //===-- glue_numeric_impl.h -----------------------------------------------===//
|
---|
3 | //
|
---|
4 | // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
|
---|
5 | // See https://llvm.org/LICENSE.txt for license information.
|
---|
6 | // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
|
---|
7 | //
|
---|
8 | //===----------------------------------------------------------------------===//
|
---|
9 |
|
---|
10 | #ifndef _PSTL_GLUE_NUMERIC_IMPL_H
|
---|
11 | #define _PSTL_GLUE_NUMERIC_IMPL_H
|
---|
12 |
|
---|
13 | #include <functional>
|
---|
14 |
|
---|
15 | #include "utils.h"
|
---|
16 | #include "numeric_fwd.h"
|
---|
17 |
|
---|
18 | namespace std
|
---|
19 | {
|
---|
20 |
|
---|
21 | // [reduce]
|
---|
22 |
|
---|
23 | template <class _ExecutionPolicy, class _ForwardIterator, class _Tp, class _BinaryOperation>
|
---|
24 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _Tp>
|
---|
25 | reduce(_ExecutionPolicy&& __exec, _ForwardIterator __first, _ForwardIterator __last, _Tp __init,
|
---|
26 | _BinaryOperation __binary_op)
|
---|
27 | {
|
---|
28 | return transform_reduce(std::forward<_ExecutionPolicy>(__exec), __first, __last, __init, __binary_op,
|
---|
29 | __pstl::__internal::__no_op());
|
---|
30 | }
|
---|
31 |
|
---|
32 | template <class _ExecutionPolicy, class _ForwardIterator, class _Tp>
|
---|
33 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _Tp>
|
---|
34 | reduce(_ExecutionPolicy&& __exec, _ForwardIterator __first, _ForwardIterator __last, _Tp __init)
|
---|
35 | {
|
---|
36 | return transform_reduce(std::forward<_ExecutionPolicy>(__exec), __first, __last, __init, std::plus<_Tp>(),
|
---|
37 | __pstl::__internal::__no_op());
|
---|
38 | }
|
---|
39 |
|
---|
40 | template <class _ExecutionPolicy, class _ForwardIterator>
|
---|
41 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy,
|
---|
42 | typename iterator_traits<_ForwardIterator>::value_type>
|
---|
43 | reduce(_ExecutionPolicy&& __exec, _ForwardIterator __first, _ForwardIterator __last)
|
---|
44 | {
|
---|
45 | typedef typename iterator_traits<_ForwardIterator>::value_type _ValueType;
|
---|
46 | return transform_reduce(std::forward<_ExecutionPolicy>(__exec), __first, __last, _ValueType{},
|
---|
47 | std::plus<_ValueType>(), __pstl::__internal::__no_op());
|
---|
48 | }
|
---|
49 |
|
---|
50 | // [transform.reduce]
|
---|
51 |
|
---|
52 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _Tp>
|
---|
53 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _Tp>
|
---|
54 | transform_reduce(_ExecutionPolicy&& __exec, _ForwardIterator1 __first1, _ForwardIterator1 __last1,
|
---|
55 | _ForwardIterator2 __first2, _Tp __init)
|
---|
56 | {
|
---|
57 | typedef typename iterator_traits<_ForwardIterator1>::value_type _InputType;
|
---|
58 | return __pstl::__internal::__pattern_transform_reduce(
|
---|
59 | std::forward<_ExecutionPolicy>(__exec), __first1, __last1, __first2, __init, std::plus<_InputType>(),
|
---|
60 | std::multiplies<_InputType>(),
|
---|
61 | __pstl::__internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
62 | __exec),
|
---|
63 | __pstl::__internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
64 | __exec));
|
---|
65 | }
|
---|
66 |
|
---|
67 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _Tp, class _BinaryOperation1,
|
---|
68 | class _BinaryOperation2>
|
---|
69 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _Tp>
|
---|
70 | transform_reduce(_ExecutionPolicy&& __exec, _ForwardIterator1 __first1, _ForwardIterator1 __last1,
|
---|
71 | _ForwardIterator2 __first2, _Tp __init, _BinaryOperation1 __binary_op1, _BinaryOperation2 __binary_op2)
|
---|
72 | {
|
---|
73 | return __pstl::__internal::__pattern_transform_reduce(
|
---|
74 | std::forward<_ExecutionPolicy>(__exec), __first1, __last1, __first2, __init, __binary_op1, __binary_op2,
|
---|
75 | __pstl::__internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
76 | __exec),
|
---|
77 | __pstl::__internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
78 | __exec));
|
---|
79 | }
|
---|
80 |
|
---|
81 | template <class _ExecutionPolicy, class _ForwardIterator, class _Tp, class _BinaryOperation, class _UnaryOperation>
|
---|
82 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _Tp>
|
---|
83 | transform_reduce(_ExecutionPolicy&& __exec, _ForwardIterator __first, _ForwardIterator __last, _Tp __init,
|
---|
84 | _BinaryOperation __binary_op, _UnaryOperation __unary_op)
|
---|
85 | {
|
---|
86 | return __pstl::__internal::__pattern_transform_reduce(
|
---|
87 | std::forward<_ExecutionPolicy>(__exec), __first, __last, __init, __binary_op, __unary_op,
|
---|
88 | __pstl::__internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator>(__exec),
|
---|
89 | __pstl::__internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator>(__exec));
|
---|
90 | }
|
---|
91 |
|
---|
92 | // [exclusive.scan]
|
---|
93 |
|
---|
94 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _Tp>
|
---|
95 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
96 | exclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
97 | _ForwardIterator2 __result, _Tp __init)
|
---|
98 | {
|
---|
99 | using namespace __pstl;
|
---|
100 | return __internal::__pattern_transform_scan(
|
---|
101 | std::forward<_ExecutionPolicy>(__exec), __first, __last, __result, __pstl::__internal::__no_op(), __init,
|
---|
102 | std::plus<_Tp>(), /*inclusive=*/std::false_type(),
|
---|
103 | __internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(__exec),
|
---|
104 | __internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(__exec));
|
---|
105 | }
|
---|
106 |
|
---|
107 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _Tp, class _BinaryOperation>
|
---|
108 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
109 | exclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
110 | _ForwardIterator2 __result, _Tp __init, _BinaryOperation __binary_op)
|
---|
111 | {
|
---|
112 | using namespace __pstl;
|
---|
113 | return __internal::__pattern_transform_scan(
|
---|
114 | std::forward<_ExecutionPolicy>(__exec), __first, __last, __result, __pstl::__internal::__no_op(), __init,
|
---|
115 | __binary_op, /*inclusive=*/std::false_type(),
|
---|
116 | __internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(__exec),
|
---|
117 | __internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(__exec));
|
---|
118 | }
|
---|
119 |
|
---|
120 | // [inclusive.scan]
|
---|
121 |
|
---|
122 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2>
|
---|
123 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
124 | inclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
125 | _ForwardIterator2 __result)
|
---|
126 | {
|
---|
127 | typedef typename iterator_traits<_ForwardIterator1>::value_type _InputType;
|
---|
128 | return transform_inclusive_scan(std::forward<_ExecutionPolicy>(__exec), __first, __last, __result,
|
---|
129 | std::plus<_InputType>(), __pstl::__internal::__no_op());
|
---|
130 | }
|
---|
131 |
|
---|
132 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _BinaryOperation>
|
---|
133 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
134 | inclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
135 | _ForwardIterator2 __result, _BinaryOperation __binary_op)
|
---|
136 | {
|
---|
137 | return transform_inclusive_scan(std::forward<_ExecutionPolicy>(__exec), __first, __last, __result, __binary_op,
|
---|
138 | __pstl::__internal::__no_op());
|
---|
139 | }
|
---|
140 |
|
---|
141 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _Tp, class _BinaryOperation>
|
---|
142 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
143 | inclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
144 | _ForwardIterator2 __result, _BinaryOperation __binary_op, _Tp __init)
|
---|
145 | {
|
---|
146 | return transform_inclusive_scan(std::forward<_ExecutionPolicy>(__exec), __first, __last, __result, __binary_op,
|
---|
147 | __pstl::__internal::__no_op(), __init);
|
---|
148 | }
|
---|
149 |
|
---|
150 | // [transform.exclusive.scan]
|
---|
151 |
|
---|
152 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _Tp, class _BinaryOperation,
|
---|
153 | class _UnaryOperation>
|
---|
154 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
155 | transform_exclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
156 | _ForwardIterator2 __result, _Tp __init, _BinaryOperation __binary_op,
|
---|
157 | _UnaryOperation __unary_op)
|
---|
158 | {
|
---|
159 | return __pstl::__internal::__pattern_transform_scan(
|
---|
160 | std::forward<_ExecutionPolicy>(__exec), __first, __last, __result, __unary_op, __init, __binary_op,
|
---|
161 | /*inclusive=*/std::false_type(),
|
---|
162 | __pstl::__internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
163 | __exec),
|
---|
164 | __pstl::__internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
165 | __exec));
|
---|
166 | }
|
---|
167 |
|
---|
168 | // [transform.inclusive.scan]
|
---|
169 |
|
---|
170 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _BinaryOperation,
|
---|
171 | class _UnaryOperation, class _Tp>
|
---|
172 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
173 | transform_inclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
174 | _ForwardIterator2 __result, _BinaryOperation __binary_op, _UnaryOperation __unary_op,
|
---|
175 | _Tp __init)
|
---|
176 | {
|
---|
177 | return __pstl::__internal::__pattern_transform_scan(
|
---|
178 | std::forward<_ExecutionPolicy>(__exec), __first, __last, __result, __unary_op, __init, __binary_op,
|
---|
179 | /*inclusive=*/std::true_type(),
|
---|
180 | __pstl::__internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
181 | __exec),
|
---|
182 | __pstl::__internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
183 | __exec));
|
---|
184 | }
|
---|
185 |
|
---|
186 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _UnaryOperation,
|
---|
187 | class _BinaryOperation>
|
---|
188 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
189 | transform_inclusive_scan(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
190 | _ForwardIterator2 __result, _BinaryOperation __binary_op, _UnaryOperation __unary_op)
|
---|
191 | {
|
---|
192 | if (__first != __last)
|
---|
193 | {
|
---|
194 | auto __tmp = __unary_op(*__first);
|
---|
195 | *__result = __tmp;
|
---|
196 | return transform_inclusive_scan(std::forward<_ExecutionPolicy>(__exec), ++__first, __last, ++__result,
|
---|
197 | __binary_op, __unary_op, __tmp);
|
---|
198 | }
|
---|
199 | else
|
---|
200 | {
|
---|
201 | return __result;
|
---|
202 | }
|
---|
203 | }
|
---|
204 |
|
---|
205 | // [adjacent.difference]
|
---|
206 |
|
---|
207 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2, class _BinaryOperation>
|
---|
208 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
209 | adjacent_difference(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
210 | _ForwardIterator2 __d_first, _BinaryOperation __op)
|
---|
211 | {
|
---|
212 |
|
---|
213 | if (__first == __last)
|
---|
214 | return __d_first;
|
---|
215 |
|
---|
216 | return __pstl::__internal::__pattern_adjacent_difference(
|
---|
217 | std::forward<_ExecutionPolicy>(__exec), __first, __last, __d_first, __op,
|
---|
218 | __pstl::__internal::__is_vectorization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
219 | __exec),
|
---|
220 | __pstl::__internal::__is_parallelization_preferred<_ExecutionPolicy, _ForwardIterator1, _ForwardIterator2>(
|
---|
221 | __exec));
|
---|
222 | }
|
---|
223 |
|
---|
224 | template <class _ExecutionPolicy, class _ForwardIterator1, class _ForwardIterator2>
|
---|
225 | __pstl::__internal::__enable_if_execution_policy<_ExecutionPolicy, _ForwardIterator2>
|
---|
226 | adjacent_difference(_ExecutionPolicy&& __exec, _ForwardIterator1 __first, _ForwardIterator1 __last,
|
---|
227 | _ForwardIterator2 __d_first)
|
---|
228 | {
|
---|
229 | typedef typename iterator_traits<_ForwardIterator1>::value_type _ValueType;
|
---|
230 | return adjacent_difference(std::forward<_ExecutionPolicy>(__exec), __first, __last, __d_first,
|
---|
231 | std::minus<_ValueType>());
|
---|
232 | }
|
---|
233 |
|
---|
234 | } // namespace std
|
---|
235 |
|
---|
236 | #endif /* _PSTL_GLUE_NUMERIC_IMPL_H_ */
|
---|