43 #ifndef PANZER_EVALUATOR_SCALAR_IMPL_HPP 44 #define PANZER_EVALUATOR_SCALAR_IMPL_HPP 46 #include "Intrepid2_FunctionSpaceTools.hpp" 49 #include "Kokkos_ViewFactory.hpp" 50 #include "Phalanx_DataLayout_MDALayout.hpp" 58 p.validateParameters(*valid_params);
64 integral = PHX::MDField<ScalarT>( p.get<std::string>(
"Integral Name"), dl_cell);
65 scalar = PHX::MDField<ScalarT,Cell,IP>( p.get<std::string>(
"Integrand Name"), ir->
dl_scalar);
68 this->addDependentField(
scalar);
71 if(p.isType<
double>(
"Multiplier"))
74 if (p.isType<
Teuchos::RCP<
const std::vector<std::string> > >(
"Field Multipliers")) {
75 const std::vector<std::string>& field_multiplier_names =
78 for (std::vector<std::string>::const_iterator name =
79 field_multiplier_names.begin();
80 name != field_multiplier_names.end(); ++name) {
88 this->addDependentField(*
field);
90 std::string n =
"Integrator_Scalar: " +
integral.fieldTag().name();
97 this->utils.setFieldData(
integral,fm);
98 this->utils.setFieldData(
scalar,fm);
102 this->utils.setFieldData(*
field,fm);
119 for (index_t cell = 0; cell < workset.num_cells; ++cell) {
120 for (std::size_t qp = 0; qp <
num_qp; ++qp) {
124 tmp(cell,qp) =
tmp(cell,qp) * (*field)(cell,qp);
149 for(index_t cl = 0; cl < workset.num_cells; cl++) {
158 template<
typename EvalT,
typename TRAITS>
163 p->
set<std::string>(
"Integral Name",
"?");
164 p->
set<std::string>(
"Integrand Name",
"?");
168 p->
set<
double>(
"Multiplier", 1.0);
171 p->
set(
"Field Multipliers", fms);
PHX::MDField< ScalarT, Cell > tmp
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
PHX::MDField< ScalarT > vector
Kokkos::DynRankView< typename InputArray::value_type, PHX::Device > createDynRankView(const InputArray &a, const std::string &name, const DimensionPack... dims)
Wrapper to simplify Panzer use of Sacado ViewFactory.
Teuchos::RCP< Teuchos::ParameterList > getValidParameters() const
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
Teuchos::RCP< PHX::DataLayout > dl_scalar
Data layout for scalar fields.
PHX::MDField< ScalarT > integral
PHX_EVALUATOR_CTOR(BasisValues_Evaluator, p)
PHX::MDField< const ScalarT, Cell, IP > field
Array_CellIP weighted_measure
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_POST_REGISTRATION_SETUP(BasisValues_Evaluator, sd, fm)