3607 lines
142 KiB
C++
3607 lines
142 KiB
C++
/* reducer_min_max.h -*- C++ -*-
|
||
*
|
||
* @copyright
|
||
* Copyright (C) 2009-2013, Intel Corporation
|
||
* All rights reserved.
|
||
*
|
||
* @copyright
|
||
* Redistribution and use in source and binary forms, with or without
|
||
* modification, are permitted provided that the following conditions
|
||
* are met:
|
||
*
|
||
* * Redistributions of source code must retain the above copyright
|
||
* notice, this list of conditions and the following disclaimer.
|
||
* * 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.
|
||
* * Neither the name of Intel Corporation nor the names of its
|
||
* contributors may be used to endorse or promote products derived
|
||
* from this software without specific prior written permission.
|
||
*
|
||
* @copyright
|
||
* 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.
|
||
*/
|
||
|
||
/** @file reducer_min_max.h
|
||
*
|
||
* @brief Defines classes for doing parallel minimum and maximum reductions.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*
|
||
* @see ReducersMinMax
|
||
*/
|
||
|
||
#ifndef REDUCER_MIN_MAX_H_INCLUDED
|
||
#define REDUCER_MIN_MAX_H_INCLUDED
|
||
|
||
#include <cilk/reducer.h>
|
||
|
||
#ifdef __cplusplus
|
||
|
||
#include <algorithm>
|
||
#include <limits>
|
||
|
||
/** @defgroup ReducersMinMax Minimum and Maximum Reducers
|
||
*
|
||
* Minimum and maximum reducers allow the computation of the minimum or
|
||
* maximum of a set of values in parallel.
|
||
*
|
||
* @ingroup Reducers
|
||
*
|
||
* You should be familiar with @ref pagereducers "Cilk reducers", described in
|
||
* file `reducers.md`, and particularly with @ref reducers_using, before trying
|
||
* to use the information in this file.
|
||
*
|
||
* @section redminmax_usage Usage Examples
|
||
*
|
||
* cilk::reducer< cilk::op_max<int> > rm;
|
||
* cilk_for (int i = 0; i < ARRAY_SIZE; ++i)
|
||
* {
|
||
* rm->calc_max(a[i]); // or *rm = cilk::max_of(*max, a[i])
|
||
* }
|
||
* std::cout << "maximum value is " << rm.get_value() << std::endl;
|
||
*
|
||
* and
|
||
*
|
||
* cilk::reducer< cilk::op_min_index<int, double> > rmi;
|
||
* cilk_for (int i = 0; i < ARRAY_SIZE; ++i)
|
||
* {
|
||
* rmi->calc_min(i, a[i]) // or *rmi = cilk::min_of(*rmi, i, a[i]);
|
||
* }
|
||
* std::cout << "minimum value a[" << rmi.get_value().first << "] = "
|
||
* << rmi.get_value().second << std::endl;
|
||
*
|
||
* @section redminmax_monoid The Monoid
|
||
*
|
||
* @subsection redminmax_monoid_values Value Set
|
||
*
|
||
* The value set of a minimum or maximum reducer is the set of values of
|
||
* `Type`, possibly augmented with a special identity value which is greater
|
||
* than (less than) any value of `Type`.
|
||
*
|
||
* @subsection redminmax_monoid_operator Operator
|
||
*
|
||
* In the most common case, the operator of a minimum reducer is defined as
|
||
*
|
||
* x MIN y == (x < y) ? x : y
|
||
*
|
||
* Thus, `a1 MIN a2 MIN … an` is the first `ai` which is not greater than any
|
||
* other `ai`.
|
||
*
|
||
* The operator of a maximum reducer is defined as
|
||
*
|
||
* x MAX y == (x > y) ? x : y
|
||
*
|
||
* Thus, `a1 MAX a2 MAX … an` is the first `ai` which is not less than any
|
||
* other `ai`.
|
||
*
|
||
* @subsection redminmax_monoid_comparators Comparators
|
||
*
|
||
* Min/max reducers are not limited to finding the minimum or maximum value
|
||
* determined by the `<` or `>` operator. In fact, all min/max reducers use a
|
||
* _comparator_, which is either a function or an object of a function class
|
||
* that defines a [strict weak ordering]
|
||
* (http://en.wikipedia.org/wiki/Strict_weak_ordering#Strict_weak_orderings)
|
||
* on a set of values. (This is exactly the same as the requirement for the
|
||
* comparison predicate for STL associative containers and sorting
|
||
* algorithms.)
|
||
*
|
||
* Just as with STL algorithms and containers, the comparator type parameter
|
||
* for min/max reducers is optional. If it is omitted, it defaults to
|
||
* `std::less`, which gives the behavior described in the previous section.
|
||
* Using non-default comparators (anything other than `std::less`) with
|
||
* min/max reducers is just like using them with STL containers and
|
||
* algorithms.
|
||
*
|
||
* Taking comparator objects into account, the reduction operation `MIN` for a
|
||
* minimum reducer is defined as
|
||
*
|
||
* x MIN y == compare(x, y) ? x : y
|
||
*
|
||
* where `compare()` is the reducer’s comparator. Similarly, the reduction
|
||
* operation MAX for a maximum reducer is defined as
|
||
*
|
||
* x MAX y == compare(y, x) ? x : y
|
||
*
|
||
* (If `compare(x, y) == x < y`, then `compare(y, x) == x > y`.)
|
||
*
|
||
* @subsection redminmax_monoid_identity Identity
|
||
*
|
||
* The identity value of the reducer is the value which is greater than (less
|
||
* than) any other value in the value set of the reducer. This is the
|
||
* [“special identity value”](#redminmax_monoid_values) if the reducer has
|
||
* one, or the largest (smallest) value in the value set otherwise.
|
||
*
|
||
* @section redminmax_index Value and Index Reducers
|
||
*
|
||
* Min/max reducers come in two families. The _value_ reducers, using `op_min`
|
||
* and `op_max` monoids, simply find the smallest or largest value from a set
|
||
* of values. The _index_ reducers, using `op_min_index` and `op_max_index`
|
||
* monoids, also record an index value associated with the first occurrence of
|
||
* the smallest or largest value.
|
||
*
|
||
* In the `%op_min_index` usage example [above](#redminmax_usage), the values
|
||
* are taken from an array, and the index of a value is the index of the array
|
||
* element it comes from. More generally, though, an index can be any sort of
|
||
* key which identifies a particular value in a collection of values. For
|
||
* example, if the values were taken from the nodes of a tree, then the
|
||
* “index” of a value might be a pointer to the node containing that value.
|
||
*
|
||
* A min/max index reducer is essentially the same as a min/max value reducer
|
||
* whose value type is an (index, value) pair, and whose comparator ignores
|
||
* the index part of the pair. (index, value) pairs are represented by
|
||
* `std::pair<Index, Type>` objects. This has the consequence that wherever
|
||
* the interface of a min/max value reducer has a `Type`, the interface of the
|
||
* corresponding min/max index reducer has a `std::pair<Index, Type>`. (There
|
||
* are convenience variants of the `reducer(Type)` constructor and the
|
||
* `calc_min()`, `calc_max()`, `%min_of()`, and `%max_of()` functions that
|
||
* take an index argument and a value argument instead of an index/value
|
||
* pair.)
|
||
*
|
||
* @section redminmax_operations Operations
|
||
*
|
||
* @subsection redminmax_constructors Constructors
|
||
*
|
||
* @subsubsection redminmax_constructors_value Min/Max Value Reducers
|
||
*
|
||
* reducer() // identity
|
||
* reducer(const Compare& compare) // identity
|
||
* reducer(const Type& value)
|
||
* reducer(move_in(Type& variable))
|
||
* reducer(const Type& value, const Compare& compare)
|
||
* reducer(move_in(Type& variable), const Compare& compare)
|
||
*
|
||
* @subsubsection redminmax_constructors_index Min/Max Index Reducers
|
||
*
|
||
* reducer() // identity
|
||
* reducer(const Compare& compare) // identity
|
||
* reducer(const std::pair<Index, Type>& pair)
|
||
* reducer(const Index& index, const Type& value)
|
||
* reducer(move_in(std::pair<Index, Type>& variable))
|
||
* reducer(const std::pair<Index, Type>& pair, const Compare& compare)
|
||
* reducer(const Index& index, const Type& value, const Compare& compare)
|
||
* reducer(move_in(std::pair<Index, Type>& variable), const Compare& compare)
|
||
*
|
||
* @subsection redminmax_get_set Set and Get
|
||
*
|
||
* r.set_value(const Type& value)
|
||
* Type = r.get_value() const
|
||
* r.move_in(Type& variable)
|
||
* r.move_out(Type& variable)
|
||
*
|
||
* Note that for an index reducer, the `Type` in these operations is actually a
|
||
* `std::pair<Index, Type>`. (See @ref redminmax_index.) There is _not_ a
|
||
* `set_value(value, index)` operation.
|
||
*
|
||
* @subsection redminmax_initial Initial Values and is_set()
|
||
*
|
||
* A minimum or maximum reducer without a specified initial value, before any
|
||
* MIN or MAX operation has been performed on it, represents the [identity
|
||
* value](#redminmax_monoid_identity) of its monoid. For value reducers with a
|
||
* numeric type and default comparator (`std::less`), this will be a well
|
||
* defined value. For example,
|
||
*
|
||
* reducer< op_max<unsigned> > r1;
|
||
* // r1.get_value() == 0
|
||
*
|
||
* reducer< op_min<float> > r2;
|
||
* // r2.get_value() == std::numeric_limits<float>::infinity
|
||
*
|
||
* In other cases, though (index reducers, non-numeric types, or non-default
|
||
* comparators), the actual identity value for the monoid may be unknown, or
|
||
* it may not even be a value of the reducer’s type. For example, there is no
|
||
* “largest string” to serve as the initial value for a
|
||
* `reducer< op_min<std::string> >`. In these cases, the result of calling
|
||
* `get_value()` is undefined.
|
||
*
|
||
* To avoid calling `get_value()` when its result is undefined, you can call
|
||
* the view’s `is_set()` function, which will return true if the reducer
|
||
* has a well-defined value — either because a MIN or MAX operation has been
|
||
* performed, or because it had a well-defined initial value:
|
||
*
|
||
* reducer< op_max<unsigned> > r1;
|
||
* // r1->is_set() == true
|
||
* // r1.get_value() == 0
|
||
*
|
||
* reducer< op_min<std::string> > r2;
|
||
* // r2->is_set() == false
|
||
* // r2.get_value() is undefined
|
||
* r2->calc_min("xyzzy");
|
||
* // r2->is_set() == true
|
||
* // r2.get_value() == "xyzzy"
|
||
*
|
||
* > Note: For an index reducer without a specified initial value, the
|
||
* > initial value of the index is the default value of the `Index` type.
|
||
*
|
||
* @subsection redminmax_view_ops View Operations
|
||
*
|
||
* The basic reduction operation is `x = x MIN a` for a minimum reducer, or
|
||
* `x = x MAX a` for a maximum reducer. The basic syntax for these operations
|
||
* uses the `calc_min()` and `calc_max()` member functions of the view class.
|
||
* An assignment syntax is also provided, using the %cilk::min_of() and
|
||
* %cilk::max_of() global functions:
|
||
*
|
||
* Class | Modifier | Assignment
|
||
* ---------------|---------------------|-----------
|
||
* `op_min` | `r->calc_min(x)` | `*r = min_of(*r, x)` or `*r = min_of(x, *r)`
|
||
* `op_max` | `r->calc_max(x)` | `*r = max_of(*r, x)` or `*r = max_of(x, *r)`
|
||
* `op_min_index` | `r->calc_min(i, x)` | `*r = min_of(*r, i, x)` or `*r = min_of(i, x, *r)`
|
||
* `op_max_index` | `r->calc_max(i, x)` | `*r = max_of(*r, i, x)` or `*r = max_of(i, x, *r)`
|
||
*
|
||
* Wherever an “`i`, `x`” argument pair is shown in the table above, a single
|
||
* pair argument may be passed instead. For example:
|
||
*
|
||
* Index index;
|
||
* Type value;
|
||
* std::pair<Index, Type> ind_val(index, value);
|
||
* // The following statements are all equivalent.
|
||
* r->calc_min(index, value);
|
||
* r->calc_min(ind_val);
|
||
* *r = min_of(*r, index, value);
|
||
* *r = min_of(*r, ind_val);
|
||
*
|
||
* The `calc_min()` and `calc_max()` member functions return a reference to
|
||
* the view, so they can be chained:
|
||
*
|
||
* r->calc_max(x).calc_max(y).calc_max(z);
|
||
*
|
||
* In a `%min_of()` or `%max_of()` assignment, the view on the left-hand side
|
||
* of the assignment must be the same as the view argument in the call.
|
||
* Otherwise, the behavior is undefined (but an assertion error will occur if
|
||
* the code is compiled with debugging enabled).
|
||
*
|
||
* *r = max_of(*r, x); // OK
|
||
* *r1 = max_of(*r2, y); // ERROR
|
||
*
|
||
* `%min_of()` and `%max_of()` calls can be nested:
|
||
*
|
||
* *r = max_of(max_of(max_of(*r, x), y), z);
|
||
* *r = min_of(i, a[i], min_of(j, a[j], min_of(k, a[k], *r)));
|
||
*
|
||
* @section redminmax_compatibility Compatibility Issues
|
||
*
|
||
* Most Cilk library reducers provide
|
||
* * Binary compatibility between `reducer_KIND` reducers compiled with Cilk
|
||
* library version 0.9 (distributed with Intel® C++ Composer XE version
|
||
* 13.0 and earlier) and the same reducers compiled with Cilk library
|
||
* version 1.0 and later.
|
||
* * Transparent casting between references to `reducer<op_KIND>` and
|
||
* `reducer_KIND`.
|
||
*
|
||
* This compatibility is not available in all cases for min/max reducers.
|
||
* There are two areas of incompatibility.
|
||
*
|
||
* @subsection redminmax_compatibility_stateful Non-empty Comparators
|
||
*
|
||
* There is no way to provide binary compatibility between the 0.9 and 1.0
|
||
* definitions of min/max reducers that use a non-empty comparator class or a
|
||
* comparator function. (Empty comparator classes like `std::less` are not a
|
||
* problem.)
|
||
*
|
||
* To avoid run-time surprises, the legacy `reducer_{min|max}[_index]` classes
|
||
* have been coded in the 1.0 library so that they will not even compile when
|
||
* instantiated with a non-empty comparator class.
|
||
*
|
||
* @subsection redminmax_compatibility_optimized Numeric Optimization
|
||
*
|
||
* Min/max reducers with a numeric value type and the default comparator can
|
||
* be implemented slightly more efficiently than other min/max reducers.
|
||
* However, the optimization is incompatible with the 0.9 library
|
||
* implementation of min/max reducers.
|
||
*
|
||
* The default min/max reducers implementation in the 1.0 library uses this
|
||
* numeric optimization. Code using legacy reducers compiled with the 1.0
|
||
* library can be safely used in the same program as code compiled with the
|
||
* 0.9 library, but classes compiled with the different Cilk libraries will be
|
||
* defined in different namespaces.
|
||
*
|
||
* The simplest solution is just to recompile the code that was compiled with
|
||
* the older version of Cilk. However, if this is impossible, you can define
|
||
* the `CILK_LIBRARY_0_9_REDUCER_MINMAX` macro (on the compiler command line,
|
||
* or in your source code before including `reducer_min_max.h`) when compiling
|
||
* with the new library. This will cause it to generate numeric reducers that
|
||
* will be less efficient, but will be fully compatible with previously
|
||
* compiled code. (Note that this macro has no effect on [the non-empty
|
||
* comparator incompatibility] (redminmax_compatibility_stateful).)
|
||
*
|
||
* @section redminmax_types Type Requirements
|
||
*
|
||
* `Type` and `Index` must be `Copy Constructible`, `Default Constructible`,
|
||
* and `Assignable`.
|
||
*
|
||
* `Compare` must be `Copy Constructible` if the reducer is constructed with a
|
||
* `compare` argument, and `Default Constructible` otherwise.
|
||
*
|
||
* The `Compare` function must induce a strict weak ordering on the elements
|
||
* of `Type`.
|
||
*
|
||
* @section redminmax_in_c Minimum and Maximum Reducers in C
|
||
*
|
||
* These macros can be used to do minimum and maximum reductions in C:
|
||
*
|
||
* Declaration | Type | Operation
|
||
* -----------------------------|-----------------------------------|----------
|
||
* @ref CILK_C_REDUCER_MIN |@ref CILK_C_REDUCER_MIN_TYPE |@ref CILK_C_REDUCER_MIN_CALC
|
||
* @ref CILK_C_REDUCER_MAX |@ref CILK_C_REDUCER_MAX_TYPE |@ref CILK_C_REDUCER_MAX_CALC
|
||
* @ref CILK_C_REDUCER_MIN_INDEX |@ref CILK_C_REDUCER_MIN_INDEX_TYPE |@ref CILK_C_REDUCER_MIN_INDEX_CALC
|
||
* @ref CILK_C_REDUCER_MAX_INDEX |@ref CILK_C_REDUCER_MAX_INDEX_TYPE |@ref CILK_C_REDUCER_MAX_INDEX_CALC
|
||
*
|
||
* For example:
|
||
*
|
||
* CILK_C_REDUCER_MIN(r, int, INT_MAX);
|
||
* CILK_C_REGISTER_REDUCER(r);
|
||
* cilk_for(int i = 0; i != n; ++i) {
|
||
* CILK_C_REDUCER_MIN_CALC(r, a[i]);
|
||
* }
|
||
* CILK_C_UNREGISTER_REDUCER(r);
|
||
* printf("The smallest value in a is %d\n", REDUCER_VIEW(r));
|
||
*
|
||
*
|
||
* CILK_C_REDUCER_MAX_INDEX(r, uint, 0);
|
||
* CILK_C_REGISTER_REDUCER(r);
|
||
* cilk_for(int i = 0; i != n; ++i) {
|
||
* CILK_C_REDUCER_MAX_INDEX_CALC(r, i, a[i]);
|
||
* }
|
||
* CILK_C_UNREGISTER_REDUCER(r);
|
||
* printf("The largest value in a is %u at %d\n",
|
||
* REDUCER_VIEW (r).value, REDUCER_VIEW(r).index);
|
||
*
|
||
* See @ref reducers_c_predefined.
|
||
*/
|
||
|
||
namespace cilk {
|
||
|
||
/** @defgroup ReducersMinMaxBinComp Binary compatibility
|
||
*
|
||
* If the macro CILK_LIBRARY_0_9_REDUCER_MINMAX is defined, then we generate
|
||
* reducer code and data structures which are binary-compatible with code that
|
||
* was compiled with the old min/max wrapper definitions, so we want the
|
||
* mangled names of the legacy min/max reducer wrapper classes to be the
|
||
* same as the names produced by the old definitions.
|
||
*
|
||
* Conversely, if the macro is not defined, then we generate binary-
|
||
* incompatible code, so we want different mangled names, to make sure that
|
||
* the linker does not allow new and old compiled legacy wrappers to be passed
|
||
* to one another. (Global variables are a different, and probably insoluble,
|
||
* problem.)
|
||
*
|
||
* Similarly, min/max classes compiled with and without
|
||
* CILK_LIBRARY_0_9_REDUCER_MINMAX are binary-incompatible, and must get
|
||
* different mangled names.
|
||
*
|
||
* The trick is, when compiling in normal (non-compatibility) mode, wrap
|
||
* everything in an extra namespace, and then `use` it into the top-level cilk
|
||
* namespace. Then
|
||
*
|
||
* * Classes and functions compiled in normal mode will be in
|
||
* different namespaces from the same classes and functions compiled in
|
||
* compatibility mode.
|
||
* * The legacy wrapper classes and functions will be in the same namespace
|
||
* as the same classes and functions compiled with the0.9 library if and
|
||
* only if the are compiled in compatibility mode.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
|
||
#ifndef CILK_LIBRARY_0_9_REDUCER_MINMAX
|
||
/** Namespace to wrap min/max reducer definitions when not compiling in “binary
|
||
* compatibility” mode.
|
||
*
|
||
* By default, all of the min/max reducer definitions are defined in this
|
||
* namespace and then imported into namespace ::cilk, so that they do not
|
||
* clash with the legacy definitions with the same names. However, if the
|
||
* macro `CILK_LIBRARY_0_9_REDUCER_MINMAX` is defined, then the min/max
|
||
* definitions go directly into namespace ::cilk, so that, for example,
|
||
* cilk::reducer_max defined with the 1.0 library is equivalent (to the
|
||
* linker) to cilk::reducer_max defined with the 0.9 library.
|
||
*
|
||
* @ingroup ReducersMinMaxBinComp
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
namespace cilk_lib_1_0 {
|
||
#endif
|
||
|
||
/** Namespace containing internal implementation classes and functions for
|
||
* min/max reducers.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
namespace min_max_internal {
|
||
|
||
using ::cilk::internal::binary_functor;
|
||
using ::cilk::internal::typed_indirect_binary_function;
|
||
using ::cilk::internal::class_is_empty;
|
||
|
||
/** @defgroup ReducersMinMaxIsSet The “is_set optimization”
|
||
*
|
||
* The obvious definition of the identity value for a max or min reducer is as
|
||
* the smallest (or largest) value of the value type. However, for an
|
||
* arbitrary comparator and/or an arbitrary value type, the largest / smallest
|
||
* value may not be known. It may not even be defined — what is the largest
|
||
* string?
|
||
*
|
||
* Therefore, min/max reducers represent their value internally as a pair
|
||
* `(value, is_set)`. When `is_set` is true, the pair represents the known
|
||
* value `value`; when `is_set` is false, the pair represents the identity
|
||
* value.
|
||
*
|
||
* This is an effective solution, but the most common use of min/max reducers
|
||
* is probably with numeric types and the default definition of minimum or
|
||
* maximum (using `std::less`), in which case there are well-defined, knowable
|
||
* smallest and largest values. Testing `is_set` for every comparison is then
|
||
* unnecessary and wasteful.
|
||
*
|
||
* The “is_set optimization” just means generating code that doesn’t use
|
||
* `is_set` when it isn’t needed. It is implemented using two metaprogramming
|
||
* classes:
|
||
*
|
||
* - do_is_set_optimization tests whether the optimization is applicable.
|
||
* - identity_value gets the appropriate identity value for a type.
|
||
*
|
||
* The is_set optimization is the reason that min/max reducers compiled with
|
||
* Cilk library 1.0 are binary-incompatible with the same reducers compiled
|
||
* with library 0.9, and therefore the optimization is suppressed when
|
||
* compiling in
|
||
* ReducersMinMaxBinComp "binary compatibility mode".
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
|
||
/** Test whether the ReducersMinMaxIsSet "is_set optimization" is
|
||
* applicable.
|
||
*
|
||
* The @ref do_is_set_optimization class is used to test whether the is_set
|
||
* optimization should be applied for a particular reducer. It is instantiated
|
||
* with a value type and a comparator, and defines a boolean constant,
|
||
* `value`. Then `%do_is_set_optimization<Type, Comp>::%value` can be used as
|
||
* a boolean template parameter to control the specialization of another
|
||
* class.
|
||
*
|
||
* In ReducersMinMaxBinComp "binary compatibility mode", when the
|
||
* `CILK_LIBRARY_0_9_REDUCER_MINMAX` macro is defined, `value` will always
|
||
* be false.
|
||
*
|
||
* @tparam Type The value type for the reducer.
|
||
* @tparam Compare The comparator type for the reducer.
|
||
*
|
||
* @result The `value` data member will be `true` if @a Type is a numeric
|
||
* type, @a Compare is `std::less<Type>`, and
|
||
* `CILK_LIBRARY_0_9_REDUCER_MINMAX` is not defined.
|
||
*
|
||
* @see ReducersMinMaxIsSet
|
||
* @see @ref view_content
|
||
*
|
||
* @ingroup ReducersMinMaxIsSet
|
||
*/
|
||
template < typename Type,
|
||
typename Compare >
|
||
struct do_is_set_optimization
|
||
{
|
||
/// `True` if the is_set optimization should be applied to min/max reducers
|
||
/// with this value type and comparator; `false` otherwise.
|
||
static const bool value = false;
|
||
};
|
||
|
||
#ifndef CILK_LIBRARY_0_9_REDUCER_MINMAX
|
||
/// @cond
|
||
template <typename Type>
|
||
struct do_is_set_optimization<Type, std::less<Type> >
|
||
{
|
||
/// True in the special case where optimization is possible.
|
||
static const bool value = std::numeric_limits<Type>::is_specialized;
|
||
};
|
||
/// @endcond
|
||
#endif
|
||
|
||
|
||
/** Get the identity value when using the ReducersMinMaxIsSet
|
||
* "is_set optimization".
|
||
*
|
||
* This class defines a function which assigns the appropriate identity value
|
||
* to a variable when the is_set optimization is applicable.
|
||
*
|
||
* @tparam Type The value type for the reducer.
|
||
* @tparam Compare The comparator type for the reducer.
|
||
* @tparam ForMax `true` to get the identity value for a max reducer (i.e.,
|
||
* the smallest value of @a Type), `false` to get the identity
|
||
* value for a min reducer (i.e., the largest value of
|
||
* @a Type).
|
||
*
|
||
* @result If @a Type and @a Compare qualify for the is_set optimization, the
|
||
* `set_identity()' function will set its argument variable to the
|
||
* smallest or largest value of @a Type, depending on @a ForMax.
|
||
* Otherwise, `set_identity()` will be a no-op.
|
||
*
|
||
* @see ReducersMinMaxIsSet
|
||
*
|
||
* @ingroup ReducersMinMaxIsSet
|
||
* @see @ref view_content
|
||
*/
|
||
template < typename Type,
|
||
typename Compare,
|
||
bool ForMax,
|
||
bool = std::numeric_limits<Type>::is_specialized,
|
||
bool = std::numeric_limits<Type>::has_infinity >
|
||
struct identity_value {
|
||
/// Assign the identity value to the reference parameter.
|
||
static void set_identity(Type&) {}
|
||
};
|
||
|
||
/// @cond
|
||
template <typename Type>
|
||
struct identity_value<Type, std::less<Type>, true, true, true> {
|
||
/// Floating max identity is negative infinity.
|
||
static void set_identity(Type& id)
|
||
{ id = -std::numeric_limits<Type>::infinity(); }
|
||
};
|
||
|
||
template <typename Type>
|
||
struct identity_value<Type, std::less<Type>, true, true, false> {
|
||
/// Integer max identity is minimum value of type.
|
||
static void set_identity(Type& id)
|
||
{ id = std::numeric_limits<Type>::min(); }
|
||
};
|
||
|
||
template <typename Type>
|
||
struct identity_value<Type, std::less<Type>, false, true, true> {
|
||
/// Floating min identity is positive infinity.
|
||
static void set_identity(Type& id)
|
||
{ id = std::numeric_limits<Type>::infinity(); }
|
||
};
|
||
|
||
template <typename Type>
|
||
struct identity_value<Type, std::less<Type>, false, true, false> {
|
||
/// Integer min identity is maximum value of type.
|
||
static void set_identity(Type& id)
|
||
{ id = std::numeric_limits<Type>::max(); }
|
||
};
|
||
|
||
/// @endcond
|
||
|
||
|
||
/** Adapter class to reverse the arguments of a predicate.
|
||
*
|
||
* Observe that:
|
||
*
|
||
* (x < y) == (y > x)
|
||
* max(x, y) == (x < y) ? y : x
|
||
* min(x, y) == (y < x) ? y : x == (x > y) ? y : x
|
||
*
|
||
* More generally, if `c` is a predicate defining a `Strict Weak Ordering`,
|
||
* and `c*(x, y) == c(y, x)`, then
|
||
*
|
||
* max(x, y, c) == c(x, y) ? y : x
|
||
* min(x, y, c) == c(y, x) ? y : x == c*(x, y) ? y : x == max(x, y, c*)
|
||
*
|
||
* For any predicate `C` with argument type `T`, the template class
|
||
* `%reverse_predicate<C, T>` defines a predicate which is identical to `C`,
|
||
* except that its arguments are reversed. Thus, for example, we could
|
||
* implement `%op_min_view<Type, Compare>` as
|
||
* `%op_max_view<Type, %reverse_predicate<Compare, Type> >`.
|
||
* (Actually, op_min_view and op_max_view are both implemented as subclasses
|
||
* of a common base class, view_base.)
|
||
*
|
||
* @note If `C` is an empty functor class, then `reverse_predicate(C)` will
|
||
* also be an empty functor class.
|
||
*
|
||
* @tparam Predicate The predicate whose arguments are to be reversed.
|
||
* @tparam Argument @a Predicate’s argument type.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
template <typename Predicate,
|
||
typename Argument = typename Predicate::first_argument_type>
|
||
class reverse_predicate : private binary_functor<Predicate>::type {
|
||
typedef typename binary_functor<Predicate>::type base;
|
||
public:
|
||
/// Default constructor
|
||
reverse_predicate() : base() {}
|
||
/// Constructor with predicate object
|
||
reverse_predicate(const Predicate& p) : base(p) {}
|
||
/// The reversed predicate operation
|
||
bool operator()(const Argument& x, const Argument& y) const
|
||
{ return base::operator()(y, x); }
|
||
};
|
||
|
||
|
||
/** Class to represent the comparator for a min/max view class.
|
||
*
|
||
* This class is intended to accomplish two objectives in the implementation
|
||
* of min/max views.
|
||
*
|
||
* 1. To minimize data bloat, when we have a reducer with a non-stateless
|
||
* comparator, we want to keep a single instance of the comparator object
|
||
* in the monoid, and just call it from the views.
|
||
* 2. In ReducersMinMaxBinComp "binary compatibility mode", views for
|
||
* reducers with a stateless comparator must have the same content as in
|
||
* Cilk library 0.9 — that is, they must contain only `value` and
|
||
* `is_set` data members.
|
||
*
|
||
* To achieve the first objective, we use the
|
||
* @ref internal::typed_indirect_binary_function class defined in
|
||
* metaprogramming.h to wrap a pointer to the actual comparator. If no
|
||
* pointer is needed because the actual comparator is stateless, the
|
||
* `typed_indirect_binary_function` class will be empty, too.
|
||
*
|
||
* To achieve the second objective, we make the
|
||
* `typed_indirect_binary_function` class a base class of the view rather than
|
||
* a data member, so the “empty base class” rule will ensure no that no
|
||
* additional space is allocated in the view unless it is needed.
|
||
*
|
||
* We could simply use typed_indirect_binary_function as the base class of the
|
||
* view, but this would mean writing comparisons as `(*this)(x, y)`, which is
|
||
* just weird. So, instead, we comparator_base as a subclass of
|
||
* typed_indirect_binary_function which provides function `compare()`
|
||
* as a synonym for `operator()`.
|
||
*
|
||
* @tparam Type The value type of the comparator class.
|
||
* @tparam Compare A predicate class.
|
||
*
|
||
* @see internal::typed_indirect_binary_function
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
class comparator_base : private typed_indirect_binary_function<Compare, Type, Type, bool>
|
||
{
|
||
typedef typed_indirect_binary_function<Compare, Type, Type, bool> base;
|
||
protected:
|
||
comparator_base(const Compare* f) : base(f) {} ///< Constructor.
|
||
|
||
/// Comparison function.
|
||
bool compare(const Type& a, const Type& b) const
|
||
{
|
||
return base::operator()(a, b);
|
||
}
|
||
|
||
/// Get the comparator pointer.
|
||
const Compare* compare_pointer() const { return base::pointer(); }
|
||
};
|
||
|
||
|
||
/** @defgroup ReducersMinMaxViewContent Content classes for min/max views
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*
|
||
* Minimum and maximum reducer view classes inherit from a “view content”
|
||
* class. The content class defines the actual data members for the view,
|
||
* and provides typedefs and member functions for accessing the data members
|
||
* as needed to support the view functionality.
|
||
*
|
||
* There are two content classes, which encapsulate the differences between
|
||
* simple min/max reducers and min/max with index reducers:
|
||
*
|
||
* - view_content
|
||
* - index_view_content
|
||
*
|
||
* @note An obvious, and arguably simpler, encapsulation strategy would be
|
||
* to just let the `Type` of a min/max view be an (index, value) pair
|
||
* structure for min_index and max_index reducers. Then all views
|
||
* would just have a `Type` data member and an `is_set` data member,
|
||
* and the comparator for min_index and max_index views could be
|
||
* customized to consider only the value component of the (index,
|
||
* value) `Type` pair. Unfortunately, this would break binary
|
||
* compatibility with reducer_max_index and reducer_min_index in
|
||
* Cilk library 0.9, because the memory layout of an (index, value)
|
||
* pair followed by a `bool` is different from the memory layout of an
|
||
* index data member followed by a value data member followed by a
|
||
* `bool` data member. The content class is designed to exactly
|
||
* replicate the layout of the views in library 0.9 reducers.
|
||
*
|
||
* A content class `C`, and its objects `c`, must define the following:
|
||
*
|
||
* Definition | Meaning
|
||
* ------------------------------------|--------
|
||
* `C::value_type` | A typedef for `Type` of the view. (A `std::pair<Index, Type>` for min_index and max_index views).
|
||
* `C::comp_value_type` | A typedef for the type of value compared by the view’s `compare()` function.
|
||
* `C()` | Constructs the content with the identity value.
|
||
* `C(const value_type&)` | Constructs the content with a specified value.
|
||
* `c.is_set()` | Returns true if the content has a known value.
|
||
* `c.value()` | Returns the content’s value.
|
||
* `c.set_value(const value_type&)` | Sets the content’s value. (The value becomes known.)
|
||
* `c.comp_value()` | Returns a const reference to the value or component of the value that is to be compared by the view’s comparator.
|
||
* `C::comp_value(const value_type&)` | Returns a const reference to a value or component of a value that is to be compared by the view’s comparator.
|
||
*
|
||
* @see view_base
|
||
*/
|
||
|
||
/** Content class for op_min_view and op_max_view.
|
||
*
|
||
* @tparam Type The value type of the op_min_view or op_max_view.
|
||
* @tparam Compare The comparator class specified for the op_min_view or
|
||
* op_max_view. (_Not_ the derived comparator class actually
|
||
* used by the view_base. For example, the view_content of an
|
||
* `op_min_view<int>` will have `Compare = std::less<int>`,
|
||
* but its comparator_base will have
|
||
* `Compare = reverse_predicate< std::less<int> >`.)
|
||
* @tparam ForMax `true` if this is the content class for an op_max_view,
|
||
* `false` if it is for an op_min_view.
|
||
*
|
||
* @note The general implementation of view_content uses an `is_set` data
|
||
* member. There is also a specialization which implements the
|
||
* ReducersMinMaxIsSet "is_set optimization". View classes that
|
||
* inherit from view_content do not need to know anything about the
|
||
* difference, though; the details are abstracted away in the
|
||
* view_content interface.
|
||
*
|
||
* @see ReducersMinMaxViewContent
|
||
*
|
||
* @ingroup ReducersMinMaxViewContent
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
template < typename Type
|
||
, typename Compare
|
||
, bool ForMax
|
||
, bool = do_is_set_optimization<Type, Compare>::value
|
||
>
|
||
class view_content {
|
||
Type m_value;
|
||
bool m_is_set;
|
||
public:
|
||
/// The value type of the view.
|
||
typedef Type value_type;
|
||
|
||
/// The type compared by the view’s `compare()` function (which is the same
|
||
/// as the value type for view_content).
|
||
typedef Type comp_value_type;
|
||
|
||
/// Construct with the identity value.
|
||
view_content() : m_value(), m_is_set(false) {}
|
||
|
||
/// Construct with a defined value.
|
||
view_content(const value_type& value) : m_value(value), m_is_set(true) {}
|
||
|
||
/// Get the value.
|
||
value_type value() const { return m_value; }
|
||
|
||
/// Set the value.
|
||
void set_value(const value_type& value)
|
||
{
|
||
m_value = value;
|
||
m_is_set = true;
|
||
}
|
||
|
||
/// Get the comparison value (which is the same as the value for
|
||
/// view_content).
|
||
const comp_value_type& comp_value() const { return m_value; }
|
||
|
||
/// Given an arbitrary value, get the corresponding comparison value (which
|
||
/// is the same as the value for view_content).
|
||
static const comp_value_type& comp_value(const value_type& value)
|
||
{
|
||
return value;
|
||
}
|
||
|
||
/// Get a const reference to value part of the value (which is the same as
|
||
/// the value for view_content).
|
||
const Type& get_reference() const { return m_value; }
|
||
|
||
/// Get a const reference to the index part of the value (which is
|
||
/// meaningless for non-index reducers, but required for view_base.
|
||
const Type& get_index_reference() const { return m_value; }
|
||
|
||
/// Test if the value is defined.
|
||
bool is_set() const { return m_is_set; }
|
||
};
|
||
|
||
/// @cond
|
||
|
||
/* This is the specialization of the view_content class for cases where
|
||
* `AssumeIsSet` is true (i.e., where the is_set optimization is applicable).
|
||
*/
|
||
template < typename Type
|
||
, typename Compare
|
||
, bool ForMax
|
||
>
|
||
class view_content<Type, Compare, ForMax, true> {
|
||
typedef identity_value<Type, Compare, ForMax> Identity;
|
||
Type m_value;
|
||
public:
|
||
typedef Type value_type;
|
||
typedef Type comp_value_type;
|
||
|
||
/// Construct with identity value.
|
||
view_content() { Identity::set_identity(m_value); }
|
||
|
||
view_content(const value_type& value) : m_value(value) {}
|
||
|
||
value_type value() const { return m_value; }
|
||
|
||
void set_value(const value_type& value)
|
||
{
|
||
m_value = value;
|
||
}
|
||
|
||
const comp_value_type& comp_value() const { return m_value; }
|
||
|
||
static const comp_value_type& comp_value(const value_type& value)
|
||
{
|
||
return value;
|
||
}
|
||
|
||
const Type& get_reference() const { return m_value; }
|
||
|
||
const Type& get_index_reference() const { return m_value; }
|
||
|
||
/// Test if the value is defined.
|
||
bool is_set() const { return true; }
|
||
};
|
||
|
||
/// @endcond
|
||
|
||
|
||
/** Content class for op_min_index_view and op_max_index_view.
|
||
*
|
||
* @tparam Index The index type of the op_min_index_view or
|
||
op_max_index_view.
|
||
* @tparam Type The value type of the op_min_view or op_max_view. (_Not_
|
||
* the value type of the view, which will be
|
||
* `std::pair<Index, Type>`.)
|
||
* @tparam Compare The comparator class specified for the op_min_index_view or
|
||
* op_max_index_view. (_Not_ the derived comparator class
|
||
* actually used by the view_base. For example, the
|
||
* index_view_content of an `op_min_index_view<int>` will have
|
||
* `Compare = std::less<int>`, but its comparator_base will
|
||
* have `Compare = reverse_predicate< std::less<int> >`.)
|
||
* @tparam ForMax `true` if this is the content class for an
|
||
* op_max_index_view, `false` if it is for an
|
||
* op_min_index_view.
|
||
*
|
||
* @see ReducersMinMaxViewContent
|
||
*
|
||
* @ingroup ReducersMinMaxViewContent
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
template < typename Index
|
||
, typename Type
|
||
, typename Compare
|
||
, bool ForMax
|
||
>
|
||
class index_view_content {
|
||
typedef identity_value<Type, Compare, ForMax> Identity;
|
||
|
||
Index m_index;
|
||
Type m_value;
|
||
bool m_is_set;
|
||
public:
|
||
/// The value type of the view (which is an <index, value> pair for
|
||
/// index_view_content).
|
||
typedef std::pair<Index, Type> value_type;
|
||
|
||
/// The type compared by the view’s `compare()` function (which is the data
|
||
/// value type for index_view_content).
|
||
typedef Type comp_value_type;
|
||
|
||
/// Construct with the identity value.
|
||
index_view_content() : m_index(), m_value(), m_is_set(false) {}
|
||
|
||
/// Construct with an index/value pair.
|
||
index_view_content(const value_type& value) :
|
||
m_index(value.first), m_value(value.second), m_is_set(true) {}
|
||
|
||
/// Construct with an index and a value.
|
||
index_view_content(const Index& index, const Type& value) :
|
||
m_index(index), m_value(value), m_is_set(true) {}
|
||
|
||
/// Construct with just an index.
|
||
index_view_content(const Index& index) :
|
||
m_index(index), m_value(), m_is_set(false) {}
|
||
|
||
/// Get the value.
|
||
value_type value() const { return value_type(m_index, m_value); }
|
||
|
||
/// Set value.
|
||
void set_value(const value_type& value)
|
||
{
|
||
m_index = value.first;
|
||
m_value = value.second;
|
||
m_is_set = true;
|
||
}
|
||
|
||
/// Get the comparison value (which is the value component of the
|
||
/// index/value pair for index_view_content).
|
||
const comp_value_type& comp_value() const { return m_value; }
|
||
|
||
/// Given an arbitrary value (i.e., index/value pair), get the
|
||
/// corresponding comparison value (which is the value component of the
|
||
/// index/value pair for index_view_content).
|
||
static const comp_value_type& comp_value(const value_type& value)
|
||
{ return value.second; }
|
||
|
||
/// Get a const reference to value part of the value.
|
||
const Type& get_reference() const { return m_value; }
|
||
|
||
/// Get a const reference to the index part of the value.
|
||
const Index& get_index_reference() const { return m_index; }
|
||
|
||
/// Test if the value is defined.
|
||
bool is_set() const { return m_is_set; }
|
||
};
|
||
|
||
|
||
template <typename View> class rhs_proxy;
|
||
|
||
/** Create an rhs_proxy.
|
||
*/
|
||
template <typename View>
|
||
inline rhs_proxy<View>
|
||
make_proxy(const typename View::value_type& value, const View& view);
|
||
|
||
template <typename Content, typename Less, typename Compare> class view_base;
|
||
|
||
|
||
/** Class to represent the right-hand side of
|
||
* `*reducer = {min|max}_of(*reducer, value)`.
|
||
*
|
||
* The only assignment operator for a min/max view class takes a rhs_proxy as
|
||
* its operand. This results in the syntactic restriction that the only
|
||
* expressions that can be assigned to a min/max view are ones which generate
|
||
* an rhs_proxy — that is, expressions of the form `max_of(view, value)` and
|
||
* `min_of(view, value)`.
|
||
*
|
||
* @warning
|
||
* The lhs and rhs views in such an assignment must be the same; otherwise,
|
||
* the behavior will be undefined. (I.e., `*r1 = min_of(*r1, x)` is legal;
|
||
* `*r1 = min_of(*r2, x)` is illegal.) This condition will be checked with a
|
||
* runtime assertion when compiled in debug mode.
|
||
*
|
||
* @tparam View The view class (op_{min|max}[_index]_view) that this proxy
|
||
* was created from.
|
||
*
|
||
* @see view_base
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
template <typename View>
|
||
class rhs_proxy {
|
||
typedef typename View::less_type less_type;
|
||
typedef typename View::compare_type compare_type;
|
||
typedef typename View::value_type value_type;
|
||
typedef typename View::content_type content_type;
|
||
typedef typename content_type::comp_value_type comp_value_type;
|
||
|
||
friend class view_base<content_type, less_type, compare_type>;
|
||
friend rhs_proxy make_proxy<View>(
|
||
const typename View::value_type& value,
|
||
const View& view);
|
||
|
||
typed_indirect_binary_function<
|
||
compare_type, comp_value_type, comp_value_type, bool>
|
||
m_comp;
|
||
const View* m_view;
|
||
value_type m_value;
|
||
|
||
rhs_proxy& operator=(const rhs_proxy&); // Disable assignment operator
|
||
rhs_proxy(); // Disable default constructor
|
||
|
||
// Constructor (called from view_base::make_proxy).
|
||
rhs_proxy(const View* view,
|
||
const value_type& value,
|
||
const compare_type* compare) :
|
||
m_view(view), m_value(value), m_comp(compare) {}
|
||
|
||
// Check matching view, then return value (called from view_base::assign).
|
||
value_type value(const typename View::base* view) const
|
||
{
|
||
__CILKRTS_ASSERT(view == m_view);
|
||
return m_value;
|
||
}
|
||
|
||
public:
|
||
|
||
/** Support max_of(max_of(view, value), value) and the like.
|
||
*/
|
||
rhs_proxy calc(const value_type& x) const
|
||
{
|
||
return rhs_proxy(
|
||
m_view,
|
||
m_comp( content_type::comp_value(m_value),
|
||
content_type::comp_value(x)
|
||
) ? x : m_value,
|
||
m_comp.pointer());
|
||
}
|
||
};
|
||
|
||
|
||
template <typename View>
|
||
inline rhs_proxy<View>
|
||
make_proxy(const typename View::value_type& value, const View& view)
|
||
{
|
||
return rhs_proxy<View>(&view, value, view.compare_pointer());
|
||
}
|
||
|
||
//@}
|
||
|
||
/** Base class for min and max view classes.
|
||
*
|
||
* This class accumulates the minimum or maximum of a set of values which have
|
||
* occurred as arguments to the `calc()` function, as determined by a
|
||
* comparator. The accumulated value will be the first `calc()` argument value
|
||
* `x` such that `compare(x, y)` is false for every `calc()` argument value
|
||
* `y`.
|
||
*
|
||
* If the comparator is `std::less`, then the accumulated value is the first
|
||
* argument value which is not less than any other argument value, i.e., the
|
||
* maximum. Similarly, if the comparator is `reverse_predicate<std::less>`,
|
||
* which is equivalent to `std::greater`, then the accumulated value is the
|
||
* first argument value which is not greater than any other argument value,
|
||
* i.e., the minimum.
|
||
*
|
||
* @note This class provides the definitions that are required for a class
|
||
* that will be used as the parameter of a
|
||
* min_max_internal::monoid_base specialization.
|
||
*
|
||
* @tparam Content A content class that provides the value types and data
|
||
* members for the view.
|
||
* @tparam Less A “less than” binary predicate that defines the min or
|
||
* max function.
|
||
* @tparam Compare A binary predicate to be used to compare the values.
|
||
* (The same as @a Less for max reducers; its reversal for
|
||
* min reducers.)
|
||
*
|
||
* @see ReducersMinMaxViewContent
|
||
* @see op_max_view
|
||
* @see op_min_view
|
||
* @see op_max_index_view
|
||
* @see op_min_index_view
|
||
* @see monoid_base
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
template <typename Content, typename Less, typename Compare>
|
||
class view_base :
|
||
// comparator_base comes first to ensure that it will get empty base class
|
||
// treatment
|
||
private comparator_base<typename Content::comp_value_type, Compare>,
|
||
private Content
|
||
{
|
||
typedef comparator_base<typename Content::comp_value_type, Compare> base;
|
||
using base::compare;
|
||
using Content::value;
|
||
using Content::set_value;
|
||
using Content::comp_value;
|
||
typedef Content content_type;
|
||
|
||
template <typename View> friend class rhs_proxy;
|
||
template <typename View>
|
||
friend rhs_proxy<View> make_proxy(const typename View::value_type& value, const View& view);
|
||
|
||
public:
|
||
|
||
/** @name Monoid support.
|
||
*/
|
||
//@{
|
||
|
||
/** Value type. Required by @ref monoid_with_view.
|
||
*/
|
||
typedef typename Content::value_type value_type;
|
||
|
||
/** The type of the comparator specified by the user, that defines the
|
||
* ordering on @a Type. Required by min_max::monoid_base.
|
||
*/
|
||
typedef Less less_type;
|
||
|
||
/** The type of the comparator actually used by the view. Required by
|
||
* min_max::monoid_base. (This is the same as the @ref less_type for a
|
||
* max reducer, or `reverse_predicate<less_type>` for a min reducer.)
|
||
*/
|
||
typedef Compare compare_type;
|
||
|
||
/** Reduce operation. Required by @ref monoid_with_view.
|
||
*/
|
||
void reduce(view_base* other)
|
||
{
|
||
if ( other->is_set() &&
|
||
( !this->is_set() ||
|
||
compare(this->comp_value(), other->comp_value()) ) )
|
||
{
|
||
this->set_value(other->value());
|
||
}
|
||
}
|
||
|
||
//@}
|
||
|
||
/** Default constructor. Initializes to identity value.
|
||
*/
|
||
explicit view_base(const compare_type* compare) :
|
||
base(compare), Content() {}
|
||
|
||
/** Value constructor.
|
||
*/
|
||
template <typename T1>
|
||
view_base(const T1& x1, const compare_type* compare) :
|
||
base(compare), Content(x1) {}
|
||
|
||
/** Value constructor.
|
||
*/
|
||
template <typename T1, typename T2>
|
||
view_base(const T1& x1, const T2& x2, const compare_type* compare) :
|
||
base(compare), Content(x1, x2) {}
|
||
|
||
|
||
/** Move-in constructor.
|
||
*/
|
||
explicit view_base(move_in_wrapper<value_type> w, const compare_type* compare) :
|
||
base(compare), Content(w.value()) {}
|
||
|
||
/** @name Reducer support.
|
||
*/
|
||
//@{
|
||
|
||
void view_move_in(value_type& v) { set_value(v); }
|
||
void view_move_out(value_type& v) { v = value(); }
|
||
void view_set_value(const value_type& v) { set_value(v); }
|
||
value_type view_get_value() const { return value(); }
|
||
// view_get_reference() NOT SUPPORTED
|
||
|
||
//@}
|
||
|
||
/** Is the value defined?
|
||
*/
|
||
using Content::is_set;
|
||
|
||
/** Reference to contained value data member.
|
||
* @deprecated For legacy reducers only.
|
||
*/
|
||
using Content::get_reference;
|
||
|
||
/** Reference to contained index data member.
|
||
* (Meaningless for non-index reducers.)
|
||
* @deprecated For legacy reducers only.
|
||
*/
|
||
using Content::get_index_reference;
|
||
|
||
protected:
|
||
|
||
/** Update the min/max value.
|
||
*/
|
||
void calc(const value_type& x)
|
||
{
|
||
if (!is_set() || compare(comp_value(), comp_value(x))) set_value(x);
|
||
}
|
||
|
||
/** Assign the result of a `{min|max}_of(view, value)` expression to the
|
||
* view.
|
||
*
|
||
* @see rhs_proxy
|
||
*/
|
||
template <typename View>
|
||
void assign(const rhs_proxy<View>& rhs)
|
||
{
|
||
calc(rhs.value(this));
|
||
}
|
||
|
||
};
|
||
|
||
|
||
/** Base class for min and max monoid classes.
|
||
*
|
||
* The unique characteristic of minimum and maximum reducers is that they
|
||
* incorporate a comparator functor that defines what “minimum” or “maximum”
|
||
* means. The monoid for a reducer contains the comparator that will be used
|
||
* for the reduction. If the comparator is a function or a class with state,
|
||
* then each view will have a pointer to the comparator.
|
||
*
|
||
* This means that the `construct()` functions first construct the monoid
|
||
* (possibly with an explicit comparator argument), and then construct the
|
||
* view with a pointer to the monoid’s comparator.
|
||
*
|
||
* @tparam View The view class.
|
||
* @tparam Align If true, reducers instantiated on this monoid will be
|
||
* aligned. By default, library reducers (unlike legacy
|
||
* library reducer _wrappers_) are unaligned.
|
||
*
|
||
* @see view_base
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
template <typename View, bool Align = false>
|
||
class monoid_base : public monoid_with_view<View, Align>
|
||
{
|
||
typedef typename View::compare_type compare_type;
|
||
typedef typename View::less_type less_type;
|
||
const compare_type m_compare;
|
||
|
||
const compare_type* compare_pointer() const { return &m_compare; }
|
||
|
||
using cilk::monoid_base<typename View::value_type, View>::provisional;
|
||
|
||
public:
|
||
|
||
/** Default constructor uses default comparator.
|
||
*/
|
||
monoid_base() : m_compare() {}
|
||
|
||
/** Constructor.
|
||
*
|
||
* @param compare The comparator to use.
|
||
*/
|
||
monoid_base(const compare_type& compare) : m_compare(compare) {}
|
||
|
||
/** Create an identity view.
|
||
*
|
||
* List view identity constructors take the list allocator as an argument.
|
||
*
|
||
* @param v The address of the uninitialized memory in which the view
|
||
* will be constructed.
|
||
*/
|
||
void identity(View *v) const { ::new((void*) v) View(compare_pointer()); }
|
||
|
||
/** @name construct functions
|
||
*
|
||
* Min/max monoid `construct()` functions optionally take one or two value
|
||
* arguments, a @ref move_in argument, and/or a comparator argument.
|
||
*/
|
||
//@{
|
||
|
||
template <typename Monoid>
|
||
static void construct(Monoid* monoid, View* view)
|
||
{ provisional( new ((void*)monoid) Monoid() ).confirm_if(
|
||
new ((void*)view) View(monoid->compare_pointer()) ); }
|
||
|
||
template <typename Monoid, typename T1>
|
||
static void construct(Monoid* monoid, View* view, const T1& x1)
|
||
{ provisional( new ((void*)monoid) Monoid() ).confirm_if(
|
||
new ((void*)view) View(x1, monoid->compare_pointer()) ); }
|
||
|
||
template <typename Monoid, typename T1, typename T2>
|
||
static void construct(Monoid* monoid, View* view, const T1& x1, const T2& x2)
|
||
{ provisional( new ((void*)monoid) Monoid() ).confirm_if(
|
||
new ((void*)view) View(x1, x2, monoid->compare_pointer()) ); }
|
||
|
||
template <typename Monoid>
|
||
static void construct(Monoid* monoid, View* view, const less_type& compare)
|
||
{ provisional( new ((void*)monoid) Monoid(compare) ).confirm_if(
|
||
new ((void*)view) View(monoid->compare_pointer()) ); }
|
||
|
||
template <typename Monoid, typename T1>
|
||
static void construct(Monoid* monoid, View* view, const T1& x1, const less_type& compare)
|
||
{ provisional( new ((void*)monoid) Monoid(compare) ).confirm_if(
|
||
new ((void*)view) View(x1, monoid->compare_pointer()) ); }
|
||
|
||
template <typename Monoid, typename T1, typename T2>
|
||
static void construct(Monoid* monoid, View* view, const T1& x1, const T2& x2, const less_type& compare)
|
||
{ provisional( new ((void*)monoid) Monoid(compare) ).confirm_if(
|
||
new ((void*)view) View(x1, x2, monoid->compare_pointer()) ); }
|
||
|
||
//@}
|
||
};
|
||
|
||
} //namespace min_max_internal
|
||
|
||
|
||
/** @defgroup ReducersMinMaxMaxValue Maximum reducers (value only)
|
||
*
|
||
* These reducers will find the largest value from a set of values.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
//@{
|
||
|
||
/** The maximum reducer view class.
|
||
*
|
||
* This is the view class for reducers created with
|
||
* `cilk::reducer< cilk::op_max<Type, Compare> >`. It accumulates the maximum,
|
||
* as determined by a comparator, of a set of values which have occurred as
|
||
* arguments to the `calc_max()` function. The accumulated value will be the
|
||
* first argument `x` such that `compare(x, y)` is false for every argument
|
||
* `y`.
|
||
*
|
||
* If the comparator is `std::less`, then the accumulated value is the first
|
||
* argument value which is not less than any other argument value, i.e., the
|
||
* maximum.
|
||
*
|
||
* @note The reducer “dereference” operation (`reducer::operator *()`)
|
||
* yields a reference to the view. Thus, for example, the view class’s
|
||
* `calc_max()` function would be used in an expression like
|
||
* `r->calc_max(a)` where `r` is an op_max reducer variable.
|
||
*
|
||
* @tparam Type The type of the values compared by the reducer. This will
|
||
* be the value type of a monoid_with_view that is
|
||
* instantiated with this view.
|
||
* @tparam Compare A `Strict Weak Ordering` whose argument type is @a Type. It
|
||
* defines the “less than” relation used to compute the
|
||
* maximum.
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_max
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
class op_max_view : public min_max_internal::view_base<
|
||
min_max_internal::view_content<Type, Compare, true>,
|
||
Compare,
|
||
Compare>
|
||
{
|
||
typedef min_max_internal::view_base<
|
||
min_max_internal::view_content<Type, Compare, true>,
|
||
Compare,
|
||
Compare> base;
|
||
using base::calc;
|
||
using base::assign;
|
||
friend class min_max_internal::rhs_proxy<op_max_view>;
|
||
|
||
public:
|
||
|
||
/** @name Constructors.
|
||
*
|
||
* All op_max_view constructors simply pass their arguments on to the
|
||
* @ref view_base base class.
|
||
*/
|
||
//@{
|
||
|
||
op_max_view() : base() {}
|
||
|
||
template <typename T1>
|
||
op_max_view(const T1& x1) : base(x1) {}
|
||
|
||
template <typename T1, typename T2>
|
||
op_max_view(const T1& x1, const T2& x2) : base(x1, x2) {}
|
||
|
||
//@}
|
||
|
||
/** @name View modifier operations.
|
||
*/
|
||
//@{
|
||
|
||
/** Maximize with a value.
|
||
*
|
||
* If @a x is greater than the current value of the view (as defined by
|
||
* the reducer’s comparator), or if the view was created without an
|
||
* initial value and its value has never been updated (with `calc_max()`
|
||
* or `= max_of()`), then the value of the view is set to @a x.
|
||
*
|
||
* @param x The value to maximize the view’s value with.
|
||
*
|
||
* @return A reference to the view. (Allows chaining
|
||
* `view.comp_max(a).comp_max(b)…`.)
|
||
*/
|
||
op_max_view& calc_max(const Type& x) { calc(x); return *this; }
|
||
|
||
/** Assign the result of a `max_of(view, value)` expression to the view.
|
||
*
|
||
* @param rhs An rhs_proxy value created by a `max_of(view, value)`
|
||
* expression.
|
||
*
|
||
* @return A reference to the view.
|
||
*
|
||
* @see min_max_internal::view_base::rhs_proxy
|
||
*/
|
||
op_max_view& operator=(const min_max_internal::rhs_proxy<op_max_view>& rhs)
|
||
{ assign(rhs); return *this; }
|
||
|
||
//@}
|
||
};
|
||
|
||
|
||
/** Compute the maximum of the value in an op_max_view and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or used in another max_of() call. For example,
|
||
*
|
||
* *reducer = max_of(*reducer, x);
|
||
* *reducer = max_of(x, *reducer);
|
||
*
|
||
* @see min_max_internal::rhs_proxy
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_view<Type, Compare> >
|
||
max_of(const op_max_view<Type, Compare>& view, const Type& value)
|
||
{
|
||
return min_max_internal::make_proxy(value, view);
|
||
}
|
||
|
||
/// @copydoc max_of(const op_max_view<Type, Compare>&, const Type&)
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_view<Type, Compare> >
|
||
max_of(const Type& value, const op_max_view<Type, Compare>& view)
|
||
{
|
||
return min_max_internal::make_proxy(value, view);
|
||
}
|
||
|
||
/** Nested maximum computation.
|
||
*
|
||
* Compute the maximum of the result of a max_of() call and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or wrapper, or used in another max_of() call. For example,
|
||
*
|
||
* *reducer = max_of(x, max_of(y, *reducer));
|
||
* wrapper = max_of(max_of(wrapper, x), y);
|
||
*
|
||
* @see min_max_internal::rhs_proxy
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_view<Type, Compare> >
|
||
max_of(const min_max_internal::rhs_proxy< op_max_view<Type, Compare> >& proxy,
|
||
const Type& value)
|
||
{
|
||
return proxy.calc(value);
|
||
}
|
||
|
||
/// @copydoc max_of(const min_max_internal::rhs_proxy< op_max_view<Type, Compare> >&, const Type&)
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_view<Type, Compare> >
|
||
max_of(const Type& value,
|
||
const min_max_internal::rhs_proxy< op_max_view<Type, Compare> >& proxy)
|
||
{
|
||
return proxy.calc(value);
|
||
}
|
||
|
||
|
||
/** Monoid class for maximum reductions. Instantiate the cilk::reducer template
|
||
* class with an op_max monoid to create a maximum reducer class. For example,
|
||
* to compute the maximum of a set of `int` values:
|
||
*
|
||
* cilk::reducer< cilk::op_max<int> > r;
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_max_view
|
||
*/
|
||
template <typename Type, typename Compare=std::less<Type>, bool Align = false>
|
||
class op_max :
|
||
public min_max_internal::monoid_base<op_max_view<Type, Compare>, Align>
|
||
{
|
||
typedef min_max_internal::monoid_base<op_max_view<Type, Compare>, Align>
|
||
base;
|
||
public:
|
||
/// Construct with default comparator.
|
||
op_max() {}
|
||
/// Construct with specified comparator.
|
||
op_max(const Compare& compare) : base(compare) {}
|
||
};
|
||
|
||
//@}
|
||
|
||
|
||
/** @defgroup ReducersMinMaxMinValue Minimum reducers (value only)
|
||
*
|
||
* These reducers will find the smallest value from a set of values.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
//@{
|
||
|
||
/** The minimum reducer view class.
|
||
*
|
||
* This is the view class for reducers created with
|
||
* `cilk::reducer< cilk::op_min<Type, Compare> >`. It accumulates the minimum,
|
||
* as determined by a comparator, of a set of values which have occurred as
|
||
* arguments to the `calc_min()` function. The accumulated value will be the
|
||
* first argument `x` such that `compare(y, x)` is false for every argument
|
||
* `y`.
|
||
*
|
||
* If the comparator is `std::less`, then the accumulated value is the first
|
||
* argument value which no other argument value is less than, i.e., the
|
||
* minimum.
|
||
*
|
||
* @note The reducer “dereference” operation (`reducer::operator *()`)
|
||
* yields a reference to the view. Thus, for example, the view class’s
|
||
* `calc_min()` function would be used in an expression like
|
||
* `r->calc_min(a)` where `r` is an op_min reducer variable.
|
||
*
|
||
* @tparam Type The type of the values compared by the reducer. This will
|
||
* be the value type of a monoid_with_view that is
|
||
* instantiated with this view.
|
||
* @tparam Compare A `Strict Weak Ordering` whose argument type is @a Type. It
|
||
* defines the “less than” relation used to compute the
|
||
* minimum.
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_min
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
class op_min_view : public min_max_internal::view_base<
|
||
min_max_internal::view_content<Type, Compare, false>,
|
||
Compare,
|
||
min_max_internal::reverse_predicate<Compare, Type> >
|
||
{
|
||
typedef min_max_internal::view_base<
|
||
min_max_internal::view_content<Type, Compare, false>,
|
||
Compare,
|
||
min_max_internal::reverse_predicate<Compare, Type> > base;
|
||
using base::calc;
|
||
using base::assign;
|
||
friend class min_max_internal::rhs_proxy<op_min_view>;
|
||
|
||
public:
|
||
/** @name Constructors.
|
||
*
|
||
* All op_min_view constructors simply pass their arguments on to the
|
||
* @ref view_base base class.
|
||
*/
|
||
//@{
|
||
|
||
op_min_view() : base() {}
|
||
|
||
template <typename T1>
|
||
op_min_view(const T1& x1) : base(x1) {}
|
||
|
||
template <typename T1, typename T2>
|
||
op_min_view(const T1& x1, const T2& x2) : base(x1, x2) {}
|
||
|
||
//@}
|
||
|
||
/** @name View modifier operations.
|
||
*/
|
||
//@{
|
||
|
||
/** Minimize with a value.
|
||
*
|
||
* If @a x is less than the current value of the view (as defined by the
|
||
* reducer’s comparator), or if the view was created without an initial
|
||
* value and its value has never been updated (with `calc_min()` or
|
||
* `= min_of()`), then the value of the view is set to @a x.
|
||
*
|
||
* @param x The value to minimize the view’s value with.
|
||
*
|
||
* @return A reference to the view. (Allows chaining
|
||
* `view.comp_min(a).comp_min(b)…`.)
|
||
*/
|
||
op_min_view& calc_min(const Type& x) { calc(x); return *this; }
|
||
|
||
/** Assign the result of a `min_of(view, value)` expression to the view.
|
||
*
|
||
* @param rhs An rhs_proxy value created by a `min_of(view, value)`
|
||
* expression.
|
||
*
|
||
* @return A reference to the view.
|
||
*
|
||
* @see min_max_internal::view_base::rhs_proxy
|
||
*/
|
||
op_min_view& operator=(const min_max_internal::rhs_proxy<op_min_view>& rhs)
|
||
{ assign(rhs); return *this; }
|
||
};
|
||
|
||
|
||
/** Compute the minimum of the value in a view and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or used in another min_of() call. For example,
|
||
*
|
||
* *reducer = min_of(*reducer, x);
|
||
* *reducer = min_of(x, *reducer);
|
||
*
|
||
* @see min_max_internal::view_base::rhs_proxy
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_view<Type, Compare> >
|
||
min_of(const op_min_view<Type, Compare>& view, const Type& value)
|
||
{
|
||
return min_max_internal::make_proxy(value, view);
|
||
}
|
||
|
||
/// @copydoc min_of(const op_min_view<Type, Compare>&, const Type&)
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_view<Type, Compare> >
|
||
min_of(const Type& value, const op_min_view<Type, Compare>& view)
|
||
{
|
||
return min_max_internal::make_proxy(value, view);
|
||
}
|
||
|
||
/** Nested minimum computation.
|
||
*
|
||
* Compute the minimum of the result of a min_of() call and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or wrapper, or used in another min_of() call. For example,
|
||
*
|
||
* *reducer = min_of(x, min_of(y, *reducer));
|
||
* wrapper = min_of(min_of(wrapper, x), y);
|
||
*
|
||
* @see min_max_internal::rhs_proxy
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_view<Type, Compare> >
|
||
min_of(const min_max_internal::rhs_proxy< op_min_view<Type, Compare> >& proxy,
|
||
const Type& value)
|
||
{
|
||
return proxy.calc(value);
|
||
}
|
||
|
||
/// @copydoc min_of(const min_max_internal::rhs_proxy< op_min_view<Type, Compare> >&, const Type&)
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_view<Type, Compare> >
|
||
min_of(const Type& value,
|
||
const min_max_internal::rhs_proxy< op_min_view<Type, Compare> >& proxy)
|
||
{
|
||
return proxy.calc(value);
|
||
}
|
||
|
||
|
||
/** Monoid class for minimum reductions. Instantiate the cilk::reducer template
|
||
* class with an op_min monoid to create a minimum reducer class. For example,
|
||
* to compute the minimum of a set of `int` values:
|
||
*
|
||
* cilk::reducer< cilk::op_min<int> > r;
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_min_view
|
||
*/
|
||
template <typename Type, typename Compare=std::less<Type>, bool Align = false>
|
||
class op_min : public min_max_internal::monoid_base<op_min_view<Type, Compare>, Align> {
|
||
typedef min_max_internal::monoid_base<op_min_view<Type, Compare>, Align> base;
|
||
public:
|
||
/// Construct with default comparator.
|
||
op_min() {}
|
||
/// Construct with specified comparator.
|
||
op_min(const Compare& compare) : base(compare) {}
|
||
};
|
||
|
||
//@}
|
||
|
||
|
||
/** @defgroup ReducersMinMaxMaxIndex Maximum reducers (value and index)
|
||
*
|
||
* These reducers will find the largest value from a set of values, and its
|
||
* index in the set.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
//@{
|
||
|
||
/** The maximum index reducer view class.
|
||
*
|
||
* This is the view class for reducers created with
|
||
* `cilk::reducer< cilk::op_max_index<Index, Type, Compare> >`. It accumulates
|
||
* the maximum, as determined by a comparator, of a set of values which have
|
||
* occurred as arguments to the `calc_max()` function, and records the index
|
||
* of the maximum value. The accumulated value will be the first argument `x`
|
||
* such that `compare(x, y)` is false for every argument `y`.
|
||
*
|
||
* If the comparator is `std::less`, then the accumulated value is the first
|
||
* argument value which is not less than any other argument value, i.e., the
|
||
* maximum.
|
||
*
|
||
* @note The reducer “dereference” operation (`reducer::operator *()`)
|
||
* yields a reference to the view. Thus, for example, the view class’s
|
||
* `calc_max()` function would be used in an expression like
|
||
* `r->calc_max(i, a)`where `r` is an op_max_index reducer
|
||
* variable.
|
||
*
|
||
* @note The word “index” suggests an integer index into an array, but there
|
||
* is no restriction on the index type or how it should be used. In
|
||
* general, it may be convenient to use it for any kind of key that
|
||
* can be used to locate the maximum value in the collection that it
|
||
* came from — for example:
|
||
* - An index into an array.
|
||
* - A key into an STL map.
|
||
* - An iterator into any STL container.
|
||
*
|
||
* @note A max_index reducer is essentially a max reducer whose value type
|
||
* is a `std::pair<Index, Type>`. This fact is camouflaged in the view
|
||
* `calc_max` function, the global `max_of` functions, and the reducer
|
||
* value constructor, which can all take an index argument and a value
|
||
* argument as an alternative to a single `std::pair` argument.
|
||
* However, the reducer `set_value()`, `get_value()`, `move_in()`, and
|
||
* `move_out()` functions work only with pairs, not with individual
|
||
* value and/or index arguments.
|
||
*
|
||
* @tparam Index The type of the indices associated with the values.
|
||
* @tparam Type The type of the values compared by the reducer. This will
|
||
* be the value type of a monoid_with_view that is
|
||
* instantiated with this view.
|
||
* @tparam Compare Used to compare the values. It must be a binary predicate.
|
||
* If it is omitted, then the view computes the conventional
|
||
* arithmetic maximum.
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_max_index
|
||
*/
|
||
template <typename Index, typename Type, typename Compare>
|
||
class op_max_index_view : public min_max_internal::view_base<
|
||
min_max_internal::index_view_content<Index, Type, Compare, true>,
|
||
Compare,
|
||
Compare>
|
||
{
|
||
typedef min_max_internal::view_base<
|
||
min_max_internal::index_view_content<Index, Type, Compare, true>,
|
||
Compare,
|
||
Compare> base;
|
||
using base::calc;
|
||
using base::assign;
|
||
typedef std::pair<Index, Type> pair_type;
|
||
friend class min_max_internal::rhs_proxy<op_max_index_view>;
|
||
|
||
public:
|
||
/** @name Constructors.
|
||
*
|
||
* All op_max_index_view constructors simply pass their arguments on to the
|
||
* @ref view_base base class, except for the `(index, value [, compare])`
|
||
* constructors, which create a `std::pair` containing the index and value.
|
||
*/
|
||
//@{
|
||
|
||
op_max_index_view() : base() {}
|
||
|
||
template <typename T1>
|
||
op_max_index_view(const T1& x1) : base(x1) {}
|
||
|
||
template <typename T1, typename T2>
|
||
op_max_index_view(const T1& x1, const T2& x2) : base(x1, x2) {}
|
||
|
||
template <typename T1, typename T2, typename T3>
|
||
op_max_index_view(const T1& x1, const T2& x2, const T3& x3) : base(x1, x2, x3) {}
|
||
|
||
op_max_index_view(const Index& i, const Type& v) : base(pair_type(i, v)) {}
|
||
|
||
op_max_index_view(const Index& i, const Type& v, const typename base::compare_type* c) :
|
||
base(pair_type(i, v), c) {}
|
||
|
||
//@}
|
||
|
||
/** Maximize with a value and index.
|
||
*
|
||
* If @a x is greater than the current value of the view (as defined by
|
||
* the reducer’s comparator), or if the view was created without an
|
||
* initial value and its value has never been updated (with `calc_max()`
|
||
* or `= max_of()`), then the value of the view is set to @a x, and the
|
||
* index is set to @a i..
|
||
*
|
||
* @param i The index of the value @a x.
|
||
* @param x The value to maximize the view’s value with.
|
||
*
|
||
* @return A reference to the view. (Allows
|
||
* `view.comp_max(i, a).comp_max(j, b)…`.)
|
||
*/
|
||
op_max_index_view& calc_max(const Index& i, const Type& x)
|
||
{ calc(pair_type(i, x)); return *this; }
|
||
|
||
/** Maximize with an index/value pair.
|
||
*
|
||
* If @a pair.second is greater than the current value of the view (as
|
||
* defined by the reducer’s comparator), or if the view was created
|
||
* without an initial value and its value has never been updated (with
|
||
* `calc_max()` or `= max_of()`), then the value of the view is set to
|
||
* @a pair.second, and the index is set to @a pair.first.
|
||
*
|
||
* @param pair A pair containing a value to maximize the view’s value
|
||
* with and its associated index.
|
||
*
|
||
* @return A reference to the view. (Allows
|
||
* `view.comp_max(p1).comp_max(p2)…`.)
|
||
*/
|
||
op_max_index_view& calc_max(const pair_type& pair)
|
||
{ calc(pair); return *this; }
|
||
|
||
/** Assign the result of a `max_of(view, index, value)` expression to the
|
||
* view.
|
||
*
|
||
* @param rhs An rhs_proxy value created by a `max_of(view, index, value)`
|
||
* expression.
|
||
*
|
||
* @return A reference to the view.
|
||
*
|
||
* @see min_max_internal::view_base::rhs_proxy
|
||
*/
|
||
op_max_index_view& operator=(const min_max_internal::rhs_proxy<op_max_index_view>& rhs)
|
||
{ assign(rhs); return *this; }
|
||
};
|
||
|
||
|
||
/** Compute the maximum of the value in a view and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or used in another max_of() call. For example,
|
||
*
|
||
* *reducer = max_of(*reducer, i, x);
|
||
* *reducer = max_of(i, x, *reducer);
|
||
*
|
||
* @see min_max_internal::rhs_proxy
|
||
*/
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const op_max_index_view<Index, Type, Compare>& view,
|
||
const Index& index, const Type& value)
|
||
{
|
||
return min_max_internal::make_proxy(std::pair<Index, Type>(index, value), view);
|
||
}
|
||
|
||
/// @copydoc max_of(const op_max_index_view<Index, Type, Compare>&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const Index& index, const Type& value,
|
||
const op_max_index_view<Index, Type, Compare>& view)
|
||
{
|
||
return min_max_internal::make_proxy(std::pair<Index, Type>(index, value), view);
|
||
}
|
||
|
||
/// @copydoc max_of(const op_max_index_view<Index, Type, Compare>&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const op_max_index_view<Index, Type, Compare>& view,
|
||
const std::pair<Index, Type>& pair)
|
||
{
|
||
return min_max_internal::make_proxy(pair, view);
|
||
}
|
||
|
||
/// @copydoc max_of(const op_max_index_view<Index, Type, Compare>&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const std::pair<Index, Type>& pair,
|
||
const op_max_index_view<Index, Type, Compare>& view)
|
||
{
|
||
return min_max_internal::make_proxy(pair, view);
|
||
}
|
||
|
||
/** Nested computation of the maximum of the value in a view and other values.
|
||
*
|
||
* Compute the maximum of the result of a max_of() call and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or used in another max_of() call. For example,
|
||
*
|
||
* *reducer = max_of(x, max_of(y, *reducer));
|
||
* *reducer = max_of(max_of(*reducer, x), y);
|
||
*
|
||
* @see min_max_internal::rhs_proxy
|
||
*/
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >& proxy,
|
||
const Index& index, const Type& value)
|
||
{
|
||
return proxy.calc(std::pair<Index, Type>(index, value));
|
||
}
|
||
|
||
/// @copydoc max_of(const min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const Index& index, const Type& value,
|
||
const min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >& proxy)
|
||
{
|
||
return proxy.calc(std::pair<Index, Type>(index, value));
|
||
}
|
||
|
||
/// @copydoc max_of(const min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >& proxy,
|
||
const std::pair<Index, Type>& pair)
|
||
{
|
||
return proxy.calc(pair);
|
||
}
|
||
|
||
/// @copydoc max_of(const min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >
|
||
max_of(const std::pair<Index, Type>& pair,
|
||
const min_max_internal::rhs_proxy< op_max_index_view<Index, Type, Compare> >& proxy)
|
||
{
|
||
return proxy.calc(pair);
|
||
}
|
||
|
||
|
||
/** Monoid class for maximum reductions with index. Instantiate the
|
||
* cilk::reducer template class with an op_max_index monoid to create a
|
||
* max_index reducer class. For example, to compute the maximum of an array of
|
||
* `double` values and the array index of the max value:
|
||
*
|
||
* cilk::reducer< cilk::op_max_index<unsigned, double> > r;
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_max_index_view
|
||
*/
|
||
template < typename Index
|
||
, typename Type
|
||
, typename Compare=std::less<Type>
|
||
, bool Align = false
|
||
>
|
||
class op_max_index : public min_max_internal::monoid_base<op_max_index_view<Index, Type, Compare>, Align>
|
||
{
|
||
typedef min_max_internal::monoid_base<
|
||
op_max_index_view<Index, Type, Compare>, Align> base;
|
||
public:
|
||
/// Construct with default comparator.
|
||
op_max_index() {}
|
||
/// Construct with specified comparator.
|
||
op_max_index(const Compare& compare) : base(compare) {}
|
||
};
|
||
|
||
//@}
|
||
|
||
|
||
|
||
/** @defgroup ReducersMinMaxMinIndex Minimum reducers (value and index)
|
||
*
|
||
* These reducers will find the smallest value from a set of values, and its
|
||
* index in the set.
|
||
*
|
||
* @ingroup ReducersMinMax
|
||
*/
|
||
//@{
|
||
|
||
/** The minimum index reducer view class.
|
||
*
|
||
* This is the view class for reducers created with
|
||
* `cilk::reducer<cilk::op_min_index<Index, Type, Compare> >`. It accumulates
|
||
* the minimum, as determined by a comparator, of a set of values which have
|
||
* occurred as arguments to the `calc_min()` function, and records the index
|
||
* of the minimum value. The accumulated value will be the first argument `x`
|
||
* such that `compare(y, x)` is false for every argument `y`.
|
||
*
|
||
* If the comparator is `std::less`, then the accumulated value is the first
|
||
* argument value which no other argument value is less than, i.e., the
|
||
* minimum.
|
||
*
|
||
* @note The reducer “dereference” operation (`reducer::operator *()`)
|
||
* yields a reference to the view. Thus, for example, the view class’s
|
||
* `calc_min()` function would be
|
||
* used in an expression like `r->calc_min(i, a)`where `r` is an
|
||
* op_min_index reducer variable.
|
||
*
|
||
* @note The word “index” suggests an integer index into an array, but there
|
||
* is no restriction on the index type or how it should be used. In
|
||
* general, it may be convenient to use it for any kind of key that
|
||
* can be used to locate the minimum value in the collection that it
|
||
* came from — for example:
|
||
* - An index into an array.
|
||
* - A key into an STL map.
|
||
* - An iterator into any STL container.
|
||
*
|
||
* @note A min_index reducer is essentially a min reducer whose value type
|
||
* is a `std::pair<Index, Type>`. This fact is camouflaged in the view
|
||
* `calc_min` function, the global `min_of` functions, and the reducer
|
||
* value constructor, which can all take an index argument and a value
|
||
* argument as an alternative to a single `std::pair` argument.
|
||
* However, the reducer `set_value()`, `get_value()`, `move_in()`, and
|
||
* `move_out()` functions work only with pairs, not with individual
|
||
* value and/or index arguments.
|
||
*
|
||
* @tparam Index The type of the indices associated with the values.
|
||
* @tparam Type The type of the values compared by the reducer. This will
|
||
* be the value type of a monoid_with_view that is
|
||
* instantiated with this view.
|
||
* @tparam Compare Used to compare the values. It must be a binary predicate.
|
||
* If it is omitted, then the view computes the conventional
|
||
* arithmetic minimum.
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_min_index
|
||
*/
|
||
template <typename Index, typename Type, typename Compare>
|
||
class op_min_index_view : public min_max_internal::view_base<
|
||
min_max_internal::index_view_content<Index, Type, Compare, false>,
|
||
Compare,
|
||
min_max_internal::reverse_predicate<Compare, Type> >
|
||
{
|
||
typedef min_max_internal::view_base<
|
||
min_max_internal::index_view_content<Index, Type, Compare, false>,
|
||
Compare,
|
||
min_max_internal::reverse_predicate<Compare, Type> > base;
|
||
using base::calc;
|
||
using base::assign;
|
||
typedef std::pair<Index, Type> pair_type;
|
||
friend class min_max_internal::rhs_proxy<op_min_index_view>;
|
||
|
||
public:
|
||
/** @name Constructors.
|
||
*
|
||
* All op_min_index_view constructors simply pass their arguments on to the
|
||
* @ref view_base base class, except for the `(index, value [, compare])`
|
||
* constructors, which create a `std::pair` containing the index and value.
|
||
*/
|
||
//@{
|
||
|
||
op_min_index_view() : base() {}
|
||
|
||
template <typename T1>
|
||
op_min_index_view(const T1& x1) : base(x1) {}
|
||
|
||
template <typename T1, typename T2>
|
||
op_min_index_view(const T1& x1, const T2& x2) : base(x1, x2) {}
|
||
|
||
template <typename T1, typename T2, typename T3>
|
||
op_min_index_view(const T1& x1, const T2& x2, const T3& x3) : base(x1, x2, x3) {}
|
||
|
||
op_min_index_view(const Index& i, const Type& v) : base(pair_type(i, v)) {}
|
||
|
||
op_min_index_view(const Index& i, const Type& v, const typename base::compare_type* c) :
|
||
base(pair_type(i, v), c) {}
|
||
|
||
//@}
|
||
|
||
/** Minimize with a value and index.
|
||
*
|
||
* If @a x is greater than the current value of the view (as defined by
|
||
* the reducer’s comparator), or if the view was created without an
|
||
* initial value and its value has never been updated (with `calc_min()`
|
||
* or `= min_of()`), then the value of the view is set to @a x, and the
|
||
* index is set to @a i..
|
||
*
|
||
* @param i The index of the value @a x.
|
||
* @param x The value to minimize the view’s value with.
|
||
*
|
||
* @return A reference to the view. (Allows
|
||
* `view.comp_min(i, a).comp_min(j, b)…`.)
|
||
*/
|
||
op_min_index_view& calc_min(const Index& i, const Type& x)
|
||
{ calc(pair_type(i, x)); return *this; }
|
||
|
||
/** Maximize with an index/value pair.
|
||
*
|
||
* If @a pair.second is less than the current value of the view (as
|
||
* defined by the reducer’s comparator), or if the view was created
|
||
* without an initial value and its value has never been updated (with
|
||
* `calc_min()` or `= min_of()`), then the value of the view is set to
|
||
* @a pair.second, and the index is set to @a pair.first.
|
||
*
|
||
* @param pair A pair containing a value to minimize the view’s value
|
||
* with and its associated index.
|
||
*
|
||
* @return A reference to the view. (Allows
|
||
* `view.comp_min(p1).comp_min(p2)…`.)
|
||
*/
|
||
op_min_index_view& calc_min(const pair_type& pair)
|
||
{ calc(pair); return *this; }
|
||
|
||
/** Assign the result of a `min_of(view, index, value)` expression to the
|
||
* view.
|
||
*
|
||
* @param rhs An rhs_proxy value created by a `min_of(view, index, value)`
|
||
* expression.
|
||
*
|
||
* @return A reference to the view.
|
||
*
|
||
* @see min_max_internal::view_base::rhs_proxy
|
||
*/
|
||
op_min_index_view& operator=(const min_max_internal::rhs_proxy<op_min_index_view>& rhs)
|
||
{ assign(rhs); return *this; }
|
||
};
|
||
|
||
|
||
/** Compute the minimum of the value in a view and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or used in another min_of() call. For example,
|
||
*
|
||
* *reducer = min_of(*reducer, i, x);
|
||
* *reducer = min_of(i, x, *reducer);
|
||
*
|
||
* @see min_max_internal::min_min_view_base::rhs_proxy
|
||
*/
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const op_min_index_view<Index, Type, Compare>& view,
|
||
const Index& index, const Type& value)
|
||
{
|
||
return min_max_internal::make_proxy(std::pair<Index, Type>(index, value), view);
|
||
}
|
||
|
||
/// @copydoc min_of(const op_min_index_view<Index, Type, Compare>&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const Index& index, const Type& value,
|
||
const op_min_index_view<Index, Type, Compare>& view)
|
||
{
|
||
return min_max_internal::make_proxy(std::pair<Index, Type>(index, value), view);
|
||
}
|
||
|
||
/// @copydoc min_of(const op_min_index_view<Index, Type, Compare>&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const op_min_index_view<Index, Type, Compare>& view,
|
||
const std::pair<Index, Type>& pair)
|
||
{
|
||
return min_max_internal::make_proxy(pair, view);
|
||
}
|
||
|
||
/// @copydoc min_of(const op_min_index_view<Index, Type, Compare>&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const std::pair<Index, Type>& pair,
|
||
const op_min_index_view<Index, Type, Compare>& view)
|
||
{
|
||
return min_max_internal::make_proxy(pair, view);
|
||
}
|
||
|
||
/** Nested computation of the minimum of the value in a view and other values.
|
||
*
|
||
* Compute the minimum of the result of a min_of() call and another value.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* view or used in another min_of() call. For example,
|
||
*
|
||
* *reducer = min_of(x, min_of(y, *reducer));
|
||
* *reducer = min_of(min_of(*reducer, x), y);
|
||
*
|
||
* @see min_max_internal::min_min_view_base::rhs_proxy
|
||
*/
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >& proxy,
|
||
const Index& index, const Type& value)
|
||
{
|
||
return proxy.calc(std::pair<Index, Type>(index, value));
|
||
}
|
||
|
||
/// @copydoc min_of(const min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const Index& index, const Type& value,
|
||
const min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >& proxy)
|
||
{
|
||
return proxy.calc(std::pair<Index, Type>(index, value));
|
||
}
|
||
|
||
/// @copydoc min_of(const min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >& proxy,
|
||
const std::pair<Index, Type>& pair)
|
||
{
|
||
return proxy.calc(pair);
|
||
}
|
||
|
||
/// @copydoc min_of(const min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >&, const Index&, const Type&)
|
||
template <typename Index, typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >
|
||
min_of(const std::pair<Index, Type>& pair,
|
||
const min_max_internal::rhs_proxy< op_min_index_view<Index, Type, Compare> >& proxy)
|
||
{
|
||
return proxy.calc(pair);
|
||
}
|
||
|
||
|
||
/** Monoid class for minimum reductions with index. Instantiate the
|
||
* cilk::reducer template class with an op_min_index monoid to create a
|
||
* min_index reducer class. For example, to compute the minimum of an array of
|
||
* `double` values and the array index of the min value:
|
||
*
|
||
* cilk::reducer< cilk::op_min_index<unsigned, double> > r;
|
||
*
|
||
* @see ReducersMinMax
|
||
* @see op_min_index_view
|
||
*/
|
||
template < typename Index
|
||
, typename Type
|
||
, typename Compare=std::less<Type>
|
||
, bool Align = false
|
||
>
|
||
class op_min_index : public min_max_internal::monoid_base<op_min_index_view<Index, Type, Compare>, Align>
|
||
{
|
||
typedef min_max_internal::monoid_base<
|
||
op_min_index_view<Index, Type, Compare>, Align> base;
|
||
public:
|
||
/// Construct with default comparator.
|
||
op_min_index() {}
|
||
/// Construct with specified comparator.
|
||
op_min_index(const Compare& compare) : base(compare) {}
|
||
};
|
||
|
||
//@}
|
||
|
||
|
||
/** Deprecated maximum reducer wrapper class.
|
||
*
|
||
* reducer_max is the same as @ref reducer<@ref op_max>, except that
|
||
* reducer_max is a proxy for the contained view, so that accumulator
|
||
* variable update operations can be applied directly to the reducer. For
|
||
* example, a value is maximized with a `reducer<%op_max>` with
|
||
* `r->calc_max(a)`, but a value can be maximized with a `%reducer_max` with
|
||
* `r.calc_max(a)`.
|
||
*
|
||
*
|
||
* @deprecated Users are strongly encouraged to use `reducer<monoid>`
|
||
* reducers rather than the old wrappers like reducer_max.
|
||
* The `reducer<monoid>` reducers show the reducer/monoid/view
|
||
* architecture more clearly, are more consistent in their
|
||
* implementation, and present a simpler model for new
|
||
* user-implemented reducers.
|
||
*
|
||
* @note Implicit conversions are provided between `%reducer_max`
|
||
* and `reducer<%op_max>`. This allows incremental code
|
||
* conversion: old code that used `%reducer_max` can pass a
|
||
* `%reducer_max` to a converted function that now expects a
|
||
* pointer or reference to a `reducer<%op_max>`, and vice
|
||
* versa. **But see @ref redminmax_compatibility.**
|
||
*
|
||
* @tparam Type The value type of the reducer.
|
||
* @tparam Compare The “less than” comparator type for the reducer.
|
||
*
|
||
* @see op_max
|
||
* @see op_max_view
|
||
* @see reducer
|
||
* @see ReducersMinMax
|
||
* @ingroup ReducersMinMaxMaxValue
|
||
*/
|
||
template <typename Type, typename Compare=std::less<Type> >
|
||
class reducer_max : public reducer< op_max<Type, Compare, true> >
|
||
{
|
||
__CILKRTS_STATIC_ASSERT(
|
||
::cilk::internal::class_is_empty<
|
||
typename ::cilk::internal::binary_functor<Compare>::type >::value,
|
||
"cilk::reducer_max<Type, Compare> only works with "
|
||
"an empty Compare class");
|
||
typedef reducer< op_max<Type, Compare, true> > base;
|
||
public:
|
||
|
||
/// Type of data in a reducer_max.
|
||
typedef Type basic_value_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type view_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type View;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type monoid_type;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type Monoid;
|
||
|
||
/// The view’s rhs proxy type.
|
||
typedef min_max_internal::rhs_proxy<View> rhs_proxy;
|
||
|
||
using base::view;
|
||
|
||
/** @name Constructors
|
||
*/
|
||
//@{
|
||
|
||
/// Construct the wrapper in its identity state (either `!is_set()`, or
|
||
/// `value() == identity value`).
|
||
reducer_max() : base() {}
|
||
|
||
/// Construct the wrapper with a specified initial value.
|
||
explicit reducer_max(const Type& initial_value) : base(initial_value) {}
|
||
|
||
/// Construct the wrapper in its identity state with a specified
|
||
/// comparator.
|
||
explicit reducer_max(const Compare& comp) : base(comp) {}
|
||
|
||
/// Construct the wrapper with a specified initial value and a specified
|
||
/// comparator.
|
||
reducer_max(const Type& initial_value, const Compare& comp)
|
||
: base(initial_value, comp) {}
|
||
|
||
//@}
|
||
|
||
/** @name Forwarded functions
|
||
* @details Functions that update the contained accumulator variable are
|
||
* simply forwarded to the contained @ref op_max_view. */
|
||
//@{
|
||
|
||
/// @copydoc cilk_lib_1_0::min_max_internal::view_content::is_set() const
|
||
bool is_set() const { return view().is_set(); }
|
||
|
||
/// @copydoc op_max_view::calc_max(const Type&)
|
||
reducer_max& calc_max(const Type& x)
|
||
{ view().calc_max(x); return *this; }
|
||
|
||
/// @copydoc op_max_view::operator=(const min_max_internal::rhs_proxy<op_max_view>&)
|
||
reducer_max& operator=(const rhs_proxy& rhs)
|
||
{ view() = rhs; return *this; }
|
||
|
||
//@}
|
||
|
||
/** Allow read-only access to the value within the current view.
|
||
*
|
||
* @returns A const reference to the value within the current view.
|
||
*/
|
||
const Type& get_reference() const { return view().get_reference(); }
|
||
|
||
/// @name Dereference
|
||
/** Dereferencing a wrapper is a no-op. It simply returns the wrapper.
|
||
* Combined with the rule that a wrapper forwards view operations to the
|
||
* view, this means that view operations can be written the same way on
|
||
* reducers and wrappers, which is convenient for incrementally
|
||
* converting code using wrappers to code using reducers. That is:
|
||
*
|
||
* reducer< op_max<int> > r;
|
||
* r->calc_max(a); // *r returns the view
|
||
* // calc_max is a view member function
|
||
*
|
||
* reducer_max<int> w;
|
||
* w->calc_max(a); // *w returns the wrapper
|
||
* // calc_max is a wrapper member function that
|
||
* // calls the corresponding view function
|
||
*/
|
||
//@{
|
||
reducer_max& operator*() { return *this; }
|
||
reducer_max const& operator*() const { return *this; }
|
||
|
||
reducer_max* operator->() { return this; }
|
||
reducer_max const* operator->() const { return this; }
|
||
//@}
|
||
|
||
/** @name Upcast
|
||
* @details In Cilk library 0.9, reducers were always cache-aligned. In
|
||
* library 1.0, reducer cache alignment is optional. By default, reducers
|
||
* are unaligned (i.e., just naturally aligned), but legacy wrappers
|
||
* inherit from cache-aligned reducers for binary compatibility.
|
||
*
|
||
* This means that a wrapper will automatically be upcast to its aligned
|
||
* reducer base class. The following conversion operators provide
|
||
* pseudo-upcasts to the corresponding unaligned reducer class.
|
||
*/
|
||
//@{
|
||
operator reducer< op_max<Type, Compare, false> >& ()
|
||
{
|
||
return *reinterpret_cast< reducer< op_max<Type, Compare, false> >* >(this);
|
||
}
|
||
|
||
operator const reducer< op_max<Type, Compare, false> >& () const
|
||
{
|
||
return *reinterpret_cast< const reducer< op_max<Type, Compare, false> >* >(this);
|
||
}
|
||
//@}
|
||
};
|
||
|
||
|
||
/// @cond internal
|
||
// The legacy definition of max_of(reducer_max, value) has different
|
||
// behavior and a different return type than this definition. We add an
|
||
// unused third argument to this version of the function to give it a different
|
||
// signature, so that they won’t end up sharing a single object file entry.
|
||
struct max_of_1_0_t {};
|
||
const max_of_1_0_t max_of_1_0 = {};
|
||
/// @endcond
|
||
|
||
/** Compute the maximum of the value in a reducer_max and another value.
|
||
*
|
||
* @deprecated Because reducer_max is deprecated.
|
||
*
|
||
* The result of this computation can only be assigned back to the original
|
||
* reducer or used in another max_of() call. For example,
|
||
*
|
||
* reducer = max_of(reducer, x);
|
||
* reducer = max_of(x, reducer);
|
||
*
|
||
* @see min_max_internal::rhs_proxy
|
||
*
|
||
* @ingroup ReducersMinMaxMaxValue
|
||
*/
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_view<Type, Compare> >
|
||
max_of(const reducer_max<Type, Compare>& r, const Type& value,
|
||
const max_of_1_0_t& = max_of_1_0)
|
||
{
|
||
return min_max_internal::make_proxy(value, r.view());
|
||
}
|
||
|
||
/// @copydoc max_of(const reducer_max<Type, Compare>&, const Type&, const max_of_1_0_t&)
|
||
/// @ingroup ReducersMinMaxMaxValue
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_max_view<Type, Compare> >
|
||
max_of(const Type& value, const reducer_max<Type, Compare>& r,
|
||
const max_of_1_0_t& = max_of_1_0)
|
||
{
|
||
return min_max_internal::make_proxy(value, r.view());
|
||
}
|
||
|
||
|
||
/** Deprecated minimum reducer wrapper class.
|
||
*
|
||
* reducer_min is the same as @ref reducer<@ref op_min>, except that
|
||
* reducer_min is a proxy for the contained view, so that accumulator
|
||
* variable update operations can be applied directly to the reducer. For
|
||
* example, a value is minimized with a `reducer<%op_min>` with
|
||
* `r->calc_min(a)`, but a value can be minimized with a `%reducer_min` with
|
||
* `r.calc_min(a)`.
|
||
*
|
||
*
|
||
* @deprecated Users are strongly encouraged to use `reducer<monoid>`
|
||
* reducers rather than the old wrappers like reducer_min.
|
||
* The `reducer<monoid>` reducers show the reducer/monoid/view
|
||
* architecture more clearly, are more consistent in their
|
||
* implementation, and present a simpler model for new
|
||
* user-implemented reducers.
|
||
*
|
||
* @note Implicit conversions are provided between `%reducer_min`
|
||
* and `reducer<%op_min>`. This allows incremental code
|
||
* conversion: old code that used `%reducer_min` can pass a
|
||
* `%reducer_min` to a converted function that now expects a
|
||
* pointer or reference to a `reducer<%op_min>`, and vice
|
||
* versa. **But see @ref redminmax_compatibility.**
|
||
*
|
||
* @tparam Type The value type of the reducer.
|
||
* @tparam Compare The “less than” comparator type for the reducer.
|
||
*
|
||
* @see op_min
|
||
* @see op_min_view
|
||
* @see reducer
|
||
* @see ReducersMinMax
|
||
* @ingroup ReducersMinMaxMinValue
|
||
*/
|
||
template <typename Type, typename Compare=std::less<Type> >
|
||
class reducer_min : public reducer< op_min<Type, Compare, true> >
|
||
{
|
||
__CILKRTS_STATIC_ASSERT(
|
||
::cilk::internal::class_is_empty<
|
||
typename ::cilk::internal::binary_functor<Compare>::type >::value,
|
||
"cilk::reducer_min<Type, Compare> only works with "
|
||
"an empty Compare class");
|
||
typedef reducer< op_min<Type, Compare, true> > base;
|
||
public:
|
||
|
||
/// Type of data in a reducer_min.
|
||
typedef Type basic_value_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type view_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type View;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type monoid_type;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type Monoid;
|
||
|
||
/// The view’s rhs proxy type.
|
||
typedef min_max_internal::rhs_proxy<View> rhs_proxy;
|
||
|
||
using base::view;
|
||
|
||
/** @name Constructors
|
||
*/
|
||
//@{
|
||
|
||
/// Construct the wrapper in its identity state (either `!is_set()`, or
|
||
/// `value() == identity value`).
|
||
reducer_min() : base() {}
|
||
|
||
/// Construct the wrapper with a specified initial value.
|
||
explicit reducer_min(const Type& initial_value) : base(initial_value) {}
|
||
|
||
/// Construct the wrapper in its identity state with a specified
|
||
/// comparator.
|
||
explicit reducer_min(const Compare& comp) : base(comp) {}
|
||
|
||
/// Construct the wrapper with a specified initial value and a specified
|
||
/// comparator.
|
||
reducer_min(const Type& initial_value, const Compare& comp)
|
||
: base(initial_value, comp) {}
|
||
|
||
//@}
|
||
|
||
/** @name Forwarded functions
|
||
* @details Functions that update the contained accumulator variable are
|
||
* simply forwarded to the contained @ref op_min_view. */
|
||
//@{
|
||
|
||
/// @copydoc cilk_lib_1_0::min_max_internal::view_content::is_set() const
|
||
bool is_set() const { return view().is_set(); }
|
||
|
||
/// @copydoc op_min_view::calc_min(const Type&)
|
||
reducer_min& calc_min(const Type& x)
|
||
{ view().calc_min(x); return *this; }
|
||
|
||
/// @copydoc op_min_view::operator=(const min_max_internal::rhs_proxy<op_min_view>&)
|
||
reducer_min& operator=(const rhs_proxy& rhs)
|
||
{ view() = rhs; return *this; }
|
||
|
||
//@}
|
||
|
||
/** Allow read-only access to the value within the current view.
|
||
*
|
||
* @returns A const reference to the value within the current view.
|
||
*/
|
||
const Type& get_reference() const { return view().get_reference(); }
|
||
|
||
/// @name Dereference
|
||
/** Dereferencing a wrapper is a no-op. It simply returns the wrapper.
|
||
* Combined with the rule that a wrapper forwards view operations to the
|
||
* view, this means that view operations can be written the same way on
|
||
* reducers and wrappers, which is convenient for incrementally
|
||
* converting code using wrappers to code using reducers. That is:
|
||
*
|
||
* reducer< op_min<int> > r;
|
||
* r->calc_min(a); // *r returns the view
|
||
* // calc_min is a view member function
|
||
*
|
||
* reducer_min<int> w;
|
||
* w->calc_min(a); // *w returns the wrapper
|
||
* // calc_min is a wrapper member function that
|
||
* // calls the corresponding view function
|
||
*/
|
||
//@{
|
||
reducer_min& operator*() { return *this; }
|
||
reducer_min const& operator*() const { return *this; }
|
||
|
||
reducer_min* operator->() { return this; }
|
||
reducer_min const* operator->() const { return this; }
|
||
//@}
|
||
|
||
/** @name Upcast
|
||
* @details In Cilk library 0.9, reducers were always cache-aligned. In
|
||
* library 1.0, reducer cache alignment is optional. By default, reducers
|
||
* are unaligned (i.e., just naturally aligned), but legacy wrappers
|
||
* inherit from cache-aligned reducers for binary compatibility.
|
||
*
|
||
* This means that a wrapper will automatically be upcast to its aligned
|
||
* reducer base class. The following conversion operators provide
|
||
* pseudo-upcasts to the corresponding unaligned reducer class.
|
||
*/
|
||
//@{
|
||
operator reducer< op_min<Type, Compare, false> >& ()
|
||
{
|
||
return *reinterpret_cast< reducer< op_min<Type, Compare, false> >* >(this);
|
||
}
|
||
|
||
operator const reducer< op_min<Type, Compare, false> >& () const
|
||
{
|
||
return *reinterpret_cast< const reducer< op_min<Type, Compare, false> >* >(this);
|
||
}
|
||
//@}
|
||
};
|
||
|
||
|
||
/** Compute the minimum of a reducer and a value.
|
||
*
|
||
* @deprecated Because reducer_min is deprecated.
|
||
*/
|
||
//@{
|
||
// The legacy definition of min_of(reducer_min, value) has different
|
||
// behavior and a different return type than this definition. We add an
|
||
// unused third argument to this version of the function to give it a different
|
||
// signature, so that they won’t end up sharing a single object file entry.
|
||
struct min_of_1_0_t {};
|
||
const min_of_1_0_t min_of_1_0 = {};
|
||
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_view<Type, Compare> >
|
||
min_of(const reducer_min<Type, Compare>& r, const Type& value,
|
||
const min_of_1_0_t& = min_of_1_0)
|
||
{
|
||
return min_max_internal::make_proxy(value, r.view());
|
||
}
|
||
|
||
template <typename Type, typename Compare>
|
||
inline min_max_internal::rhs_proxy< op_min_view<Type, Compare> >
|
||
min_of(const Type& value, const reducer_min<Type, Compare>& r,
|
||
const min_of_1_0_t& = min_of_1_0)
|
||
{
|
||
return min_max_internal::make_proxy(value, r.view());
|
||
}
|
||
//@}
|
||
|
||
|
||
/** Deprecated maximum with index reducer wrapper class.
|
||
*
|
||
* reducer_max_index is the same as @ref reducer<@ref op_max_index>, except
|
||
* that reducer_max_index is a proxy for the contained view, so that
|
||
* accumulator variable update operations can be applied directly to the
|
||
* reducer. For example, a value is maximized with a `reducer<%op_max_index>`
|
||
* with `r->calc_max(i, a)`, but a value can be maximized with a
|
||
* `%reducer_max` with `r.calc_max(i, aa)`.
|
||
*
|
||
*
|
||
* @deprecated Users are strongly encouraged to use `reducer<monoid>`
|
||
* reducers rather than the old wrappers like reducer_max.
|
||
* The `reducer<monoid>` reducers show the reducer/monoid/view
|
||
* architecture more clearly, are more consistent in their
|
||
* implementation, and present a simpler model for new
|
||
* user-implemented reducers.
|
||
*
|
||
* @note Implicit conversions are provided between `%reducer_max_index`
|
||
* and `reducer<%op_max_index>`. This allows incremental code
|
||
* conversion: old code that used `%reducer_max_index` can pass a
|
||
* `%reducer_max_index` to a converted function that now expects a
|
||
* pointer or reference to a `reducer<%op_max_index>`, and vice
|
||
* versa. **But see @ref redminmax_compatibility.**
|
||
*
|
||
* @tparam Index The index type of the reducer.
|
||
* @tparam Type The value type of the reducer.
|
||
* @tparam Compare The “less than” comparator type for the reducer.
|
||
*
|
||
* @see op_max_index
|
||
* @see op_max_index_view
|
||
* @see reducer
|
||
* @see ReducersMinMax
|
||
* @ingroup ReducersMinMaxMaxIndex
|
||
*/
|
||
template < typename Index
|
||
, typename Type
|
||
, typename Compare = std::less<Type>
|
||
>
|
||
class reducer_max_index :
|
||
public reducer< op_max_index<Index, Type, Compare, true> >
|
||
{
|
||
__CILKRTS_STATIC_ASSERT(
|
||
::cilk::internal::class_is_empty<
|
||
typename ::cilk::internal::binary_functor<Compare>::type >::value,
|
||
"cilk::reducer_max_index<Type, Compare> only works with "
|
||
"an empty Compare class");
|
||
typedef reducer< op_max_index<Index, Type, Compare, true> > base;
|
||
public:
|
||
|
||
/// Type of data in a reducer_max_index.
|
||
typedef Type basic_value_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type view_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type View;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type monoid_type;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type Monoid;
|
||
|
||
/// The view’s rhs proxy type.
|
||
typedef min_max_internal::rhs_proxy<View> rhs_proxy;
|
||
|
||
using base::view;
|
||
|
||
/** @name Constructors
|
||
*/
|
||
//@{
|
||
|
||
/// Construct the wrapper in its identity state (`!is_set()`).
|
||
reducer_max_index() : base() {}
|
||
|
||
/// Construct with a specified initial index and value.
|
||
reducer_max_index(const Index& initial_index,
|
||
const Type& initial_value)
|
||
: base(initial_index, initial_value) {}
|
||
|
||
/// Construct the wrapper with a specified comparator.
|
||
explicit reducer_max_index(const Compare& comp) : base(comp) {}
|
||
|
||
/// Construct the wrapper with a specified initial index, value,
|
||
/// and comparator.
|
||
reducer_max_index(const Index& initial_index,
|
||
const Type& initial_value,
|
||
const Compare& comp)
|
||
: base(initial_index, initial_value, comp) {}
|
||
|
||
//@}
|
||
|
||
/** @name Set / Get
|
||
*/
|
||
//@{
|
||
|
||
/// Set the index and value of this object.
|
||
void set_value(const Index& index, const Type& value)
|
||
{ base::set_value(std::make_pair(index, value)); }
|
||
|
||
/// Return the maximum value.
|
||
const Type& get_value() const
|
||
{ return view().get_reference(); }
|
||
|
||
/// Return the maximum index.
|
||
const Index& get_index() const
|
||
{ return view().get_index_reference(); }
|
||
|
||
/// Return a const reference to value data member in the view.
|
||
const Type& get_reference() const
|
||
{ return view().get_reference(); }
|
||
|
||
/// Return a const reference to index data member in the view.
|
||
const Index& get_index_reference() const
|
||
{ return view().get_index_reference(); }
|
||
|
||
//@}
|
||
|
||
/** @name Forwarded functions
|
||
* @details Functions that update the contained accumulator variable are
|
||
* simply forwarded to the contained @ref op_max_view. */
|
||
//@{
|
||
|
||
/// @copydoc cilk_lib_1_0::min_max_internal::view_content::is_set() const
|
||
bool is_set() const { return view().is_set(); }
|
||
|
||
/// @copydoc op_max_index_view::calc_max(const Index&, const Type&)
|
||
reducer_max_index& calc_max(const Index& i, const Type& x)
|
||
{ view().calc_max(i, x); return *this; }
|
||
|
||
/// @copydoc op_max_view::operator=(const min_max_internal::rhs_proxy<op_max_view>&)
|
||
reducer_max_index& operator=(const rhs_proxy& rhs)
|
||
{ view() = rhs; return *this; }
|
||
|
||
//@}
|
||
|
||
/// @name Dereference
|
||
/** Dereferencing a wrapper is a no-op. It simply returns the wrapper.
|
||
* Combined with the rule that a wrapper forwards view operations to the
|
||
* view, this means that view operations can be written the same way on
|
||
* reducers and wrappers, which is convenient for incrementally
|
||
* converting code using wrappers to code using reducers. That is:
|
||
*
|
||
* reducer< op_max_index<int, int> > r;
|
||
* r->calc_max(i, a); // *r returns the view
|
||
* // calc_max is a view member function
|
||
*
|
||
* reducer_max_index<int, int> w;
|
||
* w->calc_max(i, a); // *w returns the wrapper
|
||
* // calc_max is a wrapper member function that
|
||
* // calls the corresponding view function
|
||
*/
|
||
//@{
|
||
reducer_max_index& operator*() { return *this; }
|
||
reducer_max_index const& operator*() const { return *this; }
|
||
|
||
reducer_max_index* operator->() { return this; }
|
||
reducer_max_index const* operator->() const { return this; }
|
||
//@}
|
||
|
||
/** @name Upcast
|
||
* @details In Cilk library 0.9, reducers were always cache-aligned. In
|
||
* library 1.0, reducer cache alignment is optional. By default, reducers
|
||
* are unaligned (i.e., just naturally aligned), but legacy wrappers
|
||
* inherit from cache-aligned reducers for binary compatibility.
|
||
*
|
||
* This means that a wrapper will automatically be upcast to its aligned
|
||
* reducer base class. The following conversion operators provide
|
||
* pseudo-upcasts to the corresponding unaligned reducer class.
|
||
*/
|
||
//@{
|
||
operator reducer< op_max_index<Index, Type, Compare, false> >& ()
|
||
{
|
||
return *reinterpret_cast< reducer< op_max_index<Index, Type, Compare, false> >* >(this);
|
||
}
|
||
|
||
operator const reducer< op_max_index<Index, Type, Compare, false> >& () const
|
||
{
|
||
return *reinterpret_cast< const reducer< op_max_index<Index, Type, Compare, false> >* >(this);
|
||
}
|
||
//@}
|
||
|
||
};
|
||
|
||
|
||
/** Deprecated minimum with index reducer wrapper class.
|
||
*
|
||
* reducer_min_index is the same as @ref reducer<@ref op_min_index>, except
|
||
* that reducer_min_index is a proxy for the contained view, so that
|
||
* accumulator variable update operations can be applied directly to the
|
||
* reducer. For example, a value is minimized with a `reducer<%op_min_index>`
|
||
* with `r->calc_min(i, a)`, but a value can be minimized with a
|
||
* `%reducer_min` with `r.calc_min(i, aa)`.
|
||
*
|
||
*
|
||
* @deprecated Users are strongly encouraged to use `reducer<monoid>`
|
||
* reducers rather than the old wrappers like reducer_min.
|
||
* The `reducer<monoid>` reducers show the reducer/monoid/view
|
||
* architecture more clearly, are more consistent in their
|
||
* implementation, and present a simpler model for new
|
||
* user-implemented reducers.
|
||
*
|
||
* @note Implicit conversions are provided between `%reducer_min_index`
|
||
* and `reducer<%op_min_index>`. This allows incremental code
|
||
* conversion: old code that used `%reducer_min_index` can pass a
|
||
* `%reducer_min_index` to a converted function that now expects a
|
||
* pointer or reference to a `reducer<%op_min_index>`, and vice
|
||
* versa. **But see @ref redminmax_compatibility.**
|
||
*
|
||
* @tparam Index The index type of the reducer.
|
||
* @tparam Type The value type of the reducer.
|
||
* @tparam Compare The “less than” comparator type for the reducer.
|
||
*
|
||
* @see op_min_index
|
||
* @see op_min_index_view
|
||
* @see reducer
|
||
* @see ReducersMinMax
|
||
* @ingroup ReducersMinMaxMinIndex
|
||
*/
|
||
template < typename Index
|
||
, typename Type
|
||
, typename Compare = std::less<Type>
|
||
>
|
||
class reducer_min_index :
|
||
public reducer< op_min_index<Index, Type, Compare, true> >
|
||
{
|
||
__CILKRTS_STATIC_ASSERT(
|
||
::cilk::internal::class_is_empty<
|
||
typename ::cilk::internal::binary_functor<Compare>::type >::value,
|
||
"cilk::reducer_min_index<Type, Compare> only works with "
|
||
"an empty Compare class");
|
||
typedef reducer< op_min_index<Index, Type, Compare, true> > base;
|
||
public:
|
||
|
||
/// Type of data in a reducer_min_index.
|
||
typedef Type basic_value_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type view_type;
|
||
|
||
/// The view type for the reducer.
|
||
typedef typename base::view_type View;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type monoid_type;
|
||
|
||
/// The monoid type for the reducer.
|
||
typedef typename base::monoid_type Monoid;
|
||
|
||
/// The view’s rhs proxy type.
|
||
typedef min_max_internal::rhs_proxy<View> rhs_proxy;
|
||
|
||
using base::view;
|
||
|
||
/** @name Constructors
|
||
*/
|
||
//@{
|
||
|
||
/// Construct the wrapper in its identity state (`!is_set()`).
|
||
reducer_min_index() : base() {}
|
||
|
||
/// Construct with a specified initial index and value.
|
||
reducer_min_index(const Index& initial_index,
|
||
const Type& initial_value)
|
||
: base(initial_index, initial_value) {}
|
||
|
||
/// Construct the wrapper with a specified comparator.
|
||
explicit reducer_min_index(const Compare& comp) : base(comp) {}
|
||
|
||
/// Construct the wrapper with a specified initial index, value,
|
||
/// and comparator.
|
||
reducer_min_index(const Index& initial_index,
|
||
const Type& initial_value,
|
||
const Compare& comp)
|
||
: base(initial_index, initial_value, comp) {}
|
||
|
||
//@}
|
||
|
||
/** @name Set / Get
|
||
*/
|
||
//@{
|
||
|
||
/// Set the index and value of this object.
|
||
void set_value(const Index& index, const Type& value)
|
||
{ base::set_value(std::make_pair(index, value)); }
|
||
|
||
/// Return the minimum value.
|
||
const Type& get_value() const
|
||
{ return view().get_reference(); }
|
||
|
||
/// Return the minimum index.
|
||
const Index& get_index() const
|
||
{ return view().get_index_reference(); }
|
||
|
||
/// Return a const reference to value data member in the view.
|
||
const Type& get_reference() const
|
||
{ return view().get_reference(); }
|
||
|
||
/// Return a const reference to index data member in the view.
|
||
const Index& get_index_reference() const
|
||
{ return view().get_index_reference(); }
|
||
|
||
//@}
|
||
|
||
/** @name Forwarded functions
|
||
* @details Functions that update the contained accumulator variable are
|
||
* simply forwarded to the contained @ref op_min_view. */
|
||
//@{
|
||
|
||
/// @copydoc cilk_lib_1_0::min_max_internal::view_content::is_set() const
|
||
bool is_set() const { return view().is_set(); }
|
||
|
||
/// @copydoc op_min_index_view::calc_min(const Index&, const Type&)
|
||
reducer_min_index& calc_min(const Index& i, const Type& x)
|
||
{ view().calc_min(i, x); return *this; }
|
||
|
||
/// @copydoc op_min_view::operator=(const min_max_internal::rhs_proxy<op_min_view>&)
|
||
reducer_min_index& operator=(const rhs_proxy& rhs)
|
||
{ view() = rhs; return *this; }
|
||
|
||
//@}
|
||
|
||
/// @name Dereference
|
||
/** Dereferencing a wrapper is a no-op. It simply returns the wrapper.
|
||
* Combined with the rule that a wrapper forwards view operations to the
|
||
* view, this means that view operations can be written the same way on
|
||
* reducers and wrappers, which is convenient for incrementally
|
||
* converting code using wrappers to code using reducers. That is:
|
||
*
|
||
* reducer< op_min_index<int, int> > r;
|
||
* r->calc_min(i, a); // *r returns the view
|
||
* // calc_min is a view member function
|
||
*
|
||
* reducer_min_index<int, int> w;
|
||
* w->calc_min(i, a); // *w returns the wrapper
|
||
* // calc_min is a wrapper member function that
|
||
* // calls the corresponding view function
|
||
*/
|
||
//@{
|
||
reducer_min_index& operator*() { return *this; }
|
||
reducer_min_index const& operator*() const { return *this; }
|
||
|
||
reducer_min_index* operator->() { return this; }
|
||
reducer_min_index const* operator->() const { return this; }
|
||
//@}
|
||
|
||
/** @name Upcast
|
||
* @details In Cilk library 0.9, reducers were always cache-aligned. In
|
||
* library 1.0, reducer cache alignment is optional. By default, reducers
|
||
* are unaligned (i.e., just naturally aligned), but legacy wrappers
|
||
* inherit from cache-aligned reducers for binary compatibility.
|
||
*
|
||
* This means that a wrapper will automatically be upcast to its aligned
|
||
* reducer base class. The following conversion operators provide
|
||
* pseudo-upcasts to the corresponding unaligned reducer class.
|
||
*/
|
||
//@{
|
||
operator reducer< op_min_index<Index, Type, Compare, false> >& ()
|
||
{
|
||
return *reinterpret_cast< reducer< op_min_index<Index, Type, Compare, false> >* >(this);
|
||
}
|
||
|
||
operator const reducer< op_min_index<Index, Type, Compare, false> >& () const
|
||
{
|
||
return *reinterpret_cast< const reducer< op_min_index<Index, Type, Compare, false> >* >(this);
|
||
}
|
||
//@}
|
||
|
||
};
|
||
|
||
|
||
#ifndef CILK_LIBRARY_0_9_REDUCER_MINMAX
|
||
} // namespace cilk_lib_1_0
|
||
using namespace cilk_lib_1_0;
|
||
#endif
|
||
|
||
|
||
/// @cond internal
|
||
/** Metafunction specialization for reducer conversion.
|
||
*
|
||
* These specializations of the @ref legacy_reducer_downcast template class
|
||
* defined in reducer.h causes each `reducer< op_xxxx<Type> >` classes to have
|
||
* an `operator reducer_xxxx<Type>& ()` conversion operator that statically
|
||
* downcasts the `reducer<op_xxxx>` to the corresponding `reducer_xxxx` type.
|
||
* (The reverse conversion, from `reducer_xxxx` to `reducer<op_xxxx>`, is just
|
||
* an upcast, which is provided for free by the language.)
|
||
*/
|
||
template <typename Type, typename Compare, bool Align>
|
||
struct legacy_reducer_downcast< reducer< op_max<Type, Compare, Align> > >
|
||
{
|
||
typedef reducer_max<Type> type;
|
||
};
|
||
|
||
template <typename Type, typename Compare, bool Align>
|
||
struct legacy_reducer_downcast< reducer< op_min<Type, Compare, Align> > >
|
||
{
|
||
typedef reducer_min<Type> type;
|
||
};
|
||
|
||
template <typename Index, typename Type, typename Compare, bool Align>
|
||
struct legacy_reducer_downcast< reducer< op_max_index<Index, Type, Compare, Align> > >
|
||
{
|
||
typedef reducer_max_index<Index, Type> type;
|
||
};
|
||
|
||
template <typename Index, typename Type, typename Compare, bool Align>
|
||
struct legacy_reducer_downcast< reducer< op_min_index<Index, Type, Compare, Align> > >
|
||
{
|
||
typedef reducer_min_index<Index, Type> type;
|
||
};
|
||
/// @endcond
|
||
|
||
} // namespace cilk
|
||
|
||
#endif // __cplusplus
|
||
|
||
|
||
/** @name C language reducer macros
|
||
*
|
||
* These macros are used to declare and work with numeric minimum and maximum reducers in C
|
||
* code.
|
||
*
|
||
* @see @ref page_reducers_in_c
|
||
*/
|
||
//@{
|
||
|
||
|
||
#ifdef CILK_C_DEFINE_REDUCERS
|
||
|
||
/* Integer min/max constants */
|
||
#include <limits.h>
|
||
|
||
/* Wchar_t min/max constants */
|
||
#if defined(_MSC_VER) || defined(__ANDROID__)
|
||
# include <wchar.h>
|
||
#else
|
||
# include <stdint.h>
|
||
#endif
|
||
|
||
/* Floating-point min/max constants */
|
||
#include <math.h>
|
||
#ifndef HUGE_VALF
|
||
static const unsigned int __huge_valf[] = {0x7f800000};
|
||
# define HUGE_VALF (*((const float *)__huge_valf))
|
||
#endif
|
||
|
||
#ifndef HUGE_VALL
|
||
static const unsigned int __huge_vall[] = {0, 0, 0x00007f80, 0};
|
||
# define HUGE_VALL (*((const long double *)__huge_vall))
|
||
#endif
|
||
|
||
#endif
|
||
|
||
/** Max reducer type name.
|
||
*
|
||
* This macro expands into the identifier which is the name of the max reducer
|
||
* type for a specified numeric type.
|
||
*
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_TYPE(tn) \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_,tn)
|
||
|
||
/** Declare a max reducer object.
|
||
*
|
||
* This macro expands into a declaration of a max reducer object for a specified numeric
|
||
* type. For example:
|
||
*
|
||
* CILK_C_REDUCER_MAX(my_reducer, double, -DBL_MAX);
|
||
*
|
||
* @param obj The variable name to be used for the declared reducer object.
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
* @param v The initial value for the reducer. (A value which can be assigned to the
|
||
* numeric type represented by @a tn.)
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MAX(obj,tn,v) \
|
||
CILK_C_REDUCER_MAX_TYPE(tn) obj = \
|
||
CILK_C_INIT_REDUCER(_Typeof(obj.value), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_reduce_,tn), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_identity_,tn), \
|
||
__cilkrts_hyperobject_noop_destroy, v)
|
||
|
||
/** Maximize with a value.
|
||
*
|
||
* `CILK_C_REDUCER_MAX_CALC(reducer, v)` sets the current view of the
|
||
* reducer to the max of its previous value and a specified new value.
|
||
* This is equivalent to
|
||
*
|
||
* REDUCER_VIEW(reducer) = max(REDUCER_VIEW(reducer), v)
|
||
*
|
||
* @param reducer The reducer whose contained value is to be updated.
|
||
* @param v The value that it is to be maximized with.
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_CALC(reducer, v) do { \
|
||
_Typeof((reducer).value)* view = &(REDUCER_VIEW(reducer)); \
|
||
_Typeof(v) __value = (v); \
|
||
if (*view < __value) { \
|
||
*view = __value; \
|
||
} } while (0)
|
||
|
||
/// @cond internal
|
||
|
||
/** Declare the max reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into external function declarations for functions which implement
|
||
* the reducer functionality for the max reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_DECLARATION(t,tn,id) \
|
||
typedef CILK_C_DECLARE_REDUCER(t) CILK_C_REDUCER_MAX_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_max,tn,l,r); \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_max,tn);
|
||
|
||
/** Define the max reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into function definitions for functions which implement the
|
||
* reducer functionality for the max reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_DEFINITION(t,tn,id) \
|
||
typedef CILK_C_DECLARE_REDUCER(t) CILK_C_REDUCER_MAX_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_max,tn,l,r) \
|
||
{ if (*(t*)l < *(t*)r) *(t*)l = *(t*)r; } \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_max,tn) \
|
||
{ *(t*)v = id; }
|
||
|
||
//@{
|
||
/** @def CILK_C_REDUCER_MAX_INSTANCE
|
||
* @brief Declare or define implementation functions for a reducer type.
|
||
*
|
||
* In the runtime source file c_reducers.c, the macro `CILK_C_DEFINE_REDUCERS` will be defined, and
|
||
* this macro will generate reducer implementation functions. Everywhere else, `CILK_C_DEFINE_REDUCERS`
|
||
* will be undefined, and this macro will expand into external declarations for the functions.
|
||
*/
|
||
#ifdef CILK_C_DEFINE_REDUCERS
|
||
# define CILK_C_REDUCER_MAX_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MAX_DEFINITION(t,tn,id)
|
||
#else
|
||
# define CILK_C_REDUCER_MAX_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MAX_DECLARATION(t,tn,id)
|
||
#endif
|
||
//@}
|
||
|
||
/* Declare or define an instance of the reducer type and its functions for each
|
||
* numeric type.
|
||
*/
|
||
__CILKRTS_BEGIN_EXTERN_C
|
||
CILK_C_REDUCER_MAX_INSTANCE(char, char, CHAR_MIN)
|
||
CILK_C_REDUCER_MAX_INSTANCE(unsigned char, uchar, 0)
|
||
CILK_C_REDUCER_MAX_INSTANCE(signed char, schar, SCHAR_MIN)
|
||
CILK_C_REDUCER_MAX_INSTANCE(wchar_t, wchar_t, WCHAR_MIN)
|
||
CILK_C_REDUCER_MAX_INSTANCE(short, short, SHRT_MIN)
|
||
CILK_C_REDUCER_MAX_INSTANCE(unsigned short, ushort, 0)
|
||
CILK_C_REDUCER_MAX_INSTANCE(int, int, INT_MIN)
|
||
CILK_C_REDUCER_MAX_INSTANCE(unsigned int, uint, 0)
|
||
CILK_C_REDUCER_MAX_INSTANCE(unsigned int, unsigned, 0) // alternate name
|
||
CILK_C_REDUCER_MAX_INSTANCE(long, long, LONG_MIN)
|
||
CILK_C_REDUCER_MAX_INSTANCE(unsigned long, ulong, 0)
|
||
CILK_C_REDUCER_MAX_INSTANCE(long long, longlong, LLONG_MIN)
|
||
CILK_C_REDUCER_MAX_INSTANCE(unsigned long long, ulonglong, 0)
|
||
CILK_C_REDUCER_MAX_INSTANCE(float, float, -HUGE_VALF)
|
||
CILK_C_REDUCER_MAX_INSTANCE(double, double, -HUGE_VAL)
|
||
CILK_C_REDUCER_MAX_INSTANCE(long double, longdouble, -HUGE_VALL)
|
||
__CILKRTS_END_EXTERN_C
|
||
|
||
/// @endcond
|
||
|
||
/** Max_index reducer type name.
|
||
*
|
||
* This macro expands into the identifier which is the name of the max_index reducer
|
||
* type for a specified numeric type.
|
||
*
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_INDEX_TYPE(tn) \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_index_,tn)
|
||
|
||
/** Declare an op_max_index reducer object.
|
||
*
|
||
* This macro expands into a declaration of a max_index reducer object for a specified
|
||
* numeric type. For example:
|
||
*
|
||
* CILK_C_REDUCER_MAX_INDEX(my_reducer, double, -DBL_MAX_INDEX);
|
||
*
|
||
* @param obj The variable name to be used for the declared reducer object.
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
* @param v The initial value for the reducer. (A value which can be assigned to the
|
||
* numeric type represented by @a tn.)
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_INDEX(obj,tn,v) \
|
||
CILK_C_REDUCER_MAX_INDEX_TYPE(tn) obj = \
|
||
CILK_C_INIT_REDUCER(_Typeof(obj.value), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_index_reduce_,tn), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_index_identity_,tn), \
|
||
__cilkrts_hyperobject_noop_destroy, {0, v})
|
||
|
||
/** Maximize with a value.
|
||
*
|
||
* `CILK_C_REDUCER_MAX_INDEX_CALC(reducer, i, v)` sets the current view of the
|
||
* reducer to the max of its previous value and a specified new value.
|
||
* This is equivalent to
|
||
*
|
||
* REDUCER_VIEW(reducer) = max_index(REDUCER_VIEW(reducer), v)
|
||
*
|
||
* If the value of the reducer is changed to @a v, then the index of the reducer is
|
||
* changed to @a i.
|
||
*
|
||
* @param reducer The reducer whose contained value and index are to be updated.
|
||
* @param i The index associated with the new value.
|
||
* @param v The value that it is to be maximized with.
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_INDEX_CALC(reducer, i, v) do { \
|
||
_Typeof((reducer).value)* view = &(REDUCER_VIEW(reducer)); \
|
||
_Typeof(v) __value = (v); \
|
||
if (view->value < __value) { \
|
||
view->index = (i); \
|
||
view->value = __value; \
|
||
} } while (0)
|
||
|
||
/// @cond internal
|
||
|
||
/** Declare the max_index view type.
|
||
*
|
||
* The view of a max_index reducer is a structure containing both the
|
||
* maximum value for the reducer and the index that was associated with
|
||
* that value in the sequence of input values.
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_INDEX_VIEW(t,tn) \
|
||
typedef struct { \
|
||
__STDNS ptrdiff_t index; \
|
||
t value; \
|
||
} __CILKRTS_MKIDENT(cilk_c_reducer_max_index_view_,tn)
|
||
|
||
/** Declare the max_index reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into external function declarations for functions which implement
|
||
* the reducer functionality for the max_index reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_INDEX_DECLARATION(t,tn,id) \
|
||
CILK_C_REDUCER_MAX_INDEX_VIEW(t,tn); \
|
||
typedef CILK_C_DECLARE_REDUCER( \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_index_view_,tn)) \
|
||
CILK_C_REDUCER_MAX_INDEX_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_max_index,tn,l,r); \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_max_index,tn);
|
||
|
||
/** Define the max_index reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into function definitions for functions which implement the
|
||
* reducer functionality for the max_index reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MAX_INDEX_DEFINITION(t,tn,id) \
|
||
CILK_C_REDUCER_MAX_INDEX_VIEW(t,tn); \
|
||
typedef CILK_C_DECLARE_REDUCER( \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_max_index_view_,tn)) \
|
||
CILK_C_REDUCER_MAX_INDEX_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_max_index,tn,l,r) \
|
||
{ typedef __CILKRTS_MKIDENT(cilk_c_reducer_max_index_view_,tn) view_t; \
|
||
if (((view_t*)l)->value < ((view_t*)r)->value) \
|
||
*(view_t*)l = *(view_t*)r; } \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_max_index,tn) \
|
||
{ typedef __CILKRTS_MKIDENT(cilk_c_reducer_max_index_view_,tn) view_t; \
|
||
((view_t*)v)->index = 0; ((view_t*)v)->value = id; }
|
||
|
||
//@{
|
||
/** @def CILK_C_REDUCER_MAX_INDEX_INSTANCE
|
||
* @brief Declare or define implementation functions for a reducer type.
|
||
*
|
||
* In the runtime source file c_reducers.c, the macro `CILK_C_DEFINE_REDUCERS` will be defined, and
|
||
* this macro will generate reducer implementation functions. Everywhere else, `CILK_C_DEFINE_REDUCERS`
|
||
* will be undefined, and this macro will expand into external declarations for the functions.
|
||
*/
|
||
#ifdef CILK_C_DEFINE_REDUCERS
|
||
# define CILK_C_REDUCER_MAX_INDEX_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MAX_INDEX_DEFINITION(t,tn,id)
|
||
#else
|
||
# define CILK_C_REDUCER_MAX_INDEX_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MAX_INDEX_DECLARATION(t,tn,id)
|
||
#endif
|
||
//@}
|
||
|
||
/* Declare or define an instance of the reducer type and its functions for each
|
||
* numeric type.
|
||
*/
|
||
__CILKRTS_BEGIN_EXTERN_C
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(char, char, CHAR_MIN)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(unsigned char, uchar, 0)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(signed char, schar, SCHAR_MIN)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(wchar_t, wchar_t, WCHAR_MIN)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(short, short, SHRT_MIN)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(unsigned short, ushort, 0)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(int, int, INT_MIN)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(unsigned int, uint, 0)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(unsigned int, unsigned, 0) // alternate name
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(long, long, LONG_MIN)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(unsigned long, ulong, 0)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(long long, longlong, LLONG_MIN)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(unsigned long long, ulonglong, 0)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(float, float, -HUGE_VALF)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(double, double, -HUGE_VAL)
|
||
CILK_C_REDUCER_MAX_INDEX_INSTANCE(long double, longdouble, -HUGE_VALL)
|
||
__CILKRTS_END_EXTERN_C
|
||
|
||
/// @endcond
|
||
|
||
/** Min reducer type name.
|
||
*
|
||
* This macro expands into the identifier which is the name of the min reducer
|
||
* type for a specified numeric type.
|
||
*
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_TYPE(tn) \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_,tn)
|
||
|
||
/** Declare a min reducer object.
|
||
*
|
||
* This macro expands into a declaration of a min reducer object for a specified numeric
|
||
* type. For example:
|
||
*
|
||
* CILK_C_REDUCER_MIN(my_reducer, double, DBL_MAX);
|
||
*
|
||
* @param obj The variable name to be used for the declared reducer object.
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
* @param v The initial value for the reducer. (A value which can be assigned to the
|
||
* numeric type represented by @a tn.)
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MIN(obj,tn,v) \
|
||
CILK_C_REDUCER_MIN_TYPE(tn) obj = \
|
||
CILK_C_INIT_REDUCER(_Typeof(obj.value), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_reduce_,tn), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_identity_,tn), \
|
||
__cilkrts_hyperobject_noop_destroy, v)
|
||
|
||
/** Minimize with a value.
|
||
*
|
||
* `CILK_C_REDUCER_MIN_CALC(reducer, v)` sets the current view of the
|
||
* reducer to the min of its previous value and a specified new value.
|
||
* This is equivalent to
|
||
*
|
||
* REDUCER_VIEW(reducer) = min(REDUCER_VIEW(reducer), v)
|
||
*
|
||
* @param reducer The reducer whose contained value is to be updated.
|
||
* @param v The value that it is to be minimized with.
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_CALC(reducer, v) do { \
|
||
_Typeof((reducer).value)* view = &(REDUCER_VIEW(reducer)); \
|
||
_Typeof(v) __value = (v); \
|
||
if (*view > __value) { \
|
||
*view = __value; \
|
||
} } while (0)
|
||
|
||
/// @cond internal
|
||
|
||
/** Declare the min reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into external function declarations for functions which implement
|
||
* the reducer functionality for the min reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_DECLARATION(t,tn,id) \
|
||
typedef CILK_C_DECLARE_REDUCER(t) CILK_C_REDUCER_MIN_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_min,tn,l,r); \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_min,tn);
|
||
|
||
/** Define the min reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into function definitions for functions which implement the
|
||
* reducer functionality for the min reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_DEFINITION(t,tn,id) \
|
||
typedef CILK_C_DECLARE_REDUCER(t) CILK_C_REDUCER_MIN_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_min,tn,l,r) \
|
||
{ if (*(t*)l > *(t*)r) *(t*)l = *(t*)r; } \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_min,tn) \
|
||
{ *(t*)v = id; }
|
||
|
||
//@{
|
||
/** @def CILK_C_REDUCER_MIN_INSTANCE
|
||
* @brief Declare or define implementation functions for a reducer type.
|
||
*
|
||
* In the runtime source file c_reducers.c, the macro `CILK_C_DEFINE_REDUCERS` will be defined, and
|
||
* this macro will generate reducer implementation functions. Everywhere else, `CILK_C_DEFINE_REDUCERS`
|
||
* will be undefined, and this macro will expand into external declarations for the functions.
|
||
*/
|
||
#ifdef CILK_C_DEFINE_REDUCERS
|
||
# define CILK_C_REDUCER_MIN_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MIN_DEFINITION(t,tn,id)
|
||
#else
|
||
# define CILK_C_REDUCER_MIN_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MIN_DECLARATION(t,tn,id)
|
||
#endif
|
||
//@}
|
||
|
||
/* Declare or define an instance of the reducer type and its functions for each
|
||
* numeric type.
|
||
*/
|
||
__CILKRTS_BEGIN_EXTERN_C
|
||
CILK_C_REDUCER_MIN_INSTANCE(char, char, CHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(unsigned char, uchar, CHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(signed char, schar, SCHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(wchar_t, wchar_t, WCHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(short, short, SHRT_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(unsigned short, ushort, USHRT_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(int, int, INT_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(unsigned int, uint, UINT_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(unsigned int, unsigned, UINT_MAX) // alternate name
|
||
CILK_C_REDUCER_MIN_INSTANCE(long, long, LONG_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(unsigned long, ulong, ULONG_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(long long, longlong, LLONG_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(unsigned long long, ulonglong, ULLONG_MAX)
|
||
CILK_C_REDUCER_MIN_INSTANCE(float, float, HUGE_VALF)
|
||
CILK_C_REDUCER_MIN_INSTANCE(double, double, HUGE_VAL)
|
||
CILK_C_REDUCER_MIN_INSTANCE(long double, longdouble, HUGE_VALL)
|
||
__CILKRTS_END_EXTERN_C
|
||
|
||
/// @endcond
|
||
|
||
/** Min_index reducer type name.
|
||
*
|
||
* This macro expands into the identifier which is the name of the min_index reducer
|
||
* type for a specified numeric type.
|
||
*
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_INDEX_TYPE(tn) \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_index_,tn)
|
||
|
||
/** Declare an op_min_index reducer object.
|
||
*
|
||
* This macro expands into a declaration of a min_index reducer object for a specified
|
||
* numeric type. For example:
|
||
*
|
||
* CILK_C_REDUCER_MIN_INDEX(my_reducer, double, -DBL_MIN_INDEX);
|
||
*
|
||
* @param obj The variable name to be used for the declared reducer object.
|
||
* @param tn The @ref reducers_c_type_names "numeric type name" specifying the type of the
|
||
* reducer.
|
||
* @param v The initial value for the reducer. (A value which can be assigned to the
|
||
* numeric type represented by @a tn.)
|
||
*
|
||
* @see @ref reducers_c_predefined
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_INDEX(obj,tn,v) \
|
||
CILK_C_REDUCER_MIN_INDEX_TYPE(tn) obj = \
|
||
CILK_C_INIT_REDUCER(_Typeof(obj.value), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_index_reduce_,tn), \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_index_identity_,tn), \
|
||
__cilkrts_hyperobject_noop_destroy, {0, v})
|
||
|
||
/** Minimize with a value.
|
||
*
|
||
* `CILK_C_REDUCER_MIN_INDEX_CALC(reducer, i, v)` sets the current view of the
|
||
* reducer to the min of its previous value and a specified new value.
|
||
* This is equivalent to
|
||
*
|
||
* REDUCER_VIEW(reducer) = min_index(REDUCER_VIEW(reducer), v)
|
||
*
|
||
* If the value of the reducer is changed to @a v, then the index of the reducer is
|
||
* changed to @a i.
|
||
*
|
||
* @param reducer The reducer whose contained value and index are to be updated.
|
||
* @param i The index associated with the new value.
|
||
* @param v The value that it is to be minimized with.
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_INDEX_CALC(reducer, i, v) do { \
|
||
_Typeof((reducer).value)* view = &(REDUCER_VIEW(reducer)); \
|
||
_Typeof(v) __value = (v); \
|
||
if (view->value > __value) { \
|
||
view->index = (i); \
|
||
view->value = __value; \
|
||
} } while (0)
|
||
|
||
/// @cond internal
|
||
|
||
/** Declare the min_index view type.
|
||
*
|
||
* The view of a min_index reducer is a structure containing both the
|
||
* minimum value for the reducer and the index that was associated with
|
||
* that value in the sequence of input values.
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_INDEX_VIEW(t,tn) \
|
||
typedef struct { \
|
||
__STDNS ptrdiff_t index; \
|
||
t value; \
|
||
} __CILKRTS_MKIDENT(cilk_c_reducer_min_index_view_,tn)
|
||
|
||
/** Declare the min_index reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into external function declarations for functions which implement
|
||
* the reducer functionality for the min_index reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_INDEX_DECLARATION(t,tn,id) \
|
||
CILK_C_REDUCER_MIN_INDEX_VIEW(t,tn); \
|
||
typedef CILK_C_DECLARE_REDUCER( \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_index_view_,tn)) \
|
||
CILK_C_REDUCER_MIN_INDEX_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_min_index,tn,l,r); \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_min_index,tn);
|
||
|
||
/** Define the min_index reducer functions for a numeric type.
|
||
*
|
||
* This macro expands into function definitions for functions which implement the
|
||
* reducer functionality for the min_index reducer type for a specified numeric type.
|
||
*
|
||
* @param t The value type of the reducer.
|
||
* @param tn The value “type name” identifier, used to construct the reducer type name,
|
||
* function names, etc.
|
||
*/
|
||
#define CILK_C_REDUCER_MIN_INDEX_DEFINITION(t,tn,id) \
|
||
CILK_C_REDUCER_MIN_INDEX_VIEW(t,tn); \
|
||
typedef CILK_C_DECLARE_REDUCER( \
|
||
__CILKRTS_MKIDENT(cilk_c_reducer_min_index_view_,tn)) \
|
||
CILK_C_REDUCER_MIN_INDEX_TYPE(tn); \
|
||
__CILKRTS_DECLARE_REDUCER_REDUCE(cilk_c_reducer_min_index,tn,l,r) \
|
||
{ typedef __CILKRTS_MKIDENT(cilk_c_reducer_min_index_view_,tn) view_t; \
|
||
if (((view_t*)l)->value > ((view_t*)r)->value) \
|
||
*(view_t*)l = *(view_t*)r; } \
|
||
__CILKRTS_DECLARE_REDUCER_IDENTITY(cilk_c_reducer_min_index,tn) \
|
||
{ typedef __CILKRTS_MKIDENT(cilk_c_reducer_min_index_view_,tn) view_t; \
|
||
((view_t*)v)->index = 0; ((view_t*)v)->value = id; }
|
||
|
||
//@{
|
||
/** @def CILK_C_REDUCER_MIN_INDEX_INSTANCE
|
||
* @brief Declare or define implementation functions for a reducer type.
|
||
*
|
||
* In the runtime source file c_reducers.c, the macro `CILK_C_DEFINE_REDUCERS` will be defined, and
|
||
* this macro will generate reducer implementation functions. Everywhere else, `CILK_C_DEFINE_REDUCERS`
|
||
* will be undefined, and this macro will expand into external declarations for the functions.
|
||
*/
|
||
#ifdef CILK_C_DEFINE_REDUCERS
|
||
# define CILK_C_REDUCER_MIN_INDEX_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MIN_INDEX_DEFINITION(t,tn,id)
|
||
#else
|
||
# define CILK_C_REDUCER_MIN_INDEX_INSTANCE(t,tn,id) \
|
||
CILK_C_REDUCER_MIN_INDEX_DECLARATION(t,tn,id)
|
||
#endif
|
||
//@}
|
||
|
||
/* Declare or define an instance of the reducer type and its functions for each
|
||
* numeric type.
|
||
*/
|
||
__CILKRTS_BEGIN_EXTERN_C
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(char, char, CHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(unsigned char, uchar, CHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(signed char, schar, SCHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(wchar_t, wchar_t, WCHAR_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(short, short, SHRT_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(unsigned short, ushort, USHRT_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(int, int, INT_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(unsigned int, uint, UINT_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(unsigned int, unsigned, UINT_MAX) // alternate name
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(long, long, LONG_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(unsigned long, ulong, ULONG_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(long long, longlong, LLONG_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(unsigned long long, ulonglong, ULLONG_MAX)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(float, float, HUGE_VALF)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(double, double, HUGE_VAL)
|
||
CILK_C_REDUCER_MIN_INDEX_INSTANCE(long double, longdouble, HUGE_VALL)
|
||
__CILKRTS_END_EXTERN_C
|
||
|
||
/// @endcond
|
||
|
||
//@}
|
||
|
||
#endif // defined REDUCER_MAX_H_INCLUDED
|