def testExampleReadFieldOnAllEntity1(self):
from MEDLoaderDataForTest import MEDLoaderDataForTest
#! [PySnippetReadFieldOnAllEntity1_1]
fname = "PyExamples1.med"
meshName = "mesh"
fieldName = "FieldOnAll"
iteration = 3
order = 4
#! [PySnippetReadFieldOnAllEntity1_1]
#! [PySnippetWriteFieldOnAllEntity1_2]
m = MEDLoaderDataForTest.build2DMesh_3()
m = m[:10]
m.setName(meshName)
f = m.getMeasureField(False)
f = f.buildNewTimeReprFromThis(ONE_TIME, False)
f.setTime(5.5, iteration, order)
f.setName(fieldName)
# MEDCoupling finished, MEDLoader advanced API specific part starting from here
mm = MEDFileUMesh.New()
mm.setMeshAtLevel(0, m)
ff = MEDFileField1TS.New()
ff.setFieldNoProfileSBT(f)
mm.write(fname, 2)
ff.write(fname, 0)
#! [PySnippetWriteFieldOnAllEntity1_2]
#! [PySnippetReadFieldOnAllEntity1_3]
medfileField1TS = MEDFileField1TS.New(fname, fieldName, iteration,
order)
mm = MEDFileMesh.New(fname)
fread = medfileField1TS.getFieldOnMeshAtLevel(ON_CELLS, 0, mm)
fread2 = medfileField1TS.getFieldAtLevel(ON_CELLS, 0)
self.assertTrue(fread.isEqual(f, 1e-12, 1e-12))
self.assertTrue(fread2.isEqual(f, 1e-12, 1e-12))
#! [PySnippetReadFieldOnAllEntity1_3]
#! [PySnippetReadFieldOnAllEntity1_4]
medfileFieldMTS = MEDFileFieldMultiTS.New(fname, fieldName)
mm = MEDFileMesh.New(fname)
fread = medfileFieldMTS.getFieldOnMeshAtLevel(ON_CELLS, iteration,
order, 0, mm)
fread2 = medfileFieldMTS.getFieldAtLevel(ON_CELLS, iteration, order, 0)
self.assertTrue(fread.isEqual(f, 1e-12, 1e-12))
self.assertTrue(fread2.isEqual(f, 1e-12, 1e-12))
#! [PySnippetReadFieldOnAllEntity1_4]
#! [PySnippetReadFieldOnAllEntity1_5]
medfileFieldMTS = MEDFileFieldMultiTS.New(fname, fieldName)
for medfileField1TS in medfileFieldMTS:
if medfileField1TS.getTime()[:2] == [iteration, order]:
fread = medfileField1TS.getFieldOnMeshAtLevel(ON_CELLS, 0, mm)
fread2 = medfileField1TS.getFieldAtLevel(ON_CELLS, 0)
self.assertTrue(fread.isEqual(f, 1e-12, 1e-12))
self.assertTrue(fread2.isEqual(f, 1e-12, 1e-12))
pass
pass
#! [PySnippetReadFieldOnAllEntity1_5]
pass
def testExampleReadFieldOnAllEntity1(self):
from MEDLoaderDataForTest import MEDLoaderDataForTest
#! [PySnippetReadFieldOnAllEntity1_1]
fname="PyExamples1.med"
meshName="mesh"
fieldName="FieldOnAll"
iteration=3
order=4
#! [PySnippetReadFieldOnAllEntity1_1]
#! [PySnippetWriteFieldOnAllEntity1_2]
m=MEDLoaderDataForTest.build2DMesh_3()
m=m[:10]
m.setName(meshName)
f=m.getMeasureField(False)
f=f.buildNewTimeReprFromThis(ONE_TIME,False)
f.setTime(5.5,iteration,order)
f.setName(fieldName)
# MEDCoupling finished, MEDLoader advanced API specific part starting from here
mm=MEDFileUMesh.New()
mm.setMeshAtLevel(0,m)
ff=MEDFileField1TS.New()
ff.setFieldNoProfileSBT(f)
mm.write(fname,2)
ff.write(fname,0)
#! [PySnippetWriteFieldOnAllEntity1_2]
#! [PySnippetReadFieldOnAllEntity1_3]
medfileField1TS=MEDFileField1TS.New(fname,fieldName,iteration,order)
mm=MEDFileMesh.New(fname)
fread=medfileField1TS.getFieldOnMeshAtLevel(ON_CELLS,0,mm)
fread2=medfileField1TS.getFieldAtLevel(ON_CELLS,0)
self.assertTrue(fread.isEqual(f,1e-12,1e-12))
self.assertTrue(fread2.isEqual(f,1e-12,1e-12))
#! [PySnippetReadFieldOnAllEntity1_3]
#! [PySnippetReadFieldOnAllEntity1_4]
medfileFieldMTS=MEDFileFieldMultiTS.New(fname,fieldName)
mm=MEDFileMesh.New(fname)
fread=medfileFieldMTS.getFieldOnMeshAtLevel(ON_CELLS,iteration,order,0,mm)
fread2=medfileFieldMTS.getFieldAtLevel(ON_CELLS,iteration,order,0)
self.assertTrue(fread.isEqual(f,1e-12,1e-12))
self.assertTrue(fread2.isEqual(f,1e-12,1e-12))
#! [PySnippetReadFieldOnAllEntity1_4]
#! [PySnippetReadFieldOnAllEntity1_5]
medfileFieldMTS=MEDFileFieldMultiTS.New(fname,fieldName)
for medfileField1TS in medfileFieldMTS:
if medfileField1TS.getTime()[:2]==[iteration,order]:
fread=medfileField1TS.getFieldOnMeshAtLevel(ON_CELLS,0,mm)
fread2=medfileField1TS.getFieldAtLevel(ON_CELLS,0)
self.assertTrue(fread.isEqual(f,1e-12,1e-12))
self.assertTrue(fread2.isEqual(f,1e-12,1e-12))
pass
pass
#! [PySnippetReadFieldOnAllEntity1_5]
pass
m = fs[0].getMesh()
WriteMesh("file5.med", m, True)
for f in fs:
assert (f.getMesh().getHiddenCppPointer() == m.getHiddenCppPointer())
# extra line to insist on the fact that
WriteFieldUsingAlreadyWrittenMesh("file5.med", f)
#! [UG_ReadField_6]
from MEDLoaderDataForTest import MEDLoaderDataForTest
fname = "PyExamples1.med"
meshName = "mesh"
fieldName = "FieldOnAll"
iteration = 3
order = 4
m = MEDLoaderDataForTest.build2DMesh_3()
m.setName(meshName)
f = m.getMeasureField(False)
f = f.buildNewTimeReprFromThis(ONE_TIME, False)
f.setTime(5.5, iteration, order)
f.setName(fieldName)
mesh = m
field = f
level = 0
#! [UG_RWFieldAdv_0]
from MEDLoader import MEDFileUMesh, MEDFileField1TS
mm = MEDFileUMesh.New()
mm.setMeshAtLevel(0, mesh)
ff = MEDFileField1TS.New()
ff.setFieldNoProfileSBT(field)
