Dear Steffen,
> The problem is GridTools::find_active_cell_around_point:
> 1. it find the closest vertex
> 2. and then loops through all cells to identify the correct cell that is
> connected to this vertex
> This method works fine for all uniformly refined grids.
>
> However, if you have a locally refined grid with hanging nodes, than the
> closest vertex might be a hanging node, and your point is in the
> "hanging-element". The loop through the cells will fail, because the point
> is in a cell, which has this hanging node not as vertex.
I understand your description, but I can not reproduce this here. If
you look into find_cells_adjacent_to_vertex(), there is a section
// Now we need to make
sure that the
// vertex is not a locally refined
// vertex not being part of the
// neighboring cells. So
we loop over
// all faces to which this vertex
// belongs, check the level of
// the neighbor, and if
it is coarser,
// then check whether the vertex is
// part of that neighbor or not.
...
I would really like to fix this bug, so it would be great if you can
help me reproduce it.
Best,
Timo
--
Timo Heister
http://www.math.tamu.edu/~heister/
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
