ff1.write(fname,0)
f1=ff1.getFieldOnMeshAtLevel(ON_GAUSS_NE,m1,0)
f2,p1=ff1.getFieldWithProfile(ON_GAUSS_NE,0,mm1) ; f2.setName("")
- self.assertTrue(p1.isIdentity())
- self.assertEqual(5,p1.getNumberOfTuples())
+ self.assertTrue(p1.isIdentity2(5))
self.assertTrue(f1.getArray().isEqual(f2,1e-12))
pass
# Test for getFieldAtTopLevel method
for elt in [[0,1,2,3,4,5],[1,2,3,4,5,6],[2,3,4,5,6,7],[3,4,5,6,7,8]]:#4
m0.insertNextCell(NORM_PENTA6,elt)
pass
- m0.checkCoherency2()
+ m0.checkCoherency1()
m1=MEDCouplingUMesh(); m1.setName("mesh")
m1.setMeshDimension(2);
m1.allocateCells(5);
m.changeSpaceDimension(3,0.)
infos=["aa [b]","cc [de]","gg [klm]"]
m.getCoords().setInfoOnComponents(infos)
- m.checkCoherency2()
+ m.checkCoherency1()
mm=MEDFileUMesh()
mm.setMeshAtLevel(0,m)
m1=MEDCouplingCMesh() ; m1.setCoords(arr) ; m1.setName("Mesh")
m.changeSpaceDimension(3,0.)
infos=["aa [b]","cc [de]","gg [klm]"]
m.getCoords().setInfoOnComponents(infos)
- m.checkCoherency2()
+ m.checkCoherency1()
f=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f.setMesh(m)
f.setName("Field")
arr=DataArrayDouble(25,2) ; arr.setInfoOnComponents(compos)
m.changeSpaceDimension(3,0.)
infos=["aa [b]","cc [de]","gg [klm]"]
m.getCoords().setInfoOnComponents(infos)
- m.checkCoherency2()
+ m.checkCoherency1()
f=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f.setMesh(m)
f.setName("Field")
arr=DataArrayDouble(25,2) ; arr.setInfoOnComponents(compos)
m3D=m.buildExtrudedMesh(m1D,0)
m3D.sortCellsInMEDFileFrmt()
m3D.setName(meshName)
- m2D=m ; m2D.setCoords(m3D.getCoords()) ; m2D.shiftNodeNumbersInConn(delta) ; m2D.setName(meshName) ; m2D.checkCoherency2()
+ m2D=m ; m2D.setCoords(m3D.getCoords()) ; m2D.shiftNodeNumbersInConn(delta) ; m2D.setName(meshName) ; m2D.checkCoherency1()
m1D=m2D.computeSkin() ; m1D.setName(meshName)
m0D=MEDCouplingUMesh.Build0DMeshFromCoords(m3D.getCoords()) ; m0D.setName(meshName) ; m0D=m0D[[2,4,10]]
#
""" This test is advanced to be sure that no unnecessary copies had been made during cartesianization process. """
# UMesh non cart
arr=DataArrayDouble(4) ; arr.iota() ; m=MEDCouplingCMesh() ; m.setCoords(arr,arr) ; m=m.buildUnstructured()
- mm=MEDFileUMesh() ; mm[0]=m ; mm.forceComputationOfParts() ; mm.setAxType(AX_CYL) #<- important
+ mm=MEDFileUMesh() ; mm[0]=m ; mm.forceComputationOfParts()
d0=DataArrayInt(16) ; d0[:]=0
d1=DataArrayInt(9) ; d1[:]=0
mm.setFamilyFieldArr(0,d1) ; mm.setFamilyFieldArr(1,d0)
ref2=mm[0].getNodalConnectivityIndex().getHiddenCppPointer()
ref3=mm.getDirectUndergroundSingleGeoTypeMesh(NORM_QUAD4).getNodalConnectivity().getHiddenCppPointer()
self.assertEqual(ref0,mm.getDirectUndergroundSingleGeoTypeMesh(NORM_QUAD4).getCoords().getHiddenCppPointer())
+ mm.setAxType(AX_CYL) #<- important
mm2=mm.cartesianize() # the trigger
+ self.assertEqual(mm2.getAxType(),AX_CART)
+ mm.setAxType(AX_CART) # this is here only to avoid complaints
self.assertTrue(isinstance(mm2,MEDFileUMesh))
self.assertTrue(mm.getHiddenCppPointer()!=mm2.getHiddenCppPointer())
self.assertTrue(ref0==mm.getCoords().getHiddenCppPointer()) # <- here important
# UMesh cart
mm.setAxType(AX_CART)
mm2=mm.cartesianize() # the trigger
+ self.assertEqual(mm2.getAxType(),AX_CART)
self.assertTrue(isinstance(mm2,MEDFileUMesh))
self.assertTrue(mm.getHiddenCppPointer()==mm2.getHiddenCppPointer()) # optimization
# CurveLinearMesh non cart
mm.setName("a") ; mm.setDescription("b") ; mm.setTime(3,4,5.) ; mm.addFamily("c",-4) ; mm.setFamiliesOnGroup("d",["c"]) ; mm.setTimeUnit("ms")
ref0=mm.getMesh().getCoords().getHiddenCppPointer()
mm2=mm.cartesianize() # the trigger
+ self.assertEqual(mm2.getAxType(),AX_CART)
self.assertTrue(isinstance(mm2,MEDFileCurveLinearMesh))
self.assertTrue(mm.getHiddenCppPointer()!=mm2.getHiddenCppPointer())
self.assertTrue(ref0==mm.getMesh().getCoords().getHiddenCppPointer()) # <- here important
# CurveLinearMesh cart
mm.setAxType(AX_CART)
mm2=mm.cartesianize() # the trigger
+ self.assertEqual(mm2.getAxType(),AX_CART)
self.assertTrue(isinstance(mm2,MEDFileCurveLinearMesh))
self.assertTrue(mm.getHiddenCppPointer()==mm2.getHiddenCppPointer()) # optimization
# CMesh non cart
mm.setFamilyFieldArr(0,d1) ; mm.setFamilyFieldArr(1,d0)
mm.setName("a") ; mm.setDescription("b") ; mm.setTime(3,4,5.) ; mm.addFamily("c",-4) ; mm.setFamiliesOnGroup("d",["c"]) ; mm.setTimeUnit("ms")
mm2=mm.cartesianize() # the trigger
+ self.assertEqual(mm2.getAxType(),AX_CART)
self.assertTrue(isinstance(mm2,MEDFileCurveLinearMesh))
self.assertEqual(mm2.getMesh().getNodeGridStructure(),mm.getMesh().getNodeGridStructure())
self.assertEqual(mm2.getName(),mm.getName())
# CMesh cart
mm.setAxType(AX_CART)
mm2=mm.cartesianize() # the trigger
+ self.assertEqual(mm2.getAxType(),AX_CART)
self.assertTrue(isinstance(mm2,MEDFileCMesh))
self.assertTrue(mm.getHiddenCppPointer()==mm2.getHiddenCppPointer()) # optimization
pass