43 #ifndef PANZER_RESPONSE_SCATTER_EVALUATOR_FUNCTIONAL_HPP
44 #define PANZER_RESPONSE_SCATTER_EVALUATOR_FUNCTIONAL_HPP
49 #include "PanzerDiscFE_config.hpp"
56 #include "Phalanx_Evaluator_Macros.hpp"
57 #include "Phalanx_MDField.hpp"
72 #ifdef Panzer_BUILD_HESSIAN_SUPPORT
80 template <
typename LO,
typename GO>
85 if(globalIndexer!=Teuchos::null)
86 ugis_.push_back(globalIndexer);
97 #ifdef Panzer_BUILD_HESSIAN_SUPPORT
106 std::vector<Teuchos::RCP<const panzer::GlobalIndexer> >
ugis_;
112 template<
typename EvalT,
typename Traits>
138 template <
typename LO,
typename GO>
145 std::string blockId = wda(workset).block_id;
147 std::vector<int> blockOffsets;
152 auto cellIntegral_h = Kokkos::create_mirror_view(cellIntegral.get_view());
153 Kokkos::deep_copy(cellIntegral_h, cellIntegral.get_view());
155 const std::vector<std::size_t> & localCellIds = wda(workset).cell_local_ids;
157 for(std::size_t b=0;b<ugis_.size();b++) {
158 int start = blockOffsets[b];
160 auto LIDs = ugis_[b]->getLIDs();
161 auto LIDs_h = Kokkos::create_mirror_view(LIDs);
162 Kokkos::deep_copy(LIDs_h, LIDs);
167 for(std::size_t worksetCellIndex=0;worksetCellIndex<localCellIds.size();++worksetCellIndex) {
168 std::size_t cellLocalId = localCellIds[worksetCellIndex];
171 for(std::size_t i=0;i<LIDs_h.extent(1);i++) {
172 dgdx_b[LIDs_h(cellLocalId, i)] += cellIntegral_h(worksetCellIndex).dx(start+i);
178 #ifdef Panzer_BUILD_HESSIAN_SUPPORT
179 template <
typename LO,
typename GO>
185 PHX::View<const LO*> LIDs;
188 std::string blockId = wda(workset).block_id;
190 std::vector<int> blockOffsets;
196 const std::vector<std::size_t> & localCellIds = wda(workset).cell_local_ids;
197 for(std::size_t worksetCellIndex=0;worksetCellIndex<localCellIds.size();++worksetCellIndex) {
198 std::size_t cellLocalId = localCellIds[worksetCellIndex];
200 for(std::size_t b=0;b<ugis_.size();b++) {
201 int start = blockOffsets[b];
203 LIDs = ugis_[b]->getElementLIDs(cellLocalId);
208 for(std::size_t i=0;i<LIDs.size();i++) {
209 d2gdx2_b[LIDs[i]] += cellIntegral(worksetCellIndex).dx(start+i).dx(0);
FunctionalScatter(const std::vector< Teuchos::RCP< const panzer::GlobalIndexer > > &ugis)
FunctionalScatter(const Teuchos::RCP< const panzer::GlobalIndexer > &globalIndexer)
virtual ~FunctionalScatterBase()
ResponseScatterEvaluator_Functional(const std::string &name, const CellData &cd, const Teuchos::RCP< FunctionalScatterBase > &functionalScatter)
A constructor with concrete arguments instead of a parameter list.
void scatterDerivative(const PHX::MDField< const panzer::Traits::Jacobian::ScalarT, panzer::Cell > &cellIntegral, panzer::Traits::EvalData workset, WorksetDetailsAccessor &wda, const std::vector< Teuchos::ArrayRCP< double > > &dgdx) const
Teuchos::RCP< FunctionalScatterBase > scatterObj_
virtual void scatterDerivative(const PHX::MDField< const panzer::Traits::Jacobian::ScalarT, panzer::Cell > &cellIntegral, panzer::Traits::EvalData workset, WorksetDetailsAccessor &wda, const std::vector< Teuchos::ArrayRCP< double > > &dgdx) const =0
void preEvaluate(typename Traits::PreEvalData d)
void computeBlockOffsets(const std::string &blockId, const std::vector< Teuchos::RCP< GlobalIndexer >> &ugis, std::vector< int > &blockOffsets)
std::string responseName_
PHX::MDField< const ScalarT, panzer::Cell > cellIntegral_
Data for determining cell topology and dimensionality.
Wrapper to PHX::EvaluatorWithBaseImpl that implements Panzer-specific helpers.
void evaluateFields(typename Traits::EvalData d)
Teuchos::RCP< PHX::FieldTag > scatterHolder_
void scatterHessian(const PHX::MDField< const panzer::Traits::Hessian::ScalarT, panzer::Cell > &cellIntegral, panzer::Traits::EvalData workset, WorksetDetailsAccessor &wda, const std::vector< Teuchos::ArrayRCP< double > > &d2gdx2) const
std::vector< Teuchos::RCP< const panzer::GlobalIndexer > > ugis_
#define TEUCHOS_ASSERT(assertion_test)
Teuchos::RCP< Response_Functional< EvalT > > responseObj_
virtual void scatterHessian(const PHX::MDField< const panzer::Traits::Hessian::ScalarT, panzer::Cell > &cellIntegral, panzer::Traits::EvalData workset, WorksetDetailsAccessor &wda, const std::vector< Teuchos::ArrayRCP< double > > &d2gdx2) const =0