43#include "PanzerDiscFE_config.hpp"
51Teuchos::RCP<std::vector<panzer::Workset> >
53 const std::string & elementBlock,
54 const std::vector<std::size_t>& local_cell_ids,
55 const Kokkos::DynRankView<double,PHX::Device>& vertex_coordinates);
58Teuchos::RCP<std::map<unsigned,panzer::Workset> >
60 const std::string& elementBlock,
61 const std::vector<std::size_t>& local_cell_ids,
62 const std::vector<std::size_t>& local_side_ids,
63 const Kokkos::DynRankView<double,PHX::Device>& vertex_coordinates,
64 const bool populate_value_arrays);
67Teuchos::RCP<std::map<unsigned,panzer::Workset> >
69 const std::string & blockid_a,
70 const std::vector<std::size_t>& local_cell_ids_a,
71 const std::vector<std::size_t>& local_side_ids_a,
72 const Kokkos::DynRankView<double,PHX::Device> & vertex_coordinates_a,
74 const std::string & blockid_b,
75 const std::vector<std::size_t>& local_cell_ids_b,
76 const std::vector<std::size_t>& local_side_ids_b,
77 const Kokkos::DynRankView<double,PHX::Device> & vertex_coordinates_b);
80Teuchos::RCP<std::vector<panzer::Workset> >
82 const std::string & eblock_a,
83 const std::vector<std::size_t>& local_cell_ids_a,
84 const std::vector<std::size_t>& local_side_ids_a,
85 const Kokkos::DynRankView<double,PHX::Device> & vertex_coordinates_a,
87 const std::string & eblock_b,
88 const std::vector<std::size_t>& local_cell_ids_b,
89 const std::vector<std::size_t>& local_side_ids_b,
90 const Kokkos::DynRankView<double,PHX::Device> & vertex_coordinates_b);
103 if(
details.numOwnedCells()==-1)
104 details.setNumberOfCells(num_cells,0,0);
105 if(other_details!=Teuchos::null and other_details->numOwnedCells()==-1)
106 other_details->setNumberOfCells(num_cells,0,0);
112 std::vector<RCP<const panzer::PureBasis> > bases;
113 std::vector<RCP<const panzer::IntegrationRule> > int_rules;
119 for(std::size_t i=0;i<needs.
int_rules.size();i++)
122 for(std::size_t i=0;i<needs.
bases.size();i++)
123 bases.push_back(rcp(
new PureBasis(needs.
bases[i]->type(),needs.
bases[i]->order(),side_cell_data)));
130 details.ir_degrees = rcp(
new std::vector<int>(0));
131 details.basis_names = rcp(
new std::vector<std::string>(0));
133 for(std::size_t i=0;i<int_rules.size();i++) {
135 details.ir_degrees->push_back(int_rules[i]->cubature_degree);
137 RCP<panzer::IntegrationValues2<double> > iv2 =
139 iv2->setupArrays(int_rules[i]);
140 if (Teuchos::nonnull(other_details))
141 iv2->evaluateValues(
details.cell_vertex_coordinates, other_details->int_rules[i]->ip_coordinates,num_cells);
143 iv2->evaluateValues(
details.cell_vertex_coordinates,num_cells);
145 details.int_rules.push_back(iv2);
148 for(std::size_t b=0;b<bases.size();b++) {
150 details.basis_names->push_back(b_layout->name());
152 std::size_t int_degree_index = std::distance(
details.ir_degrees->begin(),
153 std::find(
details.ir_degrees->begin(),
155 int_rules[i]->order()));
156 RCP<panzer::BasisValues2<double> > bv2 =
158 bv2->setupArrays(b_layout);
159 bv2->evaluateValues(
details.int_rules[int_degree_index]->cub_points,
160 details.int_rules[int_degree_index]->jac,
161 details.int_rules[int_degree_index]->jac_det,
162 details.int_rules[int_degree_index]->jac_inv,
163 details.int_rules[int_degree_index]->weighted_measure,
164 details.cell_vertex_coordinates,
Data for determining cell topology and dimensionality.
Teuchos::RCP< const shards::CellTopology > getCellTopology() const
Get CellTopology for the base cell.
Description and data layouts associated with a particular basis.
void populateValueArrays(std::size_t num_cells, bool isSide, const WorksetNeeds &needs, WorksetDetails &details, const Teuchos::RCP< WorksetDetails > other_details)
Teuchos::RCP< std::vector< Workset > > buildWorksets(const WorksetNeeds &needs, const std::string &elementBlock, const std::vector< std::size_t > &local_cell_ids, const ArrayT &vertex_coordinates)
Teuchos::RCP< std::map< unsigned, Workset > > buildBCWorkset(const WorksetNeeds &needs, const std::string &elementBlock, const std::vector< std::size_t > &local_cell_ids, const std::vector< std::size_t > &local_side_ids, const ArrayT &vertex_coordinates, const bool populate_value_arrays=true)
Teuchos::RCP< std::vector< Workset > > buildEdgeWorksets(const WorksetNeeds &needs_a, const std::string &eblock_a, const std::vector< std::size_t > &local_cell_ids_a, const std::vector< std::size_t > &local_side_ids_a, const ArrayT &vertex_coordinates_a, const WorksetNeeds &needs_b, const std::string &eblock_b, const std::vector< std::size_t > &local_cell_ids_b, const std::vector< std::size_t > &local_side_ids_b, const ArrayT &vertex_coordinates_b)
std::vector< Teuchos::RCP< const IntegrationRule > > int_rules
std::vector< Teuchos::RCP< const PureBasis > > bases