Exemplo n.º 1
0
 def testFieldRW2(self):
     fileName = "Pyfile8.med"
     VAL1 = 12345.67890314
     VAL2 = -1111111111111.
     f1 = MEDLoaderDataForTest.buildVecFieldOnCells_1()
     MEDLoader.WriteFieldDep(fileName, f1, False)
     f1.setTime(10., 8, 9)
     f1.getArray().setIJ(0, 0, VAL1)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     f1.setTime(10.14, 18, 19)
     f1.getArray().setIJ(0, 0, VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     #retrieving time steps...
     f2 = MEDLoader.ReadFieldCell(fileName,
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  8, 9)
     f1.setTime(10., 8, 9)
     f1.getArray().setIJ(0, 0, VAL1)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldCell(fileName,
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  0, 1)
     f3 = MEDLoaderDataForTest.buildVecFieldOnCells_1()
     self.assertTrue(f3.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldCell(fileName,
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  18, 19)
     f1.setTime(10.14, 18, 19)
     f1.getArray().setIJ(0, 0, VAL2)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     #ON NODES
     f1 = MEDLoaderDataForTest.buildVecFieldOnNodes_1()
     fileName2 = "Pyfile9.med"
     MEDLoader.WriteFieldDep(fileName2, f1, False)
     f1.setTime(110., 108, 109)
     tmp = f1.getArray().getPointer()
     f1.getArray().setIJ(0, 3, VAL1)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName2, f1)
     f1.setTime(210., 208, 209)
     f1.getArray().setIJ(0, 3, VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName2, f1)
     f2 = MEDLoader.ReadFieldNode(fileName2,
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  108, 109)
     f1.setTime(110., 108, 109)
     f1.getArray().setIJ(0, 3, VAL1)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldNode(fileName2,
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  2, 3)
     f3 = MEDLoaderDataForTest.buildVecFieldOnNodes_1()
     self.assertTrue(f3.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldNode(fileName2,
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  208, 209)
     f1.setTime(210., 208, 209)
     f1.getArray().setIJ(0, 3, VAL2)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     pass
Exemplo n.º 2
0
 def testFieldRW2(self):
     fileName="Pyfile8.med";
     VAL1=12345.67890314;
     VAL2=-1111111111111.;
     f1=MEDLoaderDataForTest.buildVecFieldOnCells_1();
     MEDLoader.WriteField(fileName,f1,True);
     f1.setTime(10.,8,9);
     f1.getArray().setIJ(0,0,VAL1);
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     f1.setTime(10.14,18,19);
     f1.getArray().setIJ(0,0,VAL2);
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     #retrieving time steps...
     f2=MEDLoader.ReadFieldCell(fileName,f1.getMesh().getName(),0,f1.getName(),8,9);
     f1.setTime(10.,8,9);
     f1.getArray().setIJ(0,0,VAL1);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     f2=MEDLoader.ReadFieldCell(fileName,f1.getMesh().getName(),0,f1.getName(),0,1);
     f3=MEDLoaderDataForTest.buildVecFieldOnCells_1();
     self.assertTrue(f3.isEqual(f2,1e-12,1e-12));
     f2=MEDLoader.ReadFieldCell(fileName,f1.getMesh().getName(),0,f1.getName(),18,19);
     f1.setTime(10.14,18,19);
     f1.getArray().setIJ(0,0,VAL2);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     #test of throw on invalid (dt,it)
     self.assertRaises(Exception,MEDLoader.ReadFieldCell,fileName,f1.getMesh().getName(),0,f1.getName(),28,19);
     #ON NODES
     f1=MEDLoaderDataForTest.buildVecFieldOnNodes_1();
     fileName2="Pyfile9.med";
     MEDLoader.WriteField(fileName2,f1,True);
     f1.setTime(110.,108,109);
     tmp=f1.getArray().getPointer();
     f1.getArray().setIJ(0,3,VAL1);
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName2,f1);
     f1.setTime(210.,208,209);
     f1.getArray().setIJ(0,3,VAL2);
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName2,f1);
     f2=MEDLoader.ReadFieldNode(fileName2,f1.getMesh().getName(),0,f1.getName(),108,109);
     f1.setTime(110.,108,109);
     f1.getArray().setIJ(0,3,VAL1);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     f2=MEDLoader.ReadFieldNode(fileName2,f1.getMesh().getName(),0,f1.getName(),2,3);
     f3=MEDLoaderDataForTest.buildVecFieldOnNodes_1();
     self.assertTrue(f3.isEqual(f2,1e-12,1e-12));
     f2=MEDLoader.ReadFieldNode(fileName2,f1.getMesh().getName(),0,f1.getName(),208,209);
     f1.setTime(210.,208,209);
     f1.getArray().setIJ(0,3,VAL2);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     pass
Exemplo n.º 3
0
 def testFieldRW2(self):
     fileName = "Pyfile8.med"
     VAL1 = 12345.67890314
     VAL2 = -1111111111111.0
     f1 = MEDLoaderDataForTest.buildVecFieldOnCells_1()
     MEDLoader.WriteFieldDep(fileName, f1, False)
     f1.setTime(10.0, 8, 9)
     f1.getArray().setIJ(0, 0, VAL1)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     f1.setTime(10.14, 18, 19)
     f1.getArray().setIJ(0, 0, VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     # retrieving time steps...
     f2 = MEDLoader.ReadFieldCell(fileName, f1.getMesh().getName(), 0, f1.getName(), 8, 9)
     f1.setTime(10.0, 8, 9)
     f1.getArray().setIJ(0, 0, VAL1)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldCell(fileName, f1.getMesh().getName(), 0, f1.getName(), 0, 1)
     f3 = MEDLoaderDataForTest.buildVecFieldOnCells_1()
     self.assertTrue(f3.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldCell(fileName, f1.getMesh().getName(), 0, f1.getName(), 18, 19)
     f1.setTime(10.14, 18, 19)
     f1.getArray().setIJ(0, 0, VAL2)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     # ON NODES
     f1 = MEDLoaderDataForTest.buildVecFieldOnNodes_1()
     fileName2 = "Pyfile9.med"
     MEDLoader.WriteFieldDep(fileName2, f1, False)
     f1.setTime(110.0, 108, 109)
     tmp = f1.getArray().getPointer()
     f1.getArray().setIJ(0, 3, VAL1)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName2, f1)
     f1.setTime(210.0, 208, 209)
     f1.getArray().setIJ(0, 3, VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName2, f1)
     f2 = MEDLoader.ReadFieldNode(fileName2, f1.getMesh().getName(), 0, f1.getName(), 108, 109)
     f1.setTime(110.0, 108, 109)
     f1.getArray().setIJ(0, 3, VAL1)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldNode(fileName2, f1.getMesh().getName(), 0, f1.getName(), 2, 3)
     f3 = MEDLoaderDataForTest.buildVecFieldOnNodes_1()
     self.assertTrue(f3.isEqual(f2, 1e-12, 1e-12))
     f2 = MEDLoader.ReadFieldNode(fileName2, f1.getMesh().getName(), 0, f1.getName(), 208, 209)
     f1.setTime(210.0, 208, 209)
     f1.getArray().setIJ(0, 3, VAL2)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     pass
Exemplo n.º 4
0
 def testFieldRW1(self):
     f1=MEDLoaderDataForTest.buildVecFieldOnCells_1();
     WriteFieldDep("Pyfile6.med",f1,False);
     f2=ReadFieldCell("Pyfile6.med",f1.getMesh().getName(),0,f1.getName(),0,1);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     #
     f1=MEDLoaderDataForTest.buildVecFieldOnNodes_1();
     WriteFieldDep("Pyfile7.med",f1,False);
     f2=ReadFieldNode("Pyfile7.med",f1.getMesh().getName(),0,f1.getName(),2,3);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     pass
Exemplo n.º 5
0
 def testFieldRW1(self):
     f1=MEDLoaderDataForTest.buildVecFieldOnCells_1();
     MEDLoader.WriteField("Pyfile6.med",f1,True);
     f2=MEDLoader.ReadFieldCell("Pyfile6.med",f1.getMesh().getName(),0,f1.getName(),0,1);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     #
     f1=MEDLoaderDataForTest.buildVecFieldOnNodes_1();
     MEDLoader.WriteField("Pyfile7.med",f1,True);
     f2=MEDLoader.ReadFieldNode("Pyfile7.med",f1.getMesh().getName(),0,f1.getName(),2,3);
     self.assertTrue(f1.isEqual(f2,1e-12,1e-12));
     self.assertRaises(Exception,MEDLoader.ReadFieldCell,"Pyfile7.med",f1.getMesh().getName(),0,f1.getName(),2,3);
     pass
Exemplo n.º 6
0
 def testFieldRW1(self):
     f1 = MEDLoaderDataForTest.buildVecFieldOnCells_1()
     MEDLoader.WriteField("Pyfile6.med", f1, True)
     f2 = MEDLoader.ReadFieldCell("Pyfile6.med",
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  0, 1)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     #
     f1 = MEDLoaderDataForTest.buildVecFieldOnNodes_1()
     MEDLoader.WriteField("Pyfile7.med", f1, True)
     f2 = MEDLoader.ReadFieldNode("Pyfile7.med",
                                  f1.getMesh().getName(), 0, f1.getName(),
                                  2, 3)
     self.assertTrue(f1.isEqual(f2, 1e-12, 1e-12))
     self.assertRaises(Exception, MEDLoader.ReadFieldCell, "Pyfile7.med",
                       f1.getMesh().getName(), 0, f1.getName(), 2, 3)
     pass
Exemplo n.º 7
0
 def testFieldRW3(self):
     fileName = "Pyfile11.med"
     VAL1 = 12345.67890314
     VAL2 = -1111111111111.
     name1 = "AField"
     name3 = "AMesh1"
     f1 = MEDLoaderDataForTest.buildVecFieldOnCells_1()
     f1.getMesh().setName(name3)
     f1.setName(name1)
     f1.setTime(10., 8, 9)
     tmp = f1.getArray().getPointer()
     f1.getArray().setIJ(0, 0, VAL1)
     MEDLoader.WriteFieldDep(fileName, f1, False)
     f1.setTime(10.14, 18, 19)
     f1.getArray().setIJ(0, 0, VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     f1.getMesh().setName(name3)
     f1.setTime(10.55, 28, 29)
     f1.getArray().setIJ(0, 0, 3 * VAL1)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     f1.setTime(10.66, 38, 39)
     f1.getArray().setIJ(0, 0, 3 * VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     f1.setTime(10.77, 48, 49)
     f1.getArray().setIJ(0, 0, 4 * VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     #ON NODES
     f1 = MEDLoaderDataForTest.buildVecFieldOnNodes_1()
     f1.setName(name1)
     f1.getMesh().setName(name3)
     f1.setTime(110., 8, 9)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     f1.setTime(110., 108, 109)
     tmp = f1.getArray().getPointer()
     f1.getArray().setIJ(0, 3, VAL1)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     f1.setTime(210., 208, 209)
     f1.getArray().setIJ(0, 3, VAL2)
     MEDLoader.WriteFieldUsingAlreadyWrittenMesh(fileName, f1)
     #
     it1 = MEDLoader.GetCellFieldIterations(fileName, name3, name1)
     self.assertEqual(5, len(it1))
     self.assertEqual(8, it1[0][0])
     self.assertEqual(9, it1[0][1])
     self.assertEqual(18, it1[1][0])
     self.assertEqual(19, it1[1][1])
     self.assertEqual(28, it1[2][0])
     self.assertEqual(29, it1[2][1])
     self.assertEqual(38, it1[3][0])
     self.assertEqual(39, it1[3][1])
     self.assertEqual(48, it1[4][0])
     self.assertEqual(49, it1[4][1])
     it3 = MEDLoader.GetNodeFieldIterations(fileName, name3, name1)
     self.assertEqual(3, len(it3))
     self.assertEqual(8, it3[0][0])
     self.assertEqual(9, it3[0][1])
     self.assertEqual(108, it3[1][0])
     self.assertEqual(109, it3[1][1])
     self.assertEqual(208, it3[2][0])
     self.assertEqual(209, it3[2][1])
     #
     #
     f1 = MEDLoader.ReadFieldCell(fileName, name3, 0, name1, 8, 9)
     self.assertAlmostEqual(VAL1,
                            f1.getArray().getIJ(0, 0), 13)
     f1 = MEDLoader.ReadFieldCell(fileName, name3, 0, name1, 18, 19)
     self.assertAlmostEqual(VAL2,
                            f1.getArray().getIJ(0, 0), 13)
     f1 = MEDLoader.ReadFieldCell(fileName, name3, 0, name1, 28, 29)
     self.assertAlmostEqual(3 * VAL1,
                            f1.getArray().getIJ(0, 0), 13)
     f1 = MEDLoader.ReadFieldCell(fileName, name3, 0, name1, 38, 39)
     self.assertAlmostEqual(3 * VAL2,
                            f1.getArray().getIJ(0, 0), 13)
     f1 = MEDLoader.ReadFieldCell(fileName, name3, 0, name1, 48, 49)
     self.assertAlmostEqual(4 * VAL2,
                            f1.getArray().getIJ(0, 0), 13)
     #
     f1 = MEDLoader.ReadFieldNode(fileName, name3, 0, name1, 8, 9)
     self.assertAlmostEqual(71.,
                            f1.getArray().getIJ(0, 3), 13)
     f1 = MEDLoader.ReadFieldNode(fileName, name3, 0, name1, 108, 109)
     self.assertAlmostEqual(VAL1,
                            f1.getArray().getIJ(0, 3), 13)
     f1 = MEDLoader.ReadFieldNode(fileName, name3, 0, name1, 208, 209)
     self.assertAlmostEqual(VAL2,
                            f1.getArray().getIJ(0, 3), 13)
     pass
Exemplo n.º 8
0
 def testFieldRW3(self):
     fileName="Pyfile11.med";
     VAL1=12345.67890314;
     VAL2=-1111111111111.;
     name1="AField";
     name3="AMesh1";
     f1=MEDLoaderDataForTest.buildVecFieldOnCells_1();
     f1.getMesh().setName(name3);
     f1.setName(name1);
     f1.setTime(10.,8,9);
     tmp=f1.getArray().getPointer();
     f1.getArray().setIJ(0,0,VAL1);
     WriteFieldDep(fileName,f1,False);
     f1.setTime(10.14,18,19);
     f1.getArray().setIJ(0,0,VAL2);
     WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     f1.getMesh().setName(name3);
     f1.setTime(10.55,28,29);
     f1.getArray().setIJ(0,0,3*VAL1);
     WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     f1.setTime(10.66,38,39);
     f1.getArray().setIJ(0,0,3*VAL2);
     WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     f1.setTime(10.77,48,49);
     f1.getArray().setIJ(0,0,4*VAL2);
     WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     #ON NODES
     f1=MEDLoaderDataForTest.buildVecFieldOnNodes_1();
     f1.setName(name1);
     f1.getMesh().setName(name3);
     f1.setTime(110.,8,9);
     WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     f1.setTime(110.,108,109);
     tmp=f1.getArray().getPointer();
     f1.getArray().setIJ(0,3,VAL1);
     WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     f1.setTime(210.,208,209);
     f1.getArray().setIJ(0,3,VAL2);
     WriteFieldUsingAlreadyWrittenMesh(fileName,f1);
     #
     it1=GetCellFieldIterations(fileName,name3,name1);
     self.assertEqual(5,len(it1));
     self.assertEqual(8,it1[0][0]); self.assertEqual(9,it1[0][1]);
     self.assertEqual(18,it1[1][0]); self.assertEqual(19,it1[1][1]);
     self.assertEqual(28,it1[2][0]); self.assertEqual(29,it1[2][1]);
     self.assertEqual(38,it1[3][0]); self.assertEqual(39,it1[3][1]);
     self.assertEqual(48,it1[4][0]); self.assertEqual(49,it1[4][1]);
     it3=GetNodeFieldIterations(fileName,name3,name1);
     self.assertEqual(3,len(it3));
     self.assertEqual(8,it3[0][0]); self.assertEqual(9,it3[0][1]);
     self.assertEqual(108,it3[1][0]); self.assertEqual(109,it3[1][1]);
     self.assertEqual(208,it3[2][0]); self.assertEqual(209,it3[2][1]);
     #
     #
     f1=ReadFieldCell(fileName,name3,0,name1,8,9);
     self.assertAlmostEqual(VAL1,f1.getArray().getIJ(0,0),13);
     f1=ReadFieldCell(fileName,name3,0,name1,18,19);
     self.assertAlmostEqual(VAL2,f1.getArray().getIJ(0,0),13);
     f1=ReadFieldCell(fileName,name3,0,name1,28,29);
     self.assertAlmostEqual(3*VAL1,f1.getArray().getIJ(0,0),13);
     f1=ReadFieldCell(fileName,name3,0,name1,38,39);
     self.assertAlmostEqual(3*VAL2,f1.getArray().getIJ(0,0),13);
     f1=ReadFieldCell(fileName,name3,0,name1,48,49);
     self.assertAlmostEqual(4*VAL2,f1.getArray().getIJ(0,0),13);
     #
     f1=ReadFieldNode(fileName,name3,0,name1,8,9);
     self.assertAlmostEqual(71.,f1.getArray().getIJ(0,3),13);
     f1=ReadFieldNode(fileName,name3,0,name1,108,109);
     self.assertAlmostEqual(VAL1,f1.getArray().getIJ(0,3),13);
     f1=ReadFieldNode(fileName,name3,0,name1,208,209);
     self.assertAlmostEqual(VAL2,f1.getArray().getIJ(0,3),13);
     pass