[GL] Add eps to the bounding box search.

GeoLib/AABB.h

@@ -127,11 +127,14 @@ public:
      * check if point is in the axis aligned bounding box
      */
     template <typename T>
-    bool containsPoint(T const & pnt) const
+    bool containsPoint(T const & pnt, double eps) const
     {
-        if (pnt[0] < _min_pnt[0] || _max_pnt[0] <= pnt[0]) return false;
-        if (pnt[1] < _min_pnt[1] || _max_pnt[1] <= pnt[1]) return false;
-        if (pnt[2] < _min_pnt[2] || _max_pnt[2] <= pnt[2]) return false;
+        if (pnt[0] < _min_pnt[0] - eps || _max_pnt[0] + eps <= pnt[0])
+            return false;
+        if (pnt[1] < _min_pnt[1] - eps || _max_pnt[1] + eps <= pnt[1])
+            return false;
+        if (pnt[2] < _min_pnt[2] - eps || _max_pnt[2] + eps <= pnt[2])
+            return false;
         return true;
     }
 
@@ -166,7 +169,8 @@ public:
      */
     bool containsAABB(AABB const& other_aabb) const
     {
-        return containsPoint(other_aabb.getMinPoint()) && containsPoint(other_aabb.getMaxPoint());
+        return containsPoint(other_aabb.getMinPoint(), 0) &&
+               containsPoint(other_aabb.getMaxPoint(), 0);
     }
 
 protected:

GeoLib/Surface.cpp

@@ -94,9 +94,10 @@ const Triangle* Surface::operator[](std::size_t i) const
     return _sfc_triangles[i];
 }
 
-bool Surface::isPntInBoundingVolume(MathLib::Point3d const& pnt) const
+bool Surface::isPntInBoundingVolume(MathLib::Point3d const& pnt,
+                                    double eps) const
 {
-    return _bounding_volume->containsPoint(pnt);
+    return _bounding_volume->containsPoint(pnt, eps);
 }
 
 bool Surface::isPntInSfc(MathLib::Point3d const& pnt, double eps) const

GeoLib/Surface.h

@@ -62,7 +62,7 @@ public:
     /**
      * is the given point in the bounding volume of the surface
      */
-    bool isPntInBoundingVolume(MathLib::Point3d const& pnt) const;
+    bool isPntInBoundingVolume(MathLib::Point3d const& pnt, double eps) const;
 
     /**
      * is the given point pnt located in the surface

MeshGeoToolsLib/MeshNodesAlongSurface.cpp

@@ -35,7 +35,7 @@ MeshNodesAlongSurface::MeshNodesAlongSurface(MeshLib::Mesh const& mesh,
     // loop over all nodes
     for (std::size_t i = 0; i < n_nodes; i++) {
         auto* node = mesh_nodes[i];
-        if (!sfc.isPntInBoundingVolume(*node))
+        if (!sfc.isPntInBoundingVolume(*node, epsilon_radius))
             continue;
         if (sfc.isPntInSfc(*node, epsilon_radius)) {
             _msh_node_ids.push_back(node->getID());

MeshLib/MeshSearch/ElementSearch.cpp

@@ -59,7 +59,7 @@ std::size_t ElementSearch::searchByBoundingBox(
         [&aabb](MeshLib::Element* e) {
             std::size_t const nElemNodes (e->getNumberOfBaseNodes());
             for (std::size_t n=0; n < nElemNodes; ++n)
-                if (aabb.containsPoint(*e->getNode(n)))
+                if (aabb.containsPoint(*e->getNode(n), 0))
                     return true;    // any node of element is in aabb.
             return false;    // no nodes of element are in aabb.
         });

NumLib/Function/LinearInterpolationOnSurface.cpp

@@ -42,7 +42,7 @@ LinearInterpolationOnSurface::LinearInterpolationOnSurface(
 
 double LinearInterpolationOnSurface::operator()(const MathLib::Point3d& pnt) const
 {
-    if (!_sfc.isPntInBoundingVolume(pnt))
+    if (!_sfc.isPntInBoundingVolume(pnt, 0))
         return _default_value;
     auto* tri = _sfc.findTriangle(pnt);
     if (tri == nullptr)

Tests/GeoLib/TestAABB.cpp

@@ -284,9 +284,9 @@ TEST(GeoLib, AABBSinglePoint)
                 else if (k==0) p[2] = pnts.front()[2];
                 else p[2] = std::nextafter(pnts.front()[2], to_max);
                 if (i == 0 && j == 0 && k == 0)
-                    ASSERT_TRUE(aabb.containsPoint(p));
+                    ASSERT_TRUE(aabb.containsPoint(p, 0));
                 else
-                    ASSERT_FALSE(aabb.containsPoint(p));
+                    ASSERT_FALSE(aabb.containsPoint(p, 0));
             }
         }
     }