42 #ifndef EPETRA_EXT_MODEL_EVALUATOR_HPP 43 #define EPETRA_EXT_MODEL_EVALUATOR_HPP 47 #include "Teuchos_RCP.hpp" 48 #include "Teuchos_Describable.hpp" 49 #include "Teuchos_Polynomial.hpp" 50 #include "Teuchos_Array.hpp" 51 #include "Teuchos_CompilerCodeTweakMacros.hpp" 53 #ifdef HAVE_PYTRILINOS 56 typedef _object PyObject;
62 class Epetra_Operator;
66 class EpetraVectorOrthogPoly;
67 class EpetraMultiVectorOrthogPoly;
68 class EpetraOperatorOrthogPoly;
70 template <
typename ordinal_type,
typename scalar_type>
class Quadrature;
74 class ProductEpetraVector;
75 class ProductEpetraMultiVector;
76 class ProductEpetraOperator;
94 typedef Teuchos::RCP<Stokhos::ProductEpetraVector>
mp_vector_t;
148 void set_x_dot(
const Teuchos::RCP<const Epetra_Vector> &x_dot );
150 void set_x_dotdot(
const Teuchos::RCP<const Epetra_Vector> &x_dotdot );
152 Teuchos::RCP<const Epetra_Vector>
get_x_dot()
const;
156 void set_x(
const Teuchos::RCP<const Epetra_Vector> &x );
158 Teuchos::RCP<const Epetra_Vector>
get_x()
const;
160 const Teuchos::RCP<
const Teuchos::Polynomial<Epetra_Vector> > &x_poly
163 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
get_x_poly()
const;
166 const Teuchos::RCP<
const Teuchos::Polynomial<Epetra_Vector> > &x_dot_poly
169 const Teuchos::RCP<
const Teuchos::Polynomial<Epetra_Vector> > &x_dotdot_poly
172 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
get_x_dot_poly()
const;
173 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
get_x_dotdot_poly()
const;
195 void set_p(
int l,
const Teuchos::RCP<const Epetra_Vector> &p_l );
197 Teuchos::RCP<const Epetra_Vector>
get_p(
int l)
const;
207 void set_t(
double t );
221 double get_t()
const;
227 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
get_sg_basis()
const;
235 Teuchos::RCP<Stokhos::OrthogPolyExpansion<int,double,Stokhos::StandardStorage<int,double> > >
get_sg_expansion()
const;
257 typedef Teuchos::Array<Teuchos::RCP<const Epetra_Vector> >
p_t;
258 typedef Teuchos::Array<sg_const_vector_t >
p_sg_t;
259 typedef Teuchos::Array<mp_const_vector_t >
p_mp_t;
263 Teuchos::RCP<const Epetra_Vector>
x_dot_;
265 Teuchos::RCP<const Epetra_Vector>
x_;
266 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
x_dot_poly_;
268 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
x_poly_;
284 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis_;
285 Teuchos::RCP<const Stokhos::Quadrature<int,double> >
sg_quad_;
286 Teuchos::RCP<Stokhos::OrthogPolyExpansion<int,double,Stokhos::StandardStorage<int,double> > >
sg_exp_;
305 template<
class ObjType>
320 { this->operator=(obj);
evalType_ = evalType; }
373 switch(mvOrientation) {
376 default: TEUCHOS_TEST_FOR_EXCEPT(
true);
389 switch(mvOrientation) {
392 default: TEUCHOS_TEST_FOR_EXCEPT(
true);
394 TEUCHOS_UNREACHABLE_RETURN(
false);
442 const Teuchos::RCP<Epetra_MultiVector> &mv
444 ,
const Teuchos::Array<int> ¶mIndexes = Teuchos::Array<int>()
459 Teuchos::RCP<Epetra_MultiVector>
mv_;
476 const Teuchos::RCP<Epetra_MultiVector> &mv
478 ) :
dmv_(mv,orientation) {}
498 Teuchos::RCP<Epetra_Operator>
lo_;
510 bool isAlreadyInverted_ )
530 const Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > &mv
532 ,
const Teuchos::Array<int> ¶mIndexes = Teuchos::Array<int>()
547 Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly >
mv_;
560 SGDerivative(
const Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > &lo )
564 const Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > &mv
566 ) :
dmv_(mv,orientation) {}
571 Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
getLinearOp()
const 586 Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
lo_;
602 ,
const Teuchos::Array<int> ¶mIndexes = Teuchos::Array<int>()
637 ) :
dmv_(mv,orientation) {}
764 typedef Teuchos::RCP<Stokhos::EpetraVectorOrthogPoly>
sg_vector_t;
839 void set_W(
const Teuchos::RCP<Epetra_Operator> &W );
840 void set_WPrec(
const Teuchos::RCP<Epetra_Operator> &WPrec );
842 Teuchos::RCP<Epetra_Operator>
get_W()
const;
843 Teuchos::RCP<Epetra_Operator>
get_WPrec()
const;
938 void set_f_poly(
const Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> > &f_poly );
940 Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> >
get_f_poly()
const;
1031 typedef Teuchos::Array<Evaluation<Epetra_Vector> >
g_t;
1065 Teuchos::RCP<Epetra_Operator>
W_;
1079 Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> >
f_poly_;
1141 virtual Teuchos::RCP<const Epetra_Map>
get_x_map()
const = 0;
1144 virtual Teuchos::RCP<const Epetra_Map>
get_f_map()
const = 0;
1147 virtual Teuchos::RCP<const Epetra_Map>
get_p_map(
int l)
const;
1163 virtual Teuchos::RCP<const Teuchos::Array<std::string> >
get_p_names(
int l)
const;
1166 virtual Teuchos::RCP<const Epetra_Map>
get_g_map(
int j)
const;
1182 virtual Teuchos::ArrayView<const std::string>
get_g_names(
int j)
const;
1190 virtual Teuchos::RCP<const Epetra_Vector>
get_x_init()
const;
1193 virtual Teuchos::RCP<const Epetra_Vector>
get_x_dot_init()
const;
1199 virtual Teuchos::RCP<const Epetra_Vector>
get_p_init(
int l)
const;
1244 virtual Teuchos::RCP<Epetra_Operator>
create_W()
const;
1245 virtual Teuchos::RCP<EpetraExt::ModelEvaluator::Preconditioner>
create_WPrec()
const;
1248 virtual Teuchos::RCP<Epetra_Operator>
create_DfDp_op(
int l)
const;
1257 virtual Teuchos::RCP<Epetra_Operator>
create_DgDx_op(
int j)
const;
1260 virtual Teuchos::RCP<Epetra_Operator>
create_DgDp_op(
int j,
int l )
const;
1276 #ifdef HAVE_PYTRILINOS 1278 friend InArgs convertInArgsFromPython(PyObject * source);
1281 friend OutArgs convertOutArgsFromPython(PyObject * source);
1394 Teuchos::RCP<Epetra_Operator>
1396 const std::string &modelEvalDescription,
1398 const std::string &derivName
1402 Teuchos::RCP<Epetra_MultiVector>
1404 const std::string &modelEvalDescription,
1406 const std::string &derivName,
1411 Teuchos::RCP<Epetra_Operator>
1418 Teuchos::RCP<Epetra_MultiVector>
1425 Teuchos::RCP<Epetra_MultiVector>
1433 Teuchos::RCP<Epetra_MultiVector>
1441 Teuchos::RCP<Epetra_MultiVector>
1449 Teuchos::RCP<Epetra_MultiVector>
1470 {
return p_.size(); }
1490 { assert_supports(
IN_ARG_x); x_ = x; }
1494 { assert_supports(
IN_ARG_x);
return x_; }
1505 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
1510 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
1519 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
1579 { assert_l(l); p_[l] = p_l; }
1583 { assert_l(l);
return p_[l]; }
1588 { assert_supports(
IN_ARG_p_sg, l); p_sg_[l] = p_sg_l; }
1593 { assert_supports(
IN_ARG_p_sg, l);
return p_sg_[l]; }
1598 { assert_supports(
IN_ARG_p_mp, l); p_mp_[l] = p_mp_l; }
1603 { assert_supports(
IN_ARG_p_mp, l);
return p_mp_[l]; }
1607 { assert_supports(
IN_ARG_t); t_ = t; }
1611 { assert_supports(
IN_ARG_t);
return t_; }
1658 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
1667 Teuchos::RCP<const Stokhos::Quadrature<int,double> >
1676 Teuchos::RCP<Stokhos::OrthogPolyExpansion<int,double,Stokhos::StandardStorage<int,double> > >
1683 modelEvalDescription_ = new_modelEvalDescription;
1690 p_sg_.resize(new_Np);
1691 p_mp_.resize(new_Np);
1692 supports_p_sg_.resize(new_Np);
1693 supports_p_mp_.resize(new_Np);
1702 {
return modelEvalDescription_; }
1707 return DfDp_.size();
1780 {
return W_properties_; }
1783 {
return WPrec_properties_; }
1805 return DfDp_properties_[l];
1812 DfDp_sg_[l] = DfDp_sg_l;
1828 return DfDp_sg_properties_[l];
1835 DfDp_mp_[l] = DfDp_mp_l;
1851 return DfDp_mp_properties_[l];
1858 DgDx_dot_[j] = DgDx_dot_j;
1866 return DgDx_dot_[j];
1874 return DgDx_dot_properties_[j];
1881 DgDx_dot_sg_[j] = DgDx_dot_sg_j;
1889 return DgDx_dot_sg_[j];
1897 return DgDx_dot_sg_properties_[j];
1904 DgDx_dot_mp_[j] = DgDx_dot_mp_j;
1912 return DgDx_dot_mp_[j];
1920 return DgDx_dot_mp_properties_[j];
1927 DgDx_dotdot_[j] = DgDx_dotdot_j;
1935 return DgDx_dotdot_[j];
1943 return DgDx_dotdot_properties_[j];
1950 DgDx_dotdot_sg_[j] = DgDx_dotdot_sg_j;
1958 return DgDx_dotdot_sg_[j];
1966 return DgDx_dotdot_sg_properties_[j];
1973 DgDx_dotdot_mp_[j] = DgDx_dotdot_mp_j;
1981 return DgDx_dotdot_mp_[j];
1989 return DgDx_dotdot_mp_properties_[j];
2012 return DgDx_properties_[j];
2019 DgDx_sg_[j] = DgDx_sg_j;
2035 return DgDx_sg_properties_[j];
2042 DgDx_mp_[j] = DgDx_mp_j;
2058 return DgDx_mp_properties_[j];
2065 DgDp_[ j*Np() + l ] = DgDp_j_l;
2073 return DgDp_[ j*Np() + l ];
2081 return DgDp_properties_[ j*Np() + l ];
2088 DgDp_sg_[ j*Np() + l ] = DgDp_sg_j_l;
2096 return DgDp_sg_[ j*Np() + l ];
2104 return DgDp_sg_properties_[ j*Np() + l ];
2111 DgDp_mp_[ j*Np() + l ] = DgDp_mp_j_l;
2119 return DgDp_mp_[ j*Np() + l ];
2127 return DgDp_mp_properties_[ j*Np() + l ];
2132 { f_poly_ = f_poly; }
2135 Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> >
2176 this->_setModelEvalDescription(new_modelEvalDescription);
2181 { this->_set_Np(new_Np); }
2185 { this->_setSupports(arg,new_supports); }
2189 { this->_setSupports(arg,l,new_supports); }
2193 { this->_setSupports(arg,l,new_supports); }
2202 this->_setModelEvalDescription(new_modelEvalDescription);
2207 { this->_set_Np_Ng(new_Np,new_Ng); }
2211 { this->_setSupports(arg,new_supports); }
2215 { this->_setSupports(arg,l,new_supports); }
2219 { this->_setSupports(arg,j,new_supports); }
2223 { this->_setSupports(arg,j,new_supports); }
2227 { this->_setSupports(arg,j,new_supports); }
2231 { this->_setSupports(arg,j,l,new_supports); }
2235 { this->_setSupports(arg,j,new_supports); }
2239 { this->_setSupports(arg,l,new_supports); }
2243 { this->_setSupports(arg,j,new_supports); }
2247 { this->_setSupports(arg,j,new_supports); }
2251 { this->_setSupports(arg,j,new_supports); }
2255 { this->_setSupports(arg,j,l,new_supports); }
2259 { this->_setSupports(arg,j,new_supports); }
2263 { this->_setSupports(arg,l,new_supports); }
2267 { this->_setSupports(arg,j,new_supports); }
2271 { this->_setSupports(arg,j,new_supports); }
2275 { this->_setSupports(arg,j,new_supports); }
2279 { this->_setSupports(arg,j,l,new_supports); }
2283 { this->_set_W_properties(properties); }
2286 { this->_set_WPrec_properties(properties); }
2291 this->_set_DfDp_properties(l,properties);
2297 this->_set_DgDx_dot_properties(j,properties);
2303 this->_set_DgDx_dotdot_properties(j,properties);
2309 this->_set_DgDx_properties(j,properties);
2315 this->_set_DgDp_properties(j,l,properties);
2321 this->_set_DfDp_sg_properties(l,properties);
2327 this->_set_DgDx_dot_sg_properties(j,properties);
2333 this->_set_DgDx_dotdot_sg_properties(j,properties);
2339 this->_set_DgDx_sg_properties(j,properties);
2345 this->_set_DgDp_sg_properties(j,l,properties);
2351 this->_set_DfDp_mp_properties(l,properties);
2357 this->_set_DgDx_dot_mp_properties(j,properties);
2363 this->_set_DgDx_dotdot_mp_properties(j,properties);
2369 this->_set_DgDx_mp_properties(j,properties);
2375 this->_set_DgDp_mp_properties(j,l,properties);
2380 #endif // EPETRA_EXT_MODEL_EVALUATOR_HPP virtual Teuchos::RCP< const Teuchos::Array< std::string > > get_p_names(int l) const
Get the names of the parameters associated with parameter subvector l if available.
Simple aggregate class that stores a derivative object as a general linear operator or as a multi-vec...
EEvalType getType() const
void set_stage_number(int stage_number)
virtual double getInfBound() const
Return the value of an infinite bound.
deriv_properties_t DgDp_mp_properties_
Derivative get_DgDp(int j, int l) const
void setFailed() const
Set that the evaluation as a whole failed.
supports_t supports_DgDx_dotdot_sg_
EDerivativeMultiVectorOrientation getMultiVectorOrientation() const
Residual vector Taylor polynomial.
int get_stage_number() const
supports_t supports_DgDx_dot_
MPDerivative(const mp_operator_t &lo)
void set_WPrec_properties(const DerivativeProperties &properties)
Teuchos::RCP< const Epetra_Vector > x_dotdot_
MPDerivativeMultiVector(const mp_multivector_t &mv, const EDerivativeMultiVectorOrientation orientation=DERIV_MV_BY_COL, const Teuchos::Array< int > ¶mIndexes=Teuchos::Array< int >())
SGDerivative get_DgDp_sg(int j, int l) const
void set_W_properties(const DerivativeProperties &properties)
virtual InArgs createInArgs() const =0
EDerivativeMultiVectorOrientation getOrientation() const
void assert_j(int j) const
MPDerivativeMultiVector dmv_
void set_step_size(double step_size)
mp_operator_t get_W_mp() const
Get multi-point W.
Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > x_poly_
Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > get_x_dot_poly() const
Get time derivative vector Taylor polynomial.
EDerivativeMultiVectorOrientation getOrientation() const
SGDerivative get_DgDx_sg(int j) const
void set_DgDx(int j, const Derivative &DgDx_j)
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > getLinearOp() const
void setSupports(EOutArgsMembers arg, bool supports=true)
std::string modelEvalDescription_
supports_g_sg_t supports_g_mp_
A very approximate derivative (i.e. for a preconditioner)
mp_const_vector_t get_x_dot_mp() const
Get multi-point time derivative vector.
void set_x_dotdot_poly(const Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > &x_dotdot_poly)
DerivativeProperties get_DgDp_sg_properties(int j, int l) const
Teuchos::RCP< Epetra_MultiVector > getMultiVector() const
DerivativeMultiVector dmv_
Preconditioner()
Default constructor of null Operatir.
DerivativeProperties get_DgDx_dotdot_properties(int j) const
mp_const_vector_t get_p_mp(int l) const
Get multi-point parameter vector.
Time second derivative vector Taylor polynomial.
std::string modelEvalDescription_
void _setSupports(EOutArgsMembers arg, bool supports)
Teuchos::Array< DerivativeSupport > supports_t
void _set_DfDp_properties(int l, const DerivativeProperties &properties)
Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > get_x_poly() const
Get solution vector Taylor polynomial.
double get_step_size() const
Exact function evaluation.
void set_DgDx_sg_properties(int j, const DerivativeProperties &properties)
virtual Teuchos::RCP< const Epetra_Map > get_f_map() const =0
.
virtual Teuchos::RCP< Epetra_Operator > create_DgDx_dot_op(int j) const
supports_t supports_DgDp_
virtual Teuchos::RCP< const Epetra_Map > get_x_map() const =0
.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > sg_operator_t
Short-hand for stochastic Galerkin operator type.
MPDerivativeMultiVector()
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > get_sg_basis() const
DerivativeProperties get_DgDx_dot_properties(int j) const
Teuchos::Array< mp_vector_t > g_mp_t
DerivativeProperties WPrec_properties_
void set_beta(double beta)
virtual ~ModelEvaluator()
static const int NUM_E_IN_ARGS_MEMBERS
deriv_properties_t DfDp_mp_properties_
DerivativeProperties get_DfDp_sg_properties(int l) const
bool isAlreadyInverted
Bool flag.
Preconditioner(const Teuchos::RCP< Epetra_Operator > &PrecOp_, bool isAlreadyInverted_)
Usable constructor to set the (Epetra_Operator,bool) pair.
void set_x(const Teuchos::RCP< const Epetra_Vector > &x)
void set_x_dot(const Teuchos::RCP< const Epetra_Vector > &x_dot)
void set_x_dotdot_mp(const mp_const_vector_t &x_dotdot_mp)
Teuchos::RCP< const Stokhos::Quadrature< int, double > > get_sg_quadrature() const
Teuchos::RCP< Epetra_Operator > getLinearOp(const std::string &modelEvalDescription, const ModelEvaluator::Derivative &deriv, const std::string &derivName)
Multi-point solution vector.
Teuchos::RCP< const Epetra_Vector > get_x_dotdot() const
Teuchos::RCP< Epetra_MultiVector > get_DgDp_mv(const int j, const int l, const ModelEvaluator::OutArgs &outArgs, const ModelEvaluator::EDerivativeMultiVectorOrientation mvOrientation)
Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > sg_basis_
MPDerivative get_DfDp_mp(int l) const
SGDerivativeMultiVector(const Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > &mv, const EDerivativeMultiVectorOrientation orientation=DERIV_MV_BY_COL, const Teuchos::Array< int > ¶mIndexes=Teuchos::Array< int >())
virtual Teuchos::RCP< const Epetra_Vector > get_x_upper_bounds() const
supports_t supports_DfDp_sg_
mp_operator_t getLinearOp() const
EDerivativeMultiVectorOrientation orientation_
Teuchos::Array< MPDerivative > mp_deriv_t
EDerivativeMultiVectorOrientation orientation_
Teuchos::RCP< Epetra_Operator > get_WPrec() const
Teuchos::RCP< Epetra_MultiVector > get_DgDx_dot_mv(const int j, const ModelEvaluator::OutArgs &outArgs, const ModelEvaluator::EDerivativeMultiVectorOrientation mvOrientation)
void set_DfDp(int l, const Derivative &DfDp_l)
void set_WPrec(const Teuchos::RCP< Epetra_Operator > &WPrec)
Teuchos::RCP< Stokhos::ProductEpetraMultiVector > mp_multivector_t
Stochastic Galerkin quadrature.
Stochastic Galerkin expansion.
void _set_DgDp_mp_properties(int j, int l, const DerivativeProperties &properties)
Teuchos::RCP< Epetra_Operator > W_
Teuchos::RCP< Stokhos::ProductEpetraOperator > mp_operator_t
EDerivativeMultiVectorOrientation getOrientation() const
deriv_properties_t DfDp_sg_properties_
void set_DgDx_dotdot(int j, const Derivative &DgDx_dotdot_j)
SGDerivativeMultiVector getDerivativeMultiVector() const
void set_g_sg(int j, const sg_vector_t &g_sg_j)
Set stochastic Galerkin vector polynomial response.
Teuchos::RCP< const Stokhos::ProductEpetraMultiVector > mp_const_multivector_t
void set_f(const Evaluation< Epetra_Vector > &f)
void setSupports(EInArgsMembers arg, bool supports=true)
Teuchos::RCP< Epetra_Operator > PrecOp
Accessor for the Epetra_Operator.
virtual Teuchos::RCP< const Epetra_Vector > get_p_upper_bounds(int l) const
DerivativeSupport(EDerivativeMultiVectorOrientation mvOrientation)
DerivativeProperties get_DgDp_mp_properties(int j, int l) const
SGDerivative get_DgDx_dot_sg(int j) const
Teuchos::RCP< const Epetra_Vector > get_x_dot() const
virtual Teuchos::RCP< Epetra_Operator > create_DfDp_op(int l) const
Teuchos::Array< sg_vector_t > g_sg_t
Teuchos::Array< int > paramIndexes_
supports_t supports_DgDx_sg_
virtual double get_t_init() const
void set_DfDp_sg(int l, const SGDerivative &DfDp_sg_l)
virtual Teuchos::RCP< Epetra_Operator > create_DgDx_op(int j) const
Simple aggregate class that stores a derivative object as a general linear operator or as a multi-vec...
Evaluation< Epetra_Vector > f_
DerivativeProperties get_DgDx_dot_sg_properties(int j) const
Teuchos::RCP< const Epetra_Vector > x_dot_
EDerivativeMultiVectorOrientation orientation_
supports_t supports_DgDx_mp_
void assert_l(int l) const
SGDerivative get_DgDx_dotdot_sg(int j) const
Teuchos::RCP< Epetra_MultiVector > getMultiVector(const std::string &modelEvalDescription, const ModelEvaluator::Derivative &deriv, const std::string &derivName, const ModelEvaluator::EDerivativeMultiVectorOrientation mvOrientation)
mp_vector_t get_g_mp(int j) const
Get multi-point response.
void changeOrientation(const EDerivativeMultiVectorOrientation orientation)
DerivativeProperties get_DgDx_dot_mp_properties(int j) const
deriv_properties_t DgDx_dotdot_mp_properties_
EpetraExt::BlockCrsMatrix: A class for constructing a distributed block matrix.
Simple aggregate class that stores a derivative object as a general linear operator or as a multi-vec...
deriv_properties_t DgDx_dot_mp_properties_
supports_t supports_DgDx_dot_sg_
void set_sg_basis(const Teuchos::RCP< const Stokhos::OrthogPolyBasis< int, double > > &basis)
mp_const_vector_t get_x_mp() const
Get multi-point solution vector.
EDerivativeLinearity linearity
sg_vector_t get_g_sg(int j) const
Get stochastic Galerkin vector polynomial response.
Teuchos::RCP< Stokhos::OrthogPolyExpansion< int, double, Stokhos::StandardStorage< int, double > > > sg_exp_
Multi-point "W" operator.
deriv_properties_t DgDx_dot_properties_
DerivativeSupport(EDerivativeLinearOp)
Simple aggregate class for a derivative object represented as a column-wise multi-vector or its trans...
void set_x_dotdot(const Teuchos::RCP< const Epetra_Vector > &x_dotdot)
bool supports(EDerivativeMultiVectorOrientation mvOrientation) const
Stochastic Galerkin solution vector polynomial.
Teuchos::RCP< const Epetra_Vector > get_x() const
Set solution vector Taylor polynomial.
Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > x_dot_poly_
sg_const_vector_t get_x_dotdot_sg() const
SGDerivativeMultiVector dmv_
void set_W_sg(const sg_operator_t &W_sg)
Set stochastic Galerkin W operator polynomial.
Teuchos::RCP< Stokhos::EpetraVectorOrthogPoly > sg_vector_t
Short-hand for stochastic Galerkin vector type.
void set_DgDx_dot_sg(int j, const SGDerivative &DgDx_dot_j)
virtual Teuchos::RCP< EpetraExt::ModelEvaluator::Preconditioner > create_WPrec() const
void set_DgDx_dot(int j, const Derivative &DgDx_dot_j)
bool funcOrDerivesAreSet(EOutArgsMembers arg) const
Return true if the function or its derivatives are set.
MPDerivative get_DgDx_dotdot_mp(int j) const
An approximate derivative (i.e. for a Jacobian)
supports_p_sg_t supports_p_mp_
void set_DgDx_dot_mp_properties(int j, const DerivativeProperties &properties)
void set_DgDp_properties(int j, int l, const DerivativeProperties &properties)
void set_DgDx_dotdot_mp(int j, const MPDerivative &DgDx_dotdot_j)
Evaluation(const Teuchos::RCP< ObjType > &obj)
DerivativeProperties get_DgDx_dotdot_sg_properties(int j) const
void _set_DgDx_dotdot_mp_properties(int j, const DerivativeProperties &properties)
void set_DgDx_sg(int j, const SGDerivative &DgDx_j)
Derivative get_DfDp(int l) const
sg_vector_t get_f_sg() const
Get stochastic Galerkin residual vector polynomial.
virtual Teuchos::RCP< const Epetra_Vector > get_x_init() const
Teuchos::RCP< Teuchos::Polynomial< Epetra_Vector > > get_f_poly() const
Get residual vector Taylor polynomial.
Teuchos::RCP< const Stokhos::EpetraVectorOrthogPoly > sg_const_vector_t
Short-hand for stochastic Galerkin vector type.
virtual Teuchos::RCP< const Epetra_Vector > get_x_dotdot_init() const
void set_p_mp(int l, const mp_const_vector_t &p_mp_l)
Set multi-point parameter vector.
void set_DgDx_mp(int j, const MPDerivative &DgDx_j)
void set_DgDp_sg(int j, int l, const SGDerivative &DgDp_sg_j_l)
void set_DfDp_properties(int l, const DerivativeProperties &properties)
virtual Teuchos::ArrayView< const std::string > get_g_names(int j) const
Get the names of the response functions associated with response subvector j if available.
void setModelEvalDescription(const std::string &modelEvalDescription)
sg_const_vector_t get_x_sg() const
Get stochastic Galerkin solution vector polynomial.
Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > lo_
void set_x_dotdot_sg(const sg_const_vector_t &x_dotdot_sg)
mp_multivector_t getMultiVector() const
MPDerivativeMultiVector getDerivativeMultiVector() const
Teuchos::RCP< Epetra_MultiVector > get_DfDp_mv(const int l, const ModelEvaluator::OutArgs &outArgs)
EDerivativeMultiVectorOrientation getMultiVectorOrientation() const
supports_t supports_DgDx_dotdot_mp_
Teuchos::RCP< Epetra_Operator > get_W() const
Derivative(const DerivativeMultiVector &dmv)
void set_p(int l, const Teuchos::RCP< const Epetra_Vector > &p_l)
void _set_DgDx_mp_properties(int j, const DerivativeProperties &properties)
void set_DgDp_sg_properties(int j, int l, const DerivativeProperties &properties)
supports_t supports_DgDp_sg_
virtual Teuchos::RCP< Epetra_Operator > create_DgDx_dotdot_op(int j) const
Simple aggregate class for a derivative object represented as a column-wise multi-vector or its trans...
void set_DgDx_dot_properties(int j, const DerivativeProperties &properties)
SGDerivativeMultiVector()
MPDerivative get_DgDx_mp(int j) const
Teuchos::RCP< Epetra_Operator > get_DfDp_op(const int l, const ModelEvaluator::OutArgs &outArgs)
void set_x_dot_sg(const sg_const_vector_t &x_dot_sg)
Set stochastic Galerkin time derivative vector polynomial.
Teuchos::RCP< Epetra_MultiVector > mv_
MPDerivative(const MPDerivativeMultiVector &dmv)
DerivativeSupport & plus(EDerivativeMultiVectorOrientation mvOrientation)
Multi-point time second derivative vector.
mp_const_vector_t get_x_dotdot_mp() const
Preconditioner operator (approx Jacobian)
MPDerivative(const mp_multivector_t &mv, const EDerivativeMultiVectorOrientation orientation=DERIV_MV_BY_COL)
Teuchos::RCP< const Stokhos::ProductEpetraVector > mp_const_vector_t
deriv_properties_t DgDx_sg_properties_
bool isFailed() const
Return if the evaluation failed or not.
Teuchos::RCP< Epetra_MultiVector > get_DgDx_mv(const int j, const ModelEvaluator::OutArgs &outArgs, const ModelEvaluator::EDerivativeMultiVectorOrientation mvOrientation)
Teuchos::Array< int > paramIndexes_
const Teuchos::Array< int > & getParamIndexes() const
const Teuchos::Array< int > & getParamIndexes() const
void _set_WPrec_properties(const DerivativeProperties &WPrec_properties)
sg_const_vector_t get_p_sg(int l) const
Get stochastic Galerkin vector polynomial parameter.
void set_W(const Teuchos::RCP< Epetra_Operator > &W)
void _set_DgDx_dot_mp_properties(int j, const DerivativeProperties &properties)
Derivative(const Teuchos::RCP< Epetra_Operator > &lo)
void _set_DgDx_dot_properties(int j, const DerivativeProperties &properties)
void assert_supports(EInArgsMembers arg) const
EDerivativeMultiVectorOrientation
void assert_supports(EOutArgsMembers arg) const
Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > mv_
Multi-point residual vector.
virtual double get_t_upper_bound() const
virtual Teuchos::RCP< Epetra_Operator > create_W() const
If supported, create a Epetra_Operator object for W to be evaluated.
sg_const_vector_t get_x_dot_sg() const
Get stochastic Galerkin time derivative vector polynomial.
deriv_properties_t DgDx_dot_sg_properties_
virtual Teuchos::RCP< Epetra_Operator > create_DgDp_op(int j, int l) const
void set_sg_expansion(const Teuchos::RCP< Stokhos::OrthogPolyExpansion< int, double, Stokhos::StandardStorage< int, double > > > &exp)
Teuchos::Array< int > paramIndexes_
void setModelEvalDescription(const std::string &modelEvalDescription)
virtual Teuchos::RCP< const Epetra_Map > get_g_map(int j) const
.
void set_f_sg(const sg_vector_t &f_sg)
Set stochastic Galerkin residual vector polynomial.
DerivativeProperties(EDerivativeLinearity in_linearity, ERankStatus in_rank, bool in_supportsAdjoint)
SGDerivative(const Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > &lo)
DerivativeMultiVector getDerivativeMultiVector() const
mp_const_vector_t x_dot_mp_
DerivativeProperties get_DgDx_mp_properties(int j) const
sg_operator_t get_W_sg() const
Get stochastic Galerkin W operator polynomial.
Teuchos::RCP< const Stokhos::Quadrature< int, double > > sg_quad_
void set_DgDp(int j, int l, const Derivative &DgDp_j_l)
void set_W_mp(const mp_operator_t &W_sg)
Set multi-point W.
Teuchos::Array< Evaluation< Epetra_Vector > > g_t
Derivative get_DgDx_dot(int j) const
void set_f_poly(const Teuchos::RCP< Teuchos::Polynomial< Epetra_Vector > > &f_poly)
Set residual vector Taylor polynomial.
void set_DgDx_dotdot_sg_properties(int j, const DerivativeProperties &properties)
void _set_DgDx_dot_sg_properties(int j, const DerivativeProperties &properties)
Teuchos::RCP< const Epetra_Vector > get_p(int l) const
Multi-point time derivative vector.
Derivative get_DgDx(int j) const
supports_t supports_DfDp_
Teuchos::Array< DerivativeProperties > deriv_properties_t
deriv_properties_t DgDx_properties_
virtual Teuchos::RCP< const Epetra_Vector > get_x_lower_bounds() const
void set_x_dot_mp(const mp_const_vector_t &x_dot_mp)
Set multi-point time derivative vector.
void set_DfDp_mp_properties(int l, const DerivativeProperties &properties)
void set_DgDx_dot_mp(int j, const MPDerivative &DgDx_dot_j)
Derivative(const Teuchos::RCP< Epetra_MultiVector > &mv, const EDerivativeMultiVectorOrientation orientation=DERIV_MV_BY_COL)
void set_x_poly(const Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > &x_poly)
DerivativeMultiVector(const Teuchos::RCP< Epetra_MultiVector > &mv, const EDerivativeMultiVectorOrientation orientation=DERIV_MV_BY_COL, const Teuchos::Array< int > ¶mIndexes=Teuchos::Array< int >())
bool supports(EDerivativeLinearOp) const
void changeOrientation(const EDerivativeMultiVectorOrientation orientation)
void reset(const Teuchos::RCP< ObjType > &obj, EEvalType evalType)
Simple aggregate struct that stores a preconditioner as an Epetra_Operator and a bool, about whether it is inverted or not.
void _set_DfDp_mp_properties(int l, const DerivativeProperties &properties)
Solution vector Taylor polynomial.
Stochastic Galerkin time derivative vector polynomial.
EDerivativeMultiVectorOrientation getMultiVectorOrientation() const
Teuchos::Array< Teuchos::RCP< const Epetra_Vector > > p_t
void set_omega(double omega)
supports_t supports_DgDx_dot_mp_
void _set_DgDp_properties(int j, int l, const DerivativeProperties &properties)
Stochastic Galerkin residual vector polynomial.
DerivativeProperties get_WPrec_properties() const
bool supports(EOutArgsMembers arg) const
virtual Teuchos::RCP< const Epetra_Vector > get_x_dot_init() const
Stochastic Galerkin "W" operator polyomial.
deriv_properties_t DgDp_properties_
Teuchos::Array< sg_const_vector_t > p_sg_t
void set_alpha(double alpha)
Teuchos::RCP< Epetra_MultiVector > get_DgDx_dotdot_mv(const int j, const ModelEvaluator::OutArgs &outArgs, const ModelEvaluator::EDerivativeMultiVectorOrientation mvOrientation)
deriv_properties_t DgDx_dotdot_properties_
SGDerivative(const Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > &mv, const EDerivativeMultiVectorOrientation orientation=DERIV_MV_BY_COL)
Teuchos::RCP< Epetra_MultiVector > getMultiVector() const
Teuchos::RCP< const Epetra_Vector > x_
DerivativeSupport(EDerivativeLinearOp, EDerivativeMultiVectorOrientation mvOrientation)
void set_DfDp_sg_properties(int l, const DerivativeProperties &properties)
DerivativeProperties get_DfDp_properties(int l) const
Teuchos::Array< mp_const_vector_t > p_mp_t
std::string modelEvalDescription() const
virtual void evalModel(const InArgs &inArgs, const OutArgs &outArgs) const =0
supports_t supports_DgDx_
supports_t supports_DgDx_dotdot_
supports_t supports_DgDp_mp_
Teuchos::RCP< Teuchos::Polynomial< Epetra_Vector > > f_poly_
Teuchos::RCP< Stokhos::ProductEpetraVector > mp_vector_t
void set_DfDp_mp(int l, const MPDerivative &DfDp_mp_l)
DerivativeProperties get_DgDx_properties(int j) const
void set_DgDx_dotdot_mp_properties(int j, const DerivativeProperties &properties)
sg_deriv_t DgDx_dotdot_sg_
virtual Teuchos::RCP< const Epetra_Vector > get_p_init(int l) const
Time derivative vector Taylor polynomial.
Simple aggregate class for a derivative object represented as a column-wise multi-vector or its trans...
virtual Teuchos::RCP< const Epetra_Map > get_p_map(int l) const
.
bool supportsTransMVByRow_
void _set_DgDx_properties(int j, const DerivativeProperties &properties)
Evaluation< Epetra_Vector > get_g(int j) const
Get g(j) where 0 <= j && j < this->Ng().
DerivativeProperties get_W_properties() const
DerivativeProperties get_DgDp_properties(int j, int l) const
void set_DgDx_properties(int j, const DerivativeProperties &properties)
deriv_properties_t DgDp_sg_properties_
Teuchos::Array< Derivative > deriv_t
SGDerivative get_DfDp_sg(int l) const
Teuchos::Array< bool > supports_g_sg_t
sg_const_vector_t x_dot_sg_
void set_DgDp_mp(int j, int l, const MPDerivative &DgDp_mp_j_l)
void set_DgDx_dot_sg_properties(int j, const DerivativeProperties &properties)
Teuchos::RCP< Epetra_Operator > lo_
void _set_Np_Ng(int Np, int Ng)
const Teuchos::Array< int > & getParamIndexes() const
SGDerivative(const SGDerivativeMultiVector &dmv)
std::string modelEvalDescription() const
supports_g_sg_t supports_g_sg_
static const int NUM_E_OUT_ARGS_MEMBERS
void set_x_dot_poly(const Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > &x_dot_poly)
Set time derivative vector Taylor polynomial.
void _set_DgDp_sg_properties(int j, int l, const DerivativeProperties &properties)
void assert_l(int l) const
void set_g_mp(int j, const mp_vector_t &g_mp_j)
Set multi-point response.
void set_x_mp(const mp_const_vector_t &x_mp)
Set multi-point solution vector.
mp_vector_t get_f_mp() const
Get multi-point residual vector.
Evaluation< Epetra_Vector > get_f() const
Teuchos::RCP< Epetra_Operator > getLinearOp() const
Stochastic Galerkin basis.
supports_t supports_DfDp_mp_
Stochastic Galerkin time second derivative vector polynomial.
DerivativeProperties get_DfDp_mp_properties(int l) const
void set_sg_quadrature(const Teuchos::RCP< const Stokhos::Quadrature< int, double > > &quad)
DerivativeProperties get_DgDx_dotdot_mp_properties(int j) const
void set_x_sg(const sg_const_vector_t &x_sg)
Set stochastic Galerkin solution vector polynomial.
void set_DgDx_dotdot_sg(int j, const SGDerivative &DgDx_dotdot_j)
DerivativeProperties get_DgDx_sg_properties(int j) const
DerivativeProperties W_properties_
void _set_DgDx_dotdot_properties(int j, const DerivativeProperties &properties)
void _setSupports(EInArgsMembers arg, bool supports)
bool supports_[NUM_E_OUT_ARGS_MEMBERS]
Teuchos::RCP< Epetra_Operator > WPrec_
void set_f_mp(const mp_vector_t &f_mp)
Set multi-point residual vector.
void _set_DgDx_sg_properties(int j, const DerivativeProperties &properties)
virtual Teuchos::RCP< const Epetra_Vector > get_p_lower_bounds(int l) const
void _set_W_properties(const DerivativeProperties &W_properties)
void changeOrientation(const EDerivativeMultiVectorOrientation orientation)
mp_deriv_t DgDx_dotdot_mp_
sg_const_vector_t x_dotdot_sg_
Derivative get_DgDx_dotdot(int j) const
Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > getMultiVector() const
deriv_properties_t DfDp_properties_
DerivativeSupport(EDerivativeMultiVectorOrientation mvOrientation1, EDerivativeMultiVectorOrientation mvOrientation2)
DerivativeSupport & plus(EDerivativeLinearOp)
virtual OutArgs createOutArgs() const =0
void _set_DgDx_dotdot_sg_properties(int j, const DerivativeProperties &properties)
virtual double get_t_lower_bound() const
Teuchos::Array< SGDerivative > sg_deriv_t
void set_p_sg(int l, const sg_const_vector_t &p_sg_l)
Set stochastic Galerkin vector polynomial parameter.
Base interface for evaluating a stateless "model".
deriv_properties_t DgDx_dotdot_sg_properties_
void set_DgDx_dotdot_properties(int j, const DerivativeProperties &properties)
supports_p_sg_t supports_p_sg_
void _setModelEvalDescription(const std::string &modelEvalDescription)
Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > x_dotdot_poly_
void set_Np_Ng(int Np, int Ng)
MPDerivative get_DgDp_mp(int j, int l) const
void _set_DfDp_sg_properties(int l, const DerivativeProperties &properties)
std::string toString(const int &x)
bool supports(EInArgsMembers arg) const
Teuchos::RCP< const Teuchos::Polynomial< Epetra_Vector > > get_x_dotdot_poly() const
Teuchos::Array< bool > supports_p_sg_t
mp_const_vector_t x_dotdot_mp_
Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > getMultiVector() const
void _setModelEvalDescription(const std::string &modelEvalDescription)
Teuchos::RCP< const Stokhos::ProductEpetraOperator > mp_const_operator_t
void set_g(int j, const Evaluation< Epetra_Vector > &g_j)
Set g(j) where 0 <= j && j < this->Ng().
bool supports_[NUM_E_IN_ARGS_MEMBERS]
mp_multivector_t getMultiVector() const
void set_DgDx_mp_properties(int j, const DerivativeProperties &properties)
MPDerivative get_DgDx_dot_mp(int j) const
void set_DgDp_mp_properties(int j, int l, const DerivativeProperties &properties)
deriv_properties_t DgDx_mp_properties_
Evaluation(const Teuchos::RCP< ObjType > &obj, EEvalType evalType)
Teuchos::RCP< Stokhos::OrthogPolyExpansion< int, double, Stokhos::StandardStorage< int, double > > > get_sg_expansion() const