From: Michael Jermyn on
Is there an alternative to tsearchn that works for non-convex domains (tetrahedron mesh in 3d, triangle mesh in 2d)? I wrote my own, but it is very slow (it needs to be called a huge number of times in the package I work on). My version just finds the closest points and then searches all the nearby tetrahedrons/triangles.