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"
56template<
typename EvalT,
typename Traits>
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;
83template<
typename EvalT,
typename Traits>
90 num_qp = normals.extent(1);
91 num_dim = normals.extent(2);
97template<
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);
Normals(const Teuchos::ParameterList &p)
PHX::MDField< ScalarT, Cell, Point, Dim > normals
void postRegistrationSetup(typename Traits::SetupData d, PHX::FieldManager< Traits > &fm)
void evaluateFields(typename Traits::EvalData d)
int num_cells
DEPRECATED - use: numCells()
std::vector< int >::size_type getIntegrationRuleIndex(int ir_degree, const panzer::Workset &workset, WorksetDetailsAccessor &wda)
Teuchos::RCP< const std::vector< panzer::Workset > > worksets_