43 #ifndef PANZER_DIRICHLET_RESIDUAL_IMPL_HPP
44 #define PANZER_DIRICHLET_RESIDUAL_IMPL_HPP
53 template<
typename EvalT,
typename Traits>
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");
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()
T & get(const std::string &name, T def_value)
DirichletResidual(const Teuchos::ParameterList &p)
void evaluateFields(typename Traits::EvalData d)
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)