/usr/include/trilinos/ConstrainedOptPack_DecompositionSystemOrthogonal.hpp is in libtrilinos-dev 10.4.0.dfsg-1ubuntu2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 | // @HEADER
// ***********************************************************************
//
// Moocho: Multi-functional Object-Oriented arCHitecture for Optimization
// Copyright (2003) Sandia Corporation
//
// Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
// license for use of this work by or on behalf of the U.S. Government.
//
// This library is free software; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2.1 of the
// License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA
// Questions? Contact Roscoe A. Bartlett (rabartl@sandia.gov)
//
// ***********************************************************************
// @HEADER
#ifndef DECOMPOSITION_SYSTEM_ORTHOGONAL_H
#define DECOMPOSITION_SYSTEM_ORTHOGONAL_H
#include "ConstrainedOptPack_DecompositionSystemVarReductImp.hpp"
#include "Teuchos_StandardCompositionMacros.hpp"
namespace ConstrainedOptPack {
/** \brief Orthogonal variable reduction subclass.
*
* This is the interface for the coordinate variable reduction decomposition
* where:
\verbatim
Y = [ I ] = [ I ] (class MatrixIdentConcatStd)
[ N'*inv(C') ] [ -D' ]
R = Gc(:,con_decomp)'*Y = [ C N ] * [ I ] = (C + N*N'*inv(C'))
[ N'*inv(C') ]
= C*(I + inv(C)*N*N'*inv(C'))
= C*(I + D*D')
Uy = Gc(:,con_undecomp)'*Y = [ E F ] * [ I ]
[ -D' ]
= E - F*D'
\endverbatim
* See the matrix subclass <tt>MatrixDecompRangeOrthog</tt> for
* details on how linear systems with <tt>R</tt> are solved for.
*
* For now the copy constructor and the assignment operator are not defined.
*/
class DecompositionSystemOrthogonal : public DecompositionSystemVarReductImp {
public:
/** @name Constructors / initializers */
//@{
/** \brief . */
DecompositionSystemOrthogonal(
const VectorSpace::space_ptr_t &space_x = Teuchos::null
,const VectorSpace::space_ptr_t &space_c = Teuchos::null
,const basis_sys_ptr_t &basis_sys = Teuchos::null
,const basis_sys_tester_ptr_t &basis_sys_tester = Teuchos::null
,EExplicitImplicit D_imp = MAT_IMP_EXPLICIT
,EExplicitImplicit Uz_imp = MAT_IMP_EXPLICIT
);
//@}
/** @name Overridden from DecompositionSystem */
//@{
/** \brief . */
const mat_fcty_ptr_t factory_Y() const;
/** \brief . */
const mat_nonsing_fcty_ptr_t factory_R() const;
/** \brief . */
const mat_fcty_ptr_t factory_Uy() const;
//@}
protected:
/** @name Overridden from DecompositionSystemVarReductImp */
//@{
/** \brief . */
void update_D_imp_used(EExplicitImplicit *D_imp_used) const;
/** \brief . */
mat_nonsing_fcty_ptr_t::element_type::obj_ptr_t uninitialize_matrices(
std::ostream *out
,EOutputLevel olevel
,MatrixOp *Y
,MatrixOpNonsing *R
,MatrixOp *Uy
) const;
/** \brief . */
void initialize_matrices(
std::ostream *out
,EOutputLevel olevel
,const mat_nonsing_fcty_ptr_t::element_type::obj_ptr_t &C
,const mat_fcty_ptr_t::element_type::obj_ptr_t &D
,MatrixOp *Y
,MatrixOpNonsing *R
,MatrixOp *Uy
,EMatRelations mat_rel
) const;
/** \brief . */
void print_update_matrices(
std::ostream& out, const std::string& leading_str ) const;
//@}
private:
// ////////////////////////
// Private types
typedef Teuchos::RCP<MatrixSymOpNonsing> S_ptr_t;
// ////////////////////////
// Private data members
mutable S_ptr_t S_ptr_;
// We just hold on to this between calls of uninitialize_matrices() and initialize_matrices()
// ////////////////////////
// Private member functions
// Not defined and not to be called
DecompositionSystemOrthogonal(const DecompositionSystemOrthogonal&);
DecompositionSystemOrthogonal& operator=(const DecompositionSystemOrthogonal&);
}; // end class DecompositionSystemOrthogonal
} // end namespace ConstrainedOptPack
#endif // DECOMPOSITION_SYSTEM_ORTHOGONAL_H
|