def testMesh3DRW(self):
mesh=MEDLoaderDataForTest.build3DMesh_1();
mesh.checkConsistencyLight();
WriteUMeshDep("Pyfile5.med",mesh,False);
mesh_rw=ReadUMeshFromFile("Pyfile5.med",mesh.getName(),0);
self.assertTrue(mesh.isEqual(mesh_rw,1e-12));
pass
def testMesh3DRW(self):
mesh = MEDLoaderDataForTest.build3DMesh_1()
mesh.checkCoherency()
MEDLoader.WriteUMeshDep("Pyfile5.med", mesh, False)
mesh_rw = MEDLoader.ReadUMeshFromFile("Pyfile5.med", mesh.getName(), 0)
self.assertTrue(mesh.isEqual(mesh_rw, 1e-12))
pass
def testFieldProfilRW1(self):
fileName="Pyfile12.med";
mesh1=MEDLoaderDataForTest.build3DMesh_1();
da,b,newNbOfNodes=mesh1.mergeNodes(1e-12);
MEDLoader.WriteUMesh(fileName,mesh1,True);
part1=[1,2,4,13,15]
mesh2=mesh1.buildPartOfMySelf(part1,True);
mesh2.setName(mesh1.getName());#<- important for the test
#
nbOfCells=mesh2.getNumberOfCells();
self.assertEqual(5,nbOfCells);
f1=MEDLoader.MEDCouplingFieldDouble.New(MEDLoader.ON_CELLS,MEDLoader.ONE_TIME);
f1.setName("VectorFieldOnCells");
f1.setMesh(mesh2);
array=MEDLoader.DataArrayDouble.New();
array.alloc(nbOfCells,2);
f1.setArray(array);
arr1=[71.,171.,10.,110.,20.,120.,30.,130.,40.,140.]
array.setValues(arr1,nbOfCells,2);
f1.setTime(3.14,2,7);
f1.checkConsistencyLight();
#
MEDLoader.WriteField(fileName,f1,False);#<- False important for the test
#
f2=MEDLoader.ReadFieldCell(fileName,f1.getMesh().getName(),0,f1.getName(),2,7);
tt=MEDLoader.GetTypesOfField(fileName,f1.getMesh().getName(),f1.getName());
self.assertEqual(tt,[MEDLoader.ON_CELLS]);
f2.checkConsistencyLight();
self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
#
pass
def testMesh3DRW(self):
mesh = MEDLoaderDataForTest.build3DMesh_1()
mesh.checkCoherency()
MEDLoader.WriteUMeshDep("Pyfile5.med", mesh, False)
mesh_rw = MEDLoader.ReadUMeshFromFile("Pyfile5.med", mesh.getName(), 0)
self.assertTrue(mesh.isEqual(mesh_rw, 1e-12))
pass
def testMultiMeshRW1(self):
fileName = "Pyfile10.med"
mesh1 = MEDLoaderDataForTest.build3DMesh_1()
part1 = [1, 2, 4, 13, 15]
mesh2 = mesh1.buildPartOfMySelf(part1, True)
mesh2.setName("mesh2")
part2 = [3, 4, 13, 14]
mesh3 = mesh1.buildPartOfMySelf(part2, True)
mesh3.setName("mesh3")
mesh4 = MEDCouplingUMesh.New()
mesh4.setName("mesh4")
mesh4.setMeshDimension(3)
mesh4.allocateCells(1)
conn = [0, 11, 1, 3]
mesh4.insertNextCell(NORM_TETRA4, 4, conn[0:4])
mesh4.finishInsertingCells()
mesh4.setCoords(mesh1.getCoords())
meshes = [mesh1, mesh2, mesh3, mesh4]
mnane = "3DToto"
MEDLoader.WriteUMeshesPartitionDep(fileName, mnane, meshes, False)
#
mesh5 = MEDLoader.ReadUMeshFromFile(fileName, mnane)
mesh1.setName(mnane)
part3 = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17]
mesh6 = mesh5.buildPartOfMySelf(part3, True)
mesh6.setName(mnane)
self.assertTrue(mesh6.isEqual(mesh1, 1e-12))
grps = MEDLoader.GetMeshGroupsNames(fileName, mnane)
self.assertEqual(4, len(grps))
grps.index("mesh2")
grps.index("mesh3")
grps.index("mesh4")
grps.index("3DMesh_1")
#
vec = ["mesh2"]
mesh2_2 = MEDLoader.ReadUMeshFromGroups(fileName, mnane, 0, vec)
self.assertTrue(mesh2_2.isEqual(mesh2, 1e-12))
vec = ["mesh3"]
mesh3_2 = MEDLoader.ReadUMeshFromGroups(fileName, mnane, 0, vec)
self.assertTrue(mesh3_2.isEqual(mesh3, 1e-12))
vec = ["mesh4"]
mesh4_2 = MEDLoader.ReadUMeshFromGroups(fileName, mnane, 0, vec)
self.assertTrue(mesh4_2.isEqual(mesh4, 1e-12))
vec = ["3DMesh_1"]
mesh1_2 = MEDLoader.ReadUMeshFromGroups(fileName, mnane, 0, vec)
mesh1.setName("3DMesh_1")
self.assertTrue(mesh1_2.isEqual(mesh1, 1e-12))
#
vec = ["Family_-5", "Family_-3"]
mesh2_2 = MEDLoader.ReadUMeshFromFamilies(fileName, mnane, 0, vec)
mesh2_2.setName("mesh2")
self.assertTrue(mesh2_2.isEqual(mesh2, 1e-12))
pass
def testMultiMeshRW1(self):
fileName="Pyfile10.med";
mesh1=MEDLoaderDataForTest.build3DMesh_1();
part1=[1,2,4,13,15]
mesh2=mesh1.buildPartOfMySelf(part1,True);
mesh2.setName("mesh2");
part2=[3,4,13,14]
mesh3=mesh1.buildPartOfMySelf(part2,True);
mesh3.setName("mesh3");
mesh4=MEDCouplingUMesh.New();
mesh4.setName("mesh4");
mesh4.setMeshDimension(3);
mesh4.allocateCells(1);
conn=[0,11,1,3]
mesh4.insertNextCell(NORM_TETRA4,4,conn[0:4])
mesh4.finishInsertingCells();
mesh4.setCoords(mesh1.getCoords());
meshes=[mesh1,mesh2,mesh3,mesh4]
mnane="3DToto";
WriteUMeshesPartitionDep(fileName,mnane,meshes,False);
#
mesh5=ReadUMeshFromFile(fileName,mnane);
mesh1.setName(mnane);
part3=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17]
mesh6=mesh5.buildPartOfMySelf(part3,True);
mesh6.setName(mnane);
self.assertTrue(mesh6.isEqual(mesh1,1e-12));
grps=GetMeshGroupsNames(fileName,mnane);
self.assertEqual(4,len(grps));
grps.index("mesh2");
grps.index("mesh3");
grps.index("mesh4");
grps.index("3DMesh_1");
#
vec=["mesh2"];
mesh2_2=ReadUMeshFromGroups(fileName,mnane,0,vec);
self.assertTrue(mesh2_2.isEqual(mesh2,1e-12));
vec=["mesh3"];
mesh3_2=ReadUMeshFromGroups(fileName,mnane,0,vec);
self.assertTrue(mesh3_2.isEqual(mesh3,1e-12));
vec=["mesh4"];
mesh4_2=ReadUMeshFromGroups(fileName,mnane,0,vec);
self.assertTrue(mesh4_2.isEqual(mesh4,1e-12));
vec=["3DMesh_1"];
mesh1_2=ReadUMeshFromGroups(fileName,mnane,0,vec);
mesh1.setName("3DMesh_1");
self.assertTrue(mesh1_2.isEqual(mesh1,1e-12));
#
vec=["Family_-5","Family_-3"];
mesh2_2=ReadUMeshFromFamilies(fileName,mnane,0,vec);
mesh2_2.setName("mesh2");
self.assertTrue(mesh2_2.isEqual(mesh2,1e-12));
pass
def testFieldProfilRW1(self):
fileName = "Pyfile12.med"
mesh1 = MEDLoaderDataForTest.build3DMesh_1()
da, b, newNbOfNodes = mesh1.mergeNodes(1e-12)
MEDLoader.MEDLoader.WriteUMesh(fileName, mesh1, True)
part1 = [1, 2, 4, 13, 15]
mesh2 = mesh1.buildPartOfMySelf(part1, True)
mesh2.setName(mesh1.getName())
#<- important for the test
#
nbOfCells = mesh2.getNumberOfCells()
self.assertEqual(5, nbOfCells)
f1 = MEDLoader.MEDCouplingFieldDouble.New(MEDLoader.ON_CELLS,
MEDLoader.ONE_TIME)
f1.setName("VectorFieldOnCells")
f1.setMesh(mesh2)
array = MEDLoader.DataArrayDouble.New()
array.alloc(nbOfCells, 2)
f1.setArray(array)
arr1 = [71., 171., 10., 110., 20., 120., 30., 130., 40., 140.]
array.setValues(arr1, nbOfCells, 2)
f1.setTime(3.14, 2, 7)
f1.checkCoherency()
#
MEDLoader.MEDLoader.WriteField(fileName, f1, False)
#<- False important for the test
#
f2 = MEDLoader.MEDLoader.ReadFieldCell(fileName,
f1.getMesh().getName(), 0,
f1.getName(), 2, 7)
tt = MEDLoader.MEDLoader.GetTypesOfField(fileName,
f1.getMesh().getName(),
f1.getName())
self.assertEqual(tt, [MEDLoader.ON_CELLS])
f2.checkCoherency()
self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
#
pass
#! [UG_ReadMeshFromFile_0]
from MEDLoader import WriteMesh
WriteMesh("file2.med", m, True)
#! [UG_ReadMeshFromFile_0]
#! [UG_ReadMeshFromFile_1]
from MEDLoader import ReadMeshFromFile
m = ReadMeshFromFile("file2.med")
#! [UG_ReadMeshFromFile_1]
#! [UG_ReadMeshFromFile_2]
m = ReadMeshFromFile("file2.med", "mesh2")
assert (m.getName() == "mesh2")
#! [UG_ReadMeshFromFile_2]
mesh3D = MEDLoaderDataForTest.build3DMesh_1()
mesh2D = mesh3D.computeSkin()
mesh1D = mesh2D.computeSkin()
#! [UG_ReadMeshFromFile_4]
from MEDLoader import MEDFileUMesh
mm = MEDFileUMesh.New()
mm.setMeshAtLevel(0, mesh3D)
mm.setMeshAtLevel(-1, mesh2D)
#! [UG_ReadMeshFromFile_4]
otherCoordArray = mesh3D.getCoords()
#! [UG_ReadMeshFromFile_5]
mm.setCoords(otherCoordArray)
#! [UG_ReadMeshFromFile_5]
#! [UG_ReadMeshFromFile_6]
groupNodes = DataArrayInt([1, 3, 4, 5])
