ConstrainedOptPack_MeritFuncCalcNLP.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 "ConstrainedOptPack_MeritFuncCalcNLP.hpp"
00030 #include "DenseLinAlgPack_DVectorClass.hpp"
00031 
00032 namespace ConstrainedOptPack {
00033 
00034 MeritFuncCalcNLP::MeritFuncCalcNLP( const MeritFuncNLP* phi, const NLP* nlp )
00035   : phi_(phi), nlp_(nlp)
00036 {}
00037 
00038 value_type MeritFuncCalcNLP::operator()(const Vector& x) const
00039 {
00040   const size_type
00041     m  = nlp().m(),
00042     ns = nlp().ns();
00043   nlp().calc_f(x);
00044   if(m)  nlp().calc_c(x,false);
00045   return phi().value(
00046     nlp().f()
00047     ,m  ? &nlp().c()  : NULL
00048     ,NULL  // h
00049     ,NULL  // hl
00050     ,NULL  // hu
00051     );
00052 /* RAB: 20020112: ToDo: Get this working
00053   if(m)  nlp().calc_c_breve(x,false);
00054   if(ns) nlp().calc_h_breve(x,false);
00055   return phi().value(
00056     nlp().f()
00057     ,m  ? &nlp().c_breve()  : NULL
00058     ,ns ? &nlp().h_breve()  : NULL
00059     ,ns ? &nlp().hl_breve() : NULL
00060     ,ns ? &nlp().hu_breve() : NULL
00061     );
00062 */
00063 }
00064 
00065 value_type MeritFuncCalcNLP::deriv() const {
00066   return phi().deriv();
00067 }
00068 
00069 void MeritFuncCalcNLP::print_merit_func(
00070   std::ostream& out, const std::string& L
00071   ) const
00072 {
00073   out << L << "*** MeritFuncCalcNLP\n"
00074     << L << "f = f(x), c = c_breve(x_breve), h = h_breve(x_breve)\n";
00075   phi().print_merit_func(out,L);
00076 }
00077 
00078 } // end namespace ConstrainedOptPack

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