58 const Teuchos::ParameterList& p)
59 : evaluateOnClosure_(false)
62 p.validateParameters(*valid_params);
64 const std::string inName = p.get<std::string>(
"Field Name");
65 const std::string outName = p.get<std::string>(
"Sum Name");
66 Teuchos::RCP<const PureBasis> basis = p.get< Teuchos::RCP<const PureBasis> >(
"Basis");
68 if(p.isType<
bool>(
"Evaluate On Closure"))
71 inField = PHX::MDField<const ScalarT,Cell,BASIS>( inName, basis->functional);
72 outField = PHX::MDField<ScalarT,Cell>( outName, basis->cell_data);
74 this->addDependentField(
inField);
80 std::string n =
"SubcellSum: " +
outField.fieldTag().name();
91 std::vector<int> indices;
96 if(evaluateOnClosure_)
97 fieldPattern_->getSubcellClosureIndices(workset.
subcell_dim,this->wda(workset).subcell_index,indices);
99 indices = fieldPattern_->getSubcellIndices(workset.
subcell_dim,this->wda(workset).subcell_index);
101 auto outField_h = Kokkos::create_mirror_view(outField.get_static_view());
102 auto inField_h = Kokkos::create_mirror_view(inField.get_static_view());
103 Kokkos::deep_copy(inField_h, inField.get_static_view());
104 for(index_t c=0;c<workset.
num_cells;c++) {
108 for(std::size_t i=0;i<indices.size();i++)
109 outField_h(c) += inField_h(c,indices[i]);
114 Kokkos::deep_copy(outField.get_static_view(), outField_h);