Bug correction on semi partial field read.

diff --git a/src/MEDLoader/MEDFileField.cxx b/src/MEDLoader/MEDFileField.cxx
index dbdff70f15cb6ecf2b1e00ac215b5649e850208a..5a959dff765f6db2e9fac4705ee58f3dcc29b081 100644 (file)
--- a/src/MEDLoader/MEDFileField.cxx
+++ b/src/MEDLoader/MEDFileField.cxx
@@ -1689,7 +1689,7 @@ MEDCouplingFieldDouble *MEDFileFieldPerMesh::finishField2(TypeOfField type, cons
   if(da->isIdentity())
     {
       int nbOfTuples=da->getNumberOfTuples();
-      if(nbOfTuples==ComputeNbOfElems(glob,type,geoTypes,dads,locs))
+      if(nbOfTuples==mesh->getNumberOfCells())
         return finishField(type,glob,dads,locs,mesh,isPfl);
     }
   MEDCouplingAutoRefCountObjectPtr<MEDCouplingFieldDouble> ret=finishField(type,glob,dads,locs,mesh,isPfl);

diff --git a/src/MEDLoader/Swig/MEDLoaderTest3.py b/src/MEDLoader/Swig/MEDLoaderTest3.py
index ebf4b695e32cbae48b75d4b59ec5967c44c3f4de..37813d61ee561de3f26b12cacc586ff1cfa43d09 100644 (file)
--- a/src/MEDLoader/Swig/MEDLoaderTest3.py
+++ b/src/MEDLoader/Swig/MEDLoaderTest3.py
@@ -1316,6 +1316,37 @@ class MEDLoaderTest(unittest.TestCase):
         m=MEDCouplingMultiFields([])
         m=MEDCouplingFieldOverTime([])
         pass
+
+    # This is a non regression test. When a field lies partially on a mesh but fully on one of its geometric type.
+    def testBugSemiPartialField(self):
+        fname="Pyfile46.med"
+        m=MEDLoaderDataForTest.build2DMesh_3()
+        m=m[:10] ; m.setName("mesh")
+        f=m.getMeasureField(ON_CELLS)
+        f=f.buildNewTimeReprFromThis(ONE_TIME,False)
+        f.setTime(5.5,3,4)
+        f.setName("SemiPartialField")
+        #
+        f1=f[:6] ; f1.getMesh().setName(m.getName())
+        f2=f[6:] ; f2.getMesh().setName(m.getName())
+        #
+        mm=MEDFileUMesh.New()
+        mm.setMeshAtLevel(0,m)
+        ff=MEDFileField1TS.New()
+        ff.setFieldProfile(f1,mm,0,DataArrayInt.Range(0,6,1)) # no name on profile -> normally it is an error but in this special case
+        mm.write(fname,2)
+        ff.write(fname,0)
+        #
+        ff2=MEDFileField1TS.New(fname,f.getName(),f.getTime()[1],f.getTime()[2])
+        fread=ff2.getFieldOnMeshAtLevel(ON_CELLS,0,mm)
+        fread2=ff2.getFieldAtLevel(ON_CELLS,0)
+        #
+        fread.checkCoherency()
+        fread2.checkCoherency()
+        self.assertTrue(fread.isEqual(f1,1e-12,1e-12))
+        self.assertTrue(fread2.isEqual(f1,1e-12,1e-12))
+        pass
+
     pass
 
 unittest.main()