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.
Between Timo and me I think we've fixed the GridTools::find_active_cell_around_point function. Can you try again with the current development version and see if that also fixed the problem you had?
Best W. ------------------------------------------------------------------------ Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
