ConstrainedOptPack_ComputeMinMult.cpp

Go to the documentation of this file.
00001 // @HEADER
00002 // ***********************************************************************
00003 // 
00004 // Moocho: Multi-functional Object-Oriented arCHitecture for Optimization
00005 //                  Copyright (2003) Sandia Corporation
00006 // 
00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00008 // license for use of this work by or on behalf of the U.S. Government.
00009 // 
00010 // This library is free software; you can redistribute it and/or modify
00011 // it under the terms of the GNU Lesser General Public License as
00012 // published by the Free Software Foundation; either version 2.1 of the
00013 // License, or (at your option) any later version.
00014 //  
00015 // This library is distributed in the hope that it will be useful, but
00016 // WITHOUT ANY WARRANTY; without even the implied warranty of
00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00018 // Lesser General Public License for more details.
00019 //  
00020 // You should have received a copy of the GNU Lesser General Public
00021 // License along with this library; if not, write to the Free Software
00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00023 // USA
00024 // Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov) 
00025 // 
00026 // ***********************************************************************
00027 // @HEADER
00028 
00029 #include <limits>
00030 
00031 #include "ConstrainedOptPack_ComputeMinMult.hpp"
00032 #include "AbstractLinAlgPack_SpVectorClass.hpp"
00033 #include "DenseLinAlgPack_DVectorClass.hpp"
00034 
00035 namespace {
00036 template< class T >
00037 inline
00038 T my_min( const T& v1, const T& v2 ) { return v1 < v2 ? v1 : v2; }
00039 } // end namespace
00040 
00041 ConstrainedOptPack::value_type
00042 ConstrainedOptPack ::min_abs( const DVectorSlice& mu )
00043 {
00044   if( !mu.dim() )
00045     return 0.0;
00046   value_type min = std::fabs(mu(1));
00047   for( DVectorSlice::const_iterator itr = mu.begin() + 1; itr != mu.end(); )
00048     min = my_min( min, std::fabs(*itr++) );
00049   return min;
00050 }
00051 
00052 ConstrainedOptPack::value_type
00053 ConstrainedOptPack ::min_abs( const SpVectorSlice& mu )
00054 {
00055   if( !mu.dim() )
00056     return 0.0;
00057   if( !mu.nz() )
00058     return 0.0;
00059   value_type min = std::fabs(mu.begin()->value());
00060   for( SpVectorSlice::const_iterator itr = mu.begin() + 1; itr != mu.end(); ++itr )
00061     min = my_min( min, std::fabs(itr->value()) );
00062   return min;
00063 }

Generated on Tue Jul 13 09:30:50 2010 for MOOCHO (Single Doxygen Collection) by  doxygen 1.4.7