45#include "Teuchos_ScalarTraits.hpp"
53DiagonalQuadraticResponseOnlyModelEvaluator
54::DiagonalQuadraticResponseOnlyModelEvaluator(
55 const Teuchos::RCP<Epetra_Comm> &comm,
56 const int localDim,
const double &pt,
const double &p0,
const double &scale
58 :epetra_comm_(comm), scale_(scale)
80Teuchos::RCP<const Epetra_Map>
87Teuchos::RCP<const Epetra_Map>
94Teuchos::RCP<const Epetra_Map>
97 TEUCHOS_TEST_FOR_EXCEPT(l!=0);
102Teuchos::RCP<const Epetra_Map>
105 TEUCHOS_TEST_FOR_EXCEPT(j!=0);
110Teuchos::RCP<const Epetra_Vector>
113 TEUCHOS_TEST_FOR_EXCEPT(l!=0);
133 outArgs.set_Np_Ng(1, 1);
135 outArgs.set_DgDp_properties(
152 using Teuchos::dyn_cast;
153 using Teuchos::rcp_dynamic_cast;
165 const RCP<Epetra_Vector> g_out = outArgs.get_g(0);
167 const RCP<Epetra_MultiVector> DgDp_trans_out =
174 if (nonnull(g_out) || nonnull(DgDp_trans_out)) {
179 p_minus_pt.Update(-1.0, *
pt_, 1.0);
181 if (nonnull(g_out)) {
183 p_minus_pt.Dot(p_minus_pt, dot);
184 (*g_out)[0] =
scale_ * 0.5 * dot[0];
187 if (nonnull(DgDp_trans_out)) {
188 (*DgDp_trans_out) = p_minus_pt;
189 DgDp_trans_out->Scale(
scale_);
void evalModel(const InArgs &inArgs, const OutArgs &outArgs) const
InArgs createInArgs() const
Teuchos::RCP< const Epetra_Map > get_p_map(int l) const
\breif .
OutArgs createOutArgs() const
Teuchos::RCP< const Epetra_Comm > epetra_comm_
Teuchos::RCP< const Epetra_Vector > get_p_init(int l) const
Teuchos::RCP< const Epetra_Map > get_g_map(int j) const
\breif .
Teuchos::RCP< const Epetra_Map > map_g_
Teuchos::RCP< const Epetra_Map > get_f_map() const
Teuchos::RCP< const Epetra_Map > get_x_map() const
Teuchos::RCP< Epetra_Vector > p0_
Teuchos::RCP< Epetra_Vector > pt_
Teuchos::RCP< const Epetra_Map > map_p_
void setModelEvalDescription(const std::string &modelEvalDescription)
void setModelEvalDescription(const std::string &modelEvalDescription)
@ DERIV_LINEARITY_NONCONST
EpetraExt::BlockCrsMatrix: A class for constructing a distributed block matrix.
Teuchos::RCP< Epetra_MultiVector > get_DgDp_mv(const int j, const int l, const ModelEvaluator::OutArgs &outArgs, const ModelEvaluator::EDerivativeMultiVectorOrientation mvOrientation)