reduce.h 6.36 KB
Newer Older
1
/*
Marcus Winter committed
2
 * Copyright (c) 2014-2016, Siemens AG. All rights reserved.
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
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions are met:
 *
 * 1. Redistributions of source code must retain the above copyright notice,
 * this list of conditions and the following disclaimer.
 *
 * 2. Redistributions in binary form must reproduce the above copyright notice,
 * this list of conditions and the following disclaimer in the documentation
 * and/or other materials provided with the distribution.
 *
 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
 * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
 * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
 * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 * POSSIBILITY OF SUCH DAMAGE.
 */

#ifndef EMBB_ALGORITHMS_REDUCE_H_
#define EMBB_ALGORITHMS_REDUCE_H_

30
#include <embb/tasks/execution_policy.h>
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
#include <embb/algorithms/identity.h>

namespace embb {
namespace algorithms {

/**
 * \defgroup CPP_ALGORITHMS_REDUCTION Reduction
 * Parallel reduction computation
 * \ingroup CPP_ALGORITHMS
 * \{
 */

#ifdef DOXYGEN

/**
 * Performs a parallel reduction operation on a range of elements.
 *
 * The range consists of the elements from \c first to \c last, excluding the
 * last element. The type of the result (\c ReturnType) is deduced from the
 * \c neutral element.
 *
 * \return
 * <tt>reduction(transformation(*first), ..., transformation(*(last-1)))</tt>
 * where the reduction function is applied pairwise.
 * \throws embb::base::ErrorException if not enough MTAPI tasks can be created
 *         to satisfy the requirements of the algorithm.
 * \threadsafe if the elements in the range are not modified by another thread
 *             while the algorithm is executed.
 * \note No guarantee is given on the order in which the functions \c reduction
 *       and \c transformation are applied to the elements.\n
 *       For all \c x of type \c ReturnType it must hold that
 *       <tt>reduction(x, neutral) == x</tt>. \n
 *       The reduction operation need not be commutative but must be
 *       associative, i.e., <tt>reduction(x, reduction(y, z)) ==
 *       reduction(reduction(x, y), z))</tt> for all \c x, \c y, \c z of type
 *       \c ReturnType.
67
 * \see embb::mtapi::ExecutionPolicy, ZipIterator, Identity
68 69 70
 * \tparam RAI Random access iterator
 * \tparam ReturnType Type of result of reduction operation, deduced from
 *         \c neutral
71
 * \tparam ReductionFunction Binary reduction function object with signature
72
 *         <tt>ReturnType ReductionFunction(ReturnType, ReturnType)</tt>.
73 74
 * \tparam TransformationFunction Unary transformation function object with
 *         signature <tt>ReturnType TransformationFunction(typename
75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91
 *         std::iterator_traits<RAI>::value_type)</tt>
 */
template<typename RAI, typename ReturnType, typename ReductionFunction,
         typename TransformationFunction>
ReturnType Reduce(
  RAI first,
  /**< [IN] Random access iterator pointing to the first element of the range */
  RAI last,
  /**< [IN] Random access iterator pointing to the last plus one element of the
            range */
  ReturnType neutral,
  /**< [IN] Neutral element of the reduction operation. */
  ReductionFunction reduction,
  /**< [IN] Reduction operation to be applied to the elements of the range */
  TransformationFunction transformation = Identity(),
  /**< [IN] Transforms the elements of the range before the reduction operation
            is applied */
92
  const embb::mtapi::ExecutionPolicy& policy = embb::mtapi::ExecutionPolicy(),
93
  /**< [IN] embb::mtapi::ExecutionPolicy for the reduction computation */
94 95 96 97 98 99 100 101 102 103 104 105
  size_t block_size = 0
  /**< [IN] Lower bound for partitioning the range of elements into blocks that
            are treated in parallel. Partitioning of a block stops if its size
            is less than or equal to \c block_size. The default value 0 means
            that the minimum block size is determined automatically depending on
            the number of elements in the range divided by the number of
            available cores. */
  );

#else // DOXYGEN

/**
106 107 108 109 110 111 112 113 114 115
 * Overload of above described Doxygen dummy.
 */
template<typename RAI, typename ReturnType, typename ReductionFunction,
         typename TransformationFunction>
ReturnType Reduce(
  RAI first,
  RAI last,
  ReturnType neutral,
  ReductionFunction reduction,
  TransformationFunction transformation,
116
  const embb::tasks::ExecutionPolicy& policy,
117 118 119 120
  size_t block_size
  );

/**
121 122 123 124 125 126 127 128 129
 * Overload of above described Doxygen dummy with less arguments.
 */
template<typename RAI, typename ReturnType, typename ReductionFunction>
ReturnType Reduce(
  RAI first,
  RAI last,
  ReturnType neutral,
  ReductionFunction reduction
  ) {
130
  return Reduce(first, last, neutral, reduction, Identity(),
131
    embb::tasks::ExecutionPolicy(), 0);
132 133 134 135 136 137 138 139 140 141 142 143 144 145 146
}

/**
 * Overload of above described Doxygen dummy with less arguments.
 */
template<typename RAI, typename ReturnType, typename ReductionFunction,
         typename TransformationFunction>
ReturnType Reduce(
  RAI first,
  RAI last,
  ReturnType neutral,
  ReductionFunction reduction,
  TransformationFunction transformation
  ) {
  return Reduce(first, last, neutral, reduction, transformation,
147
    embb::tasks::ExecutionPolicy(), 0);
148 149 150 151 152 153 154 155 156 157 158 159 160
}

/**
 * Overload of above described Doxygen dummy with less arguments.
 */
template<typename RAI, typename ReturnType, typename ReductionFunction,
         typename TransformationFunction>
ReturnType Reduce(
  RAI first,
  RAI last,
  ReturnType neutral,
  ReductionFunction reduction,
  TransformationFunction transformation,
161
  const embb::tasks::ExecutionPolicy& policy
162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177
  ) {
  return Reduce(first, last, neutral, reduction, transformation, policy, 0);
}

#endif // else DOXYGEN

/**
 * \}
 */

}  // namespace algorithms
}  // namespace embb

#include <embb/algorithms/internal/reduce-inl.h>

#endif  // EMBB_ALGORITHMS_REDUCE_H_