43 #ifndef PANZER_NEUMANN_RESIDUAL_IMPL_HPP 44 #define PANZER_NEUMANN_RESIDUAL_IMPL_HPP 51 #include "Intrepid2_FunctionSpaceTools.hpp" 52 #include "Teuchos_RCP.hpp" 57 template<
typename EvalT,
typename Traits>
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 =
"Neumann Residual Evaluator";
91 template<
typename EvalT,
typename Traits>
98 num_ip = flux.extent(1);
99 num_dim = flux.extent(2);
101 TEUCHOS_ASSERT(flux.extent(1) == normal.extent(1));
102 TEUCHOS_ASSERT(flux.extent(2) == normal.extent(2));
108 template<
typename EvalT,
typename Traits>
114 residual.deep_copy(
ScalarT(0.0));
115 auto normal_dot_flux_v = normal_dot_flux.get_static_view();
116 auto normal_v = normal.get_static_view();
117 auto flux_v = flux.get_static_view();
119 std::size_t l_num_ip = num_ip, l_num_dim = num_dim;
120 Kokkos::parallel_for(
"NeumannResidual", workset.
num_cells, KOKKOS_LAMBDA (
const index_t cell) {
121 for (std::size_t ip = 0; ip < l_num_ip; ++ip) {
122 normal_dot_flux_v(cell,ip) =
ScalarT(0.0);
123 for (std::size_t dim = 0; dim < l_num_dim; ++dim) {
124 normal_dot_flux_v(cell,ip) += normal_v(cell,ip,dim) * flux_v(cell,ip,dim);
130 const Teuchos::RCP<const BasisValues2<double> > bv = this->wda(workset).bases[basis_index];
131 auto weighted_basis_scalar = this->wda(workset).bases[basis_index]->weighted_basis_scalar.get_static_view();
132 auto residual_v = residual.get_static_view();
133 Kokkos::parallel_for(
"NeumannResidual", workset.
num_cells, KOKKOS_LAMBDA (
const index_t cell) {
134 for (std::size_t basis = 0; basis < residual_v.extent(1); ++basis) {
135 for (std::size_t qp = 0; qp < l_num_ip; ++qp) {
136 residual_v(cell,basis) += normal_dot_flux_v(cell,qp)*weighted_basis_scalar(cell,basis,qp);
142 Intrepid2::FunctionSpaceTools<PHX::exec_space>::
143 integrate<ScalarT>(residual.get_view(),
144 normal_dot_flux.get_view(),
145 (this->wda(workset).bases[basis_index])->weighted_basis_scalar.get_view());
int num_cells
DEPRECATED - use: numCells()
void evaluateFields(typename Traits::EvalData d)
Teuchos::RCP< panzer::BasisIRLayout > basisIRLayout(std::string basis_type, const int basis_order, const PointRule &pt_rule)
Nonmember constructor.
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)
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.
typename EvalT::ScalarT ScalarT
NeumannResidual(const Teuchos::ParameterList &p)
Teuchos::RCP< const std::vector< panzer::Workset > > worksets_