X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FMEDLoader%2FSwig%2FMEDLoaderTest4.py;h=e5b5be72294fd92572b68d8c6e748b3a4942ecf7;hb=06b22e5321361b9c5fdfdbf4fd81e12d37ead050;hp=58932a087aa2ed67d47f272ec3428bd52719f945;hpb=f2494a83efd03107a558640c7ece7773cb06dfea;p=tools%2Fmedcoupling.git diff --git a/src/MEDLoader/Swig/MEDLoaderTest4.py b/src/MEDLoader/Swig/MEDLoaderTest4.py index 58932a087..e5b5be722 100644 --- a/src/MEDLoader/Swig/MEDLoaderTest4.py +++ b/src/MEDLoader/Swig/MEDLoaderTest4.py @@ -1,5 +1,5 @@ # -*- coding: iso-8859-1 -*- -# Copyright (C) 2007-2015 CEA/DEN, EDF R&D +# Copyright (C) 2007-2016 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 @@ -17,7 +17,7 @@ # # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com # -# Author : Anthony Geay (CEA/DEN) +# Author : Anthony Geay (EDF R&D) from MEDLoader import * import unittest @@ -37,13 +37,13 @@ class MEDLoaderTest4(unittest.TestCase): tri=MEDCouplingUMesh("tri",2) tri.allocateCells() ; tri.insertNextCell(NORM_TRI3,[0,1,2]) tri.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,0.)])) - tris=[tri.deepCpy() for i in xrange(4)] + tris = [tri.deepCopy() for i in range(4)] for i,elt in enumerate(tris): elt.translate([i,0]) tris=MEDCouplingUMesh.MergeUMeshes(tris) quad=MEDCouplingUMesh("quad",2) quad.allocateCells() ; quad.insertNextCell(NORM_QUAD4,[0,1,2,3]) quad.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,1.),(1.,0.)])) - quads=[quad.deepCpy() for i in xrange(5)] + quads = [quad.deepCopy() for i in range(5)] for i,elt in enumerate(quads): elt.translate([5+i,0]) quads=MEDCouplingUMesh.MergeUMeshes(quads) m=MEDCouplingUMesh.MergeUMeshes(tris,quads) @@ -64,13 +64,13 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) fs.pushBackTimeStep(f) ##### Time step 1 @@ -80,14 +80,14 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(1100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(1200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) fs.pushBackTimeStep(f) ##### Time step 2 @@ -97,14 +97,14 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(2100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(2200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) fs.pushBackTimeStep(f) ##### Time step 3 @@ -115,14 +115,14 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(3100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(3200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) # fs.pushBackTimeStep(f) @@ -134,19 +134,19 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(4100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(4200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) fs.pushBackTimeStep(f) mm.write(fname,2) fs.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,m1,mm,fs,f,fCell0,fCell1 ########## GO for reading in MEDReader, by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -181,7 +181,7 @@ class MEDLoaderTest4(unittest.TestCase): mml=fcscp.buildFromScratchDataSetSupport(0,fields) mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) - for i in xrange(1,5): + for i in range(1, 5): self.assertTrue(fcscp.isDataSetSupportEqualToThePreviousOne(i,fields)) pass ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() @@ -198,7 +198,8 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,100,101,102,103,104,105,106,107,108]))) self.assertTrue(not a9) - for i in xrange(5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(5): fsst=MEDFileField1TSStructItem.BuildItemFrom(fields[0][i],mst) fields[0][i].loadArraysIfNecessary() tmpMem=fields.getHeapMemorySize() @@ -221,13 +222,13 @@ class MEDLoaderTest4(unittest.TestCase): tri=MEDCouplingUMesh("tri",2) tri.allocateCells() ; tri.insertNextCell(NORM_TRI3,[0,1,2]) tri.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,0.)])) - tris=[tri.deepCpy() for i in xrange(4)] + tris = [tri.deepCopy() for i in range(4)] for i,elt in enumerate(tris): elt.translate([i,0]) tris=MEDCouplingUMesh.MergeUMeshes(tris) quad=MEDCouplingUMesh("quad",2) quad.allocateCells() ; quad.insertNextCell(NORM_QUAD4,[0,1,2,3]) quad.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,1.),(1.,0.)])) - quads=[quad.deepCpy() for i in xrange(5)] + quads = [quad.deepCopy() for i in range(5)] for i,elt in enumerate(quads): elt.translate([5+i,0]) quads=MEDCouplingUMesh.MergeUMeshes(quads) m=MEDCouplingUMesh.MergeUMeshes(tris,quads) @@ -244,20 +245,20 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) # fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName) ; fNode.setMesh(m1) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(300) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs.pushBackTimeStep(f) ##### Time step 1 @@ -267,21 +268,21 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(1100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(1200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) # fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName) ; fNode.setMesh(m1) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(1300) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs.pushBackTimeStep(f) ##### Time step 2 @@ -291,21 +292,21 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName) ; fNode.setMesh(m1) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(2300) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) # fCell0=MEDCouplingFieldDouble(ON_CELLS) ; fCell0.setTime(float(i),i,0) fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(2100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(2200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) fs.pushBackTimeStep(f) ##### Time step 3 @@ -316,21 +317,21 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(3100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(3200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) # fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName) ; fNode.setMesh(m1) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(3300) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) # fs.pushBackTimeStep(f) @@ -342,27 +343,27 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(4100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(4200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) # fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName) ; fNode.setMesh(m1) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(4300) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) # fs.pushBackTimeStep(f) mm.write(fname,2) fs.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,m1,mm,fs,f,fCell0,fCell1 ########## GO for reading in MEDReader, by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -398,7 +399,7 @@ class MEDLoaderTest4(unittest.TestCase): mml=fcscp.buildFromScratchDataSetSupport(0,fields) mml2=mml.prepare() assert isinstance(mml2,MEDUMeshMultiLev) - for i in xrange(1,5): + for i in range(1, 5): self.assertTrue(fcscp.isDataSetSupportEqualToThePreviousOne(i,fields)) pass ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() @@ -409,8 +410,9 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([2,0,1,2,1,2,2,2,0,2,3,4,2,4,5,2,5,3,2,6,7,2,7,8,2,8,6,2,9,10,2,10,11,2,11,9,2,12,13,2,13,14,2,14,15,2,15,12,2,16,17,2,17,18,2,18,19,2,19,16,2,20,21,2,21,22,2,22,23,2,23,20,2,24,25,2,25,26,2,26,27,2,27,24,2,28,29,2,29,30,2,30,31,2,31,28,3,0,1,2,3,3,4,5,3,6,7,8,3,9,10,11,4,12,13,14,15,4,16,17,18,19,4,20,21,22,23,4,24,25,26,27,4,28,29,30,31]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) # for cells - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst)# Second 0 is for cells f.loadArraysIfNecessary() @@ -419,7 +421,7 @@ class MEDLoaderTest4(unittest.TestCase): vExp=DataArrayDouble([200.,201.,202.,203.,204.,205.,206.,207.,208.,209.,210.,211.,212.,213.,214.,215.,216.,217.,218.,219.,220.,221.,222.,223.,224.,225.,226.,227.,228.,229.,230.,231.,232.,233.,234.,235.,236.,237.,238.,239.,240.,241.,242.,243.,244.,245.,246.,247.,248.,249.,250.,251.,252.,253.,254.,255.,256.,257.,258.,259.,260.,261.,262.,263.,100.,101.,102.,103.,104.,105.,106.,107.,108.,109.,110.,111.,112.,113.,114.,115.,116.,117.],41,2) ; vExp.setInfoOnComponents(['Comp1 [m]','Com2 [s^2]']) ; vExp+=i*1000 self.assertTrue(v.isEqual(vExp,1e-12)) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][1][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst)# Second 0 is for cells f.loadArraysIfNecessary() @@ -439,13 +441,13 @@ class MEDLoaderTest4(unittest.TestCase): tri=MEDCouplingUMesh("tri",2) tri.allocateCells() ; tri.insertNextCell(NORM_TRI3,[0,1,2]) tri.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,0.)])) - tris=[tri.deepCpy() for i in xrange(4)] + tris = [tri.deepCopy() for i in range(4)] for i,elt in enumerate(tris): elt.translate([i,0]) tris=MEDCouplingUMesh.MergeUMeshes(tris) quad=MEDCouplingUMesh("quad",2) quad.allocateCells() ; quad.insertNextCell(NORM_QUAD4,[0,1,2,3]) quad.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,1.),(1.,0.)])) - quads=[quad.deepCpy() for i in xrange(5)] + quads = [quad.deepCopy() for i in range(5)] for i,elt in enumerate(quads): elt.translate([5+i,0]) quads=MEDCouplingUMesh.MergeUMeshes(quads) m=MEDCouplingUMesh.MergeUMeshes(tris,quads) @@ -462,13 +464,13 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) fs.pushBackTimeStep(f) ##### Time step 1 on nodes @@ -478,7 +480,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName) ; fNode.setMesh(m1) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(1300) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs.pushBackTimeStep(f) ##### Time step 2 on cells @@ -488,14 +490,14 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(2100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(2200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) fs.pushBackTimeStep(f) ##### Time step 3 on nodes @@ -505,7 +507,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName) ; fNode.setMesh(m1) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(3300) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs.pushBackTimeStep(f) ##### Time step 4 @@ -515,20 +517,20 @@ class MEDLoaderTest4(unittest.TestCase): fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfCells()) ; arr.iota(4100) ; arr.rearrange(2) fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell0.checkCoherency() + fCell0.checkConsistencyLight() f.setFieldNoProfileSBT(fCell0) # fCell1=MEDCouplingFieldDouble(ON_CELLS) ; fCell1.setTime(float(i),i,0) fCell1.setName(fieldName) ; fCell1.setMesh(m1) arr=DataArrayDouble(2*m1.getNumberOfCells()) ; arr.iota(4200) ; arr.rearrange(2) fCell1.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]","Com2 [s^2]"]) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() f.setFieldNoProfileSBT(fCell1) # fs.pushBackTimeStep(f) mm.write(fname,2) fs.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,m1,mm,fs,f,fCell0,fCell1 ########## GO for reading in MEDReader, by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -562,7 +564,7 @@ class MEDLoaderTest4(unittest.TestCase): mml=fcscp.buildFromScratchDataSetSupport(0,fields) mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) - for i in xrange(1,3): + for i in range(1, 3): self.assertTrue(fcscp.isDataSetSupportEqualToThePreviousOne(i,fields)) pass ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() @@ -573,8 +575,9 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([2,0,1,2,1,2,2,2,0,2,3,4,2,4,5,2,5,3,2,6,7,2,7,8,2,8,6,2,9,10,2,10,11,2,11,9,2,12,13,2,13,14,2,14,15,2,15,12,2,16,17,2,17,18,2,18,19,2,19,16,2,20,21,2,21,22,2,22,23,2,23,20,2,24,25,2,25,26,2,26,27,2,27,24,2,28,29,2,29,30,2,30,31,2,31,28,3,0,1,2,3,3,4,5,3,6,7,8,3,9,10,11,4,12,13,14,15,4,16,17,18,19,4,20,21,22,23,4,24,25,26,27,4,28,29,30,31]))) assert a4 is None assert a5 is None + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) # for cells - for i in xrange(3): + for i in range(3): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst)# Second 0 is for cells f.loadArraysIfNecessary() @@ -592,7 +595,7 @@ class MEDLoaderTest4(unittest.TestCase): mml=fcscp.buildFromScratchDataSetSupport(0,fields) mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) - for i in xrange(1,2): + for i in range(1, 2): self.assertTrue(fcscp.isDataSetSupportEqualToThePreviousOne(i,fields)) pass ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() @@ -603,7 +606,8 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,0,1,2,3,3,4,5,3,6,7,8,3,9,10,11,4,12,13,14,15,4,16,17,18,19,4,20,21,22,23,4,24,25,26,27,4,28,29,30,31,2,0,1,2,1,2,2,2,0,2,3,4,2,4,5,2,5,3,2,6,7,2,7,8,2,8,6,2,9,10,2,10,11,2,11,9,2,12,13,2,13,14,2,14,15,2,15,12,2,16,17,2,17,18,2,18,19,2,19,16,2,20,21,2,21,22,2,22,23,2,23,20,2,24,25,2,25,26,2,26,27,2,27,24,2,28,29,2,29,30,2,30,31,2,31,28]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(2): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst)# Second 0 is for cells f.loadArraysIfNecessary() @@ -623,13 +627,13 @@ class MEDLoaderTest4(unittest.TestCase): tri=MEDCouplingUMesh("tri",2) tri.allocateCells() ; tri.insertNextCell(NORM_TRI3,[0,1,2]) tri.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,0.)])) - tris=[tri.deepCpy() for i in xrange(4)] + tris = [tri.deepCopy() for i in range(4)] for i,elt in enumerate(tris): elt.translate([i,0]) tris=MEDCouplingUMesh.MergeUMeshes(tris) quad=MEDCouplingUMesh("quad",2) quad.allocateCells() ; quad.insertNextCell(NORM_QUAD4,[0,1,2,3]) quad.setCoords(DataArrayDouble([(0.,0.),(0.,1.),(1.,1.),(1.,0.)])) - quads=[quad.deepCpy() for i in xrange(5)] + quads = [quad.deepCopy() for i in range(5)] for i,elt in enumerate(quads): elt.translate([5+i,0]) quads=MEDCouplingUMesh.MergeUMeshes(quads) m=MEDCouplingUMesh.MergeUMeshes(tris,quads) @@ -647,7 +651,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(0+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) # @@ -656,7 +660,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName2) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(100+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs2.pushBackTimeStep(f) # @@ -665,7 +669,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName3) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(200+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs3.pushBackTimeStep(f) ##### Time step 1 @@ -675,7 +679,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(0+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) # @@ -684,7 +688,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName2) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(100+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs2.pushBackTimeStep(f) # @@ -693,7 +697,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName3) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(200+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs3.pushBackTimeStep(f) ##### Time step 2 @@ -703,7 +707,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(0+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) # @@ -712,7 +716,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName2) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(100+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs2.pushBackTimeStep(f) # @@ -721,7 +725,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName3) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(200+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs3.pushBackTimeStep(f) ##### Time step 3 @@ -731,7 +735,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(0+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) # @@ -740,7 +744,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName2) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(100+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs2.pushBackTimeStep(f) # @@ -749,7 +753,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName3) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(200+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs3.pushBackTimeStep(f) ##### Time step 4 @@ -759,7 +763,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(0+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) # @@ -768,7 +772,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName2) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(100+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs2.pushBackTimeStep(f) # @@ -777,13 +781,13 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setName(fieldName3) ; fNode.setMesh(m) arr=DataArrayDouble(2*m.getNumberOfNodes()) ; arr.iota(200+1000*i) ; arr.rearrange(2) fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) - fNode.checkCoherency() + fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs3.pushBackTimeStep(f) # mm.write(fname,2) fs1.write(fname,0) ; fs2.write(fname,0) ; fs3.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs1,fs2,fs3,f,fNode ########## GO for reading in MEDReader, by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -819,7 +823,7 @@ class MEDLoaderTest4(unittest.TestCase): mml=fcscp.buildFromScratchDataSetSupport(0,fields) mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) - for i in xrange(1,5): + for i in range(1, 5): self.assertTrue(fcscp.isDataSetSupportEqualToThePreviousOne(i,fields)) pass ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() @@ -830,8 +834,9 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,0,1,2,3,3,4,5,3,6,7,8,3,9,10,11,4,12,13,14,15,4,16,17,18,19,4,20,21,22,23,4,24,25,26,27,4,28,29,30,31]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) # test all the time steps of the 1/1 time step serie, on field 1 - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -842,7 +847,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(v.isEqual(vExp,1e-12)) pass # test all the time steps of the 1/1 time step serie, on field 2 - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][1][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -853,7 +858,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(v.isEqual(vExp,1e-12)) pass # test all the time steps of the 1/1 time step serie, on field 3 - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][2][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -888,7 +893,7 @@ class MEDLoaderTest4(unittest.TestCase): fieldName3="zeField3" ; pfl3=DataArrayInt([0,1,2,3,4,5,9,10]) ; pfl3.setName("pfl3") # on cells but different support fs1=MEDFileFieldMultiTS() ; fs2=MEDFileFieldMultiTS() ; fs3=MEDFileFieldMultiTS() # - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName1) @@ -915,7 +920,7 @@ class MEDLoaderTest4(unittest.TestCase): pass mm.write(fname,2) fs1.write(fname,0) ; fs2.write(fname,0) ; fs3.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs1,fs2,fs3,f,fNode ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -950,7 +955,7 @@ class MEDLoaderTest4(unittest.TestCase): mml=fcscp.buildFromScratchDataSetSupport(0,fields) mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) - for i in xrange(1,5): + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() @@ -967,10 +972,11 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([100,101,102,103,104,105]))) self.assertTrue(not a9) - for i in xrange(5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(5): nbOfT=[6,8] fieldNames=[fieldName1,fieldName2] - for j in xrange(2): + for j in range(2): m={"i":j} f=allFMTSLeavesPerCommonSupport[0][0][j][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) @@ -987,7 +993,7 @@ class MEDLoaderTest4(unittest.TestCase): mml=fcscp.buildFromScratchDataSetSupport(0,fields) mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) - for i in xrange(1,5): + for i in range(1, 5): self.assertTrue(fcscp.isDataSetSupportEqualToThePreviousOne(i,fields)) pass ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() @@ -1004,7 +1010,8 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([100,101,102,103,104,105,109,110]))) self.assertTrue(not a9) - for i in xrange(5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(5): f=allFMTSLeavesPerCommonSupport[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1036,12 +1043,12 @@ class MEDLoaderTest4(unittest.TestCase): fieldName4="zeField4" ;# on nodes fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() ; fs2=MEDFileFieldMultiTS() ; fs3=MEDFileFieldMultiTS() ; fs4=MEDFileFieldMultiTS() # - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) arr=DataArrayDouble(2*8) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) # @@ -1076,7 +1083,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName4) ; fNode.setMesh(m) arr=DataArrayDouble(2*15) ; arr.iota(400+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_4 [m]","Com2_4 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_4 [m]","Com2_4 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs4.pushBackTimeStep(f) pass @@ -1121,7 +1128,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a.isEqual(coordsX,1e-12)) self.assertTrue(b.isEqual(coordsY,1e-12)) self.assertTrue(isinstance(mml2,MEDCMeshMultiLev)) - for i in xrange(1,5): + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass a6,a7=mml2.retrieveFamilyIdsOnCells() @@ -1133,7 +1140,8 @@ class MEDLoaderTest4(unittest.TestCase): a10,a11=mml2.retrieveNumberIdsOnNodes() self.assertTrue(a10.isEqual(DataArrayInt([200,201,202,203,204,205,206,207,208,209,210,211,212,213,214]))) self.assertTrue(a11) # True because no copy - for i in xrange(5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1170,10 +1178,11 @@ class MEDLoaderTest4(unittest.TestCase): a10,a11=mml2.retrieveNumberIdsOnNodes() self.assertTrue(a10.isEqual(DataArrayInt([202,203,204,207,208,209,212,213,214]))) self.assertTrue(not a11) # False because copy - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1213,7 +1222,8 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([102,103,105,107]))) self.assertTrue(not a9) # False because copy - for i in xrange(5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(5): f=allFMTSLeavesPerCommonSupport[2][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1245,12 +1255,12 @@ class MEDLoaderTest4(unittest.TestCase): fieldName4="zeField4" ;# on nodes fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() ; fs2=MEDFileFieldMultiTS() ; fs3=MEDFileFieldMultiTS() ; fs4=MEDFileFieldMultiTS() # - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) arr=DataArrayDouble(2*8) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) # @@ -1285,7 +1295,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName4) ; fNode.setMesh(m) arr=DataArrayDouble(2*15) ; arr.iota(400+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_4 [m]","Com2_4 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_4 [m]","Com2_4 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs4.pushBackTimeStep(f) pass @@ -1336,10 +1346,11 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([100,101,102,103,104,105,106,107]))) self.assertTrue(a9) # True because no copy - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1364,7 +1375,7 @@ class MEDLoaderTest4(unittest.TestCase): mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDCurveLinearMeshMultiLev)) # here the 2nd support is a part of CMesh that is also a CMesh -> CMesh not a UMesh a,b,c=mml2.buildVTUArrays() - self.assertTrue(not c)#False here because a is the result of a computation not the internal strucutre + self.assertTrue(not c)#False here because a is the result of a computation not the internal structure self.assertTrue(a.isEqual(a0Exp[pfl2],1e-12)) self.assertEqual(b,[3,3]) a6,a7=mml2.retrieveFamilyIdsOnCells() @@ -1373,10 +1384,11 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([102,103,106,107]))) self.assertTrue(not a9) # False because copy - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1416,7 +1428,8 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([102,103,105,107]))) self.assertTrue(not a9) # False because copy - for i in xrange(5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(5): f=allFMTSLeavesPerCommonSupport[2][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1448,12 +1461,12 @@ class MEDLoaderTest4(unittest.TestCase): fieldName2="zeField2" fieldName3="zeField3" fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() ; fs2=MEDFileFieldMultiTS() ; fs3=MEDFileFieldMultiTS() - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_GAUSS_NE) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) arr=DataArrayDouble(2*38) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) # @@ -1461,7 +1474,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*11) ; arr.iota(100+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) # @@ -1473,7 +1486,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setGaussLocalizationOnCells([6,7,8],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2],[0.8,0.05,0.1,0.04]) fNode.setGaussLocalizationOnCells([9,10],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2,0.3,0.3,0.4,0.4,0.8,0.8],[0.8,0.05,0.1,0.01,0.02,0.005,0.005]) arr=DataArrayDouble(2*(4*2+2*5+3*4+2*7)) ; arr.iota(300+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs2.pushBackTimeStep(f) # @@ -1481,7 +1494,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName3) ; fNode.setMesh(m) arr=DataArrayDouble(2*15) ; arr.iota(400+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_3 [m]","Com2_3 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_3 [m]","Com2_3 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs3.pushBackTimeStep(f) # @@ -1489,7 +1502,7 @@ class MEDLoaderTest4(unittest.TestCase): # mm.write(fname,2) fs0.write(fname,0) ; fs1.write(fname,0) ; fs2.write(fname,0) ; fs3.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs1,fs2,fs3,f,fNode ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -1531,10 +1544,11 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,2,7,3,3,7,8,3,3,3,8,4,3,8,9,4,3,13,9,8,3,13,14,9,4,0,5,6,1,4,1,6,7,2,4,5,10,11,6,4,6,11,12,7,4,12,13,8,7]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1596,7 +1610,7 @@ class MEDLoaderTest4(unittest.TestCase): pfl1=DataArrayInt([0,1,7,9,10]) ; pfl1.setName("pfl1") # on cells pfl2=DataArrayInt([1,2,3,6,7,8,11,12,13]) ; pfl2.setName("pfl2") # on nodes fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() ; fs2=MEDFileFieldMultiTS() ; fs3=MEDFileFieldMultiTS() - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_GAUSS_NE) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) @@ -1621,7 +1635,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setGaussLocalizationOnCells([2,3],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2],[0.8,0.05,0.1,0.04]) fNode.setGaussLocalizationOnCells([4],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2,0.3,0.3,0.4,0.4,0.8,0.8],[0.8,0.05,0.1,0.01,0.02,0.005,0.005]) arr=DataArrayDouble(2*(2*1+5*1+4*2+7*1)) ; arr.iota(300+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldProfile(fNode,mm,0,pfl1) fs2.pushBackTimeStep(f) # @@ -1637,7 +1651,7 @@ class MEDLoaderTest4(unittest.TestCase): # mm.write(fname,2) fs0.write(fname,0) ; fs1.write(fname,0) ; fs2.write(fname,0) ; fs3.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs1,fs2,fs3,f,fNode ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -1679,10 +1693,11 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,1,4,2,3,4,5,2,4,0,3,4,1,4,3,6,7,4,4,7,8,5,4]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1741,7 +1756,7 @@ class MEDLoaderTest4(unittest.TestCase): fieldName2="zeField2" pfl1=DataArrayInt([1,2,3,6,7,8,11,12,13]) ; pfl1.setName("pfl1") # on nodes fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() ; fs2=MEDFileFieldMultiTS() - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_NODES) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) @@ -1770,7 +1785,7 @@ class MEDLoaderTest4(unittest.TestCase): # mm.write(fname,2) fs0.write(fname,0) ; fs1.write(fname,0) ; fs2.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs1,fs2,f,fNode ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -1812,10 +1827,11 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,1,4,2,3,4,5,2,4,0,3,4,1,4,3,6,7,4,4,7,8,5,4]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1858,7 +1874,7 @@ class MEDLoaderTest4(unittest.TestCase): # fieldName0="zeField0" fs0=MEDFileFieldMultiTS() - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_GAUSS_PT) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) @@ -1867,13 +1883,13 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setGaussLocalizationOnCells([6,10,13],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2],[0.8,0.05,0.1,0.04]) fNode.setGaussLocalizationOnCells([11,12,14],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2,0.3,0.3,0.4,0.4,0.8,0.8],[0.8,0.05,0.1,0.01,0.02,0.005,0.005]) arr=DataArrayDouble(2*(2*6+5*4+4*3+7*3)) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) pass mm.write(fname,2) fs0.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs0,f,fNode ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -1915,10 +1931,11 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([4,1,0,5,6,4,3,2,7,8,4,4,3,8,9,4,6,5,10,11,4,9,8,13,14,4,19,18,23,24,4,2,1,6,7,4,7,6,11,12,4,11,10,15,16,4,12,11,16,17,4,8,7,12,13,4,13,12,17,18,4,17,16,21,22,4,14,13,18,19,4,16,15,20,21,4,18,17,22,23]))) # <- here the mesh is renumbered : the mesh is equal to m[[0,2,3,4,7,15, 1,5,8,9, 6,10,13, 11,12,14]] self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -1945,7 +1962,7 @@ class MEDLoaderTest4(unittest.TestCase): fieldName0="zeField0" fieldName1="zeField1" fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_GAUSS_PT) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) @@ -1954,7 +1971,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setGaussLocalizationOnCells([6,10,13],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2],[0.8,0.05,0.1,0.04]) fNode.setGaussLocalizationOnCells([11,12,14],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2,0.3,0.3,0.4,0.4,0.8,0.8],[0.8,0.05,0.1,0.01,0.02,0.005,0.005]) arr=DataArrayDouble(2*(2*6+5*4+4*3+7*3)) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) # @@ -1962,13 +1979,13 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*16) ; arr.iota(300+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) pass mm.write(fname,2) fs0.write(fname,0) ; fs1.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs0,fs1,f,fNode ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -2010,10 +2027,11 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([4,1,0,5,6,4,2,1,6,7,4,3,2,7,8,4,4,3,8,9,4,6,5,10,11,4,7,6,11,12,4,8,7,12,13,4,9,8,13,14,4,11,10,15,16,4,12,11,16,17,4,13,12,17,18,4,14,13,18,19,4,16,15,20,21,4,17,16,21,22,4,18,17,22,23,4,19,18,23,24]))) # <- here the mesh is NOT renumbered : the mesh is equal to m self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -2050,12 +2068,12 @@ class MEDLoaderTest4(unittest.TestCase): fieldName1="zeField1" fieldName2="zeField2" ; pfl1=DataArrayInt([2,3]) ; pfl1.setName("pfl1") fieldName3="zefield3" ; pfl2=DataArrayInt([2,3,4]) ; pfl2.setName("pfl2") - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) arr=DataArrayDouble(2*5) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) # @@ -2063,7 +2081,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*5) ; arr.iota(100+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) # @@ -2071,7 +2089,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName2) ; fNode.setMesh(m[pfl1]) arr=DataArrayDouble(2*2) ; arr.iota(200+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_2 [m]","Com2_2 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldProfile(fNode,mm,0,pfl1) fs2.pushBackTimeStep(f) # @@ -2079,13 +2097,13 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName3) ; fNode.setMesh(m[pfl2]) arr=DataArrayDouble(2*3) ; arr.iota(300+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_3 [m]","Com2_3 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_3 [m]","Com2_3 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldProfile(fNode,mm,0,pfl2) fs3.pushBackTimeStep(f) pass mm.write(fname,2) fs0.write(fname,0) ; fs1.write(fname,0) ; fs2.write(fname,0) ; fs3.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs0 ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -2129,11 +2147,12 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([8,1,0,6,7,13,12,18,19,8,2,1,7,8,14,13,19,20,8,2,3,8,9,14,15,20,21,8,3,4,9,10,15,16,21,22,8,4,5,10,11,16,17,22,23]))) self.assertTrue(a4.isEqual(DataArrayInt([-1,-1,0,31,62]))) self.assertTrue(a5.isEqual(DataArrayInt([6,4,3,2,8,9,4,15,21,20,14,4,3,15,14,2,4,2,14,20,8,4,8,20,21,9,4,9,21,15,3,6,4,4,3,9,10,4,16,22,21,15,4,4,16,15,3,4,3,15,21,9,4,9,21,22,10,4,10,22,16,4,6,4,5,4,10,11,4,17,23,22,16,4,5,17,16,4,4,4,16,22,10,4,10,22,23,11,4,11,23,17,5]))) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -2165,11 +2184,12 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([8,2,3,8,9,14,15,20,21,8,3,4,9,10,15,16,21,22]))) self.assertTrue(a4.isEqual(DataArrayInt([0,31]))) self.assertTrue(a5.isEqual(DataArrayInt([6,4,3,2,8,9,4,15,21,20,14,4,3,15,14,2,4,2,14,20,8,4,8,20,21,9,4,9,21,15,3,6,4,4,3,9,10,4,16,22,21,15,4,4,16,15,3,4,3,15,21,9,4,9,21,22,10,4,10,22,16,4]))) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -2192,11 +2212,11 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([8,2,3,8,9,14,15,20,21,8,3,4,9,10,15,16,21,22,8,4,5,10,11,16,17,22,23]))) self.assertTrue(a4.isEqual(DataArrayInt([0,31,62]))) self.assertTrue(a5.isEqual(DataArrayInt([6,4,3,2,8,9,4,15,21,20,14,4,3,15,14,2,4,2,14,20,8,4,8,20,21,9,4,9,21,15,3,6,4,4,3,9,10,4,16,22,21,15,4,4,16,15,3,4,3,15,21,9,4,9,21,22,10,4,10,22,16,4,6,4,5,4,10,11,4,17,23,22,16,4,5,17,16,4,4,4,16,22,10,4,10,22,23,11,4,11,23,17,5]))) - for i in xrange(1,5): + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[2][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -2222,12 +2242,12 @@ class MEDLoaderTest4(unittest.TestCase): fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() fieldName0="zeField0" fieldName1="zeField1" - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) arr=DataArrayDouble(2*3) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) # @@ -2235,13 +2255,13 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*3) ; arr.iota(100+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) pass mm.write(fname,2) fs0.write(fname,0) ; fs1.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs0 ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -2283,7 +2303,8 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([8,2,3,8,9,14,15,20,21,8,3,4,9,10,15,16,21,22,8,4,5,10,11,16,17,22,23]))) self.assertTrue(a4.isEqual(DataArrayInt([0,31,62]))) self.assertTrue(a5.isEqual(DataArrayInt([6,4,3,2,8,9,4,15,21,20,14,4,3,15,14,2,4,2,14,20,8,4,8,20,21,9,4,9,21,15,3,6,4,4,3,9,10,4,16,22,21,15,4,4,16,15,3,4,3,15,21,9,4,9,21,22,10,4,10,22,16,4,6,4,5,4,10,11,4,17,23,22,16,4,5,17,16,4,4,4,16,22,10,4,10,22,23,11,4,11,23,17,5]))) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass a6,a7=mml2.retrieveFamilyIdsOnCells() @@ -2292,7 +2313,7 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8 is None) self.assertTrue(a9) - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -2348,12 +2369,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m0.getNumberOfCells()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) fCell.setArray(a) fCell.setTime(*t) - fCell.checkCoherency() + fCell.checkConsistencyLight() a=DataArrayDouble(m0.getNumberOfNodes()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode.setArray(a) fNode.setTime(*t) - fNode.checkCoherency() + fNode.checkConsistencyLight() f1ts.setFieldNoProfileSBT(fCell) f1ts.setFieldNoProfileSBT(fNode) ffs.pushBackTimeStep(f1ts) @@ -2363,12 +2384,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m0.getNumberOfCells()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) fCell.setArray(a) fCell.setTime(*t) - fCell.checkCoherency() + fCell.checkConsistencyLight() a=DataArrayDouble(m0.getNumberOfNodes()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode.setArray(a) fNode.setTime(*t) - fNode.checkCoherency() + fNode.checkConsistencyLight() f1ts.setFieldNoProfileSBT(fCell) f1ts.setFieldNoProfileSBT(fNode) ffs.pushBackTimeStep(f1ts) @@ -2378,12 +2399,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m0.getNumberOfCells()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) fCell.setArray(a) fCell.setTime(*t) - fCell.checkCoherency() + fCell.checkConsistencyLight() a=DataArrayDouble(m0.getNumberOfNodes()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode.setArray(a) fNode.setTime(*t) - fNode.checkCoherency() + fNode.checkConsistencyLight() f1ts.setFieldNoProfileSBT(fCell) f1ts.setFieldNoProfileSBT(fNode) ffs.pushBackTimeStep(f1ts) @@ -2393,12 +2414,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m0.getNumberOfCells()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) fCell.setArray(a) fCell.setTime(*t) - fCell.checkCoherency() + fCell.checkConsistencyLight() a=DataArrayDouble(m0.getNumberOfNodes()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode.setArray(a) fNode.setTime(*t) - fNode.checkCoherency() + fNode.checkConsistencyLight() f1ts.setFieldNoProfileSBT(fCell) f1ts.setFieldNoProfileSBT(fNode) ffs.pushBackTimeStep(f1ts) @@ -2451,6 +2472,7 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([100,101,102,103,104,105,106,107]))) self.assertTrue(a9) # no copy here + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) pass def test16(self): @@ -2483,7 +2505,7 @@ class MEDLoaderTest4(unittest.TestCase): mm1.setFamiliesIdsOnGroup("Grp1_1",[0,1]) ; mm1.setFamiliesIdsOnGroup("Grp1_2",[2,3]) mms.pushMesh(mm1) ; del mm1 # - m1=m0.deepCpy() ; m1.translate([2.5,0.,0.]) ; m1.setName("mesh2") + m1=m0.deepCopy() ; m1.translate([2.5,0.,0.]) ; m1.setName("mesh2") # fCell2=MEDCouplingFieldDouble(ON_CELLS) fCell2.setName("zeField3_1") @@ -2515,12 +2537,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m0.getNumberOfCells()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) fCell1.setArray(a) fCell1.setTime(*t) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() a=DataArrayDouble(m0.getNumberOfNodes()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode1.setArray(a) fNode1.setTime(*t) - fNode1.checkCoherency() + fNode1.checkConsistencyLight() f1ts1.setFieldNoProfileSBT(fCell1) ; ffs1_1.pushBackTimeStep(f1ts1) f1ts2.setFieldNoProfileSBT(fNode1) ; ffs1_2.pushBackTimeStep(f1ts2) # @@ -2529,12 +2551,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m1.getNumberOfCells()) ; a.iota(1000.+off) ; a.setInfoOnComponents(["xx [m]"]) fCell2.setArray(a) fCell2.setTime(*t) - fCell2.checkCoherency() + fCell2.checkConsistencyLight() a=DataArrayDouble(m1.getNumberOfNodes()) ; a.iota(1000+off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode2.setArray(a) fNode2.setTime(*t) - fNode2.checkCoherency() + fNode2.checkConsistencyLight() f1ts1.setFieldNoProfileSBT(fCell2) ; ffs2_1.pushBackTimeStep(f1ts1) f1ts2.setFieldNoProfileSBT(fNode2) ; ffs2_2.pushBackTimeStep(f1ts2) # TimeStep 1 @@ -2544,12 +2566,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m0.getNumberOfCells()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) fCell1.setArray(a) fCell1.setTime(*t) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() a=DataArrayDouble(m0.getNumberOfNodes()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode1.setArray(a) fNode1.setTime(*t) - fNode1.checkCoherency() + fNode1.checkConsistencyLight() f1ts1.setFieldNoProfileSBT(fCell1) ; ffs1_1.pushBackTimeStep(f1ts1) f1ts2.setFieldNoProfileSBT(fNode1) ; ffs1_2.pushBackTimeStep(f1ts2) # @@ -2558,12 +2580,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m1.getNumberOfCells()) ; a.iota(1000.+off) ; a.setInfoOnComponents(["xx [m]"]) fCell2.setArray(a) fCell2.setTime(*t) - fCell2.checkCoherency() + fCell2.checkConsistencyLight() a=DataArrayDouble(m1.getNumberOfNodes()) ; a.iota(1000+off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode2.setArray(a) fNode2.setTime(*t) - fNode2.checkCoherency() + fNode2.checkConsistencyLight() f1ts1.setFieldNoProfileSBT(fCell2) ; ffs2_1.pushBackTimeStep(f1ts1) f1ts2.setFieldNoProfileSBT(fNode2) ; ffs2_2.pushBackTimeStep(f1ts2) # TimeStep 2 @@ -2573,12 +2595,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m0.getNumberOfCells()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) fCell1.setArray(a) fCell1.setTime(*t) - fCell1.checkCoherency() + fCell1.checkConsistencyLight() a=DataArrayDouble(m0.getNumberOfNodes()) ; a.iota(off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode1.setArray(a) fNode1.setTime(*t) - fNode1.checkCoherency() + fNode1.checkConsistencyLight() f1ts1.setFieldNoProfileSBT(fCell1) ; ffs1_1.pushBackTimeStep(f1ts1) f1ts2.setFieldNoProfileSBT(fNode1) ; ffs1_2.pushBackTimeStep(f1ts2) # @@ -2587,12 +2609,12 @@ class MEDLoaderTest4(unittest.TestCase): a=DataArrayDouble(m1.getNumberOfCells()) ; a.iota(1000.+off) ; a.setInfoOnComponents(["xx [m]"]) fCell2.setArray(a) fCell2.setTime(*t) - fCell2.checkCoherency() + fCell2.checkConsistencyLight() a=DataArrayDouble(m1.getNumberOfNodes()) ; a.iota(1000+off) ; a.setInfoOnComponents(["xx [m]"]) a=a.negate() fNode2.setArray(a) fNode2.setTime(*t) - fNode2.checkCoherency() + fNode2.checkConsistencyLight() f1ts1.setFieldNoProfileSBT(fCell2) ; ffs2_1.pushBackTimeStep(f1ts1) f1ts2.setFieldNoProfileSBT(fNode2) ; ffs2_2.pushBackTimeStep(f1ts2) # @@ -2642,37 +2664,37 @@ class MEDLoaderTest4(unittest.TestCase): m.allocateCells() m.insertNextCell(NORM_QUAD4,[0,3,4,1]) m.insertNextCell(NORM_QUAD4,[1,4,5,2]) - m.checkCoherency1() + m.checkConsistency() # t=(1.1,0,-1) f=MEDCouplingFieldDouble(ON_GAUSS_NE) ; f.setTime(*t) ; f.setMesh(m) f.setArray(DataArrayDouble([3.,5.,7.,6.,2.,3.,11.,8.])) f.setName(fieldName1) - f.checkCoherency() - MEDLoader.WriteField(fname,f,True) + f.checkConsistencyLight() + WriteField(fname,f,True) f2=MEDCouplingFieldDouble(ON_CELLS) ; f2.setTime(*t) ; f2.setMesh(m) f2.setArray(DataArrayDouble([7.,11.],2,1)) f2.setName(fieldName2) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f2) + WriteFieldUsingAlreadyWrittenMesh(fname,f2) f3=MEDCouplingFieldDouble(ON_NODES) ; f3.setTime(*t) ; f3.setMesh(m) f3.setArray(DataArrayDouble([1.,2.,4.,1.,2.,4.],6,1)) f3.setName(fieldName3) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f3) + WriteFieldUsingAlreadyWrittenMesh(fname,f3) # t=(2.1,1,-1) f=MEDCouplingFieldDouble(ON_GAUSS_NE) ; f.setTime(*t) ; f.setMesh(m) f.setArray(DataArrayDouble([7.,6.,3.,5.,11.,8.,2.,3.])) f.setName(fieldName1) - f.checkCoherency() - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f) + f.checkConsistencyLight() + WriteFieldUsingAlreadyWrittenMesh(fname,f) f2=MEDCouplingFieldDouble(ON_CELLS) ; f2.setTime(*t) ; f2.setMesh(m) f2.setArray(DataArrayDouble([11.,7.],2,1)) f2.setName(fieldName2) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f2) + WriteFieldUsingAlreadyWrittenMesh(fname,f2) f3=MEDCouplingFieldDouble(ON_NODES) ; f3.setTime(*t) ; f3.setMesh(m) f3.setArray(DataArrayDouble([4.,2.,1.,4.,2.,1.],6,1)) f3.setName(fieldName3) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f3) + WriteFieldUsingAlreadyWrittenMesh(fname,f3) ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() fields=MEDFileFields(fname,False) @@ -2719,12 +2741,13 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([0,1]))) self.assertTrue(a9) # no copy here - for i in xrange(1,2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 2): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass vExp0=[DataArrayDouble([7.,11.]),DataArrayDouble([11.,7.])] vExp1=[DataArrayDouble([3.,5.,7.,6.,2.,3.,11.,8.]),DataArrayDouble([7.,6.,3.,5.,11.,8.,2.,3.])] - for i in xrange(2): + for i in range(2): f=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -2756,39 +2779,39 @@ class MEDLoaderTest4(unittest.TestCase): m.allocateCells() m.insertNextCell(NORM_QUAD4,[0,3,4,1]) m.insertNextCell(NORM_QUAD4,[1,4,5,2]) - m.checkCoherency1() + m.checkConsistency() # t=(1.1,0,-1) f=MEDCouplingFieldDouble(ON_GAUSS_PT) ; f.setTime(*t) ; f.setMesh(m) f.setGaussLocalizationOnType(NORM_QUAD4,[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.2,0.2,0.8,0.8],[0.7,0.3]) f.setArray(DataArrayDouble([3.,5.,4.,6.])) ; f.getArray().setInfoOnComponents(["Smth"]) f.setName(fieldName1) - f.checkCoherency() - MEDLoader.WriteField(fname,f,True) + f.checkConsistencyLight() + WriteField(fname,f,True) f2=MEDCouplingFieldDouble(ON_CELLS) ; f2.setTime(*t) ; f2.setMesh(m) f2.setArray(DataArrayDouble([7.,11.],2,1)) f2.setName(fieldName2) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f2) + WriteFieldUsingAlreadyWrittenMesh(fname,f2) f3=MEDCouplingFieldDouble(ON_NODES) ; f3.setTime(*t) ; f3.setMesh(m) f3.setArray(DataArrayDouble([1.,2.,4.,1.,2.,4.],6,1)) f3.setName(fieldName3) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f3) + WriteFieldUsingAlreadyWrittenMesh(fname,f3) # t=(2.1,1,-1) f=MEDCouplingFieldDouble(ON_GAUSS_PT) ; f.setTime(*t) ; f.setMesh(m) f.setGaussLocalizationOnType(NORM_QUAD4,[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.2,0.2,0.8,0.8],[0.7,0.3]) f.setArray(DataArrayDouble([5.,3.,6.,4.])) ; f.getArray().setInfoOnComponents(["Smth"]) f.setName(fieldName1) - f.checkCoherency() - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f) + f.checkConsistencyLight() + WriteFieldUsingAlreadyWrittenMesh(fname,f) f2=MEDCouplingFieldDouble(ON_CELLS) ; f2.setTime(*t) ; f2.setMesh(m) f2.setArray(DataArrayDouble([11.,7.],2,1)) f2.setName(fieldName2) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f2) + WriteFieldUsingAlreadyWrittenMesh(fname,f2) f3=MEDCouplingFieldDouble(ON_NODES) ; f3.setTime(*t) ; f3.setMesh(m) f3.setArray(DataArrayDouble([4.,2.,1.,4.,2.,1.],6,1)) f3.setName(fieldName3) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f3) + WriteFieldUsingAlreadyWrittenMesh(fname,f3) ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() fields=MEDFileFields(fname,False) @@ -2836,12 +2859,13 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([0,1]))) self.assertTrue(a9) # no copy here - for i in xrange(1,2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 2): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass vExp0=[DataArrayDouble([7.,11.]),DataArrayDouble([11.,7.])] vExp1=[DataArrayDouble([3.,5.,4.,6.]),DataArrayDouble([5.,3.,6.,4.])] - for i in xrange(2): + for i in range(2): f=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -2904,12 +2928,13 @@ class MEDLoaderTest4(unittest.TestCase): a8,a9=mml2.retrieveNumberIdsOnCells() self.assertTrue(a8.isEqual(DataArrayInt([0,1]))) self.assertTrue(a9) # no copy here - for i in xrange(1,2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 2): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass vExp0=[DataArrayDouble([7.,11.]),DataArrayDouble([11.,7.])] vExp1=[DataArrayDouble([3.,5.,4.,6.]),DataArrayDouble([5.,3.,6.,4.])] - for i in xrange(2): + for i in range(2): f=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) # no load needed here v=mml.buildDataArray(fsst,fields,f.getUndergroundDataArray()) @@ -2948,7 +2973,7 @@ class MEDLoaderTest4(unittest.TestCase): f.setName(fieldName) arr=DataArrayDouble(24) ; arr.iota() ; arr.setInfoOnComponents(["AStr"]) f.setArray(arr) - f.checkCoherency() + f.checkConsistencyLight() f1ts=MEDFileField1TS() ; f1ts.setFieldNoProfileSBT(f) ff.pushBackTimeStep(f1ts) # time 1 @@ -2957,7 +2982,7 @@ class MEDLoaderTest4(unittest.TestCase): f.setName(fieldName) arr=DataArrayDouble(24) ; arr.iota() ; arr.reverse() ; arr.setInfoOnComponents(["AStr"]) f.setArray(arr) - f.checkCoherency() + f.checkConsistencyLight() f1ts=MEDFileField1TS() ; f1ts.setFieldNoProfileSBT(f) ff.pushBackTimeStep(f1ts) # @@ -3015,10 +3040,11 @@ class MEDLoaderTest4(unittest.TestCase): a12,a13=mml2.retrieveNumberIdsOnNodes() self.assertTrue(not a12) self.assertTrue(a13) # no copy here - for i in xrange(1,2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 2): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(2): + for i in range(2): f=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -3066,16 +3092,16 @@ class MEDLoaderTest4(unittest.TestCase): f0.setName(fieldName0) ; f0.setTime(*t) da=m.getCoords().magnitude() ; da.setInfoOnComponents(["zeInfo"]) f0.setArray(da) - f0.checkCoherency() + f0.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldNoProfileSBT(f0) fmts0.pushBackTimeStep(f1ts) # f1=MEDCouplingFieldDouble(ON_CELLS) ; f1.setMesh(m) f1.setName(fieldName1) ; f1.setTime(*t) - da=m.getBarycenterAndOwner().magnitude() ; da.setInfoOnComponents(["zeInfoCell"]) + da=m.computeCellCenterOfMass().magnitude() ; da.setInfoOnComponents(["zeInfoCell"]) f1.setArray(da) - f1.checkCoherency() + f1.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldNoProfileSBT(f1) fmts1.pushBackTimeStep(f1ts) @@ -3084,16 +3110,16 @@ class MEDLoaderTest4(unittest.TestCase): f2.setName(fieldName2) ; f2.setTime(*t) da=m.getCoords().magnitude()[pfl2] ; da.setInfoOnComponents(["zzzz"]) f2.setArray(da) - f2.checkCoherency() + f2.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldProfile(f2,mm,0,pfl2) fmts2.pushBackTimeStep(f1ts) # f3=MEDCouplingFieldDouble(ON_CELLS) ; mTmp=m[pfl3] ; mTmp.setName(m.getName()) ; f3.setMesh(mTmp) f3.setName(fieldName3) ; f3.setTime(*t) - da=mTmp.getBarycenterAndOwner().magnitude() ; da.setInfoOnComponents(["abcdefg"]) + da=mTmp.computeCellCenterOfMass().magnitude() ; da.setInfoOnComponents(["abcdefg"]) f3.setArray(da) - f3.checkCoherency() + f3.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldProfile(f3,mm,0,pfl3) fmts3.pushBackTimeStep(f1ts) @@ -3103,16 +3129,16 @@ class MEDLoaderTest4(unittest.TestCase): f0.setName(fieldName0) ; f0.setTime(*t) da=m.getCoords().magnitude() ; da.reverse() ; da.setInfoOnComponents(["zeInfo"]) f0.setArray(da) - f0.checkCoherency() + f0.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldNoProfileSBT(f0) fmts0.pushBackTimeStep(f1ts) # f1=MEDCouplingFieldDouble(ON_CELLS) ; f1.setMesh(m) f1.setName(fieldName1) ; f1.setTime(*t) - da=m.getBarycenterAndOwner().magnitude() ; da.reverse() ; da.setInfoOnComponents(["zeInfoCell"]) + da=m.computeCellCenterOfMass().magnitude() ; da.reverse() ; da.setInfoOnComponents(["zeInfoCell"]) f1.setArray(da) - f1.checkCoherency() + f1.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldNoProfileSBT(f1) fmts1.pushBackTimeStep(f1ts) @@ -3121,16 +3147,16 @@ class MEDLoaderTest4(unittest.TestCase): f2.setName(fieldName2) ; f2.setTime(*t) da=m.getCoords().magnitude()[pfl2] ; da.reverse() ; da.setInfoOnComponents(["zzzz"]) f2.setArray(da) - f2.checkCoherency() + f2.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldProfile(f2,mm,0,pfl2) fmts2.pushBackTimeStep(f1ts) # f3=MEDCouplingFieldDouble(ON_CELLS) ; mTmp=m[pfl3] ; mTmp.setName(m.getName()) ; f3.setMesh(mTmp) f3.setName(fieldName3) ; f3.setTime(*t) - da=mTmp.getBarycenterAndOwner().magnitude() ; da.reverse() ; da.setInfoOnComponents(["abcdefg"]) + da=mTmp.computeCellCenterOfMass().magnitude() ; da.reverse() ; da.setInfoOnComponents(["abcdefg"]) f3.setArray(da) - f3.checkCoherency() + f3.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldProfile(f3,mm,0,pfl3) fmts3.pushBackTimeStep(f1ts) @@ -3200,10 +3226,11 @@ class MEDLoaderTest4(unittest.TestCase): a12,a13=mml2.retrieveNumberIdsOnNodes() self.assertTrue(a12.isEqual(DataArrayInt([0,10,20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200,210,220,230,240]))) self.assertTrue(a13) # no copy here - for i in xrange(1,2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 2): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(2): + for i in range(2): f=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -3252,7 +3279,8 @@ class MEDLoaderTest4(unittest.TestCase): a12,a13=mml2.retrieveNumberIdsOnNodes() self.assertTrue(a12.isEqual(DataArrayInt([50,60,70,100,110,120,150,160,170,200,210,220]))) self.assertTrue(not a13) # copy here - for i in xrange(2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(2): f=allFMTSLeavesPerCommonSupport1[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -3331,7 +3359,7 @@ class MEDLoaderTest4(unittest.TestCase): f0.setName(fieldName0) ; f0.setTime(*t) da=DataArrayDouble(9) ; da.iota() ; da.setInfoOnComponents(["zeInfo"]) f0.setArray(da) - f0.checkCoherency() + f0.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldNoProfileSBT(f0) fmts0.pushBackTimeStep(f1ts) @@ -3341,7 +3369,7 @@ class MEDLoaderTest4(unittest.TestCase): f0.setName(fieldName0) ; f0.setTime(*t) da=DataArrayDouble(9) ; da.iota() ; da.reverse() ; da.setInfoOnComponents(["zeInfo"]) f0.setArray(da) - f0.checkCoherency() + f0.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldNoProfileSBT(f0) fmts0.pushBackTimeStep(f1ts) @@ -3401,6 +3429,7 @@ class MEDLoaderTest4(unittest.TestCase): a12,a13=mml2.retrieveNumberIdsOnNodes() self.assertTrue(not a12) self.assertTrue(a13) # no copy here + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) # f=allFMTSLeavesPerCommonSupport1[0][0][0][0] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) @@ -3444,7 +3473,7 @@ class MEDLoaderTest4(unittest.TestCase): f0.setName(fieldName0) ; f0.setTime(*t) da=DataArrayDouble(4) ; da.iota() ; da.setInfoOnComponents(["zeInfo"]) f0.setArray(da) - f0.checkCoherency() + f0.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldProfile(f0,mm,-1,pfl) fmts0.pushBackTimeStep(f1ts) @@ -3454,7 +3483,7 @@ class MEDLoaderTest4(unittest.TestCase): f0.setName(fieldName0) ; f0.setTime(*t) da=DataArrayDouble(4) ; da.iota() ; da.reverse() ; da.setInfoOnComponents(["zeInfo"]) f0.setArray(da) - f0.checkCoherency() + f0.checkConsistencyLight() f1ts=MEDFileField1TS() f1ts.setFieldProfile(f0,mm,-1,pfl) fmts0.pushBackTimeStep(f1ts) @@ -3512,6 +3541,7 @@ class MEDLoaderTest4(unittest.TestCase): a12,a13=mml2.retrieveNumberIdsOnNodes() self.assertTrue(not a12) self.assertTrue(a13) # no copy here + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) # f=allFMTSLeavesPerCommonSupport1[0][0][0][0] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) @@ -3602,7 +3632,8 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a.isEqual(arr0,1e-12)) self.assertTrue(b.isEqual(arr1,1e-12)) self.assertTrue(c.isEqual(arr2,1e-12)) - for i in xrange(2): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(2): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -3667,7 +3698,7 @@ class MEDLoaderTest4(unittest.TestCase): pfl1=DataArrayInt([1,3,5]) ; pfl1.setName(pflName1) tmp=m0[pfl1] ; f2.setMesh(tmp) f2.setGaussLocalizationOnType(NORM_QUAD4,[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[-0.5,-0.5,0.5,-0.5,0.5,0.5,-0.5,0.5,0.,0.],[0.1,0.1,0.1,0.1,0.6]) - f2.checkCoherency() + f2.checkConsistencyLight() f1ts.setFieldProfile(f2,m,0,pfl1) fmts.pushBackTimeStep(f1ts) # @@ -3721,6 +3752,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([2,10,12,2,12,13]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) a6,a7=mml2.retrieveFamilyIdsOnCells() self.assertTrue(a6.isEqual(DataArrayInt([-7,-8]))) self.assertTrue(not a7) # copy here because profile on cells @@ -3733,6 +3765,7 @@ class MEDLoaderTest4(unittest.TestCase): a12,a13=mml2.retrieveNumberIdsOnNodes() self.assertTrue(a12.isEqual(DataArrayInt([203,204,205,206,207,208,209,210,211,212,213,214,215,216,217]))) self.assertTrue(a13) # no copy here + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) fff0=allFMTSLeavesPerCommonSupport1[0][0][0][0] fsst=MEDFileField1TSStructItem.BuildItemFrom(fff0,mst) fff0.loadArraysIfNecessary() @@ -3758,6 +3791,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([4,1,4,5,2,4,4,7,8,5,4,7,10,11,8]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) a6,a7=mml2.retrieveFamilyIdsOnCells() self.assertTrue(a6.isEqual(DataArrayInt([-2,-4,-6]))) self.assertTrue(not a7) # copy here because profile on cells @@ -3793,6 +3827,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([2,10,12,2,12,13]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) a6,a7=mml2.retrieveFamilyIdsOnCells() self.assertTrue(a6.isEqual(DataArrayInt([-7,-8]))) self.assertTrue(not a7) # copy here because profile on cells @@ -3929,6 +3964,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,4,2,1,4,0,3,4,1,1,5]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) a6,a7=mml2.retrieveFamilyIdsOnCells() self.assertTrue(a6.isEqual(DataArrayInt([-3,-7,13]))) self.assertTrue(not a7) # copy here because profile on cells @@ -3946,6 +3982,7 @@ class MEDLoaderTest4(unittest.TestCase): fsst=MEDFileField1TSStructItem.BuildItemFrom(fff0,mst) fff0.loadArraysIfNecessary() v=mml2.buildDataArray(fsst,fields,fff0.getUndergroundDataArray()) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertEqual(fff0.getName(),"f0Node") self.assertEqual(v.getHiddenCppPointer(),fff0.getUndergroundDataArray().getHiddenCppPointer()) vExp=DataArrayDouble([0.,1.,2.,3.,4.,5.]) @@ -3955,6 +3992,7 @@ class MEDLoaderTest4(unittest.TestCase): fsst=MEDFileField1TSStructItem.BuildItemFrom(fff1,mst) fff1.loadArraysIfNecessary() v=mml2.buildDataArray(fsst,fields,fff1.getUndergroundDataArray()) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertEqual(fff1.getName(),"f1Node") self.assertEqual(v.getHiddenCppPointer(),fff1.getUndergroundDataArray().getHiddenCppPointer()) vExp=DataArrayDouble([5.,4.,3.,2.,1.,0.]) @@ -3981,7 +4019,7 @@ class MEDLoaderTest4(unittest.TestCase): f0=MEDCouplingFieldDouble(ON_NODES) ; f0.setMesh(m0) ; f0.setName("f0NoPfl") arr0=DataArrayDouble([0.,1.,2.,3.,1.,1.5,2.2,3.1,2.,2.2,3.,3.1,3.,3.1,3.5,4.]) f0.setArray(arr0) - f0.checkCoherency() + f0.checkConsistencyLight() f1ts0.setFieldNoProfileSBT(f0) self.assertEqual(f1ts0.getMeshName(),"mesh") # @@ -4047,6 +4085,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) # fff0=allFMTSLeavesPerCommonSupport1[0][0][0][0] fsst=MEDFileField1TSStructItem.BuildItemFrom(fff0,mst) @@ -4080,7 +4119,7 @@ class MEDLoaderTest4(unittest.TestCase): pass def test28(self): - """ This test defines 2 fields f0,f1,f2,f3 lying on an unstructured mesh whith cells including NORM_POINT1. + """ This test defines 2 fields f0,f1,f2,f3 lying on an unstructured mesh with cells including NORM_POINT1. Both f0 and f1 are on NODES and f2 and f3 are on cells. f1 and f2 share the same support. f0 is on a nodal support that is not matchable with any cells (including NORM_POINT1) This test is a more aggressive version of test26. @@ -4190,6 +4229,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,4,2,1,4,0,3,4,1,1,3,1,5]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) a6,a7=mml2.retrieveFamilyIdsOnCells() self.assertTrue(a6.isEqual(DataArrayInt([-3,-7,-12,13]))) self.assertTrue(not a7) # copy here because profile on cells @@ -4207,6 +4247,7 @@ class MEDLoaderTest4(unittest.TestCase): fsst=MEDFileField1TSStructItem.BuildItemFrom(fff0,mst) fff0.loadArraysIfNecessary() v=mml2.buildDataArray(fsst,fields,fff0.getUndergroundDataArray()) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertEqual(fff0.getName(),"f0Node") self.assertEqual(v.getHiddenCppPointer(),fff0.getUndergroundDataArray().getHiddenCppPointer()) vExp=DataArrayDouble([0.,1.,2.,3.,4.,5.]) @@ -4224,10 +4265,12 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,2,6,3,3,6,7,3,4,0,4,5,1,4,1,5,6,2,4,6,8,9,7,4,8,10,11,9]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) fff1=allFMTSLeavesPerCommonSupport1[0][0][0][0] fsst=MEDFileField1TSStructItem.BuildItemFrom(fff1,mst) fff1.loadArraysIfNecessary() v=mml2.buildDataArray(fsst,fields,fff1.getUndergroundDataArray()) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertEqual(fff1.getName(),"f2Cell") self.assertNotEqual(v.getHiddenCppPointer(),fff0.getUndergroundDataArray().getHiddenCppPointer()) vExp=DataArrayDouble([2,3,0,1,4,5]) @@ -4257,6 +4300,7 @@ class MEDLoaderTest4(unittest.TestCase): fsst=MEDFileField1TSStructItem.BuildItemFrom(fff3,mst) fff3.loadArraysIfNecessary() v=mml2.buildDataArray(fsst,fields,fff3.getUndergroundDataArray()) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertEqual(fff3.getName(),"f3Cell") self.assertNotEqual(v.getHiddenCppPointer(),fff0.getUndergroundDataArray().getHiddenCppPointer()) vExp=DataArrayDouble([-1.1,-3.1,5,4,3,2,1,0]) @@ -4276,7 +4320,7 @@ class MEDLoaderTest4(unittest.TestCase): fCell=MEDCouplingFieldDouble(ON_CELLS) ; fCell.setName("fCell") arrCell=DataArrayDouble([7.]) ; arrCell.setInfoOnComponent(0,"smth") ; fCell.setArray(arrCell) fCell.setMesh(m) - MEDLoader.WriteField(fname,fCell,True) + WriteField(fname,fCell,True) refCoo=[-1.,-1.,-1.,-1.,1.,-1.,1.,1.,-1.,1.,-1.,-1.,-1.,-1.,1.,-1.,1.,1.,1.,1.,1.,1.,-1.,1.,-1.,0.,-1.,0.,1.,-1.,1.,0.,-1.,0.,-1.,-1.,-1.,0.,1.,0.,1.,1.,1.,0.,1.,0.,-1.,1.,-1.,-1.,0.,-1.,1.,0.,1.,1.,0.,1.,-1.,0.,0.,0.,-1.,-1.,0.,0.,0.,1.,0.,1.,0.,0.,0.,-1.,0.,0.,0.,1.,0.,0.,0.] weights=[0.1714677640603571,0.27434842249657115,0.1714677640603571,0.27434842249657115,0.43895747599451346,0.27434842249657115,0.1714677640603571,0.27434842249657115,0.1714677640603571,0.27434842249657115,0.43895747599451346,0.27434842249657115,0.43895747599451346,0.7023319615912209,0.43895747599451346,0.27434842249657115,0.43895747599451346,0.27434842249657115,0.1714677640603571,0.27434842249657115,0.1714677640603571,0.27434842249657115,0.43895747599451346,0.27434842249657115,0.1714677640603571,0.27434842249657115,0.1714677640603571] gCoords=[-0.774596669241483,-0.774596669241483,-0.774596669241483,-0.774596669241483,-0.774596669241483,0.0,-0.774596669241483,-0.774596669241483,0.774596669241483,-0.774596669241483,0.0,-0.774596669241483,-0.774596669241483,0.0,0.0,-0.774596669241483,0.0,0.774596669241483,-0.774596669241483,0.774596669241483,-0.774596669241483,-0.774596669241483,0.774596669241483,0.0,-0.774596669241483,0.774596669241483,0.774596669241483,0.0,-0.774596669241483,-0.774596669241483,0.0,-0.774596669241483,0.0,0.0,-0.774596669241483,0.774596669241483,0.0,0.0,-0.774596669241483,0.0,0.0,0.0,0.0,0.0,0.774596669241483,0.0,0.774596669241483,-0.774596669241483,0.0,0.774596669241483,0.0,0.0,0.774596669241483,0.774596669241483,0.774596669241483,-0.774596669241483,-0.774596669241483,0.774596669241483,-0.774596669241483,0.0,0.774596669241483,-0.774596669241483,0.774596669241483,0.774596669241483,0.0,-0.774596669241483,0.774596669241483,0.0,0.0,0.774596669241483,0.0,0.774596669241483,0.774596669241483,0.774596669241483,-0.774596669241483,0.774596669241483,0.774596669241483,0.0,0.774596669241483,0.774596669241483,0.774596669241483] @@ -4285,7 +4329,7 @@ class MEDLoaderTest4(unittest.TestCase): fGauss.setGaussLocalizationOnType(NORM_HEXA27,refCoo,gCoords,weights) arrGauss=DataArrayDouble(fGauss.getNumberOfTuplesExpected()) ; arrGauss.setInfoOnComponent(0,"gaussc") ; arrGauss.iota() fGauss.setArray(arrGauss) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,fGauss) + WriteFieldUsingAlreadyWrittenMesh(fname,fGauss) ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() fields=MEDFileFields(fname,False) @@ -4319,6 +4363,7 @@ class MEDLoaderTest4(unittest.TestCase): mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertTrue(a0.isEqual(coo,1e-12)) self.assertTrue(a1.isEqual(DataArrayByte([29]))) self.assertTrue(a2.isEqual(DataArrayInt([0]))) @@ -4330,6 +4375,7 @@ class MEDLoaderTest4(unittest.TestCase): fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() v=mml2.buildDataArray(fsst,fields,ffCell.getUndergroundDataArray()) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertEqual(v.getHiddenCppPointer(),ffCell.getUndergroundDataArray().getHiddenCppPointer()) self.assertEqual(ffCell.getName(),"fCell") self.assertTrue(v.isEqual(arrCell,1e-12)) ; self.assertTrue(v.isEqualWithoutConsideringStr(DataArrayDouble([7.]),1e-12)) ; self.assertEqual(v.getInfoOnComponents(),["smth"]) @@ -4339,9 +4385,10 @@ class MEDLoaderTest4(unittest.TestCase): fsst=MEDFileField1TSStructItem.BuildItemFrom(ffGauss,mst) ffGauss.loadArraysIfNecessary() v=mml2.buildDataArray(fsst,fields,ffGauss.getUndergroundDataArray()) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertEqual(v.getHiddenCppPointer(),ffGauss.getUndergroundDataArray().getHiddenCppPointer()) self.assertEqual(ffGauss.getName(),"fGauss") - self.assertTrue(v.isEqual(arrGauss,1e-12)) ; self.assertTrue(v.isEqualWithoutConsideringStr(DataArrayDouble(range(27)),1e-12)) ; self.assertEqual(v.getInfoOnComponents(),["gaussc"]) + self.assertTrue(v.isEqual(arrGauss, 1e-12)) ; self.assertTrue(v.isEqualWithoutConsideringStr(DataArrayDouble(list(range(27))), 1e-12)) ; self.assertEqual(v.getInfoOnComponents(), ["gaussc"]) ffGauss=allFMTSLeavesPerCommonSupport1[0][0][1][0] pass @@ -4370,7 +4417,7 @@ class MEDLoaderTest4(unittest.TestCase): # fmts0=MEDFileFieldMultiTS() fmts1=MEDFileFieldMultiTS() - for i in xrange(30): + for i in range(30): f1ts=MEDFileField1TS() fFaces=MEDCouplingFieldDouble(ON_CELLS) ; fFaces.setName("FieldOnFaces") arr=DataArrayDouble(98) ; arr.iota() ; arr[i]=100. @@ -4432,7 +4479,8 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a.isEqual(arrX,1e-12)) self.assertTrue(b.isEqual(arrY,1e-12)) self.assertTrue(c.isEqual(arrZ,1e-12)) - for i in xrange(30): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(30): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -4456,7 +4504,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([4,0,12,15,3,4,12,24,27,15,4,24,36,39,27,4,36,48,51,39,4,3,15,18,6,4,15,27,30,18,4,27,39,42,30,4,39,51,54,42,4,6,18,21,9,4,18,30,33,21,4,30,42,45,33,4,42,54,57,45,4,1,13,16,4,4,13,25,28,16,4,25,37,40,28,4,37,49,52,40,4,4,16,19,7,4,16,28,31,19,4,28,40,43,31,4,40,52,55,43,4,7,19,22,10,4,19,31,34,22,4,31,43,46,34,4,43,55,58,46,4,2,14,17,5,4,14,26,29,17,4,26,38,41,29,4,38,50,53,41,4,5,17,20,8,4,17,29,32,20,4,29,41,44,32,4,41,53,56,44,4,8,20,23,11,4,20,32,35,23,4,32,44,47,35,4,44,56,59,47,4,0,12,13,1,4,12,24,25,13,4,24,36,37,25,4,36,48,49,37,4,1,13,14,2,4,13,25,26,14,4,25,37,38,26,4,37,49,50,38,4,3,15,16,4,4,15,27,28,16,4,27,39,40,28,4,39,51,52,40,4,4,16,17,5,4,16,28,29,17,4,28,40,41,29,4,40,52,53,41,4,6,18,19,7,4,18,30,31,19,4,30,42,43,31,4,42,54,55,43,4,7,19,20,8,4,19,31,32,20,4,31,43,44,32,4,43,55,56,44,4,9,21,22,10,4,21,33,34,22,4,33,45,46,34,4,45,57,58,46,4,10,22,23,11,4,22,34,35,23,4,34,46,47,35,4,46,58,59,47,4,0,1,4,3,4,3,4,7,6,4,6,7,10,9,4,1,2,5,4,4,4,5,8,7,4,7,8,11,10,4,12,13,16,15,4,15,16,19,18,4,18,19,22,21,4,13,14,17,16,4,16,17,20,19,4,19,20,23,22,4,24,25,28,27,4,27,28,31,30,4,30,31,34,33,4,25,26,29,28,4,28,29,32,31,4,31,32,35,34,4,36,37,40,39,4,39,40,43,42,4,42,43,46,45,4,37,38,41,40,4,40,41,44,43,4,43,44,47,46,4,48,49,52,51,4,51,52,55,54,4,54,55,58,57,4,49,50,53,52,4,52,53,56,55,4,55,56,59,58]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(30): + for i in range(30): ffCell=allFMTSLeavesPerCommonSupport1[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -4489,7 +4537,7 @@ class MEDLoaderTest4(unittest.TestCase): fmts0=MEDFileFieldMultiTS() fmts1=MEDFileFieldMultiTS() pfl=DataArrayInt(11) ; pfl.iota() ; pfl.setName("PflOnHECA8") - for i in xrange(30): + for i in range(30): f1ts=MEDFileField1TS() fFaces=MEDCouplingFieldDouble(ON_CELLS) ; fFaces.setName("FieldOnCells") arr=DataArrayDouble(11) ; arr.iota() ; arr[i%11]=100. @@ -4537,6 +4585,7 @@ class MEDLoaderTest4(unittest.TestCase): mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) # here UMesh is important ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertTrue(not ncc)# here ncc=False because the coordinates are not in ms neither in children. This is the most important line in the test. self.assertTrue(a0.isEqual(DataArrayDouble([0.,0.,0.,1.,0.,0.,2.,0.,0.,0.,1.,0.,1.,1.,0.,2.,1.,0.,0.,2.,0.,1.,2.,0.,2.,2.,0.,0.,3.,0.,1.,3.,0.,2.,3.,0.,0.,0.,1.,1.,0.,1.,2.,0.,1.,0.,1.,1.,1.,1.,1.,2.,1.,1.,0.,2.,1.,1.,2.,1.,2.,2.,1.,0.,3.,1.,1.,3.,1.,2.,3.,1.,0.,0.,2.,1.,0.,2.,2.,0.,2.,0.,1.,2.,1.,1.,2.,2.,1.,2.,0.,2.,2.,1.,2.,2.,2.,2.,2.,0.,3.,2.,1.,3.,2.,2.,3.,2.,0.,0.,3.,1.,0.,3.,2.,0.,3.,0.,1.,3.,1.,1.,3.,2.,1.,3.,0.,2.,3.,1.,2.,3.,2.,2.,3.,0.,3.,3.,1.,3.,3.,2.,3.,3.,0.,0.,4.,1.,0.,4.,2.,0.,4.,0.,1.,4.,1.,1.,4.,2.,1.,4.,0.,2.,4.,1.,2.,4.,2.,2.,4.,0.,3.,4.,1.,3.,4.,2.,3.,4.],60,3),1e-12)) self.assertTrue(a1.isEqual(DataArrayByte([12,12,12,12,12,12,12,12,12,12,12]))) @@ -4544,7 +4593,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([8,1,0,3,4,13,12,15,16,8,2,1,4,5,14,13,16,17,8,4,3,6,7,16,15,18,19,8,5,4,7,8,17,16,19,20,8,7,6,9,10,19,18,21,22,8,8,7,10,11,20,19,22,23,8,13,12,15,16,25,24,27,28,8,14,13,16,17,26,25,28,29,8,16,15,18,19,28,27,30,31,8,17,16,19,20,29,28,31,32,8,19,18,21,22,31,30,33,34]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(30): + for i in range(30): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -4579,7 +4628,7 @@ class MEDLoaderTest4(unittest.TestCase): fmts0=MEDFileFieldMultiTS() fmts1=MEDFileFieldMultiTS() pfl=DataArrayInt(31) ; pfl.iota() ; pfl.setName("PflOnQUAD4") - for i in xrange(30): + for i in range(30): f1ts=MEDFileField1TS() fFaces=MEDCouplingFieldDouble(ON_CELLS) ; fFaces.setName("FieldOnFaces") arr=DataArrayDouble(31) ; arr.iota() ; arr[i]=100. @@ -4642,7 +4691,8 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a.isEqual(arrX,1e-12)) self.assertTrue(b.isEqual(arrY,1e-12)) self.assertTrue(c.isEqual(arrZ,1e-12)) - for i in xrange(30): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(30): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -4658,6 +4708,7 @@ class MEDLoaderTest4(unittest.TestCase): mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) # here UMesh is important ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) self.assertTrue(ncc)# True because, the coords are computed by the implicit unstructured level -1 structured mesh self.assertTrue(a0.isEqual(DataArrayDouble([0.,0.,0.,1.,0.,0.,2.,0.,0.,0.,1.,0.,1.,1.,0.,2.,1.,0.,0.,2.,0.,1.,2.,0.,2.,2.,0.,0.,3.,0.,1.,3.,0.,2.,3.,0.,0.,0.,1.,1.,0.,1.,2.,0.,1.,0.,1.,1.,1.,1.,1.,2.,1.,1.,0.,2.,1.,1.,2.,1.,2.,2.,1.,0.,3.,1.,1.,3.,1.,2.,3.,1.,0.,0.,2.,1.,0.,2.,2.,0.,2.,0.,1.,2.,1.,1.,2.,2.,1.,2.,0.,2.,2.,1.,2.,2.,2.,2.,2.,0.,3.,2.,1.,3.,2.,2.,3.,2.,0.,0.,3.,1.,0.,3.,2.,0.,3.,0.,1.,3.,1.,1.,3.,2.,1.,3.,0.,2.,3.,1.,2.,3.,2.,2.,3.,0.,3.,3.,1.,3.,3.,2.,3.,3.,0.,0.,4.,1.,0.,4.,2.,0.,4.,0.,1.,4.,1.,1.,4.,2.,1.,4.,0.,2.,4.,1.,2.,4.,2.,2.,4.,0.,3.,4.,1.,3.,4.,2.,3.,4.],60,3),1e-12)) self.assertTrue(a1.isEqual(DataArrayByte([9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9]))) @@ -4668,7 +4719,7 @@ class MEDLoaderTest4(unittest.TestCase): a6,a7=mml2.retrieveFamilyIdsOnCells() self.assertTrue(a6 is None) self.assertTrue(a7) - for i in xrange(30): + for i in range(30): ffCell=allFMTSLeavesPerCommonSupport1[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -4700,13 +4751,13 @@ class MEDLoaderTest4(unittest.TestCase): mm=MEDFileUMesh() mm.setMeshAtLevel(0,m) mm.write(fname,2) - for i in xrange(15): + for i in range(15): fCell0=MEDCouplingFieldDouble(ON_CELLS) ; fCell0.setTime(float(i)+0.1,i,0) fCell0.setName(fieldName) ; fCell0.setMesh(m) arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota(0) ; arr[i%10]=100. fCell0.setArray(arr) ; arr.setInfoOnComponents(["Comp1 [m]"]) - fCell0.checkCoherency() - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,fCell0) + fCell0.checkConsistencyLight() + WriteFieldUsingAlreadyWrittenMesh(fname,fCell0) pass ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -4741,14 +4792,15 @@ class MEDLoaderTest4(unittest.TestCase): mml2=mml.prepare() self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() - self.assertTrue(not ncc)# false beacause 2D in MED file + self.assertTrue(not ncc)# false because 2D in MED file self.assertTrue(a0.isEqual(DataArrayDouble([(5.5,0.5,0),(5.5,-0.5,0),(6.5,0.5,0),(6.5,-0.5,0),(6.5,1.5,0),(7.5,0.5,0),(7.5,-0.5,0),(7.5,1.5,0),(7.5,2.5,0),(8.5,0.5,0),(8.5,-0.5,0),(8.5,1.5,0),(8.5,2.5,0),(8.5,3.5,0),(8.55,0.5,0),(8.55,-0.5,0),(8.55,1.5,0),(8.55,2.5,0),(8.55,3.5,0)]),1e-12)) self.assertTrue(a1.isEqual(DataArrayByte([9,9,9,9,9,9,7,7,7,7]))) self.assertTrue(a2.isEqual(DataArrayInt([0,5,10,15,20,25,30,35,40,45])))# the bug was here. self.assertTrue(a3.isEqual(DataArrayInt([4,1,0,2,3,4,3,2,5,6,4,2,4,7,5,4,6,5,9,10,4,5,7,11,9,4,7,8,12,11,4,10,9,14,15,4,9,11,16,14,4,11,12,17,16,4,12,13,18,17]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(15): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(15): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -4775,7 +4827,7 @@ class MEDLoaderTest4(unittest.TestCase): fieldName0="zeField0" fieldName1="zeField1" fs0=MEDFileFieldMultiTS() ; fs1=MEDFileFieldMultiTS() - for i in xrange(5): + for i in range(5): f=MEDFileField1TS() fNode=MEDCouplingFieldDouble(ON_GAUSS_PT) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName0) ; fNode.setMesh(m) @@ -4784,7 +4836,7 @@ class MEDLoaderTest4(unittest.TestCase): fNode.setGaussLocalizationOnCells([6,10,13],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2],[0.8,0.05,0.1,0.04]) fNode.setGaussLocalizationOnCells([11,12,14],[-1.,-1.,1.,-1.,1.,1.,-1.,1.],[0.5,0.5,0.7,0.7,0.1,0.1,0.2,0.2,0.3,0.3,0.4,0.4,0.8,0.8],[0.8,0.05,0.1,0.01,0.02,0.005,0.005]) arr=DataArrayDouble(2*(2*6+5*4+4*3+7*3)) ; arr.iota(0+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_0 [m]","Com2_0 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs0.pushBackTimeStep(f) # @@ -4792,13 +4844,13 @@ class MEDLoaderTest4(unittest.TestCase): fNode=MEDCouplingFieldDouble(ON_CELLS) ; fNode.setTime(float(i),i,0) fNode.setName(fieldName1) ; fNode.setMesh(m) arr=DataArrayDouble(2*16) ; arr.iota(300+1000*i) ; arr.rearrange(2) - fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkCoherency() + fNode.setArray(arr) ; arr.setInfoOnComponents(["Comp1_1 [m]","Com2_1 [s^2]"]) ; fNode.checkConsistencyLight() f.setFieldNoProfileSBT(fNode) fs1.pushBackTimeStep(f) pass mm.write(fname,2) fs0.write(fname,0) ; fs1.write(fname,0) - a0Exp=mm.getCoords().deepCpy() + a0Exp=mm.getCoords().deepCopy() del m,mm,fs0,fs1,f,fNode ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes(fname) ; ms.cartesianizeMe() @@ -4825,7 +4877,7 @@ class MEDLoaderTest4(unittest.TestCase): pass allFMTSLeavesToDisplay.append(allFMTSLeavesToDisplay2) pass - # Here 2 MED fields in input and at the end 5 ! 1+4 ! 4 fields have been built from zeField0 due to subspliting per dis / per geo type + # Here 2 MED fields in input and at the end 5 ! 1+4 ! 4 fields have been built from zeField0 due to subsplitting per dis / per geo type self.assertEqual(len(allFMTSLeavesToDisplay),1) self.assertEqual(len(allFMTSLeavesToDisplay[0]),5) allFMTSLeavesPerTimeSeries=MEDFileAnyTypeFieldMultiTS.SplitIntoCommonTimeSeries(sum(allFMTSLeavesToDisplay,[])) @@ -4833,7 +4885,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertEqual(len(allFMTSLeavesPerTimeSeries[0]),5) allFMTSLeavesPerCommonSupport=MEDFileAnyTypeFieldMultiTS.SplitPerCommonSupport(allFMTSLeavesPerTimeSeries[0],ms[ms.getMeshesNames()[0]]) self.assertEqual(len(allFMTSLeavesPerCommonSupport),5) - for i in xrange(5): + for i in range(5): self.assertEqual(len(allFMTSLeavesPerCommonSupport[i][0]),1) # mst=MEDFileMeshStruct.New(ms[0]) @@ -4850,10 +4902,11 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([4,1,0,5,6,4,2,1,6,7,4,3,2,7,8,4,4,3,8,9,4,6,5,10,11,4,7,6,11,12,4,8,7,12,13,4,9,8,13,14,4,11,10,15,16,4,12,11,16,17,4,13,12,17,18,4,14,13,18,19,4,16,15,20,21,4,17,16,21,22,4,18,17,22,23,4,19,18,23,24]))) # <- here the mesh is NOT renumbered : the mesh is equal to m self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1,5): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1, 5): self.assertTrue((fcscp.isDataSetSupportEqualToThePreviousOne(i,fields))) pass - for i in xrange(5): + for i in range(5): f=allFMTSLeavesPerCommonSupport[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(f,mst) f.loadArraysIfNecessary() @@ -4906,14 +4959,14 @@ class MEDLoaderTest4(unittest.TestCase): pass def test35(self): - """ Emulate MEDReader in // mode context. Here a Simple mesh having more nodes than really needed. This test focuses on that point particulary.""" + """ Emulate MEDReader in // mode context. Here a Simple mesh having more nodes than really needed. This test focuses on that point particularly.""" fname="ForMEDReader35.med" arrX=DataArrayDouble(7) ; arrX.iota() arrY=DataArrayDouble([0.,1.]) m=MEDCouplingCMesh() ; m.setCoords(arrX,arrY) ; m=m.buildUnstructured() ; m=m[[0,5,1,4,2,3]] ; m.changeSpaceDimension(3,0.) ; m.setName("Mesh") f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(m) ; f.setName("Field") ; f.setArray(DataArrayDouble([(0.1,1.1),(2.1,3.1),(4.1,5.1),(6.1,7.1),(8.1,9.1),(10.1,11.1)])) ; f.getArray().setInfoOnComponents(["aa","bbb"]) - MEDLoader.WriteUMesh(fname,m,True) - MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fname,f) + WriteUMesh(fname,m,True) + WriteFieldUsingAlreadyWrittenMesh(fname,f) ########## GO for reading in MEDReader,by not loading all. Mesh is fully loaded but not fields values ms=MEDFileMeshes() # here we reproduce what is done by ParaMEDFileMeshes.ParaNew ms.pushMesh(MEDFileUMesh.LoadPartOf(fname,"Mesh",[NORM_QUAD4],[0,2,1],-1,-1)); @@ -4967,6 +5020,7 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([4,1,0,4,5,4,3,2,6,7]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny().isEqual(DataArrayInt([0,1,5,6,7,8,12,13]))) f2=allFMTSLeavesPerCommonSupport[0][0][0][0] fsst=MEDFileField1TSStructItem.BuildItemFrom(f2,mst) f2.loadArraysIfNecessary() @@ -5000,7 +5054,7 @@ class MEDLoaderTest4(unittest.TestCase): arr=DataArrayDouble(7) ; arr.iota(2000) f1.setArray(arr) f1.setName("bbb") - f1.checkCoherency() + f1.checkConsistencyLight() f1.setTime(*zeTime) ff1.setFieldNoProfileSBT(f1) # @@ -5009,7 +5063,7 @@ class MEDLoaderTest4(unittest.TestCase): arr=DataArrayDouble(9) ; arr.iota(4000) f2.setArray(arr) f2.setName("ddd") - f2.checkCoherency() + f2.checkConsistencyLight() f2.setTime(*zeTime) ff2.setFieldNoProfileSBT(f2) # @@ -5018,18 +5072,18 @@ class MEDLoaderTest4(unittest.TestCase): f3.setGaussLocalizationOnType(NORM_TRI3,[0,0,1,0,0,1],[0.333333,0.333333],[0.5]) arr=DataArrayDouble(3) ; arr.iota(1000) f3.setArray(arr) - f3.checkCoherency() + f3.checkConsistencyLight() f3.setTime(*zeTime) f3.setName("aaa") ff3.setFieldNoProfileSBT(f3) # ff4=MEDFileField1TS() - m0d=m0.deepCpy() ; m0d.zipCoords() + m0d=m0.deepCopy() ; m0d.zipCoords() f4=MEDCouplingFieldDouble(ON_NODES) ; f4.setMesh(m0d) arr=DataArrayDouble(5) ; arr.iota(3000) f4.setArray(arr) f4.setName("ccc") - f4.checkCoherency() + f4.checkConsistencyLight() f4.setTime(*zeTime) pfl=DataArrayInt([0,1,2,3,4]) ; pfl.setName("PFL") ff4.setFieldProfile(f4,mm,0,pfl) @@ -5078,7 +5132,8 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,0,2,1,3,3,2,0,3,2,4,1]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -5105,7 +5160,8 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(a3.isEqual(DataArrayInt([3,0,2,1,3,3,2,0,3,2,4,1,2,5,0,2,0,3,2,3,2,2,2,6]))) self.assertTrue(a4 is None) self.assertTrue(a5 is None) - for i in xrange(1): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1): ffCell=allFMTSLeavesPerCommonSupport1[1][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -5126,9 +5182,9 @@ class MEDLoaderTest4(unittest.TestCase): arrZ=DataArrayDouble(5) ; arrZ.iota() ; arrZ*=1.6 ; arrZ-=8. ; arrZ.setInfoOnComponent(0,comps[2]) m=MEDCouplingCMesh() ; m.setCoords(arrX,arrY,arrZ) ; m.setName(meshName) mm=MEDFileCMesh() ; mm.setMesh(m) ; mm.setDescription(description) - mm.setAxType(AX_CYL) # the test is here ! + mm.setAxisType(AX_CYL) # the test is here ! f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(m) ; f.setName("Field") - arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkCoherency() + arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkConsistencyLight() ff=MEDFileField1TS() ; ff.setFieldNoProfileSBT(f) fmts=MEDFileFieldMultiTS() ; fmts.pushBackTimeStep(ff) # @@ -5174,7 +5230,8 @@ class MEDLoaderTest4(unittest.TestCase): ref_a.setInfoOnComponents(comps) self.assertTrue(a.isEqual(ref_a,1e-14)) self.assertEqual(b,[3,4,5]) - for i in xrange(1): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -5195,9 +5252,9 @@ class MEDLoaderTest4(unittest.TestCase): arrZ=DataArrayDouble(5) ; arrZ.iota() ; arrZ*=2*pi/(len(arrZ)-1) ; arrZ.setInfoOnComponent(0,comps[2]) m=MEDCouplingCMesh() ; m.setCoords(arrX,arrY,arrZ) ; m.setName(meshName) mm=MEDFileCMesh() ; mm.setMesh(m) ; mm.setDescription(description) - mm.setAxType(AX_SPHER) # the test is here ! + mm.setAxisType(AX_SPHER) # the test is here ! f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(m) ; f.setName("Field") - arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkCoherency() + arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkConsistencyLight() ff=MEDFileField1TS() ; ff.setFieldNoProfileSBT(f) fmts=MEDFileFieldMultiTS() ; fmts.pushBackTimeStep(ff) # @@ -5243,7 +5300,8 @@ class MEDLoaderTest4(unittest.TestCase): ref_a.setInfoOnComponents(comps) self.assertTrue(a.isEqual(ref_a,1e-14)) self.assertEqual(b,[3,4,5]) - for i in xrange(1): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -5264,9 +5322,9 @@ class MEDLoaderTest4(unittest.TestCase): arrZ=DataArrayDouble(5) ; arrZ.iota() ; arrZ*=2*pi/(len(arrZ)-1) ; arrZ.setInfoOnComponent(0,comps[2]) m=MEDCouplingCMesh() ; m.setCoords(arrX,arrY,arrZ) ; m.setName(meshName) ; m=m.buildUnstructured() mm=MEDFileUMesh() ; mm[0]=m ; mm.setDescription(description) # the test is here : UMesh ! - mm.setAxType(AX_SPHER) # the test is here ! + mm.setAxisType(AX_SPHER) # the test is here ! f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(m) ; f.setName("Field") - arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkCoherency() + arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkConsistencyLight() ff=MEDFileField1TS() ; ff.setFieldNoProfileSBT(f) fmts=MEDFileFieldMultiTS() ; fmts.pushBackTimeStep(ff) # @@ -5312,7 +5370,8 @@ class MEDLoaderTest4(unittest.TestCase): ref_a=DataArrayDouble([0.,0.,0.,0.,0.,0.8,0.,0.,1.6,0.,0.,0.,0.6928203230275509,0.,0.4,1.3856406460551018,0.,0.8,0.,0.,-0.,0.692820323027551,0.,-0.4,1.385640646055102,0.,-0.8,0.,0.,-0.,0.,0.,-0.8,0.,0.,-1.6,0.,0.,0.,0.,0.,0.8,0.,0.,1.6,0.,0.,0.,0.,0.6928203230275509,0.4,0.,1.3856406460551018,0.8,0.,0.,-0.,0.,0.692820323027551,-0.4,0.,1.385640646055102,-0.8,0.,0.,-0.,0.,0.,-0.8,0.,0.,-1.6,-0.,0.,0.,-0.,0.,0.8,-0.,0.,1.6,-0.,0.,0.,-0.6928203230275509,0.,0.4,-1.3856406460551018,0.,0.8,-0.,0.,-0.,-0.692820323027551,0.,-0.4,-1.385640646055102,0.,-0.8,-0.,0.,-0.,0.,0.,-0.8,0.,0.,-1.6,-0.,-0.,0.,-0.,-0.,0.8,-0.,-0.,1.6,-0.,-0.,0.,0.,-0.6928203230275509,0.4,0.,-1.3856406460551018,0.8,-0.,-0.,-0.,0.,-0.692820323027551,-0.4,0.,-1.385640646055102,-0.8,-0.,-0.,-0.,0.,0.,-0.8,0.,0.,-1.6,0.,-0.,0.,0.,-0.,0.8,0.,-0.,1.6,0.,-0.,0.,0.6928203230275509,0.,0.4,1.3856406460551018,0.,0.8,0.,-0.,-0.,0.692820323027551,0.,-0.4,1.385640646055102,0.,-0.8,0.,-0.,-0.,0.,0.,-0.8,0.,0.,-1.6],60,3) ref_a.setInfoOnComponents(comps) self.assertTrue(a0.isEqual(ref_a,1e-14))#<- Test is here - for i in xrange(1): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -5332,9 +5391,9 @@ class MEDLoaderTest4(unittest.TestCase): arrZ=DataArrayDouble(5) ; arrZ.iota() ; arrZ*=2*pi/(len(arrZ)-1) ; arrZ.setInfoOnComponent(0,comps[2]) m=MEDCouplingCMesh() ; m.setCoords(arrX,arrY,arrZ) ; m.setName(meshName) ; m=m.buildCurveLinear() mm=MEDFileCurveLinearMesh() ; mm.setMesh(m) ; mm.setDescription(description) # the test is here CLMesh! - mm.setAxType(AX_SPHER) # the test is here ! + mm.setAxisType(AX_SPHER) # the test is here ! f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(m) ; f.setName("Field") - arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkCoherency() + arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkConsistencyLight() ff=MEDFileField1TS() ; ff.setFieldNoProfileSBT(f) fmts=MEDFileFieldMultiTS() ; fmts.pushBackTimeStep(ff) # @@ -5377,7 +5436,8 @@ class MEDLoaderTest4(unittest.TestCase): ref_a.setInfoOnComponents(comps) self.assertTrue(a.isEqual(ref_a,1e-14))#<- Test is here self.assertEqual(b,[3,4,5]) - for i in xrange(1): + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1): ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) ffCell.loadArraysIfNecessary() @@ -5386,6 +5446,143 @@ class MEDLoaderTest4(unittest.TestCase): self.assertTrue(v.isEqual(DataArrayDouble([0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.1,2.2,2.3]),1e-14)) pass + def test41(self): + """This test focused on bug revealed with // load of multi nodes field with no profile. The error was the first node field (dataarray partdef) change the partdef for the others ! """ + fname="ForMEDReader41.med" + meshName="mesh" + nx=5 + arr=DataArrayDouble(nx) ; arr.iota() + m=MEDCouplingCMesh() ; m.setCoords(arr,arr) ; m=m.buildUnstructured() ; m.setName(meshName) + renum=DataArrayInt.Aggregate([DataArrayInt.Range(0,m.getNumberOfCells(),2),DataArrayInt.Range(1,m.getNumberOfCells(),2)]) + m=m[renum] # important think to renum if not we see nothing. The goal if to have dynamic_cast!=NULL + mm=MEDFileUMesh() ; mm[0]=m + mm.write(fname,2) + f0=MEDCouplingFieldDouble(ON_NODES) ; f0.setMesh(m) ; f0.setName("aaa") + arr0=DataArrayDouble(nx*nx) ; arr0.iota() ; f0.setArray(arr0) + ff0=MEDFileField1TS() ; ff0.setFieldNoProfileSBT(f0) + f1=MEDCouplingFieldDouble(ON_NODES) ; f1.setMesh(m) ; f1.setName("bbb") + arr1=DataArrayDouble(nx*nx) ; arr1.iota() ; arr1+=100 ; f1.setArray(arr1) + ff1=MEDFileField1TS() ; ff1.setFieldNoProfileSBT(f1) + ff0.write(fname,0) ; ff1.write(fname,0) + # + a=8 ; b=16 + ms=MEDFileMeshes() + mm=MEDFileUMesh.LoadPartOf(fname,meshName,[NORM_QUAD4],[a,b,1],-1,-1) + ms.pushMesh(mm) + ms[0].zipCoords() + ms.cartesianizeMe() + fields=MEDFileFields.LoadPartOf(fname,False,ms); + fields.removeFieldsWithoutAnyTimeStep() + fields_per_mesh=[fields.partOfThisLyingOnSpecifiedMeshName(meshName) for meshName in ms.getMeshesNames()] + allFMTSLeavesToDisplay=[] + for fields in fields_per_mesh: + allFMTSLeavesToDisplay2=[] + for fmts in fields: + tmp=fmts.splitDiscretizations() + for itmp in tmp: + if itmp.presenceOfMultiDiscPerGeoType(): + tmp2=itmp.splitMultiDiscrPerGeoTypes() + for iii,itmp2 in enumerate(tmp2): + name="%s_%i"%(itmp2.getName(),iii) + itmp2.setName(name) + allFMTSLeavesToDisplay2.append(itmp2) + pass + pass + else: + allFMTSLeavesToDisplay2.append(itmp) + pass + pass + allFMTSLeavesToDisplay.append(allFMTSLeavesToDisplay2) + pass + # GO for reading in MEDReader, by not loading all. Mesh is fully loaded but not fields values + allFMTSLeavesPerTimeSeries=MEDFileAnyTypeFieldMultiTS.SplitIntoCommonTimeSeries(sum(allFMTSLeavesToDisplay,[])) + allFMTSLeavesPerCommonSupport=MEDFileAnyTypeFieldMultiTS.SplitPerCommonSupport(allFMTSLeavesPerTimeSeries[0],ms[ms.getMeshesNames()[0]]) + mst=MEDFileMeshStruct.New(ms[0]) + fcscp=allFMTSLeavesPerCommonSupport[0][1] + mml=fcscp.buildFromScratchDataSetSupport(0,fields) + mml2=mml.prepare() + # + f2=allFMTSLeavesPerCommonSupport[0][0][0][0] + fsst=MEDFileField1TSStructItem.BuildItemFrom(f2,mst) + f2.loadArraysIfNecessary() + v0=mml.buildDataArray(fsst,fields,f2.getUndergroundDataArray()) + assert(v0.isEqual(DataArrayDouble([1,2,3,4,6,7,8,9,11,12,13,14,16,17,18,19,21,22,23,24]),1e-12)) + # + f2=allFMTSLeavesPerCommonSupport[0][0][1][0] + fsst=MEDFileField1TSStructItem.BuildItemFrom(f2,mst) + f2.loadArraysIfNecessary() + v1=mml.buildDataArray(fsst,fields,f2.getUndergroundDataArray()) + assert(v1.isEqual(DataArrayDouble([101,102,103,104,106,107,108,109,111,112,113,114,116,117,118,119,121,122,123,124]),1e-12)) + pass + + def test42(self): + """ EDF14869 - SEG4 """ + fname="ForMEDReader42.med" + meshName="mesh" + # + a0exp=DataArrayDouble([0.,1.,0.3,0.7]) + m=MEDCouplingUMesh("mesh",1) + m.setCoords(a0exp) + m.allocateCells() + m.insertNextCell(NORM_SEG4,[0,1,2,3]) + mm=MEDFileUMesh() ; mm[0]=m + # + f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(m) ; f.setName("Field") + arr=DataArrayDouble(m.getNumberOfCells()) ; arr.iota() ; arr*=0.1 ; f.setArray(arr) ; f.checkConsistencyLight() + ff=MEDFileField1TS() ; ff.setFieldNoProfileSBT(f) + fmts=MEDFileFieldMultiTS() ; fmts.pushBackTimeStep(ff) + # + ms=MEDFileMeshes() ; ms.pushMesh(mm) + fields=MEDFileFields() ; fields.pushField(fmts) + ms.write(fname,2) ; fields.write(fname,0) + # + ms=MEDFileMeshes(fname) + fields=MEDFileFields(fname,False) + fields.removeFieldsWithoutAnyTimeStep() + fields_per_mesh=[fields.partOfThisLyingOnSpecifiedMeshName(meshName) for meshName in ms.getMeshesNames()] + allFMTSLeavesToDisplay=[] + for fields in fields_per_mesh: + allFMTSLeavesToDisplay2=[] + for fmts in fields: + tmp=fmts.splitDiscretizations() + for itmp in tmp: + self.assertTrue(not itmp.presenceOfMultiDiscPerGeoType()) + pass + allFMTSLeavesToDisplay2+=tmp + pass + allFMTSLeavesToDisplay.append(allFMTSLeavesToDisplay2) + pass + # + self.assertEqual(len(allFMTSLeavesToDisplay),1) + self.assertEqual(len(allFMTSLeavesToDisplay[0]),1) + allFMTSLeavesPerTimeSeries=MEDFileAnyTypeFieldMultiTS.SplitIntoCommonTimeSeries(sum(allFMTSLeavesToDisplay,[])) + self.assertEqual(len(allFMTSLeavesPerTimeSeries),1) + allFMTSLeavesPerCommonSupport1=MEDFileAnyTypeFieldMultiTS.SplitPerCommonSupport(allFMTSLeavesToDisplay[0],ms[ms.getMeshesNames()[0]]) + self.assertEqual(len(allFMTSLeavesPerCommonSupport1),1) + # + mst=MEDFileMeshStruct.New(ms[0]) + fcscp=allFMTSLeavesPerCommonSupport1[0][1] + mml=fcscp.buildFromScratchDataSetSupport(0,fields) + mml2=mml.prepare() + self.assertTrue(isinstance(mml2,MEDUMeshMultiLev)) + ncc,a0,a1,a2,a3,a4,a5=mml2.buildVTUArrays() + self.assertTrue(not ncc) + self.assertTrue(a0.isEqual(a0exp.changeNbOfComponents(3,0.),1e-12)) + self.assertTrue(a1.isEqual(DataArrayByte([35])))# VTK_CUBIC_LINE + self.assertTrue(a2.isEqual(DataArrayInt([0]))) + self.assertTrue(a3.isEqual(DataArrayInt([4,0,1,2,3]))) + self.assertTrue(a4 is None) + self.assertTrue(a5 is None) + self.assertTrue(mml2.retrieveGlobalNodeIdsIfAny() is None) + for i in range(1): + ffCell=allFMTSLeavesPerCommonSupport1[0][0][0][i] + fsst=MEDFileField1TSStructItem.BuildItemFrom(ffCell,mst) + ffCell.loadArraysIfNecessary() + v=mml2.buildDataArray(fsst,fields,ffCell.getUndergroundDataArray()) + self.assertEqual(v.getHiddenCppPointer(),ffCell.getUndergroundDataArray().getHiddenCppPointer()) + self.assertTrue(v.isEqual(DataArrayDouble([0.0]),1e-14)) + pass + pass if __name__ == "__main__":