cl.checkCoherency2()
li4=[sqrt(2.)*elt for elt in [1.,3.,5.,7.]]
li4_1=[0.5,0.5,2.5,2.5,6.5,6.5,12.5,12.5]
+ self.assertEqual(2,cl.getSpaceDimension())
+ self.assertEqual(1,cl.getMeshDimension())
+ self.assertEqual(4,cl.getNumberOfCells())
+ self.assertEqual(5,cl.getNumberOfNodes())
self.assertTrue(cl.getMeasureField(False).getArray().isEqual(DataArrayDouble(li4),1e-14))
self.assertTrue(cl.buildUnstructured().getMeasureField(False).getArray().isEqual(DataArrayDouble(li4),1e-14))
self.assertTrue(cl.getBarycenterAndOwner().isEqual(DataArrayDouble(li4_1,4,2),1e-14))
self.assertTrue(m.getBoundingBoxForBBTree().isEqual(DataArrayDouble([-0.5,0.5,-0.5,0.5,-0.5,0.5,-0.5,-0.31819805153394637],2,4),1e-12))
pass
+ def testSwig2CartBuildUnstructuredOnExoticCases1(self):
+ """ Test focusing on traduction from cartesian to unstructured mesh when spaceDim greater than meshDim.
+ """
+ #
+ m=MEDCouplingCMesh()
+ arrX=DataArrayDouble(3) ; arrX.iota()
+ arrY=DataArrayDouble(4) ; arrY.iota()
+ arrZ=DataArrayDouble(1) ; arrZ.iota()
+ m.setCoords(arrX,arrY,arrZ)
+ self.assertEqual(2,m.getMeshDimension())
+ self.assertEqual(3,m.getSpaceDimension())
+ mu=m.buildUnstructured()
+ self.assertTrue(mu.getNodalConnectivity().isEqual(DataArrayInt([4,1,0,3,4,4,2,1,4,5,4,4,3,6,7,4,5,4,7,8,4,7,6,9,10,4,8,7,10,11])))
+ self.assertTrue(mu.getNodalConnectivityIndex().isEqual(DataArrayInt([0,5,10,15,20,25,30])))
+ coo0=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)])
+ self.assertTrue(mu.getCoords().isEqual(coo0,1e-12))
+ mu.writeVTK("tutu.vtu")
+ #
+ m=MEDCouplingCMesh()
+ arrX=DataArrayDouble(3) ; arrX.iota()
+ arrY=DataArrayDouble(1) ; arrY.iota()
+ arrZ=DataArrayDouble(4) ; arrZ.iota()
+ m.setCoords(arrX,arrY,arrZ)
+ self.assertEqual(2,m.getMeshDimension())
+ self.assertEqual(3,m.getSpaceDimension())
+ mu=m.buildUnstructured()
+ self.assertTrue(mu.getNodalConnectivity().isEqual(DataArrayInt([4,1,0,3,4,4,2,1,4,5,4,4,3,6,7,4,5,4,7,8,4,7,6,9,10,4,8,7,10,11])))
+ self.assertTrue(mu.getNodalConnectivityIndex().isEqual(DataArrayInt([0,5,10,15,20,25,30])))
+ coo1=DataArrayDouble([(0,0,0),(1,0,0),(2,0,0),(0,0,1),(1,0,1),(2,0,1),(0,0,2),(1,0,2),(2,0,2),(0,0,3),(1,0,3),(2,0,3)])
+ self.assertTrue(mu.getCoords().isEqual(coo1,1e-12))
+ #
+ m=MEDCouplingCMesh()
+ arrX=DataArrayDouble(1) ; arrX.iota() ; arrX+=9
+ arrY=DataArrayDouble(3) ; arrY.iota()
+ arrZ=DataArrayDouble(4) ; arrZ.iota()
+ m.setCoords(arrX,arrY,arrZ)
+ self.assertEqual(2,m.getMeshDimension())
+ self.assertEqual(3,m.getSpaceDimension())
+ mu=m.buildUnstructured()
+ self.assertTrue(mu.getNodalConnectivity().isEqual(DataArrayInt([4,1,0,3,4,4,2,1,4,5,4,4,3,6,7,4,5,4,7,8,4,7,6,9,10,4,8,7,10,11])))
+ self.assertTrue(mu.getNodalConnectivityIndex().isEqual(DataArrayInt([0,5,10,15,20,25,30])))
+ coo2=DataArrayDouble([(9,0,0),(9,1,0),(9,2,0),(9,0,1),(9,1,1),(9,2,1),(9,0,2),(9,1,2),(9,2,2),(9,0,3),(9,1,3),(9,2,3)])
+ self.assertTrue(mu.getCoords().isEqual(coo2,1e-12))
+ #
+ m=MEDCouplingCMesh()
+ arrX=DataArrayDouble(3) ; arrX.iota()
+ arrY=DataArrayDouble(1) ; arrY.iota(7)
+ arrZ=DataArrayDouble(1) ; arrZ.iota(8)
+ m.setCoords(arrX,arrY,arrZ)
+ self.assertEqual(1,m.getMeshDimension())
+ self.assertEqual(3,m.getSpaceDimension())
+ mu=m.buildUnstructured()
+ self.assertTrue(mu.getNodalConnectivity().isEqual(DataArrayInt([1,0,1,1,1,2])))
+ self.assertTrue(mu.getNodalConnectivityIndex().isEqual(DataArrayInt([0,3,6])))
+ coo3=DataArrayDouble([(0,7,8),(1,7,8),(2,7,8)])
+ self.assertTrue(mu.getCoords().isEqual(coo3,1e-12))
+ #
+ m=MEDCouplingCMesh()
+ arrX=DataArrayDouble(1) ; arrX.iota(7)
+ arrY=DataArrayDouble(1) ; arrY.iota(8)
+ arrZ=DataArrayDouble(3) ; arrZ.iota()
+ m.setCoords(arrX,arrY,arrZ)
+ self.assertEqual(1,m.getMeshDimension())
+ self.assertEqual(3,m.getSpaceDimension())
+ mu=m.buildUnstructured()
+ self.assertTrue(mu.getNodalConnectivity().isEqual(DataArrayInt([1,0,1,1,1,2])))
+ self.assertTrue(mu.getNodalConnectivityIndex().isEqual(DataArrayInt([0,3,6])))
+ coo4=DataArrayDouble([(7,8,0),(7,8,1),(7,8,2)])
+ self.assertTrue(mu.getCoords().isEqual(coo4,1e-12))
+ #
+ m=MEDCouplingCMesh()
+ arrX=DataArrayDouble(3) ; arrX.iota()
+ arrY=DataArrayDouble(1) ; arrY.iota(7)
+ m.setCoords(arrX,arrY)
+ self.assertEqual(1,m.getMeshDimension())
+ self.assertEqual(2,m.getSpaceDimension())
+ mu=m.buildUnstructured()
+ self.assertTrue(mu.getNodalConnectivity().isEqual(DataArrayInt([1,0,1,1,1,2])))
+ self.assertTrue(mu.getNodalConnectivityIndex().isEqual(DataArrayInt([0,3,6])))
+ coo5=DataArrayDouble([(0,7),(1,7),(2,7)])
+ self.assertTrue(mu.getCoords().isEqual(coo5,1e-12))
+ #
+ m=MEDCouplingCMesh()
+ arrX=DataArrayDouble(1) ; arrX.iota(7)
+ arrY=DataArrayDouble(3) ; arrY.iota()
+ m.setCoords(arrX,arrY)
+ self.assertEqual(1,m.getMeshDimension())
+ self.assertEqual(2,m.getSpaceDimension())
+ mu=m.buildUnstructured()
+ self.assertTrue(mu.getNodalConnectivity().isEqual(DataArrayInt([1,0,1,1,1,2])))
+ self.assertTrue(mu.getNodalConnectivityIndex().isEqual(DataArrayInt([0,3,6])))
+ coo6=DataArrayDouble([(7,0),(7,1),(7,2)])
+ self.assertTrue(mu.getCoords().isEqual(coo6,1e-12))
+ pass
+
+ def testSwig2Colinearize2D2(self):
+ """ simple non regression test but that has revealed a bug"""
+ coo=DataArrayDouble([(0,0),(0,0.5),(0,1),(1,1),(1,0),(0.5,0)])
+ m=MEDCouplingUMesh("mesh",2) ; m.setCoords(coo)
+ m.allocateCells() ; m.insertNextCell(NORM_POLYGON,[0,1,2,3,4,5])
+ m.checkCoherency2()
+ refPtr=m.getCoords().getHiddenCppPointer()
+ #
+ m.colinearize2D(1e-12)
+ m.checkCoherency2()
+ self.assertEqual(refPtr,m.getCoords().getHiddenCppPointer())
+ self.assertTrue(m.getNodalConnectivity().isEqual(DataArrayInt([5,0,2,3,4])))
+ self.assertTrue(m.getNodalConnectivityIndex().isEqual(DataArrayInt([0,5])))
+ pass
+
def setUp(self):
pass
pass
