RTOp Package Browser (Single Doxygen Collection) Version of the Day
RTOpPack_ROpMaxIndexLessThanBound.hpp
Go to the documentation of this file.
00001 // @HEADER
00002 // ***********************************************************************
00003 // 
00004 // RTOp: Interfaces and Support Software for Vector Reduction Transformation
00005 //       Operations
00006 //                Copyright (2006) Sandia Corporation
00007 // 
00008 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00009 // license for use of this work by or on behalf of the U.S. Government.
00010 // 
00011 // This library is free software; you can redistribute it and/or modify
00012 // it under the terms of the GNU Lesser General Public License as
00013 // published by the Free Software Foundation; either version 2.1 of the
00014 // License, or (at your option) any later version.
00015 //  
00016 // This library is distributed in the hope that it will be useful, but
00017 // WITHOUT ANY WARRANTY; without even the implied warranty of
00018 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00019 // Lesser General Public License for more details.
00020 //  
00021 // You should have received a copy of the GNU Lesser General Public
00022 // License along with this library; if not, write to the Free Software
00023 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00024 // USA
00025 // Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov) 
00026 // 
00027 // ***********************************************************************
00028 // @HEADER
00029 
00030 #ifndef RTOPPACK_ROP_MAX_INDEX_LESS_THAN_BOUND_HPP
00031 #define RTOPPACK_ROP_MAX_INDEX_LESS_THAN_BOUND_HPP
00032 
00033 #include "RTOpPack_ROpMaxIndex.hpp"
00034 #include "Teuchos_StandardMemberCompositionMacros.hpp"
00035 
00036 
00037 namespace RTOpPack {
00038 
00039 
00041 template<class Scalar>
00042 class ROpMaxIndexLessThanBoundEleWiseReductionOp {
00043 public:
00045   ROpMaxIndexLessThanBoundEleWiseReductionOp(
00046     const Scalar &bound = ScalarTraits<Scalar>::zero()
00047     )
00048     :bound_(bound)
00049     {}
00051   void operator()(const index_type i, const Scalar &v0,
00052     ScalarIndex<Scalar> &reduct) const
00053     {
00054       if(
00055         v0 <  bound_
00056         &&
00057         (
00058           v0 > reduct.scalar
00059           ||
00060           ( v0 == reduct.scalar && i < reduct.index )
00061           )
00062         )
00063       {
00064         reduct = ScalarIndex<Scalar>(v0, i);
00065       }
00066     }
00067 private:
00068   Scalar bound_;
00069 };
00070 
00071 
00081 template<class Scalar>
00082 class ROpMaxIndexLessThanBound
00083   : public ROp_1_CoordVariantScalarReduction<
00084       Scalar,
00085       ScalarIndex<Scalar>,
00086       ROpMaxIndexLessThanBoundEleWiseReductionOp<Scalar>,
00087       ROpMaxIndexReductObjReductionOp<Scalar> >
00088 {
00089 public:
00091   ROpMaxIndexLessThanBound(
00092     const Scalar &bound_in = Teuchos::ScalarTraits<Scalar>::zero()
00093     )
00094     {
00095       this->setOpNameBase("ROpMaxIndexLessThanBound");
00096       bound(bound_in);
00097       this->initReductObjValue(
00098         ScalarIndex<Scalar>(-ScalarTraits<Scalar>::rmax(), -1));
00099     }
00101   void bound(const Scalar& bound_in)
00102     { 
00103       this->setEleWiseReduction(
00104         ROpMaxIndexLessThanBoundEleWiseReductionOp<Scalar>(bound_in)
00105         );
00106     }
00108   ScalarIndex<Scalar> operator()(const ReductTarget& reduct_obj ) const
00109     { return this->getRawVal(reduct_obj); }
00110 };
00111 
00112 
00113 } // namespace RTOpPack
00114 
00115 
00116 #endif // RTOPPACK_ROP_MAX_INDEX_LESS_THAN_BOUND_HPP
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines