NLPInterfacePack: C++ Interfaces and Implementation for Non-Linear Programs Version of the Day
NLPInterfacePack_NLPSecondOrder.cpp
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 // Redistribution and use in source and binary forms, with or without
00011 // modification, are permitted provided that the following conditions are
00012 // met:
00013 //
00014 // 1. Redistributions of source code must retain the above copyright
00015 // notice, this list of conditions and the following disclaimer.
00016 //
00017 // 2. Redistributions in binary form must reproduce the above copyright
00018 // notice, this list of conditions and the following disclaimer in the
00019 // documentation and/or other materials provided with the distribution.
00020 //
00021 // 3. Neither the name of the Corporation nor the names of the
00022 // contributors may be used to endorse or promote products derived from
00023 // this software without specific prior written permission.
00024 //
00025 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
00026 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00027 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
00028 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
00029 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00030 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00031 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00032 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00033 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00034 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00035 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00036 //
00037 // Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov) 
00038 // 
00039 // ***********************************************************************
00040 // @HEADER
00041 
00042 #include "NLPInterfacePack_NLPSecondOrder.hpp"
00043 #include "Teuchos_Assert.hpp"
00044 
00045 namespace {
00046   const char name_HL[] = "HL";
00047 }
00048 
00049 namespace NLPInterfacePack {
00050 
00051 // constructors
00052 
00053 NLPSecondOrder::NLPSecondOrder()
00054   : HL_(NULL)
00055 {}
00056 
00057 
00058 void NLPSecondOrder::initialize(bool test_setup) {
00059   num_HL_evals_ = 0;
00060   NLPFirstOrder::initialize(test_setup);
00061 }
00062 
00063 // <<std aggr>> members for HL
00064 
00065 void NLPSecondOrder::set_HL(MatrixSymOp* HL)
00066 {
00067   HL_ = HL;
00068 }
00069 
00070 MatrixSymOp* NLPSecondOrder::get_HL()
00071 {
00072   return StandardCompositionRelationshipsPack::get_role_name(HL_, false, name_HL);
00073 }
00074 
00075 MatrixSymOp& NLPSecondOrder::HL()
00076 {
00077   return StandardCompositionRelationshipsPack::role_name(HL_, false, name_HL);
00078 }
00079 
00080 const MatrixSymOp& NLPSecondOrder::HL() const
00081 {
00082   return StandardCompositionRelationshipsPack::role_name(HL_, false, name_HL);
00083 }
00084 
00085 void NLPSecondOrder::unset_quantities()
00086 {
00087   NLPFirstOrder::unset_quantities();
00088   HL_ = NULL;
00089 }
00090 
00091 // calculations
00092 
00093 void NLPSecondOrder::calc_HL(
00094   const Vector& x, const Vector* lambda, bool newpoint
00095   ) const
00096 {
00097 #ifdef TEUCHOS_DEBUG
00098   TEUCHOS_TEST_FOR_EXCEPTION( lambda  && this->m()  == 0, std::logic_error, "" );
00099 #endif
00100   StandardCompositionRelationshipsPack::assert_role_name_set(HL_, "NLP::calc_HL()", name_HL);
00101   imp_calc_HL(x,lambda,newpoint,second_order_info());
00102 }
00103 
00104 size_type NLPSecondOrder::num_HL_evals() const
00105 {
00106   return num_HL_evals_;
00107 }
00108 
00109 } // namespace NLPInterfacePack
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends