43 #ifndef PANZER_RESPONSE_SCATTER_EVALUATOR_IPCoordinates_IMPL_HPP 44 #define PANZER_RESPONSE_SCATTER_EVALUATOR_IPCoordinates_IMPL_HPP 49 #include "PanzerDiscFE_config.hpp" 51 #include "Phalanx_Evaluator_Macros.hpp" 52 #include "Phalanx_MDField.hpp" 53 #include "Phalanx_DataLayout_MDALayout.hpp" 65 template<
typename EvalT,
typename Traits>
69 : responseName_(name), ir_order_(ir_order)
77 RCP<PHX::DataLayout> dl_dummy = rcp(
new PHX::MDALayout<panzer::Dummy>(0));
81 std::string n =
"IPCoordinates Response Scatter: " + name;
85 template<
typename EvalT,
typename Traits>
95 template<
typename EvalT,
typename Traits>
103 template<
typename EvalT,
typename Traits>
110 if (tmpCoords_.size() != Teuchos::as<std::size_t>(iv.
ip_coordinates.extent(2))) {
112 for(std::size_t dim=0;dim<tmpCoords_.size();dim++)
113 tmpCoords_[dim].clear();
116 auto ip_coordinates_h = Kokkos::create_mirror_view(PHX::as_view(iv.
ip_coordinates));
117 Kokkos::deep_copy(ip_coordinates_h, PHX::as_view(iv.
ip_coordinates));
123 for (index_t cell = 0; cell < workset.
num_cells; ++cell)
125 tmpCoords_[dim].push_back(ip_coordinates_h(static_cast<int>(cell),ip,dim));
129 template<
typename EvalT,
typename Traits>
133 std::vector<panzer::Traits::Residual::ScalarT> & coords = *responseObj_->getNonconstCoords();
136 for (std::size_t dim = 0; dim < tmpCoords_.size(); ++dim) {
137 for (
typename std::vector<ScalarT>::const_iterator x=tmpCoords_[dim].begin(); x != tmpCoords_[dim].end(); ++ x)
138 coords.push_back(Sacado::ScalarValue<ScalarT>::eval(*x));
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)
Teuchos::RCP< GlobalEvaluationDataContainer > gedc
int num_cells
DEPRECATED - use: numCells()
static std::string buildLookupName(const std::string &responseName)
ResponseScatterEvaluator_IPCoordinates(const std::string &name, int ir_order)
A constructor with concrete arguments instead of a parameter list.
std::vector< int >::size_type getIntegrationRuleIndex(int ir_degree, const panzer::Workset &workset, WorksetDetailsAccessor &wda)
void evaluateFields(typename Traits::EvalData d)
Array_CellIPDim ip_coordinates
void postEvaluate(typename Traits::PostEvalData d)
void preEvaluate(typename Traits::PreEvalData d)
Teuchos::RCP< const std::vector< panzer::Workset > > worksets_
Teuchos::RCP< PHX::FieldTag > scatterHolder_