[Intersect2D] Bug fix: residual cell construction

[tools/medcoupling.git] / src / INTERP_KERNEL / BBTreePts.txx

diff --git a/src/INTERP_KERNEL/BBTreePts.txx b/src/INTERP_KERNEL/BBTreePts.txx
index 7beb4b37bafc91e7f44b1a1fcb4075198f6ce54c..2961f63a8113509e07decf54b24f33f6eeb14639 100644 (file)
--- a/src/INTERP_KERNEL/BBTreePts.txx
+++ b/src/INTERP_KERNEL/BBTreePts.txx
@@ -1,4 +1,4 @@
-// Copyright (C) 2007-2016  CEA/DEN, EDF R&D
+// Copyright (C) 2007-2021  CEA/DEN, EDF R&D
 //
 // This library is free software; you can redistribute it and/or
 // modify it under the terms of the GNU Lesser General Public
@@ -90,9 +90,9 @@ public:
     new_elems_right.reserve(nbelems/2+1);
     double max_left = -std::numeric_limits<double>::max();
     double min_right=  std::numeric_limits<double>::max();
-    for(int i=0;i<nbelems;i++)
+    for(ConnType i=0;i<nbelems;i++)
       {
-        int elem;
+        ConnType elem;
         if(elems!=0)
           elem= elems[i];
         else
@@ -115,11 +115,11 @@ public:
     tmp=0;
     if(!new_elems_left.empty())
       tmp=&(new_elems_left[0]);
-    _left=new BBTreePts(pts, tmp, level+1, (int)new_elems_left.size(),_epsilon);
+    _left=new BBTreePts(pts, tmp, level+1, (ConnType)new_elems_left.size(),_epsilon);
     tmp=0;
     if(!new_elems_right.empty())
       tmp=&(new_elems_right[0]);
-    _right=new BBTreePts(pts, tmp, level+1, (int)new_elems_right.size(),_epsilon);
+    _right=new BBTreePts(pts, tmp, level+1, (ConnType)new_elems_right.size(),_epsilon);
   }
 
 
@@ -130,12 +130,12 @@ public:
     delete _right;
   }
 
-  /*! returns in \a elems the list of elements potentially containing the point pointed to by \a xx
+  /*! returns in \a elem an element potentially containing the point pointed to by \a xx
       Contrary to BBTreePts::getElementsAroundPoint the norm 2 is used here.
 
-    \param [in] xx pointer to query point coords
-    \param [in] threshold
-    \param elems list of elements (given in 0-indexing) intersecting the bounding box
+    \param[in] xx pointer to query point coords
+    \param[in] threshold detection precision
+    \param[out] elem the element intersecting the bounding box
     \sa BBTreePts::getElementsAroundPoint
   */
   double getElementsAroundPoint2(const double *xx, double threshold, ConnType& elem) const
@@ -164,7 +164,7 @@ public:
       return _left->getElementsAroundPoint2(xx,threshold,elem);
     if(xx[_level%dim]-s>_max_left)
       return _right->getElementsAroundPoint2(xx,threshold,elem);
-    int eleml,elemr;
+    ConnType eleml,elemr;
     double retl=_left->getElementsAroundPoint2(xx,threshold,eleml);
     double retr=_right->getElementsAroundPoint2(xx,threshold,elemr);
     if(retl<retr)
@@ -211,7 +211,7 @@ public:
     _right->getElementsAroundPoint(xx,elems);
   }
   
-  int size() const
+  ConnType size() const
   {
     if(_terminal)
       return _nbelems;