59 const Teuchos::ParameterList& p)
63 const std::string name = p.get<std::string>(
"Name");
64 side_id = p.get<
int>(
"Side ID");
65 Teuchos::RCP<panzer::IntegrationRule> quadRule
66 = p.get< Teuchos::RCP<panzer::IntegrationRule> >(
"IR");
67 if(p.isParameter(
"Normalize"))
71 Teuchos::RCP<PHX::DataLayout> vector_dl = quadRule->dl_vector;
75 normals = PHX::MDField<ScalarT,Cell,Point,Dim>(name, vector_dl);
76 this->addEvaluatedField(
normals);
78 std::string n =
"Normals: " + name;
106 Intrepid2::CellTools<PHX::exec_space>::getPhysicalSideNormals(normals.get_view(),
107 this->wda(workset).int_rules[quad_index]->jac.get_view(),
108 side_id, *this->wda(workset).int_rules[quad_index]->int_rule->topology);
113 auto local_normals = normals;
114 auto local_num_qp = num_qp;
115 auto local_num_dim = num_dim;
117 if (Sacado::IsADType<ScalarT>::value) {
118 scratch = PHX::View<ScalarT*>(
"normals:scratch",workset.
num_cells,PHX::getFadSize(local_normals.get_static_view()));
121 scratch = PHX::View<ScalarT*>(
"normals:scratch",workset.
num_cells);
125 Kokkos::parallel_for(
"normalize", workset.
num_cells, KOKKOS_LAMBDA (index_t c) {
126 for(std::size_t q=0;q<local_num_qp;q++) {
130 for(std::size_t d=0;d<local_num_dim;d++)
131 norm(c) += local_normals(c,q,d)*local_normals(c,q,d);
134 norm(c) = sqrt(norm(c));
135 for(std::size_t d=0;d<local_num_dim;d++)
136 local_normals(c,q,d) /= norm(c);