On 1/8/19 12:43 AM, 'Maxi Miller' via deal.II User Group wrote:
>
> I assume it tries to find the point, which is not available due to being
> within a cell located in another thread. I assume I have to modify the cell
> iterator in that case. Is that correct?
You mean owned by some other MPI process? I don't know whether you
interpretation is correct. The exception is found here:
// The first time around, we check for vertices in the hint_cell. If
// that does not work, we set the cell iterator to an invalid one, and
// look for a global vertex close to the point. If that does not work,
// we are in trouble, and just throw an exception.
//
// If we got here, then we did not find the point. If the
// current_cell.state() here is not IteratorState::valid, it means that
// the user did not provide a hint_cell, and at the beginning of the
// while loop we performed an actual global search on the mesh
// vertices. Not finding the point then means the point is outside the
// domain.
AssertThrow(current_cell.state() == IteratorState::valid,
ExcPointNotFound<spacedim>(p));
Do you think you can come up with a small testcase?
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
www: http://www.math.colostate.edu/~bangerth/
--
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
