glue_numeric_impl.h 12 KB

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