43 #ifndef PANZER_NORMALS_IMPL_HPP
44 #define PANZER_NORMALS_IMPL_HPP
49 #include "Intrepid2_FunctionSpaceTools.hpp"
50 #include "Intrepid2_CellTools.hpp"
51 #include "Phalanx_Scratch_Utilities.hpp"
56 template<
typename EvalT,
typename Traits>
63 const std::string name = p.
get<std::string>(
"Name");
76 this->addEvaluatedField(
normals);
78 std::string n =
"Normals: " + name;
83 template<
typename EvalT,
typename Traits>
90 num_qp = normals.extent(1);
91 num_dim = normals.extent(2);
97 template<
typename EvalT,
typename Traits>
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);
int num_cells
DEPRECATED - use: numCells()
T & get(const std::string &name, T def_value)
Normals(const Teuchos::ParameterList &p)
std::vector< int >::size_type getIntegrationRuleIndex(int ir_degree, const panzer::Workset &workset, WorksetDetailsAccessor &wda)
bool isParameter(const std::string &name) const
PHX::MDField< ScalarT, Cell, Point, Dim > normals
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)
void evaluateFields(typename Traits::EvalData d)
Teuchos::RCP< const std::vector< panzer::Workset > > worksets_