43 #ifndef PANZER_RESPONSE_SCATTER_EVALUATOR_FUNCTIONAL_IMPL_HPP 44 #define PANZER_RESPONSE_SCATTER_EVALUATOR_FUNCTIONAL_IMPL_HPP 49 #include "PanzerDiscFE_config.hpp" 51 #include "Phalanx_Evaluator_Macros.hpp" 52 #include "Phalanx_MDField.hpp" 53 #include "Phalanx_DataLayout_MDALayout.hpp" 58 #include "Thyra_DefaultProductVector.hpp" 59 #include "Thyra_SpmdVectorBase.hpp" 60 #include "Thyra_ProductVectorBase.hpp" 62 #include "Teuchos_ArrayRCP.hpp" 69 template<
typename EvalT,
typename Traits>
75 , scatterObj_(functionalScatter)
89 cellIntegral_ = PHX::MDField<const ScalarT,panzer::Cell>(name,dl_cell);
92 std::string n =
"Functional Response Scatter: " + name;
96 template<
typename EvalT,
typename Traits>
99 const std::string & responseName,
102 : responseName_(responseName)
103 , scatterObj_(functionalScatter)
117 cellIntegral_ = PHX::MDField<const ScalarT,panzer::Cell>(integrandName,dl_cell);
120 std::string n =
"Functional Response Scatter: " + responseName;
124 template<
typename EvalT,
typename Traits>
134 template<
typename EvalT,
typename Traits>
139 this->utils.setFieldData(cellIntegral_,fm);
142 template<
typename EvalT,
typename Traits>
147 responseObj_->value += cellIntegral_(i);
156 using Teuchos::rcp_dynamic_cast;
157 using Thyra::SpmdVectorBase;
165 std::vector<Teuchos::ArrayRCP<double> > local_dgdxs;
166 for(
int b=0;b<prod_dgdx->productSpace()->numBlocks();b++) {
170 dgdx->getNonconstLocalData(ptrFromRef(local_dgdx));
174 local_dgdxs.push_back(local_dgdx);
177 scatterObj_->scatterDerivative(cellIntegral_,d,this->wda,local_dgdxs);
180 #ifdef Panzer_BUILD_HESSIAN_SUPPORT 186 using Teuchos::rcp_dynamic_cast;
187 using Thyra::SpmdVectorBase;
195 std::vector<Teuchos::ArrayRCP<double> > local_dgdxs;
196 for(
int b=0;b<prod_dgdx->productSpace()->numBlocks();b++) {
200 dgdx->getNonconstLocalData(ptrFromRef(local_dgdx));
204 local_dgdxs.push_back(local_dgdx);
208 scatterObj_->scatterHessian(cellIntegral_,d,this->wda,local_dgdxs);
ResponseScatterEvaluator_Functional(const std::string &name, const CellData &cd, const Teuchos::RCP< FunctionalScatterBase > &functionalScatter)
A constructor with concrete arguments instead of a parameter list.
static std::string buildLookupName(const std::string &responseName)
void preEvaluate(typename Traits::PreEvalData d)
Teuchos::RCP< GlobalEvaluationData > getDataObject(const std::string &key) const
PHX::MDField< const ScalarT, panzer::Cell > cellIntegral_
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Data for determining cell topology and dimensionality.
void evaluateFields(typename Traits::EvalData d)
GlobalEvaluationDataContainer gedc
Teuchos::RCP< PHX::FieldTag > scatterHolder_
#define TEUCHOS_ASSERT(assertion_test)
std::size_t numCells() const