From 2be14e573fc988ee2e8c745fb24f09d5ee026507 Mon Sep 17 00:00:00 2001
From: Anthony Geay <anthony.geay@edf.fr>
Date: Wed, 24 Aug 2022 15:37:09 +0200
Subject: [PATCH] Fix problem for not locatable points in the context of NodeFE
 spatial discr

---
 src/INTERP_KERNEL/InterpKernelRootsMultiDim.hxx |  5 ++++-
 src/MEDCoupling_Swig/MEDCouplingRemapperTest.py | 11 +++++++++++
 2 files changed, 15 insertions(+), 1 deletion(-)

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)):
-- 
2.39.2