74 const Teuchos::RCP<FunctionalScatterBase> & functionalScatter)
76 , scatterObj_(functionalScatter)
84 RCP<PHX::DataLayout> dl_dummy = rcp(
new PHX::MDALayout<panzer::Dummy>(0));
89 RCP<PHX::DataLayout> dl_cell = rcp(
new PHX::MDALayout<panzer::Cell>(cd.
numCells()));
90 cellIntegral_ = PHX::MDField<const ScalarT,panzer::Cell>(name,dl_cell);
93 std::string n =
"Functional Response Scatter: " + name;
100 const std::string & responseName,
102 const Teuchos::RCP<FunctionalScatterBase> & functionalScatter)
103 : responseName_(responseName)
104 , scatterObj_(functionalScatter)
112 RCP<PHX::DataLayout> dl_dummy = rcp(
new PHX::MDALayout<panzer::Dummy>(0));
117 RCP<PHX::DataLayout> dl_cell = rcp(
new PHX::MDALayout<panzer::Cell>(cd.
numCells()));
118 cellIntegral_ = PHX::MDField<const ScalarT,panzer::Cell>(integrandName,dl_cell);
121 std::string n =
"Functional Response Scatter: " + responseName;
151 using Teuchos::rcp_dynamic_cast;
152 using Thyra::SpmdVectorBase;
155 TEUCHOS_ASSERT(scatterObj_!=Teuchos::null);
156 TEUCHOS_ASSERT(responseObj_->getGhostedVector()!=Teuchos::null);
158 RCP<ProductVectorBase<double> > prod_dgdx = Thyra::castOrCreateNonconstProductVectorBase(responseObj_->getGhostedVector());
160 std::vector<Teuchos::ArrayRCP<double> > local_dgdxs;
161 for(
int b=0;b<prod_dgdx->productSpace()->numBlocks();b++) {
163 Teuchos::ArrayRCP<double> local_dgdx;
164 RCP<SpmdVectorBase<double> > dgdx = rcp_dynamic_cast<SpmdVectorBase<double> >(prod_dgdx->getNonconstVectorBlock(b));
165 dgdx->getNonconstLocalData(ptrFromRef(local_dgdx));
167 TEUCHOS_ASSERT(!local_dgdx.is_null());
169 local_dgdxs.push_back(local_dgdx);
172 scatterObj_->scatterDerivative(cellIntegral_,d,this->wda,local_dgdxs);
181 using Teuchos::rcp_dynamic_cast;
182 using Thyra::SpmdVectorBase;
185 TEUCHOS_ASSERT(scatterObj_!=Teuchos::null);
186 TEUCHOS_ASSERT(responseObj_->getGhostedVector()!=Teuchos::null);
188 RCP<ProductVectorBase<double> > prod_dgdx = Thyra::castOrCreateNonconstProductVectorBase(responseObj_->getGhostedVector());
190 std::vector<Teuchos::ArrayRCP<double> > local_dgdxs;
191 for(
int b=0;b<prod_dgdx->productSpace()->numBlocks();b++) {
193 Teuchos::ArrayRCP<double> local_dgdx;
194 RCP<SpmdVectorBase<double> > dgdx = rcp_dynamic_cast<SpmdVectorBase<double> >(prod_dgdx->getNonconstVectorBlock(b));
195 dgdx->getNonconstLocalData(ptrFromRef(local_dgdx));
197 TEUCHOS_ASSERT(!local_dgdx.is_null());
199 local_dgdxs.push_back(local_dgdx);
203 scatterObj_->scatterHessian(cellIntegral_,d,this->wda,local_dgdxs);