58 const Teuchos::ParameterList& p) :
59 residual( p.get<
std::string>(
"Residual Name"),
61 scalar( p.get<
std::string>(
"Value Name"),
66 p.validateParameters(*valid_params);
68 Teuchos::RCP<const PureBasis> basis
69 = p.get< Teuchos::RCP<BasisIRLayout> >(
"Basis")->getBasis();
72 TEUCHOS_TEST_FOR_EXCEPTION(!basis->isScalarBasis(),std::logic_error,
73 "Integrator_TransientBasisTimesScalar: Basis of type \"" << basis->name() <<
"\" is not a "
77 this->addDependentField(
scalar);
82 if (p.isType<Teuchos::RCP<
const std::vector<std::string> > >(
"Field Multipliers")) {
83 const std::vector<std::string>& field_multiplier_names =
84 *(p.get<Teuchos::RCP<const std::vector<std::string> > >(
"Field Multipliers"));
86 for (std::vector<std::string>::const_iterator name =
87 field_multiplier_names.begin();
88 name != field_multiplier_names.end(); ++name) {
89 PHX::MDField<const ScalarT,Cell,IP> tmp_field(*name, p.get< Teuchos::RCP<panzer::IntegrationRule> >(
"IR")->dl_scalar);
94 for (
typename std::vector<PHX::MDField<const ScalarT,Cell,IP> >::iterator
field =
field_multipliers.begin();
96 this->addDependentField(*
field);
98 std::string n =
"Integrator_TransientBasisTimesScalar: " +
residual.fieldTag().name();
110 this->utils.setFieldData(residual,fm);
111 this->utils.setFieldData(scalar,fm);
113 num_nodes = residual.extent(1);
114 num_qp = scalar.extent(1);
118 tmp = Kokkos::createDynRankView(residual.get_static_view(),
"tmp",scalar.extent(0), num_qp);
133 Kokkos::deep_copy (residual.get_static_view(),
ScalarT(0.0));
135 for (index_t cell = 0; cell < workset.
num_cells; ++cell) {
136 for (std::size_t qp = 0; qp < num_qp; ++qp) {
138 for (
typename std::vector<PHX::MDField<const ScalarT,Cell,IP> >::iterator
field = field_multipliers.begin();
140 tmp(cell,qp) = tmp(cell,qp) * (*field)(cell,qp);
145 Intrepid2::FunctionSpaceTools<PHX::exec_space>::
146 integrate<ScalarT>(residual.get_view(),
148 (this->wda(workset).bases[basis_index])->weighted_basis_scalar.get_view());
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.
std::vector< std::string >::size_type getBasisIndex(std::string basis_name, const panzer::Workset &workset, WorksetDetailsAccessor &wda)
Returns the index in the workset bases for a particular BasisIRLayout name.