[1166] | 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_ */
|
---|