Skip to content

Improve performance of interpolation

This MR enhances the performance of point interpolation using an rtree for element lookup.

The linear search through elements is replaced by an rtree, populated with the bounding boxes of the elements. For a given point, the rtree is queried which returns the subset of bounding boxes containing that point. Only these bounding boxes are then searched to determine correct element. Bounding boxes for curved elements are defined based on quadrature points with a 10% margin.

Fixes #137 (closed).

Edited by Chris Cantwell

Merge request reports