43 #ifndef __Panzer_CellExtreme_impl_hpp__ 44 #define __Panzer_CellExtreme_impl_hpp__ 48 #include "Intrepid2_FunctionSpaceTools.hpp" 51 #include "Phalanx_DataLayout_MDALayout.hpp" 59 p.validateParameters(*valid_params);
63 if(p.isType<
bool>(
"Use Max"))
64 use_max = p.get<
bool>(
"Use Max");
66 Teuchos::RCP<panzer::IntegrationRule> ir = p.get< Teuchos::RCP<panzer::IntegrationRule> >(
"IR");
69 Teuchos::RCP<PHX::DataLayout> dl_cell = Teuchos::rcp(
new PHX::MDALayout<Cell>(ir->dl_scalar->dimension(0)));
70 extreme = PHX::MDField<ScalarT>( p.get<std::string>(
"Extreme Name"), dl_cell);
71 scalar = PHX::MDField<const ScalarT,Cell,IP>( p.get<std::string>(
"Field Name"), ir->dl_scalar);
73 this->addEvaluatedField(
extreme);
74 this->addDependentField(
scalar);
77 if(p.isType<
double>(
"Multiplier"))
80 if (p.isType<Teuchos::RCP<
const std::vector<std::string> > >(
"Field Multipliers")) {
81 const std::vector<std::string>& field_multiplier_names =
82 *(p.get<Teuchos::RCP<const std::vector<std::string> > >(
"Field Multipliers"));
84 for (std::vector<std::string>::const_iterator name =
85 field_multiplier_names.begin();
86 name != field_multiplier_names.end(); ++name) {
87 PHX::MDField<const ScalarT,Cell,IP> tmp_field(*name, p.get< Teuchos::RCP<panzer::IntegrationRule> >(
"IR")->dl_scalar);
94 this->addDependentField(*
field);
96 std::string n =
"CellExtreme: " +
extreme.fieldTag().name();
103 this->utils.setFieldData(
extreme,fm);
104 this->utils.setFieldData(
scalar,fm);
108 this->utils.setFieldData(*
field,fm);
118 for (index_t cell = 0; cell < workset.num_cells; ++cell) {
120 for (std::size_t qp = 0; qp <
num_qp; ++qp) {
124 current *= (*
field)(cell,qp);
140 template<
typename EvalT,
typename TRAITS>
141 Teuchos::RCP<Teuchos::ParameterList>
144 Teuchos::RCP<Teuchos::ParameterList> p = Teuchos::rcp(
new Teuchos::ParameterList);
145 p->set<std::string>(
"Extreme Name",
"?");
146 p->set<std::string>(
"Field Name",
"?");
148 Teuchos::RCP<panzer::IntegrationRule> ir;
150 p->set<
bool>(
"Use Max",
true);
151 p->set<
double>(
"Multiplier", 1.0);
153 Teuchos::RCP<const std::vector<std::string> > fms;
154 p->set(
"Field Multipliers", fms);
PHX::MDField< ScalarT > vector
Teuchos::RCP< Teuchos::ParameterList > getValidParameters() const
PHX_EVALUATOR_CTOR(BasisValues_Evaluator, p)
PHX::MDField< const ScalarT, Cell, IP > field
PHX::MDField< const ScalarT, Cell, IP > scalar
PHX_EVALUATE_FIELDS(BasisValues_Evaluator, workset)
std::vector< std::string >::size_type getIntegrationRuleIndex(int ir_degree, panzer::Workset &workset, WorksetDetailsAccessor &wda)
std::vector< PHX::MDField< const ScalarT, Cell, IP > > field_multipliers
PHX::MDField< ScalarT > extreme
PHX_POST_REGISTRATION_SETUP(BasisValues_Evaluator, sd, fm)