From: Roy S. <roy...@ic...> - 2007-07-05 13:34:04
|
On Thu, 5 Jul 2007, Tahar Amari wrote: > By the way is there somewhere in libmesh a way of determining to which > tetrahedra belongs a points with '(x,y,z) coordinates given a priori, or > do I have to make it myself ( I can figure out how to do it), but in case it > exists in libmesh I guess it will be optimized (there are several possible > algorithms but those are differently efficient !). The PointLocator class currently defaults to PointLocatorTree, which builds an octree structure for finding the element containing a point. This probably costs a bit of overhead when you first create the structure (especially if you do much adaptive mesh refinement, which invalidates old PointLocators), but it should be roughly O(log N) for each point you find with it. --- Roy |