57 const Teuchos::RCP<const FieldPattern> & geomAggPattern)
68 const Teuchos::RCP<const FieldPattern> & geomAggPattern)
70 TEUCHOS_ASSERT(patterns.size()>0);
73 if(geomAggPattern==Teuchos::null) {
74 std::vector<std::pair<FieldType,FPPtr>> patternVec;
77 auto itr = patterns.cbegin();
78 for(;itr!=patterns.cend();++itr)
79 patternVec.push_back(std::make_pair(std::get<1>(*itr),std::get<2>(*itr)));
97 TEUCHOS_TEST_FOR_EXCEPTION(geomPattern==Teuchos::null,std::logic_error,
98 "Geometric field pattern not yet set, call buildPatterns first");
100 return geomPattern->getDimension();
106 TEUCHOS_TEST_FOR_EXCEPTION(geomPattern==Teuchos::null,std::logic_error,
107 "Geometric field pattern not yet set, call buildPatterns first");
109 return geomPattern->getCellTopology();
124 TEUCHOS_TEST_FOR_EXCEPTION(
true,std::logic_error,
125 "FieldAggPattern::getSubcellClosureIndices should not be called");
132 os <<
"FieldPattern: FieldAggPattern" << std::endl;
133 os <<
"FieldPattern: |numFields| = " <<
numFields_.size() << std::endl;
134 os <<
"FieldPattern: numFields = [ ";
136 for(std::size_t i=0;i<
numFields_.size();i++) {
140 os <<
"]" << std::endl;
141 os <<
"FieldPattern: |fieldIds| = " <<
fieldIds_.size() <<
" (" << total <<
")" << std::endl;
142 os <<
"FieldPattern: fieldIds = [ ";
143 for(std::size_t i=0;i<
fieldIds_.size();i++)
145 os <<
"]" << std::endl;
146 os <<
"FieldPattern: local offsets\n";
148 std::map<int,int>::const_iterator itr;
150 int fieldId = itr->first;
152 os <<
"FieldPattern: field " << itr->first <<
" = [ ";
153 for(std::size_t i=0;i<
offsets.size();i++)
155 os <<
"]" << std::endl;
163 "FieldID = " << fieldId <<
" not defined in this pattern");
165 return std::get<2>(
patterns_[idxIter->second]);
172 "FieldID = " << fieldId <<
" not defined in this pattern");
174 return std::get<1>(
patterns_[idxIter->second]);
182 for(;itr!=
patterns_.cend();++itr,++index)
192 numFields_.resize(geomAggPattern->numberIds(),0);
212 for(
int dim=0;dim<numDimensions;dim++) {
213 int numSubcell = geomAggPattern->getSubcellCount(dim);
214 for(
int subcell=0;subcell<numSubcell;subcell++) {
219 const std::vector<int> * geomIndices =
nullptr;
225 if (geomIndices->size() > 0) {
226 const int geomIndex = (*geomIndices)[0];
231 if (std::get<1>(*itr) == fieldType) {
232 const std::size_t fieldSize = std::get<2>(*itr)->getSubcellIndices(dim,subcell).size();
235 for (std::size_t i=0;i<fieldSize;++i)
257 for(
int d=0;d<dimension+1;d++) {
258 int numSubcell = geomIdsPattern->getSubcellCount(d);
262 for(
int sc=0;sc<numSubcell;sc++) {
265 const std::vector<int> & geomIds = geomIdsPattern->getSubcellIndices(d,sc);
266 TEUCHOS_ASSERT(geomIds.size()<=1);
267 if (geomIds.size() > 0) {
268 const int geomId = geomIds[0];
280 std::map<int,std::vector<int> >::const_iterator itr =
fieldOffsets_.find(fieldId);
296 const auto hostOffsetsStdVector = this->
localOffsets(fieldId);
297 PHX::View<int*>
offsets(
"panzer::FieldAggPattern::localOffsetsKokkos",hostOffsetsStdVector.size());
298 auto hostOffsets = Kokkos::create_mirror_view(
offsets);
299 for (
size_t i=0; i < hostOffsetsStdVector.size(); ++i)
300 hostOffsets(i) = hostOffsetsStdVector[i];
301 Kokkos::deep_copy(
offsets,hostOffsets);
308 if(a[0] < b[0])
return true;
309 if(a[0] > b[0])
return false;
312 if(a[1] < b[1])
return true;
313 if(a[1] > b[1])
return false;
316 if(a[2] < b[2])
return true;
317 if(a[2] > b[2])
return false;
324const std::pair<std::vector<int>,std::vector<int> > &
328 typedef std::map<Teuchos::Tuple<int,3>, std::pair<std::vector<int>,std::vector<int> >,
LessThan> OffsetMap;
330 Teuchos::Tuple<int,3> subcellTuple = Teuchos::tuple<int>(fieldId,subcellDim,subcellId);
332 OffsetMap::const_iterator itr
338 TEUCHOS_TEST_FOR_EXCEPTION(subcellDim >=
getDimension(),std::logic_error,
339 "FieldAggPattern::localOffsets_closure precondition subcellDim<getDimension() failed");
340 TEUCHOS_TEST_FOR_EXCEPTION(subcellId < 0,std::logic_error,
341 "FieldAggPattern::localOffsets_closure precondition subcellId>=0 failed");
342 TEUCHOS_TEST_FOR_EXCEPTION(subcellId>=
getSubcellCount(subcellDim),std::logic_error,
343 "FieldAggPattern::localOffsets_closure precondition subcellId<getSubcellCount(subcellDim) failed");
347 const std::vector<int> & fieldOffsets =
localOffsets(fieldId);
350 std::vector<int> closureOffsets;
352 fieldPattern->getSubcellClosureIndices(subcellDim,subcellId,closureOffsets);
355 std::pair<std::vector<int>,std::vector<int> > & indicesPair
358 std::vector<int> & closureIndices = indicesPair.first;
359 for(std::size_t i=0;i<closureOffsets.size();i++)
360 closureIndices.push_back(fieldOffsets[closureOffsets[i]]);
362 std::vector<int> & basisIndices = indicesPair.second;
363 basisIndices.assign(closureOffsets.begin(),closureOffsets.end());
384 offsets.resize(fieldPattern->numberIds(),-111111);
388 std::vector<int> fieldIdsGeomOrder;
389 for(std::size_t i=0;i<
fieldIds_.size();++i) {
391 fieldIdsGeomOrder.push_back(i);
395 TEUCHOS_ASSERT((
int) fieldIdsGeomOrder.size()==fieldPattern->numberIds());
402 const std::vector<int> & fIndices = fieldPattern->getSubcellIndices(dim,sc);
404 for(std::size_t i=0;i<fIndices.size();i++)
405 offsets[fIndices[i]] = fieldIdsGeomOrder[cnt++];
410 for(std::size_t i=0;i<
offsets.size();i++) {
PHX::View< const int * > offsets
void buildFieldIdToPatternIdx()
void buildFieldIdsVector()
std::map< int, int > fieldIdToPatternIdx_
virtual void buildPattern(const std::vector< std::tuple< int, panzer::FieldType, Teuchos::RCP< const FieldPattern > > > &patterns, const Teuchos::RCP< const FieldPattern > &geomAggPattern=Teuchos::null)
const std::pair< std::vector< int >, std::vector< int > > & localOffsets_closure(int fieldId, int subcellDim, int subcellId) const
std::vector< int > numFields_
virtual void getSubcellClosureIndices(int, int, std::vector< int > &) const
virtual Teuchos::RCP< const FieldPattern > getFieldPattern(int fieldId) const
void mergeFieldPatterns(const FieldType &fieldType)
virtual int getDimension() const
Teuchos::RCP< const FieldPattern > geomAggPattern_
void buildFieldPatternData()
std::map< Teuchos::Tuple< int, 3 >, std::pair< std::vector< int >, std::vector< int > >, LessThan > fieldSubcellOffsets_closure_
const std::vector< int > & localOffsets(int fieldId) const
std::vector< std::tuple< int, panzer::FieldType, Teuchos::RCP< const FieldPattern > > > patterns_
virtual const std::vector< int > & getSubcellIndices(int dimension, int subcell) const
void localOffsets_build(int fieldId, std::vector< int > &offsets) const
virtual void print(std::ostream &os) const
Print this pattern.
Teuchos::RCP< const FieldPattern > FPPtr
std::map< int, PHX::View< int * > > fieldOffsetsKokkos_
Stores the Field offsets for the fieldId key. Note that the key is the fieldId, not the index into th...
std::vector< std::vector< std::vector< int > > > patternData_
virtual Teuchos::RCP< const FieldPattern > getGeometricAggFieldPattern() const
virtual int getSubcellCount(int dimension) const
std::vector< int > fieldIds_
virtual FieldType getFieldType(int fieldId) const
std::map< int, std::vector< int > > fieldOffsets_
Stores the Field offsets for the fieldId key. Note that the key is the fieldId, not the index into th...
virtual shards::CellTopology getCellTopology() const
const PHX::View< const int * > localOffsetsKokkos(int fieldId) const
virtual void print(std::ostream &os) const
virtual int getSubcellCount(int dim) const
virtual int getDimension() const
FieldType
The type of discretization to use for a field pattern.
@ DG
Continuous Galerkin Formulation.
bool operator()(const Teuchos::Tuple< int, 3 > &a, const Teuchos::Tuple< int, 3 > &b) const