MOOCHO (Single Doxygen Collection) Version of the Day
MoochoPack_EvalNewPointTailoredApproachCoordinate_Step.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 // 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 "MoochoPack_EvalNewPointTailoredApproachCoordinate_Step.hpp"
00043 #include "ConstrainedOptPack_MatrixIdentConcatStd.hpp"
00044 #include "NLPInterfacePack_NLPDirect.hpp"
00045 #include "AbstractLinAlgPack_MatrixOp.hpp"
00046 #include "AbstractLinAlgPack_MatrixZero.hpp"
00047 #include "AbstractLinAlgPack_VectorMutable.hpp"
00048 #include "Teuchos_dyn_cast.hpp"
00049 
00050 namespace MoochoPack {
00051 
00052 EvalNewPointTailoredApproachCoordinate_Step::EvalNewPointTailoredApproachCoordinate_Step(
00053   const deriv_tester_ptr_t&     deriv_tester
00054   ,const bounds_tester_ptr_t&   bounds_tester
00055   ,EFDDerivTesting              fd_deriv_testing
00056   )
00057   :EvalNewPointTailoredApproach_Step(deriv_tester,bounds_tester,fd_deriv_testing)
00058 {}
00059 
00060 // protected
00061 
00062 void EvalNewPointTailoredApproachCoordinate_Step::uninitialize_Y_Uy(
00063   MatrixOp         *Y
00064   ,MatrixOp        *Uy
00065   )
00066 {
00067   // Nothing to free
00068 }
00069 
00070 void EvalNewPointTailoredApproachCoordinate_Step::calc_py_Y_Uy(
00071   const NLPDirect       &nlp
00072   ,const D_ptr_t        &D
00073   ,VectorMutable        *py
00074   ,MatrixOp             *Y
00075   ,MatrixOp             *Uy
00076   ,EJournalOutputLevel  olevel
00077   ,std::ostream         &out
00078   )
00079 {
00080   namespace rcp = MemMngPack;
00081   using Teuchos::dyn_cast;
00082 
00083   MatrixIdentConcatStd
00084     &cY = dyn_cast<MatrixIdentConcatStd>(*Y);
00085   //
00086   // Y = [      I     ] space_xD  
00087   //     [    Zero    ] space_xI
00088   //        space_xD
00089   //
00090   VectorSpace::space_ptr_t
00091     space_x  = nlp.space_x(),
00092     space_xD = space_x->sub_space(nlp.var_dep())->clone(),
00093     space_xI = space_x->sub_space(nlp.var_indep())->clone();
00094   cY.initialize(
00095     space_x                                                // space_cols
00096     ,space_xD                                              // space_rows
00097     ,MatrixIdentConcatStd::BOTTOM                          // top_or_bottom
00098     ,1.0                                                   // alpha
00099     ,Teuchos::rcp(
00100       new MatrixZero(
00101         space_xI    // space_cols
00102         ,space_xD   // space_rows
00103         ) )                                            // D_ptr
00104     ,BLAS_Cpp::no_trans                                    // D_trans
00105     );
00106   // py is not altered here!
00107 }
00108 
00109 void EvalNewPointTailoredApproachCoordinate_Step::recalc_py(
00110   const MatrixOp          &D
00111   ,VectorMutable          *py
00112   ,EJournalOutputLevel    olevel
00113   ,std::ostream           &out
00114   )
00115 {
00116   // py is not altered here!
00117 }
00118 
00119 void EvalNewPointTailoredApproachCoordinate_Step::print_calc_py_Y_Uy(
00120   std::ostream& out, const std::string& L
00121   ) const
00122 {
00123   out
00124     << L << "*** Coordinate decomposition\n"
00125     << L << "py_k = py_k\n"
00126     << L << "Y = [ I ; 0 ] <: R^(n x m) [0 represented using MatrixZero]\n"
00127     << L << "Uy = Gc(var_dep,con_undecomp)\'\n"
00128     ;
00129 }
00130 
00131 } // end namespace MoochoPack 
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines