RTOp Package Browser (Single Doxygen Collection)  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RTOpPack_ROpMaxIndex.hpp
Go to the documentation of this file.
1 // @HEADER
2 // ***********************************************************************
3 //
4 // RTOp: Interfaces and Support Software for Vector Reduction Transformation
5 // Operations
6 // Copyright (2006) Sandia Corporation
7 //
8 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
9 // license for use of this work by or on behalf of the U.S. Government.
10 //
11 // Redistribution and use in source and binary forms, with or without
12 // modification, are permitted provided that the following conditions are
13 // met:
14 //
15 // 1. Redistributions of source code must retain the above copyright
16 // notice, this list of conditions and the following disclaimer.
17 //
18 // 2. Redistributions in binary form must reproduce the above copyright
19 // notice, this list of conditions and the following disclaimer in the
20 // documentation and/or other materials provided with the distribution.
21 //
22 // 3. Neither the name of the Corporation nor the names of the
23 // contributors may be used to endorse or promote products derived from
24 // this software without specific prior written permission.
25 //
26 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37 //
38 // Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov)
39 //
40 // ***********************************************************************
41 // @HEADER
42 
43 #ifndef RTOPPACK_ROP_MAX_INDEX_HPP
44 #define RTOPPACK_ROP_MAX_INDEX_HPP
45 
46 #include "RTOpPack_RTOpTHelpers.hpp"
47 
48 namespace RTOpPack {
49 
50 
52 template<class Scalar>
54 public:
56  void operator()(const index_type i, const Scalar &v0,
57  ScalarIndex<Scalar> &reduct) const
58  {
59  if(
60  v0 > reduct.scalar
61  ||
62  (v0 == reduct.scalar && i < reduct.index)
63  )
64  {
65  reduct = ScalarIndex<Scalar>(v0, i);
66  }
67  }
68 };
69 
70 
72 template<class Scalar>
74 public:
76  void operator()(
77  const ScalarIndex<Scalar>& in_reduct, ScalarIndex<Scalar>& inout_reduct
78  ) const
79  {
80  if(
81  in_reduct.scalar > inout_reduct.scalar
82  ||
83  (
84  in_reduct.scalar == inout_reduct.scalar
85  &&
86  in_reduct.index < inout_reduct.index
87  )
88  )
89  {
90  inout_reduct = in_reduct;
91  }
92  }
93 };
94 
95 
101 template<class Scalar>
104  Scalar,
105  ScalarIndex<Scalar>,
106  ROpMaxIndexEleWiseReductionOp<Scalar>,
107  ROpMaxIndexReductObjReductionOp<Scalar> >
108 {
109 public:
112  {
113  this->setOpNameBase("ROpMaxIndex");
114  this->initReductObjValue(
116  }
118  ScalarIndex<Scalar> operator()(const ReductTarget& reduct_obj) const
119  { return this->getRawVal(reduct_obj); }
120 };
121 
122 
123 } // namespace RTOpPack
124 
125 
126 #endif // RTOPPACK_ROP_MAX_INDEX_HPP
Base class for coordinate-variant scalar reduction RTOps with one input vector.
Teuchos_Ordinal index_type
Simple struct for a Scalar and an Ordinal object.
Abstract base class for all reduction objects.
Returns the maximum element and its index: result.scalar = x(k) and result.index = k such that x(k) &gt;...
void setOpNameBase(const std::string &op_name_base)
Just set the operator name.
void operator()(const index_type i, const Scalar &v0, ScalarIndex< Scalar > &reduct) const
ScalarIndex< Scalar > operator()(const ReductTarget &reduct_obj) const
void operator()(const ScalarIndex< Scalar > &in_reduct, ScalarIndex< Scalar > &inout_reduct) const