Added a function to find all mesh elements by coordinates.

Geo/GModel.cpp

@@ -697,6 +697,14 @@ MElement *GModel::getMeshElementByCoord(SPoint3 &p, int dim)
   }
   return _octree->find(p.x(), p.y(), p.z(), dim);
 }
+std::vector<MElement*> GModel::getMeshElementsByCoord(SPoint3 &p, int dim)
+{
+  if(!_octree){
+    Msg::Debug("Rebuilding mesh element octree");
+    _octree = new MElementOctree(this);
+  }
+  return _octree->findAll(p.x(), p.y(), p.z(), dim);
+}
 
 MVertex *GModel::getMeshVertexByTag(int n)
 {

Geo/GModel.h

@@ -297,6 +297,7 @@ class GModel
 
   // access a mesh element by coordinates (using an octree search)
   MElement *getMeshElementByCoord(SPoint3 &p, int dim = -1);
+  std::vector<MElement*> getMeshElementsByCoord(SPoint3 &p, int dim = -1);
 
   // access a mesh element by tag, using the element cache
   MElement *getMeshElementByTag(int n);