45#include "Teuchos_ScalarTraits.hpp"
55 ,
const bool showGetInvalidArg
57 :d_(d),showGetInvalidArg_(showGetInvalidArg)
72 int indices[nx] = { 0, 1 };
73 for(
int i = 0; i < nx; ++i )
74 W_graph_->InsertGlobalIndices(i,nx,indices);
84Teuchos::RCP<const Epetra_Map>
90Teuchos::RCP<const Epetra_Map>
96Teuchos::RCP<const Epetra_Vector>
102Teuchos::RCP<Epetra_Operator>
124 outArgs.set_W_properties(
136 using Teuchos::dyn_cast;
137 using Teuchos::rcp_dynamic_cast;
156 f[0] = x[0] + x[1]*x[1] - p[0];
157 f[1] =
d_ * ( x[0]*x[0] - x[1] - p[1] );
171 values[0] = 1.0; indexes[0] = 0;
172 values[1] = 2.0*x[1]; indexes[1] = 1;
173 DfDx.SumIntoGlobalValues( 0, 2, values, indexes );
175 values[0] = 2.0*
d_*x[0]; indexes[0] = 0;
176 values[1] = -
d_; indexes[1] = 1;
177 DfDx.SumIntoGlobalValues( 1, 2, values, indexes );
void setModelEvalDescription(const std::string &modelEvalDescription)
void setModelEvalDescription(const std::string &modelEvalDescription)
@ DERIV_LINEARITY_NONCONST
Teuchos::RCP< Epetra_Vector > x0_
EpetraModelEval2DSim(const double d=10.0, const double p0=2.0, const double p1=0.0, const double x00=1.0, const double x01=1.0, const bool showGetInvalidArg=false)
OutArgs createOutArgs() const
Teuchos::RCP< Epetra_Vector > p_
Teuchos::RCP< const Epetra_Comm > epetra_comm_
Teuchos::RCP< const Epetra_Vector > get_x_init() const
Teuchos::RCP< const Epetra_Map > get_x_map() const
Teuchos::RCP< Epetra_Operator > create_W() const
InArgs createInArgs() const
Teuchos::RCP< const Epetra_Map > get_f_map() const
Teuchos::RCP< const Epetra_Map > map_x_
void evalModel(const InArgs &inArgs, const OutArgs &outArgs) const
Teuchos::RCP< Epetra_CrsGraph > W_graph_