partial_sum.hpp 1.6 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546
  1. //---------------------------------------------------------------------------//
  2. // Copyright (c) 2013 Kyle Lutz <kyle.r.lutz@gmail.com>
  3. //
  4. // Distributed under the Boost Software License, Version 1.0
  5. // See accompanying file LICENSE_1_0.txt or copy at
  6. // http://www.boost.org/LICENSE_1_0.txt
  7. //
  8. // See http://boostorg.github.com/compute for more information.
  9. //---------------------------------------------------------------------------//
  10. #ifndef BOOST_COMPUTE_ALGORITHM_PARTIAL_SUM_HPP
  11. #define BOOST_COMPUTE_ALGORITHM_PARTIAL_SUM_HPP
  12. #include <boost/static_assert.hpp>
  13. #include <boost/compute/system.hpp>
  14. #include <boost/compute/command_queue.hpp>
  15. #include <boost/compute/algorithm/inclusive_scan.hpp>
  16. #include <boost/compute/type_traits/is_device_iterator.hpp>
  17. namespace boost {
  18. namespace compute {
  19. /// Calculates the cumulative sum of the elements in the range [\p first,
  20. /// \p last) and writes the resulting values to the range beginning at
  21. /// \p result.
  22. ///
  23. /// Space complexity on GPUs: \Omega(n)<br>
  24. /// Space complexity on GPUs when \p first == \p result: \Omega(2n)<br>
  25. /// Space complexity on CPUs: \Omega(1)
  26. template<class InputIterator, class OutputIterator>
  27. inline OutputIterator
  28. partial_sum(InputIterator first,
  29. InputIterator last,
  30. OutputIterator result,
  31. command_queue &queue = system::default_queue())
  32. {
  33. BOOST_STATIC_ASSERT(is_device_iterator<InputIterator>::value);
  34. BOOST_STATIC_ASSERT(is_device_iterator<OutputIterator>::value);
  35. return ::boost::compute::inclusive_scan(first, last, result, queue);
  36. }
  37. } // end compute namespace
  38. } // end boost namespace
  39. #endif // BOOST_COMPUTE_ALGORITHM_PARTIAL_SUM_HPP