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.getNodalConnectivityIndex().isEqual(DataArrayInt([0,4])))
## false quadratic
coo2=DataArrayDouble([(-5,0),(-1,0),(4,3),(7,0),(1,6),(1,0),(-3,0),(6,1),(5,0),(3,0),(2,0),(4,0),(6,0),(6.5,0.5),(5,2),(2.5,4.5),(-2,3),(-4,0),(-2,0),(0,0)])
+ coo2.setInfoOnComponents(["aa","bbbb"])
m=MEDCouplingUMesh("mesh",2) ; m.setCoords(coo2) ; m.allocateCells()
m.insertNextCell(NORM_QPOLYG,[5,9,8,3,7,2,4,0,6,1,10,11,12,13,14,15,16,17,18,19])
refPtr=m.getCoords().getHiddenCppPointer()
self.assertTrue(m.colinearize2D(1e-12).isEqual(DataArrayInt([0])))
self.assertNotEqual(refPtr,m.getCoords().getHiddenCppPointer())#not same coordinates here
+ self.assertEqual(["aa","bbbb"],m.getCoords().getInfoOnComponents())
refPtr=m.getCoords().getHiddenCppPointer()
self.assertTrue(coo2.isEqual(m.getCoords()[:20],1e-12))
- self.assertTrue(m.getCoords()[20:].isEqual(DataArrayDouble([(1.,0.),(4.,3.)]),1e-12))
+ self.assertTrue(m.getCoords()[20:].isEqualWithoutConsideringStr(DataArrayDouble([(1.,0.),(4.,3.)]),1e-12))
self.assertTrue(m.getNodalConnectivity().isEqual(DataArrayInt([32,0,3,4,20,21,16])))
self.assertTrue(m.getNodalConnectivityIndex().isEqual(DataArrayInt([0,7])))
self.assertTrue(m.colinearize2D(1e-12).isEqual(DataArrayInt([])))
self.assertTrue(m.getNodalConnectivityIndex().isEqual(DataArrayInt([0,7])))
pass
+ def testSwig2BoundingBoxForBBTree1(self):
+ """ This test appears simple but it checks that bounding box are correctly computed for quadratic polygons. It can help a lot to reduce the amount of intersections !
+ """
+ coo=DataArrayDouble([-0.5,-0.5,-0.5,0.5,0.5,0.5,0.5,-0.5,0.45,0.,0.3181980515339464,0.31819805153394637,0.,0.45,-0.31819805153394637,0.3181980515339464,-0.45,0.,-0.3181980515339465,-0.31819805153394637,0.,-0.45,0.3181980515339463,-0.3181980515339465,-0.5,0.0,0.0,0.5,0.5,0.0,0.0,-0.5,-0.4090990257669732,-0.4090990257669732,0.40909902576697316,-0.4090990257669732],18,2)
+ m=MEDCouplingUMesh("mesh",2) ; m.setCoords(coo)
+ m.allocateCells()
+ m.insertNextCell(NORM_QPOLYG,[0,1,2,3,11,5,7,9,12,13,14,17,4,6,8,16])
+ m.insertNextCell(NORM_QPOLYG,[3,0,9,11,15,16,10,17])
+ 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
