Skip to content

Fix/manifold history points

Chris Cantwell requested to merge fix/manifold-history-points into master

This MR fixes (again) the history points for the manifold case. Since the points may not truly lie in the manifold, we have to be forgiving. The algorithm now proceeds as follows

  • Each process searches for the closest element in its partition which might contain the history point.
  • Distances are recorded and the minimum is found across all processes to decide which process should 'own' the history point.
  • In the case where multiple processes satisfy the above (e.g. point lies on partition boundary), the one of highest rank is chosen.

A warning is printed if the distance from the point to the element is greater than the geometric tolerance constant. In this case, the history point will write out data from the projection of the point on the element.

Merge request reports