From: Anthony Geay <anthony.geay@edf.fr>
Date: Wed, 24 Aug 2022 13:37:09 +0000 (+0200)
Subject: Fix problem for not locatable points in the context of NodeFE spatial discr
X-Git-Tag: V9_10_0a1~8
X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=2be14e573fc988ee2e8c745fb24f09d5ee026507;p=tools%2Fmedcoupling.git

Fix problem for not locatable points in the context of NodeFE spatial discr
---

diff --git a/src/INTERP_KERNEL/InterpKernelRootsMultiDim.hxx b/src/INTERP_KERNEL/InterpKernelRootsMultiDim.hxx
index 68eb57acb..5d2faaba9 100644
--- a/src/INTERP_KERNEL/InterpKernelRootsMultiDim.hxx
+++ b/src/INTERP_KERNEL/InterpKernelRootsMultiDim.hxx
@@ -206,7 +206,10 @@ namespace INTERP_KERNEL
             if (temp > test) test=temp;
         }
         check=(test < TOLMIN);
-        return;
+        if( check )
+          return;
+        else
+          continue;
       }
       test=0.0;
       for (i=0;i<n;i++)
diff --git a/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py b/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py
index ef5df95cd..71cde7dde 100644
--- a/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py
+++ b/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py
@@ -1642,6 +1642,17 @@ class MEDCouplingBasicsTest(unittest.TestCase):
             ftest.setGaussLocalizationOnType(gt,sum([list(elt) for elt in MEDCouplingGaussLocalization.GetDefaultReferenceCoordinatesOf(gt).getValuesAsTuple()],[]),list(cooInRefElemFoundByNewton),[1])
             self.assertTrue ( float( (ftest.getLocalizationOfDiscr()-zePt).magnitude() ) < 1e-4 )
 
+        # testing with outside point
+        for zePt in outPts:
+            #safer than
+            #self.assertRaises(InterpKernelException,sd.getCooInRefElement,srcField.getMesh(),zePt.buildDADouble())
+            try:
+                sd.getCooInRefElement(srcField.getMesh(),zePt.buildDADouble())
+            except InterpKernelException as e:
+                self.assertTrue("fail to locate point" in e.what())
+            else:
+                self.assertTrue(false,"")    
+
     def checkMatrix(self,mat1,mat2,nbCols,eps):
         self.assertEqual(len(mat1),len(mat2))
         for i in range(len(mat1)):