本文整理汇总了C++中typenameAbstractMesh::rGetContainingElementIndices方法的典型用法代码示例。如果您正苦于以下问题：C++ typenameAbstractMesh::rGetContainingElementIndices方法的具体用法？C++ typenameAbstractMesh::rGetContainingElementIndices怎么用？C++ typenameAbstractMesh::rGetContainingElementIndices使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类typenameAbstractMesh的用法示例。
在下文中一共展示了typenameAbstractMesh::rGetContainingElementIndices方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: GetNumNodes
void PottsBasedCellPopulation<DIM>::WriteVtkResultsToFile(const std::string& rDirectory)
{
#ifdef CHASTE_VTK
unsigned num_timesteps = SimulationTime::Instance()->GetTimeStepsElapsed();
std::stringstream time;
time << num_timesteps;
// Create mesh writer for VTK output
VtkMeshWriter<DIM, DIM> mesh_writer(rDirectory, "results_"+time.str(), false);
// Iterate over any cell writers that are present
unsigned num_nodes = GetNumNodes();
for (typename std::set<boost::shared_ptr<AbstractCellWriter<DIM, DIM> > >::iterator cell_writer_iter = this->mCellWriters.begin();
cell_writer_iter != this->mCellWriters.end();
++cell_writer_iter)
{
// Create vector to store VTK cell data
std::vector<double> vtk_cell_data(num_nodes);
// Iterate over nodes in the mesh
for (typename AbstractMesh<DIM,DIM>::NodeIterator iter = mpPottsMesh->GetNodeIteratorBegin();
iter != mpPottsMesh->GetNodeIteratorEnd();
++iter)
{
// Get the index of this node in the mesh and those elements (i.e. cells) that contain this node
unsigned node_index = iter->GetIndex();
std::set<unsigned> element_indices = iter->rGetContainingElementIndices();
// If there are no elements associated with this node, then we set the value of any VTK cell data to be -1 at this node...
if (element_indices.empty())
{
// Populate the vector of VTK cell data
vtk_cell_data[node_index] = -1.0;
}
else
{
// ... otherwise there should be exactly one element (i.e. cell) containing this node
assert(element_indices.size() == 1);
unsigned elem_index = *(element_indices.begin());
CellPtr p_cell = this->GetCellUsingLocationIndex(elem_index);
// Populate the vector of VTK cell data
vtk_cell_data[node_index] = (*cell_writer_iter)->GetCellDataForVtkOutput(p_cell, this);
}
}
mesh_writer.AddPointData((*cell_writer_iter)->GetVtkCellDataName(), vtk_cell_data);
}
// When outputting any CellData, we assume that the first cell is representative of all cells
unsigned num_cell_data_items = this->Begin()->GetCellData()->GetNumItems();
std::vector<std::string> cell_data_names = this->Begin()->GetCellData()->GetKeys();
std::vector<std::vector<double> > cell_data;
for (unsigned var=0; var<num_cell_data_items; var++)
{
std::vector<double> cell_data_var(num_nodes);
cell_data.push_back(cell_data_var);
}
for (typename AbstractMesh<DIM,DIM>::NodeIterator iter = mpPottsMesh->GetNodeIteratorBegin();
iter != mpPottsMesh->GetNodeIteratorEnd();
++iter)
{
// Get the index of this node in the mesh and those elements (i.e. cells) that contain this node
unsigned node_index = iter->GetIndex();
std::set<unsigned> element_indices = iter->rGetContainingElementIndices();
// If there are no elements associated with this node, then we set the value of any VTK cell data to be -1 at this node...
if (element_indices.empty())
{
for (unsigned var=0; var<num_cell_data_items; var++)
{
cell_data[var][node_index] = -1.0;
}
}
else
{
// ... otherwise there should be exactly one element (i.e. cell) containing this node
assert(element_indices.size() == 1);
unsigned elem_index = *(element_indices.begin());
CellPtr p_cell = this->GetCellUsingLocationIndex(elem_index);
for (unsigned var=0; var<num_cell_data_items; var++)
{
cell_data[var][node_index] = p_cell->GetCellData()->GetItem(cell_data_names[var]);
}
}
}
for (unsigned var=0; var<cell_data.size(); var++)
{
mesh_writer.AddPointData(cell_data_names[var], cell_data[var]);
}
/*
* The current VTK writer can only write things which inherit from AbstractTetrahedralMeshWriter.
* For now, we do an explicit conversion to NodesOnlyMesh. This can be written to VTK then visualized as glyphs.
*/
NodesOnlyMesh<DIM> temp_mesh;
//.........这里部分代码省略.........