42 #ifndef STOKHOS_REDUCED_QUADRATURE_FACTORY_HPP 43 #define STOKHOS_REDUCED_QUADRATURE_FACTORY_HPP 45 #include "Teuchos_RCP.hpp" 46 #include "Teuchos_Array.hpp" 47 #include "Teuchos_BLAS.hpp" 48 #include "Teuchos_LAPACK.hpp" 49 #include "Teuchos_SerialDenseMatrix.hpp" 50 #include "Teuchos_ParameterList.hpp" 64 template <
typename ordinal_type,
typename value_type>
78 virtual Teuchos::RCP<const Stokhos::UserDefinedQuadrature<ordinal_type, value_type> >
80 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Q,
81 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>&
Q2,
82 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& F,
83 const Teuchos::Array<value_type>& weights)
const;
88 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Q,
89 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& F,
90 const Teuchos::Array<value_type>& weights,
91 Teuchos::RCP< Teuchos::Array<value_type> >& red_weights,
92 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_points,
93 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_values
97 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Q,
98 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& F,
99 const Teuchos::Array<value_type>& weights,
100 Teuchos::RCP< Teuchos::Array<value_type> >& red_weights,
101 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_points,
102 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_values
106 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Q,
107 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& F,
108 const Teuchos::Array<value_type>& weights,
109 Teuchos::RCP< Teuchos::Array<value_type> >& red_weights,
110 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_points,
111 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_values
115 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Q,
116 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>&
Q2,
117 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& F,
118 const Teuchos::Array<value_type>& weights,
119 Teuchos::RCP< Teuchos::Array<value_type> >& red_weights,
120 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_points,
121 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_values
125 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Q,
126 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>&
Q2,
127 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& F,
128 const Teuchos::Array<value_type>& weights,
129 Teuchos::RCP< Teuchos::Array<value_type> >& red_weights,
130 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_points,
131 Teuchos::RCP< Teuchos::Array< Teuchos::Array<value_type> > >& red_values
135 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& A,
136 const Teuchos::SerialDenseVector<ordinal_type, value_type>& b,
137 Teuchos::SerialDenseVector<ordinal_type, value_type>&
x,
138 Teuchos::ETransp transa, Teuchos::EUplo uplo)
const;
141 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& A,
142 const Teuchos::SerialDenseVector<ordinal_type, value_type>& b,
143 Teuchos::SerialDenseVector<ordinal_type, value_type>&
x,
144 Teuchos::ETransp transa, Teuchos::EUplo uplo)
const;
147 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& A,
148 const Teuchos::SerialDenseVector<ordinal_type, value_type>& b,
149 Teuchos::SerialDenseVector<ordinal_type, value_type>&
x,
150 Teuchos::ETransp transa, Teuchos::EUplo uplo)
const;
153 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& A,
154 const Teuchos::SerialDenseVector<ordinal_type, value_type>& b,
155 Teuchos::SerialDenseVector<ordinal_type, value_type>&
x,
156 Teuchos::ETransp transa, Teuchos::EUplo uplo)
const;
159 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& A,
160 const Teuchos::SerialDenseVector<ordinal_type, value_type>& b,
161 Teuchos::SerialDenseVector<ordinal_type, value_type>&
x,
162 Teuchos::ETransp transa, Teuchos::EUplo uplo)
const;
165 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& A,
166 const Teuchos::SerialDenseVector<ordinal_type, value_type>& b,
167 Teuchos::SerialDenseVector<ordinal_type, value_type>&
x,
168 Teuchos::ETransp transa, Teuchos::EUplo uplo)
const;
171 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& A,
172 const Teuchos::SerialDenseVector<ordinal_type, value_type>& b,
173 Teuchos::SerialDenseVector<ordinal_type, value_type>&
x,
174 Teuchos::ETransp transa, Teuchos::EUplo uplo)
const;
177 const Teuchos::SerialDenseMatrix<ordinal_type,value_type>& R,
217 Teuchos::BLAS<ordinal_type,value_type>
blas;
value_type reduction_tol
Dimension reduction tolerance.
value_type objective_value
Value used in LP-based objective function.
bool eliminate_dependent_rows
Whether to eliminate dependent rows in constraints.
bool verbose
Whether to print a bunch of stuff out.
virtual ~ReducedQuadratureFactory()
Destructor.
std::string reduction_method
Reduction method.
void reducedQuadrature_Q_Squared(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &F, const Teuchos::Array< value_type > &weights, Teuchos::RCP< Teuchos::Array< value_type > > &red_weights, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_points, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_values) const
Teuchos::ParameterList params
Algorithm parameters.
const IndexType const IndexType const IndexType const IndexType const ValueType const ValueType * x
std::string solver_method
Underdetermined solver method.
void underdetermined_solver(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const Teuchos::SerialDenseVector< ordinal_type, value_type > &b, Teuchos::SerialDenseVector< ordinal_type, value_type > &x, Teuchos::ETransp transa, Teuchos::EUplo uplo) const
Teuchos::LAPACK< ordinal_type, value_type > lapack
void solver_GLPK(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const Teuchos::SerialDenseVector< ordinal_type, value_type > &b, Teuchos::SerialDenseVector< ordinal_type, value_type > &x, Teuchos::ETransp transa, Teuchos::EUplo uplo) const
void solver_CLP_IP(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const Teuchos::SerialDenseVector< ordinal_type, value_type > &b, Teuchos::SerialDenseVector< ordinal_type, value_type > &x, Teuchos::ETransp transa, Teuchos::EUplo uplo) const
Top-level namespace for Stokhos classes and functions.
ReducedQuadratureFactory(const Teuchos::ParameterList ¶ms)
Constructor.
void solver_CLP(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const Teuchos::SerialDenseVector< ordinal_type, value_type > &b, Teuchos::SerialDenseVector< ordinal_type, value_type > &x, Teuchos::ETransp transa, Teuchos::EUplo uplo) const
void reducedQuadrature_Q2_CPQR(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q2, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &F, const Teuchos::Array< value_type > &weights, Teuchos::RCP< Teuchos::Array< value_type > > &red_weights, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_points, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_values) const
Generate a basis from a given set of PCE expansions that is orthogonal with respect to the product me...
ordinal_type n_choose_k(const ordinal_type &n, const ordinal_type &k) const
Compute bionomial coefficient (n ; k) = n!/( k! (n-k)! )
void reducedQuadrature_Q2(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q2, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &F, const Teuchos::Array< value_type > &weights, Teuchos::RCP< Teuchos::Array< value_type > > &red_weights, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_points, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_values) const
ordinal_type computeRank(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &R, const value_type tol) const
void solver_TRSM(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const Teuchos::SerialDenseVector< ordinal_type, value_type > &b, Teuchos::SerialDenseVector< ordinal_type, value_type > &x, Teuchos::ETransp transa, Teuchos::EUplo uplo) const
void solver_qpOASES(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const Teuchos::SerialDenseVector< ordinal_type, value_type > &b, Teuchos::SerialDenseVector< ordinal_type, value_type > &x, Teuchos::ETransp transa, Teuchos::EUplo uplo) const
ReducedQuadratureFactory & operator=(const ReducedQuadratureFactory &)
void reducedQuadrature_Q_Squared_CPQR(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &F, const Teuchos::Array< value_type > &weights, Teuchos::RCP< Teuchos::Array< value_type > > &red_weights, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_points, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_values) const
void solver_CompressedSensing(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &A, const Teuchos::SerialDenseVector< ordinal_type, value_type > &b, Teuchos::SerialDenseVector< ordinal_type, value_type > &x, Teuchos::ETransp transa, Teuchos::EUplo uplo) const
virtual Teuchos::RCP< const Stokhos::UserDefinedQuadrature< ordinal_type, value_type > > createReducedQuadrature(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q2, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &F, const Teuchos::Array< value_type > &weights) const
Get reduced quadrature object.
void reducedQuadrature_Q_Squared_CPQR2(const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &Q, const Teuchos::SerialDenseMatrix< ordinal_type, value_type > &F, const Teuchos::Array< value_type > &weights, Teuchos::RCP< Teuchos::Array< value_type > > &red_weights, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_points, Teuchos::RCP< Teuchos::Array< Teuchos::Array< value_type > > > &red_values) const
Teuchos::BLAS< ordinal_type, value_type > blas