45 #include <PanzerAdaptersSTK_config.hpp>
50 namespace panzer_stk {
66 PANZER_FUNC_TIME_MONITOR(
"panzer::CustomMeshFactory::buildMesh()");
72 mesh->initialize(parallelMach);
83 PANZER_FUNC_TIME_MONITOR(
"panzer::CustomMeshFactory::buildUncomittedMesh()");
87 machRank_ = stk::parallel_machine_rank(parallelMach);
88 machSize_ = stk::parallel_machine_size(parallelMach);
101 stk::ParallelMachine parallelMach)
const
103 PANZER_FUNC_TIME_MONITOR(
"panzer::CustomMeshFactory::completeMeshConstruction()");
166 if(defaultParams == Teuchos::null) {
169 defaultParams->
set<
int>(
"Dimension",3);
171 defaultParams->
set<
int>(
"NumBlocks",0);
173 defaultParams->
set<
int>(
"NumNodesPerProc",0);
174 defaultParams->
set<
int*>(
"Nodes",NULL);
176 defaultParams->
set<
double*>(
"Coords",NULL);
178 defaultParams->
set<
int>(
"NumElementsPerProc",0);
180 defaultParams->
set<
int*>(
"BlockIDs",NULL);
181 defaultParams->
set<
int*>(
"Element2Nodes",NULL);
183 defaultParams->
set<
int>(
"OffsetToGlobalElementIDs", 0);
185 defaultParams->
set<
double*>(
"ChargeDensity",NULL);
186 defaultParams->
set<
double*>(
"ElectricPotential",NULL);
189 bcs.
set<
int>(
"Count",0);
192 return defaultParams;
208 typedef shards::Hexahedron<8> HexTopo;
209 const CellTopologyData *ctd = shards::getCellTopologyData<HexTopo>();
210 const CellTopologyData *side_ctd = shards::CellTopology(ctd).getBaseCellTopologyData(2,0);
213 std::stringstream block_id;
214 block_id <<
"eblock-" << blk;
241 std::vector<double> coords(dim,0.0);
243 for (
int k=0;k<dim;++k)
249 std::vector<std::string> block_ids;
255 std::stringstream block_id;
262 std::vector<stk::mesh::EntityId> elt2nodes(8);
264 for (
int k=0;k<8;++k)
278 stk::mesh::BulkData& bulkData = *mesh.
getBulkData();
279 const stk::mesh::EntityRank sideRank = mesh.
getSideRank();
282 stk::mesh::Part *box[6];
290 stk::mesh::Part *wall = mesh.
getSideset(
"wall");
292 std::vector<stk::mesh::Entity> elements;
296 for (std::vector<stk::mesh::Entity>::const_iterator
297 itr=elements.begin();itr!=elements.end();++itr) {
298 stk::mesh::Entity element = (*itr);
299 const size_t numSides = bulkData.num_connectivity(element, sideRank);
300 stk::mesh::Entity
const* relations = bulkData.begin(element, sideRank);
303 for (std::size_t i=0;i<numSides;++i) {
304 stk::mesh::Entity side = relations[i];
305 const size_t numNeighbors = bulkData.num_elements(side);
306 stk::mesh::Entity
const* neighbors = bulkData.begin_elements(side);
308 if (numNeighbors == 1) {
312 else if (numNeighbors == 2) {
330 constexpr std::size_t dim_1 = 8;
334 std::stringstream block_id;
335 block_id <<
"eblock-" << blk;
338 std::vector<stk::mesh::Entity> elements;
342 const auto n_elements = elements.size();
345 std::vector<std::size_t> local_ids;
350 for (
const auto& elem : elements){
354 for (std::size_t k=0;k<dim_1;++k) {
355 const auto loc = q*dim_1 + k;
365 charge_density_by_local_ids, 1.0);
370 electric_potential_by_local_ids, 1.0);
unsigned entityOwnerRank(stk::mesh::Entity entity) const
void setSolutionFieldData(const std::string &fieldName, const std::string &blockId, const std::vector< std::size_t > &localElementIds, const ArrayT &solutionValues, double scaleValue=1.0)
void setParameterList(const Teuchos::RCP< Teuchos::ParameterList > ¶mList)
From ParameterListAcceptor.
static void parsePeriodicBCList(const Teuchos::RCP< Teuchos::ParameterList > &pl, std::vector< Teuchos::RCP< const PeriodicBC_MatcherBase > > &periodicBC, bool &useBBoxSearch)
void getElementBlockNames(std::vector< std::string > &names) const
T & get(const std::string &name, T def_value)
stk::mesh::Part * getElementBlockPart(const std::string &name) const
get the block part
ParameterList & set(std::string const &name, T const &value, std::string const &docString="", RCP< const ParameterEntryValidator > const &validator=null)
void addEntityToSideset(stk::mesh::Entity entity, stk::mesh::Part *sideset)
void addSolutionField(const std::string &fieldName, const std::string &blockId)
void initializeWithDefaults()
CustomMeshFactory()
Constructor.
stk::mesh::Part * getSideset(const std::string &name) const
Teuchos::RCP< const Teuchos::ParameterList > getValidParameters() const
From ParameterListAcceptor.
Teuchos::RCP< STK_Interface > buildMesh(stk::ParallelMachine parallelMach) const
Build the mesh object.
std::size_t elementLocalId(stk::mesh::Entity elmt) const
void addSideSets(STK_Interface &mesh) const
unsigned getDimension() const
get the dimension
void addElement(const Teuchos::RCP< ElementDescriptor > &ed, stk::mesh::Part *block)
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)
stk::mesh::EntityRank getSideRank() const
void initialize(stk::ParallelMachine parallelMach, bool setupIO=true, const bool buildRefinementSupport=false)
virtual void completeMeshConstruction(STK_Interface &mesh, stk::ParallelMachine parallelMach) const
Teuchos::RCP< stk::mesh::BulkData > getBulkData() const
virtual ~CustomMeshFactory()
Destructor.
stk::mesh::EntityId elementGlobalId(std::size_t lid) const
void setMyParamList(const RCP< ParameterList > ¶mList)
void validateParametersAndSetDefaults(ParameterList const &validParamList, int const depth=1000)
double * ElectricPotential_
std::string containingBlockId(stk::mesh::Entity elmt) const
void addNode(stk::mesh::EntityId gid, const std::vector< double > &coord)
void buildSubcells()
force the mesh to build subcells: edges and faces
bool isInitialized() const
Has initialize been called on this mesh object?
void fillSolutionFieldData(STK_Interface &mesh) const
void buildElements(STK_Interface &mesh) const
int OffsetToGlobalElementIDs_
void addSideset(const std::string &name, const CellTopologyData *ctData)
std::vector< Teuchos::RCP< const PeriodicBC_MatcherBase > > periodicBCVec_
void buildLocalElementIDs()
ParameterList & sublist(const std::string &name, bool mustAlreadyExist=false, const std::string &docString="")
void getMyElements(std::vector< stk::mesh::Entity > &elements) const
virtual Teuchos::RCP< STK_Interface > buildUncommitedMesh(stk::ParallelMachine parallelMach) const
void buildMetaData(STK_Interface &mesh) const
void addElementBlock(const std::string &name, const CellTopologyData *ctData)