43 #ifndef PANZER_STK_SCATTER_CELL_AVG_VECTOR_IMPL_HPP 44 #define PANZER_STK_SCATTER_CELL_AVG_VECTOR_IMPL_HPP 46 #include "Teuchos_Assert.hpp" 48 #include "Phalanx_config.hpp" 49 #include "Phalanx_Evaluator_Macros.hpp" 50 #include "Phalanx_MDField.hpp" 51 #include "Phalanx_DataLayout.hpp" 52 #include "Phalanx_DataLayout_MDALayout.hpp" 57 #include "Teuchos_FancyOStream.hpp" 58 #include "Teuchos_ArrayRCP.hpp" 69 std::string scatterName = p.get<std::string>(
"Scatter Name");
71 const std::vector<std::string> & names =
80 for (std::size_t fd = 0; fd < names.size(); ++fd)
87 PHX::Tag<ScalarT> scatterHolder(scatterName,
Teuchos::rcp(
new PHX::MDALayout<panzer::Dummy>(0)));
88 this->addEvaluatedField(scatterHolder);
90 this->setName(scatterName+
": STK-Scatter Cell Vectors");
113 const std::vector<std::size_t> & localCellIds = this->wda(workset).cell_local_ids;
114 std::string blockId = this->wda(workset).block_id;
115 std::string d_mod[3] = {
"X",
"Y",
"Z"};
118 for(std::size_t fieldIndex = 0; fieldIndex <
scatterFields_.size(); fieldIndex++)
120 PHX::MDField<const ScalarT,panzer::Cell,panzer::Point,panzer::Dim> &
field =
scatterFields_[fieldIndex];
121 std::string fieldName =
field.fieldTag().name();
122 int numCells =
field.dimension(0);
124 int numDims =
field.dimension(2);
126 for (
int dim = 0; dim < numDims; dim++)
129 PHX::MDField<double,panzer::Cell,panzer::NODE>
average = af.
buildStaticArray<double,panzer::Cell,panzer::NODE>(
"",numCells,1);
132 for(
int i = 0; i < numCells; i++)
136 average(i,0) += Sacado::ScalarValue<ScalarT>::eval(
field(i,j,dim));
141 mesh_->setCellFieldData(fieldName+d_mod[dim],blockId,localCellIds,
average);
PHX_POST_REGISTRATION_SETUP(ScatterCellAvgQuantity, d, fm)
PHX::MDField< Scalar, T0 > buildStaticArray(const std::string &str, int d0) const
std::vector< PHX::MDField< const ScalarT, panzer::Cell, panzer::Point > > scatterFields_
PHX::MDField< ScalarT, Cell > average
PHX_EVALUATE_FIELDS(ScatterCellAvgQuantity, workset)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Teuchos::RCP< STK_Interface > mesh_
PHX_EVALUATOR_CTOR(ScatterCellAvgQuantity, p)
PHX::MDField< const ScalarT, Cell, IP > field
panzer_stk::STK_Interface::SolutionFieldType VariableField
Teuchos::RCP< PHX::DataLayout > dl_vector
Data layout for vector fields.
std::vector< VariableField * > stkFields_