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>
73 const Teuchos::RCP<FunctionalScatterBase> & functionalScatter)
75 , scatterObj_(functionalScatter)
83 RCP<PHX::DataLayout> dl_dummy = rcp(
new PHX::MDALayout<panzer::Dummy>(0));
88 RCP<PHX::DataLayout> dl_cell = rcp(
new PHX::MDALayout<panzer::Cell>(cd.
numCells()));
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,
101 const Teuchos::RCP<FunctionalScatterBase> & functionalScatter)
102 : responseName_(responseName)
103 , scatterObj_(functionalScatter)
111 RCP<PHX::DataLayout> dl_dummy = rcp(
new PHX::MDALayout<panzer::Dummy>(0));
116 RCP<PHX::DataLayout> dl_cell = rcp(
new PHX::MDALayout<panzer::Cell>(cd.
numCells()));
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;
160 TEUCHOS_ASSERT(scatterObj_!=Teuchos::null);
161 TEUCHOS_ASSERT(responseObj_->getGhostedVector()!=Teuchos::null);
163 RCP<ProductVectorBase<double> > prod_dgdx = Thyra::castOrCreateNonconstProductVectorBase(responseObj_->getGhostedVector());
165 std::vector<Teuchos::ArrayRCP<double> > local_dgdxs;
166 for(
int b=0;b<prod_dgdx->productSpace()->numBlocks();b++) {
168 Teuchos::ArrayRCP<double> local_dgdx;
169 RCP<SpmdVectorBase<double> > dgdx = rcp_dynamic_cast<SpmdVectorBase<double> >(prod_dgdx->getNonconstVectorBlock(b));
170 dgdx->getNonconstLocalData(ptrFromRef(local_dgdx));
172 TEUCHOS_ASSERT(!local_dgdx.is_null());
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;
190 TEUCHOS_ASSERT(scatterObj_!=Teuchos::null);
191 TEUCHOS_ASSERT(responseObj_->getGhostedVector()!=Teuchos::null);
193 RCP<ProductVectorBase<double> > prod_dgdx = Thyra::castOrCreateNonconstProductVectorBase(responseObj_->getGhostedVector());
195 std::vector<Teuchos::ArrayRCP<double> > local_dgdxs;
196 for(
int b=0;b<prod_dgdx->productSpace()->numBlocks();b++) {
198 Teuchos::ArrayRCP<double> local_dgdx;
199 RCP<SpmdVectorBase<double> > dgdx = rcp_dynamic_cast<SpmdVectorBase<double> >(prod_dgdx->getNonconstVectorBlock(b));
200 dgdx->getNonconstLocalData(ptrFromRef(local_dgdx));
202 TEUCHOS_ASSERT(!local_dgdx.is_null());
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)
Data for determining cell topology and dimensionality.
void evaluateFields(typename Traits::EvalData d)
GlobalEvaluationDataContainer gedc
Teuchos::RCP< PHX::FieldTag > scatterHolder_
std::size_t numCells() const