Hi Rene,
I have a problem in finding the list neighbor cells. This is the code I
wrote:
*********************
int iter3=0;
for (Triangulation<3>::active_cell_iterator cell = tr.begin_active(); cell
!= tr.end(); ++cell)
{
for (unsigned int v = 0; v < GeometryInfo<3>::vertices_per_cell;
++v)
{
auto v_to_c = GridTools::vertex_to_cell_map(tr);
auto v_to_c_d =
GridTools::vertex_to_cell_centers_directions(tr, v_to_c);
auto c_and_p = GridTools::find_active_cell_around_point(
mappinggg, tr, cell->vertex(v), v_to_c, v_to_c_d);
auto pcell = c_and_p.first;
std::cout<<"a neighbor of " << cell->id() <<"is : " <<
pcell->id()<<std::endl;
}
iter3++;
}
**********************************
The problem is when I print the neighbors of each cell in the terminal, in
the neighbors list for each cell only the cells with smaller ids comparing
to the main cell are illustrated.
For instance: for the neighbors of cell 0_3:765 only the cells: 0_3:742,
0_3:743, 0_3:760, 0_3:761, 0_3:746, 0_3:747 and 0_3:744 are shown.
Do you have any idea about the reason?
Best
Shahab
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/e0ddc6e9-f2ec-495b-9787-8887922d66c2%40googlegroups.com.
