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" 52 #ifdef HAVE_PYTRILINOS 55 typedef _object PyObject;
61 class Epetra_Operator;
65 class EpetraVectorOrthogPoly;
66 class EpetraMultiVectorOrthogPoly;
67 class EpetraOperatorOrthogPoly;
69 template <
typename ordinal_type,
typename scalar_type>
class Quadrature;
73 class ProductEpetraVector;
74 class ProductEpetraMultiVector;
75 class ProductEpetraOperator;
93 typedef Teuchos::RCP<Stokhos::ProductEpetraVector>
mp_vector_t;
147 void set_x_dot(
const Teuchos::RCP<const Epetra_Vector> &x_dot );
149 void set_x_dotdot(
const Teuchos::RCP<const Epetra_Vector> &x_dotdot );
151 Teuchos::RCP<const Epetra_Vector>
get_x_dot()
const;
155 void set_x(
const Teuchos::RCP<const Epetra_Vector> &x );
157 Teuchos::RCP<const Epetra_Vector>
get_x()
const;
159 const Teuchos::RCP<
const Teuchos::Polynomial<Epetra_Vector> > &x_poly
162 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
get_x_poly()
const;
165 const Teuchos::RCP<
const Teuchos::Polynomial<Epetra_Vector> > &x_dot_poly
168 const Teuchos::RCP<
const Teuchos::Polynomial<Epetra_Vector> > &x_dotdot_poly
171 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
get_x_dot_poly()
const;
172 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
get_x_dotdot_poly()
const;
194 void set_p(
int l,
const Teuchos::RCP<const Epetra_Vector> &p_l );
196 Teuchos::RCP<const Epetra_Vector>
get_p(
int l)
const;
206 void set_t(
double t );
220 double get_t()
const;
226 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
get_sg_basis()
const;
234 Teuchos::RCP<Stokhos::OrthogPolyExpansion<int,double,Stokhos::StandardStorage<int,double> > >
get_sg_expansion()
const;
256 typedef Teuchos::Array<Teuchos::RCP<const Epetra_Vector> >
p_t;
257 typedef Teuchos::Array<sg_const_vector_t >
p_sg_t;
258 typedef Teuchos::Array<mp_const_vector_t >
p_mp_t;
262 Teuchos::RCP<const Epetra_Vector>
x_dot_;
264 Teuchos::RCP<const Epetra_Vector>
x_;
265 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
x_dot_poly_;
267 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
x_poly_;
283 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
sg_basis_;
284 Teuchos::RCP<const Stokhos::Quadrature<int,double> >
sg_quad_;
285 Teuchos::RCP<Stokhos::OrthogPolyExpansion<int,double,Stokhos::StandardStorage<int,double> > >
sg_exp_;
304 template<
class ObjType>
319 { this->operator=(obj);
evalType_ = evalType; }
372 switch(mvOrientation) {
375 default: TEUCHOS_TEST_FOR_EXCEPT(
true);
388 switch(mvOrientation) {
391 default: TEUCHOS_TEST_FOR_EXCEPT(
true);
441 const Teuchos::RCP<Epetra_MultiVector> &mv
443 ,
const Teuchos::Array<int> ¶mIndexes = Teuchos::Array<int>()
458 Teuchos::RCP<Epetra_MultiVector>
mv_;
475 const Teuchos::RCP<Epetra_MultiVector> &mv
477 ) :
dmv_(mv,orientation) {}
497 Teuchos::RCP<Epetra_Operator>
lo_;
509 bool isAlreadyInverted_ )
529 const Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > &mv
531 ,
const Teuchos::Array<int> ¶mIndexes = Teuchos::Array<int>()
546 Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly >
mv_;
559 SGDerivative(
const Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly > &lo )
563 const Teuchos::RCP< Stokhos::EpetraMultiVectorOrthogPoly > &mv
565 ) :
dmv_(mv,orientation) {}
570 Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
getLinearOp()
const 585 Teuchos::RCP< Stokhos::EpetraOperatorOrthogPoly >
lo_;
601 ,
const Teuchos::Array<int> ¶mIndexes = Teuchos::Array<int>()
636 ) :
dmv_(mv,orientation) {}
763 typedef Teuchos::RCP<Stokhos::EpetraVectorOrthogPoly>
sg_vector_t;
838 void set_W(
const Teuchos::RCP<Epetra_Operator> &W );
839 void set_WPrec(
const Teuchos::RCP<Epetra_Operator> &WPrec );
841 Teuchos::RCP<Epetra_Operator>
get_W()
const;
842 Teuchos::RCP<Epetra_Operator>
get_WPrec()
const;
937 void set_f_poly(
const Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> > &f_poly );
939 Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> >
get_f_poly()
const;
1030 typedef Teuchos::Array<Evaluation<Epetra_Vector> >
g_t;
1064 Teuchos::RCP<Epetra_Operator>
W_;
1078 Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> >
f_poly_;
1140 virtual Teuchos::RCP<const Epetra_Map>
get_x_map()
const = 0;
1143 virtual Teuchos::RCP<const Epetra_Map>
get_f_map()
const = 0;
1146 virtual Teuchos::RCP<const Epetra_Map>
get_p_map(
int l)
const;
1162 virtual Teuchos::RCP<const Teuchos::Array<std::string> >
get_p_names(
int l)
const;
1165 virtual Teuchos::RCP<const Epetra_Map>
get_g_map(
int j)
const;
1181 virtual Teuchos::ArrayView<const std::string>
get_g_names(
int j)
const;
1189 virtual Teuchos::RCP<const Epetra_Vector>
get_x_init()
const;
1192 virtual Teuchos::RCP<const Epetra_Vector>
get_x_dot_init()
const;
1198 virtual Teuchos::RCP<const Epetra_Vector>
get_p_init(
int l)
const;
1243 virtual Teuchos::RCP<Epetra_Operator>
create_W()
const;
1244 virtual Teuchos::RCP<EpetraExt::ModelEvaluator::Preconditioner>
create_WPrec()
const;
1247 virtual Teuchos::RCP<Epetra_Operator>
create_DfDp_op(
int l)
const;
1256 virtual Teuchos::RCP<Epetra_Operator>
create_DgDx_op(
int j)
const;
1259 virtual Teuchos::RCP<Epetra_Operator>
create_DgDp_op(
int j,
int l )
const;
1275 #ifdef HAVE_PYTRILINOS 1277 friend InArgs convertInArgsFromPython(PyObject * source);
1280 friend OutArgs convertOutArgsFromPython(PyObject * source);
1393 Teuchos::RCP<Epetra_Operator>
1395 const std::string &modelEvalDescription,
1397 const std::string &derivName
1401 Teuchos::RCP<Epetra_MultiVector>
1403 const std::string &modelEvalDescription,
1405 const std::string &derivName,
1410 Teuchos::RCP<Epetra_Operator>
1417 Teuchos::RCP<Epetra_MultiVector>
1424 Teuchos::RCP<Epetra_MultiVector>
1432 Teuchos::RCP<Epetra_MultiVector>
1440 Teuchos::RCP<Epetra_MultiVector>
1448 Teuchos::RCP<Epetra_MultiVector>
1469 {
return p_.size(); }
1489 { assert_supports(
IN_ARG_x); x_ = x; }
1493 { assert_supports(
IN_ARG_x);
return x_; }
1504 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
1509 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
1518 Teuchos::RCP<const Teuchos::Polynomial<Epetra_Vector> >
1578 { assert_l(l); p_[l] = p_l; }
1582 { assert_l(l);
return p_[l]; }
1587 { assert_supports(
IN_ARG_p_sg, l); p_sg_[l] = p_sg_l; }
1592 { assert_supports(
IN_ARG_p_sg, l);
return p_sg_[l]; }
1597 { assert_supports(
IN_ARG_p_mp, l); p_mp_[l] = p_mp_l; }
1602 { assert_supports(
IN_ARG_p_mp, l);
return p_mp_[l]; }
1606 { assert_supports(
IN_ARG_t); t_ = t; }
1610 { assert_supports(
IN_ARG_t);
return t_; }
1657 Teuchos::RCP<const Stokhos::OrthogPolyBasis<int,double> >
1666 Teuchos::RCP<const Stokhos::Quadrature<int,double> >
1675 Teuchos::RCP<Stokhos::OrthogPolyExpansion<int,double,Stokhos::StandardStorage<int,double> > >
1682 modelEvalDescription_ = new_modelEvalDescription;
1689 p_sg_.resize(new_Np);
1690 p_mp_.resize(new_Np);
1691 supports_p_sg_.resize(new_Np);
1692 supports_p_mp_.resize(new_Np);
1701 {
return modelEvalDescription_; }
1706 return DfDp_.size();
1779 {
return W_properties_; }
1782 {
return WPrec_properties_; }
1804 return DfDp_properties_[l];
1811 DfDp_sg_[l] = DfDp_sg_l;
1827 return DfDp_sg_properties_[l];
1834 DfDp_mp_[l] = DfDp_mp_l;
1850 return DfDp_mp_properties_[l];
1857 DgDx_dot_[j] = DgDx_dot_j;
1865 return DgDx_dot_[j];
1873 return DgDx_dot_properties_[j];
1880 DgDx_dot_sg_[j] = DgDx_dot_sg_j;
1888 return DgDx_dot_sg_[j];
1896 return DgDx_dot_sg_properties_[j];
1903 DgDx_dot_mp_[j] = DgDx_dot_mp_j;
1911 return DgDx_dot_mp_[j];
1919 return DgDx_dot_mp_properties_[j];
1926 DgDx_dotdot_[j] = DgDx_dotdot_j;
1934 return DgDx_dotdot_[j];
1942 return DgDx_dotdot_properties_[j];
1949 DgDx_dotdot_sg_[j] = DgDx_dotdot_sg_j;
1957 return DgDx_dotdot_sg_[j];
1965 return DgDx_dotdot_sg_properties_[j];
1972 DgDx_dotdot_mp_[j] = DgDx_dotdot_mp_j;
1980 return DgDx_dotdot_mp_[j];
1988 return DgDx_dotdot_mp_properties_[j];
2011 return DgDx_properties_[j];
2018 DgDx_sg_[j] = DgDx_sg_j;
2034 return DgDx_sg_properties_[j];
2041 DgDx_mp_[j] = DgDx_mp_j;
2057 return DgDx_mp_properties_[j];
2064 DgDp_[ j*Np() + l ] = DgDp_j_l;
2072 return DgDp_[ j*Np() + l ];
2080 return DgDp_properties_[ j*Np() + l ];
2087 DgDp_sg_[ j*Np() + l ] = DgDp_sg_j_l;
2095 return DgDp_sg_[ j*Np() + l ];
2103 return DgDp_sg_properties_[ j*Np() + l ];
2110 DgDp_mp_[ j*Np() + l ] = DgDp_mp_j_l;
2118 return DgDp_mp_[ j*Np() + l ];
2126 return DgDp_mp_properties_[ j*Np() + l ];
2131 { f_poly_ = f_poly; }
2134 Teuchos::RCP<Teuchos::Polynomial<Epetra_Vector> >
2175 this->_setModelEvalDescription(new_modelEvalDescription);
2180 { this->_set_Np(new_Np); }
2184 { this->_setSupports(arg,new_supports); }
2188 { this->_setSupports(arg,l,new_supports); }
2192 { this->_setSupports(arg,l,new_supports); }
2201 this->_setModelEvalDescription(new_modelEvalDescription);
2206 { this->_set_Np_Ng(new_Np,new_Ng); }
2210 { this->_setSupports(arg,new_supports); }
2214 { this->_setSupports(arg,l,new_supports); }
2218 { this->_setSupports(arg,j,new_supports); }
2222 { this->_setSupports(arg,j,new_supports); }
2226 { this->_setSupports(arg,j,new_supports); }
2230 { this->_setSupports(arg,j,l,new_supports); }
2234 { this->_setSupports(arg,j,new_supports); }
2238 { this->_setSupports(arg,l,new_supports); }
2242 { this->_setSupports(arg,j,new_supports); }
2246 { this->_setSupports(arg,j,new_supports); }
2250 { this->_setSupports(arg,j,new_supports); }
2254 { this->_setSupports(arg,j,l,new_supports); }
2258 { this->_setSupports(arg,j,new_supports); }
2262 { this->_setSupports(arg,l,new_supports); }
2266 { this->_setSupports(arg,j,new_supports); }
2270 { this->_setSupports(arg,j,new_supports); }
2274 { this->_setSupports(arg,j,new_supports); }
2278 { this->_setSupports(arg,j,l,new_supports); }
2282 { this->_set_W_properties(properties); }
2285 { this->_set_WPrec_properties(properties); }
2290 this->_set_DfDp_properties(l,properties);
2296 this->_set_DgDx_dot_properties(j,properties);
2302 this->_set_DgDx_dotdot_properties(j,properties);
2308 this->_set_DgDx_properties(j,properties);
2314 this->_set_DgDp_properties(j,l,properties);
2320 this->_set_DfDp_sg_properties(l,properties);
2326 this->_set_DgDx_dot_sg_properties(j,properties);
2332 this->_set_DgDx_dotdot_sg_properties(j,properties);
2338 this->_set_DgDx_sg_properties(j,properties);
2344 this->_set_DgDp_sg_properties(j,l,properties);
2350 this->_set_DfDp_mp_properties(l,properties);
2356 this->_set_DgDx_dot_mp_properties(j,properties);
2362 this->_set_DgDx_dotdot_mp_properties(j,properties);
2368 this->_set_DgDx_mp_properties(j,properties);
2374 this->_set_DgDp_mp_properties(j,l,properties);
2379 #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