43 #ifndef PANZER_DIRICHLET_RESIDUAL_IMPL_HPP 44 #define PANZER_DIRICHLET_RESIDUAL_IMPL_HPP 53 template<
typename EvalT,
typename Traits>
56 const Teuchos::ParameterList& p)
58 std::string residual_name = p.get<std::string>(
"Residual Name");
59 std::string dof_name = p.get<std::string>(
"DOF Name");
60 std::string value_name = p.get<std::string>(
"Value Name");
62 Teuchos::RCP<PHX::DataLayout> data_layout =
63 p.get< Teuchos::RCP<PHX::DataLayout> >(
"Data Layout");
65 residual = PHX::MDField<ScalarT>(residual_name, data_layout);
66 dof = PHX::MDField<const ScalarT>(dof_name, data_layout);
67 value = PHX::MDField<const ScalarT>(value_name, data_layout);
69 this->addEvaluatedField(residual);
70 this->addDependentField(dof);
71 this->addDependentField(value);
73 std::string n =
"Dirichlet Residual Evaluator";
78 template<
typename EvalT,
typename Traits>
85 cell_data_size = residual.fieldTag().dataLayout().extent(1);
89 template<
typename EvalT,
typename Traits>
95 auto residual_v = residual.get_static_view();
96 auto dof_v = dof.get_static_view();
97 auto value_v = value.get_static_view();
98 auto local_cell_data_size = cell_data_size;
99 Kokkos::parallel_for (workset.
num_cells, KOKKOS_LAMBDA (index_t i) {
100 for (std::size_t j = 0; j < local_cell_data_size; ++j)
101 residual_v(i,j)=dof_v(i,j)-value_v(i,j);
int num_cells
DEPRECATED - use: numCells()
DirichletResidual(const Teuchos::ParameterList &p)
void evaluateFields(typename Traits::EvalData d)
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)