RTOpPack_ROpMinIndexGreaterThanBound.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_MIN_INDEX_GREATER_THAN_BOUND_HPP
00031 #define RTOPPACK_ROP_MIN_INDEX_GREATER_THAN_BOUND_HPP
00032 
00033 #include "RTOpPack_ROpMinIndex.hpp"
00034 #include "Teuchos_StandardMemberCompositionMacros.hpp"
00035 
00036 
00037 namespace RTOpPack {
00038 
00039 
00041 template<class Scalar>
00042 class ROpMinIndexGreaterThanBoundEleWiseReductionOp {
00043 public:
00045   ROpMinIndexGreaterThanBoundEleWiseReductionOp(
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 ROpMinIndexGreaterThanBound
00083   : public ROp_1_CoordVariantScalarReduction<
00084       Scalar,
00085       ScalarIndex<Scalar>,
00086       ROpMinIndexGreaterThanBoundEleWiseReductionOp<Scalar>,
00087       ROpMinIndexReductObjReductionOp<Scalar> >
00088 {
00089 public:
00090 
00092   ROpMinIndexGreaterThanBound(
00093     const Scalar &bound_in = Teuchos::ScalarTraits<Scalar>::zero()
00094     )
00095     {
00096       this->setOpNameBase("ROpMinIndexGreaterThanBound");
00097       bound(bound_in);
00098       initReductObjValue(
00099         ScalarIndex<Scalar>(+ScalarTraits<Scalar>::rmax(), -1));
00100     }
00102   void bound(const Scalar& bound_in)
00103     { 
00104       this->setEleWiseReduction(
00105         ROpMinIndexGreaterThanBoundEleWiseReductionOp<Scalar>(bound_in)
00106         );
00107     }
00109   ScalarIndex<Scalar> operator()(const ReductTarget& reduct_obj ) const
00110     { return this->getRawVal(reduct_obj); }
00111   
00112 };
00113 
00114 
00115 } // namespace RTOpPack
00116 
00117 
00118 #endif // RTOPPACK_ROP_MIN_INDEX_GREATER_THAN_BOUND_HPP
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines
Generated on Wed Apr 13 09:59:06 2011 for RTOp Package Browser (Single Doxygen Collection) by  doxygen 1.6.3