CellTopology get_cell_topology_new( const Bucket & bucket) { const BulkData & bulk_data = BulkData::get(bucket); const FEMMetaData & fem_meta_data = FEMMetaData::get(bulk_data); const PartVector & all_parts = fem_meta_data.get_parts(); CellTopology cell_topology; const std::pair< const unsigned *, const unsigned * > supersets = bucket.superset_part_ordinals(); if (supersets.first != supersets.second) { const Part *first_found_part = NULL; for ( const unsigned * it = supersets.first ; it != supersets.second ; ++it ) { const Part & part = * all_parts[*it] ; if ( part.primary_entity_rank() == bucket.entity_rank() ) { CellTopology top = fem_meta_data.get_cell_topology( part ); if ( ! cell_topology.getCellTopologyData() ) { cell_topology = top ; if (!first_found_part) first_found_part = ∂ } else { ThrowErrorMsgIf( top.getCellTopologyData() && top != cell_topology, "Cell topology is ambiguously defined. It is defined as " << cell_topology.getName() << " on part " << first_found_part->name() << " and as " << top.getName() << " on its superset part " << part.name() ); } } } } return cell_topology ; }
inline bool is_matching_rank(const FieldBase& f, const Bucket& b) { ThrowAssert(&f.get_mesh() == &b.mesh()); return(b.entity_rank() == static_cast<unsigned>(f.entity_rank())); }
inline unsigned field_scalars_per_entity(const FieldBase& f, const Bucket& b) { ThrowAssert(f.entity_rank() == b.entity_rank()); ThrowAssert(&f.get_mesh() == &b.mesh()); const unsigned bytes_per_scalar = f.data_traits().size_of; return f.get_meta_data_for_field()[b.bucket_id()].m_bytes_per_entity/bytes_per_scalar; }
inline unsigned field_bytes_per_entity(const FieldBase& f, const Bucket& b) { ThrowAssert(f.entity_rank() == b.entity_rank()); ThrowAssert(&f.get_mesh() == &b.mesh()); return f.get_meta_data_for_field()[b.bucket_id()].m_bytes_per_entity; }