From a812b4c783fac843e17068619a9016a03422f028 Mon Sep 17 00:00:00 2001
From: vbd <vbd>
Date: Wed, 7 May 2008 12:26:34 +0000
Subject: [PATCH] NPPAL 19079 : bug correction in PointLocator algos.

---
 src/INTERP_KERNEL/PointLocatorAlgos.txx | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

diff --git a/src/INTERP_KERNEL/PointLocatorAlgos.txx b/src/INTERP_KERNEL/PointLocatorAlgos.txx
index 412f004e9..7f18751c3 100644
--- a/src/INTERP_KERNEL/PointLocatorAlgos.txx
+++ b/src/INTERP_KERNEL/PointLocatorAlgos.txx
@@ -41,7 +41,10 @@ namespace INTERP_KERNEL{
 						}
 					for (int index= conn_index[i]; index < conn_index[i+1];index++)
 						{
-							const double* coordelem = coords+OTT<ConnType,numPol>::ind2C(conn[OTT<ConnType,numPol>::ind2C(index)]);
+							//coordelem points to the coordinates of the current node of the i-th element
+							const double* coordelem = coords+OTT<ConnType,numPol>::ind2C(conn[OTT<ConnType,numPol>::ind2C(index)])*SPACEDIM;
+
+							//the bounding box is updated by checking wheher the node is at the min/max in exach dimension
 							for (int idim=0; idim<SPACEDIM;idim++)
 								{
 									_bb[2*(i*SPACEDIM+idim)]=(coordelem[idim]<_bb[2*(i*SPACEDIM+idim)])?coordelem[idim]:_bb[2*(i*SPACEDIM+idim)];
@@ -58,6 +61,8 @@ namespace INTERP_KERNEL{
 			delete _tree;
 		}
 	
+		//returns the list of elements that contains 
+		//the point pointed to by x
 		std::list<int> locates(const double* x)
 		{
 			vector<int> candidates;
@@ -65,8 +70,9 @@ namespace INTERP_KERNEL{
 			list<int> retlist;
 			for (int i=0; i< candidates.size(); i++)
 				{
-					if (elementContainsPoint(i,x))
-						retlist.push_back(OTT<ConnType,numPol>::indFC(i));
+					int ielem=candidates[i];
+					if (elementContainsPoint(ielem,x))
+						retlist.push_back(OTT<ConnType,numPol>::indFC(ielem));
 				}
 			return retlist;
 		}
-- 
2.39.2