/usr/include/trilinos/AbstractLinAlgPack_DirectSparseSolverDense.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 148 149 150 151 152 153 154 155 | // @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 DIRECT_SPARSE_SOLVER_DENSE_H
#define DIRECT_SPARSE_SOLVER_DENSE_H
#include <valarray>
#include <vector>
#include <string>
#include "AbstractLinAlgPack_DirectSparseSolverImp.hpp"
#include "DenseLinAlgPack_DMatrixClass.hpp"
#include "DenseLinAlgPack_IVector.hpp"
#include "Teuchos_StandardMemberCompositionMacros.hpp"
namespace AbstractLinAlgPack {
/** \brief Concreate sparse solver subclass that uses the dense LAPACK routines.
*
* ToDo: Finish documentation!
*/
class DirectSparseSolverDense : public DirectSparseSolverImp {
public:
/** @name Constructors/initializers */
//@{
/** \brief Default constructor */
DirectSparseSolverDense();
//@}
/** @name Overridden from DirectSparseSolver */
//@{
/** \brief . */
const basis_matrix_factory_ptr_t basis_matrix_factory() const;
/** \brief . */
void estimated_fillin_ratio( value_type estimated_fillin_ratio );
//@}
protected:
/** @name Protected types */
//@{
/** \brief Implements the BasisMatrix object for Dense.
*/
class BasisMatrixDense : public BasisMatrixImp {
public:
/** @name Overridden from BasisMatrixImp */
//@{
/** \brief . */
Teuchos::RCP<BasisMatrixImp> create_matrix() const;
/** \brief . */
void V_InvMtV(
VectorMutable* v_lhs, BLAS_Cpp::Transp trans_rhs1
,const Vector& v_rhs2) const ;
//@}
}; // end class BasisMatrixDense
/** \brief Stores the factorization structure for Dense
*/
class FactorizationStructureDense : public FactorizationStructure {
public:
friend class DirectSparseSolverDense;
friend class BasisMatrixDense;
private:
FortranTypes::f_int m_; // Number of rows in A
FortranTypes::f_int n_; // Number of columns in A
FortranTypes::f_int nz_; // Number of nonzeros in A
FortranTypes::f_int rank_; // Rank of the basis
IVector col_perm_; // First rank entries selects the basis of A
IVector inv_col_perm_; // Inverse of col_perm_
FactorizationStructureDense();
}; // end class FactorizationStructureDense
/** \brief Stores the factorization nonzeros for Dense
*/
class FactorizationNonzerosDense : public FactorizationNonzeros {
public:
typedef FortranTypes::f_int f_int;
friend class DirectSparseSolverDense;
friend class BasisMatrixDense;
private:
DMatrix LU_;
bool rect_analyze_and_factor_; // true for n > m analyze_and_factor()
std::valarray<f_int> ipiv_; // The permutation sent to xGETRS (identity if rect_analyze_and_factor_==true)
IVector basis_perm_; // Only used if rect_analyze_and_factor_==true
}; // end class FactorizationNonzerosDense
//@}
/** @name Overridden from DirectSparseSolverImp */
//@{
/** \brief . */
const Teuchos::RCP<FactorizationStructure> create_fact_struc() const;
/** \brief . */
const Teuchos::RCP<FactorizationNonzeros> create_fact_nonzeros() const;
/** \brief . */
void imp_analyze_and_factor(
const AbstractLinAlgPack::MatrixConvertToSparse &A
,FactorizationStructure *fact_struc
,FactorizationNonzeros *fact_nonzeros
,DenseLinAlgPack::IVector *row_perm
,DenseLinAlgPack::IVector *col_perm
,size_type *rank
,std::ostream *out
);
/** \brief . */
void imp_factor(
const AbstractLinAlgPack::MatrixConvertToSparse &A
,const FactorizationStructure &fact_struc
,FactorizationNonzeros *fact_nonzeros
,std::ostream *out
);
//@}
}; // end class DirectSparseSolverDense
} // end namespace AbstractLinAlgPack
#endif // DIRECT_SPARSE_SOLVER_DENSE_H
|