43 #ifndef __Panzer_CellAverage_impl_hpp__ 44 #define __Panzer_CellAverage_impl_hpp__ 48 #include "Intrepid2_FunctionSpaceTools.hpp" 51 #include "Phalanx_DataLayout_MDALayout.hpp" 56 template<
typename EvalT,
typename Traits>
59 const Teuchos::ParameterList& p) : quad_index(0)
62 p.validateParameters(*valid_params);
64 Teuchos::RCP<panzer::IntegrationRule> ir = p.get< Teuchos::RCP<panzer::IntegrationRule> >(
"IR");
67 Teuchos::RCP<PHX::DataLayout> dl_cell = Teuchos::rcp(
new PHX::MDALayout<Cell>(ir->dl_scalar->extent(0)));
68 average = PHX::MDField<ScalarT,Cell>( p.get<std::string>(
"Average Name"), dl_cell);
69 scalar = PHX::MDField<const ScalarT,Cell,IP>( p.get<std::string>(
"Field Name"), ir->dl_scalar);
71 this->addEvaluatedField(
average);
72 this->addDependentField(
scalar);
75 if(p.isType<
double>(
"Multiplier"))
78 if (p.isType<Teuchos::RCP<
const std::vector<std::string> > >(
"Field Multipliers")) {
79 const std::vector<std::string>& field_multiplier_names =
80 *(p.get<Teuchos::RCP<const std::vector<std::string> > >(
"Field Multipliers"));
82 for (std::vector<std::string>::const_iterator name =
83 field_multiplier_names.begin();
84 name != field_multiplier_names.end(); ++name) {
85 PHX::MDField<const ScalarT,Cell,IP> tmp_field(*name, p.get< Teuchos::RCP<panzer::IntegrationRule> >(
"IR")->dl_scalar);
90 for (
typename std::vector<PHX::MDField<const ScalarT,Cell,IP> >::iterator
field =
field_multipliers.begin();
92 this->addDependentField(*
field);
94 std::string n =
"CellAverage: " +
average.fieldTag().name();
99 template<
typename EvalT,
typename Traits>
106 num_qp = scalar.extent(1);
111 template<
typename EvalT,
typename Traits>
117 for (index_t cell = 0; cell < workset.
num_cells; ++cell) {
123 for (std::size_t qp = 0; qp < num_qp; ++qp) {
125 for (
typename std::vector<PHX::MDField<const ScalarT,Cell,IP> >::iterator
field = field_multipliers.begin();
127 current *= (*
field)(cell,qp);
130 average(cell) += current/num_qp;
136 template<
typename EvalT,
typename TRAITS>
137 Teuchos::RCP<Teuchos::ParameterList>
140 Teuchos::RCP<Teuchos::ParameterList> p = Teuchos::rcp(
new Teuchos::ParameterList);
141 p->set<std::string>(
"Average Name",
"?");
142 p->set<std::string>(
"Field Name",
"?");
144 Teuchos::RCP<panzer::IntegrationRule> ir;
146 p->set<
double>(
"Multiplier", 1.0);
148 Teuchos::RCP<const std::vector<std::string> > fms;
149 p->set(
"Field Multipliers", fms);
Teuchos::RCP< Teuchos::ParameterList > getValidParameters() const
int num_cells
DEPRECATED - use: numCells()
typename EvalT::ScalarT ScalarT
std::vector< int >::size_type getIntegrationRuleIndex(int ir_degree, const panzer::Workset &workset, WorksetDetailsAccessor &wda)
double multiplier
The scalar multiplier out in front of the integral ( ).
PHX::MDField< ScalarT, Cell > average
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)
void evaluateFields(typename Traits::EvalData d)
PHX::MDField< const ScalarT, Cell, IP > scalar
PHX::MDField< ScalarT, panzer::Cell, panzer::BASIS > field
A field to which we'll contribute, or in which we'll store, the result of computing this integral...
CellAverage(const Teuchos::ParameterList &p)
std::vector< PHX::MDField< const ScalarT, Cell, IP > > field_multipliers
Teuchos::RCP< const std::vector< panzer::Workset > > worksets_