60 const Teuchos::ParameterList& p)
62 std::string residual_name = p.get<std::string>(
"Residual Name");
63 std::string flux_name = p.get<std::string>(
"Flux Name");
64 std::string normal_name = p.get<std::string>(
"Normal Name");
65 std::string normal_dot_flux_name = normal_name +
" dot " + flux_name;
67 const Teuchos::RCP<const panzer::PureBasis> basis =
68 p.get< Teuchos::RCP<const panzer::PureBasis> >(
"Basis");
70 const Teuchos::RCP<const panzer::IntegrationRule> ir =
71 p.get< Teuchos::RCP<const panzer::IntegrationRule> >(
"IR");
74 residual = PHX::MDField<ScalarT>(residual_name, basis->functional);
75 normal_dot_flux = PHX::MDField<ScalarT>(normal_dot_flux_name, ir->dl_scalar);
76 flux = PHX::MDField<const ScalarT>(flux_name, ir->dl_vector);
77 normal = PHX::MDField<const ScalarT>(normal_name, ir->dl_vector);
79 this->addEvaluatedField(residual);
80 this->addEvaluatedField(normal_dot_flux);
81 this->addDependentField(normal);
82 this->addDependentField(flux);
86 std::string n =
"Interface Residual Evaluator";
114 residual.deep_copy(
ScalarT(0.0));
116 for (index_t cell = 0; cell < workset.
num_cells; ++cell) {
117 for (std::size_t ip = 0; ip < num_ip; ++ip) {
118 normal_dot_flux(cell,ip) =
ScalarT(0.0);
119 for (std::size_t dim = 0; dim < num_dim; ++dim) {
120 normal_dot_flux(cell,ip) += normal(cell,ip,dim) * flux(cell,ip,dim);
126 const Teuchos::RCP<const BasisValues2<double> > bv = this->wda(workset).bases[basis_index];
127 for (index_t cell = 0; cell < workset.
num_cells; ++cell) {
128 for (std::size_t basis = 0; basis < residual.extent(1); ++basis) {
129 for (std::size_t qp = 0; qp < num_ip; ++qp) {
130 residual(cell,basis) += normal_dot_flux(cell,qp)*bv->weighted_basis_scalar(cell,basis,qp);
136 Intrepid2::FunctionSpaceTools<PHX::exec_space>::
137 integrate(residual.get_view(),
138 normal_dot_flux.get_view(),
139 (this->wda(workset).bases[basis_index])->weighted_basis_scalar.get_view());
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.