level3(self)
gtumesh(self)
pass
+
+ def testPenta18_1(self):
+ arr=DataArrayDouble([
+ (0.,1.,1.),(0.,0.,1.),(1.,0.,1.),
+ (0.,1.,0.),(0.,0.,0.),(1.,0.,0.),
+ (0.,0.5,1.),(0.5,0.,1.),(0.5,0.5,1.),
+ (0.,0.5,0.),(0.5,0.,0.),(0.5,0.5,0.),
+ (0.,1.,0.5),(0.,0.,0.5),(1.,0.,0.5),
+ (0.,0.5,0.5),(0.5,0.,0.5),(0.5,0.5,0.5)])
+ m=MEDCouplingUMesh("mesh",3)
+ m.setCoords(arr)
+ m.allocateCells(1)
+ m.insertNextCell(NORM_PENTA18,list(range(18)))
+ m.checkConsistencyLight()
+ self.assertTrue(m.getMeasureField(True).getArray().isEqual(DataArrayDouble([0.5]),1e-12))
+ #
+ f=MEDCouplingFieldDouble(ON_NODES)
+ f.setMesh(m)
+ f.setName("FieldOnPenta18")
+ f.setArray(DataArrayDouble(list(range(18))))
+ f.checkConsistencyLight()
+ #
+ m2,d,di,rd,rdi=m.buildDescendingConnectivity()
+ self.assertTrue(m2.getNodalConnectivity().isEqual(DataArrayInt([6,0,1,2,6,7,8,6,3,5,4,11,10,9,9,0,3,4,1,12,9,13,6,15,9,1,4,5,2,13,10,14,7,16,9,2,4,5,0,14,11,12,8,17])))
+ self.assertTrue(m2.getNodalConnectivityIndex().isEqual(DataArrayInt([0,7,14,24,34,44])))
+ self.assertTrue(d.isEqual(DataArrayInt([0,1,2,3,4])))
+ self.assertTrue(di.isEqual(DataArrayInt([0,5])))
+ self.assertTrue(rd.isEqual(DataArrayInt([0,0,0,0,0])))
+ self.assertTrue(rdi.isEqual(DataArrayInt([0,1,2,3,4,5])))
+ #
+ f2=MEDCouplingFieldDouble(ON_NODES)
+ f2.setMesh(m)
+ f2.setName("FieldOnPenta18Sub")
+ f2.setArray(DataArrayDouble(list(range(18))))
+ f2.checkConsistencyLight()
+ pass
+
+ def testSKLAReplaceDeletePacks(self):
+ index=DataArrayInt([0,3,5,6,6])
+ value=DataArrayInt([1,2,3, 2,3, 3 ])
+ sla=MEDCouplingSkyLineArray(index,value)
+ idx=DataArrayInt([0,3])
+ packs=[DataArrayInt([4,5]),DataArrayInt([6,7,8])]
+ sla.replaceSimplePacks(idx,packs)
+ self.assertTrue(sla.getIndexArray().isEqual(DataArrayInt([0,2,4,5,8])))
+ self.assertTrue(sla.getValuesArray().isEqual(DataArrayInt([4,5, 2,3, 3, 6,7,8])))
+ sla.deleteSimplePacks(idx)
+ self.assertTrue(sla.getIndexArray().isEqual(DataArrayInt([0,2,3])))
+ self.assertTrue(sla.getValuesArray().isEqual(DataArrayInt([2,3, 3])))
+ sla.deleteSimplePack(1)
+ self.assertTrue(sla.getIndexArray().isEqual(DataArrayInt([0,2])))
+ self.assertTrue(sla.getValuesArray().isEqual(DataArrayInt([2,3])))
+ pass
+
+ def testDADAsArcOfCircle(self):
+ d=DataArrayDouble([3.06915124862645,2.1464466094067824,2.85355345827285,2.3620444674400574,2.637955532559882,2.1464467447661937],3,2)
+ center,radius,ang=d.asArcOfCircle()
+ self.assertTrue((d-center).magnitude().isUniform(radius,1e-10))
+ self.assertAlmostEqual(ang,-4.712389294301196,12)
+ pass
+
+ def testDAMaxAbsValue(self):
+ d=DataArrayDouble([-2,3,1.2,-2.9])
+ a,b=d.getMaxAbsValue()
+ self.assertAlmostEqual(a,3.,13)
+ self.assertEqual(b,1)
+ a,b=(-d).getMaxAbsValue()
+ self.assertAlmostEqual(a,-3.,13)
+ self.assertEqual(b,1)
+ self.assertAlmostEqual((-d).getMaxAbsValueInArray(),-3.,13)
+ pass
+
+ def testDAIFindIdForEach1(self):
+ a1=DataArrayInt([17,27,2,10,-4,3,12,27,16])
+ b1=DataArrayInt([3,16,-4,27,17])
+ ret=a1.findIdForEach(b1)
+ self.assertTrue(ret.isEqual(DataArrayInt([5,8,4,7,0])))
+ self.assertTrue(a1[ret].isEqual(b1))
+ b2=DataArrayInt([3,16,22,27,17])
+ self.assertRaises(InterpKernelException,a1.findIdForEach,b2) # 22 not in a1 !
+ a1.rearrange(3)
+ self.assertRaises(InterpKernelException,a1.findIdForEach,b1) # a1 is not single component
+ pass
+
+ def testAttractSeg3MidPtsAroundNodes1(self):
+ """ Test of MEDCouplingUMesh.attractSeg3MidPtsAroundNodes methods """
+ ptsExpToBeModified=DataArrayInt([95,96,97,98,101,103,104,106,108,110])
+ eps=1e-12
+ a=2./3.
+ b=1./3.
+ coo=DataArrayDouble([10,0,0,10,10,0,10,0,3.+b,10,0,6.+a,10,10,3.+b,10,10,6.+a,10,3.+b,0,10,6.+a,0,3.+b,0,0,6.+a,0,0,3.+b,10,0,6.+a,10,0,10,3.+b,6.+a,10,6.+a,6.+a,10,3.+b,3.+b,10,6.+a,3.+b,3.+b,0,3.+b,3.+b,0,6.+a,6.+a,0,3.+b,6.+a,0,6.+a,6.+a,10,3.+b,6.+a,10,6.+a,3.+b,10,3.+b,3.+b,10,6.+a,3.+b,3.+b,0,6.+a,3.+b,0,3.+b,6.+a,0,6.+a,6.+a,0,3.+b,3.+b,6.+a,6.+a,3.+b,6.+a,3.+b,6.+a,6.+a,6.+a,6.+a,6.+a,3.+b,3.+b,3.+b,6.+a,3.+b,3.+b,3.+b,6.+a,3.+b,6.+a,6.+a,3.+b,10,0,1.+a,10,0,5.,10,10,1.+a,10,10,5.,10,1.+a,0,10,5.,0,10,8.+b,0,5.,0,0,8.+b,0,0,5.,10,0,8.+b,10,0,10,1.+a,6.+a,10,5.,6.+a,10,8.+b,6.+a,10,1.+a,3.+b,10,3.+b,5.,10,5.,3.+b,10,6.+a,5.,10,8.+b,3.+b,10,3.+b,1.+a,10,6.+a,1.+a,3.+b,0,1.+a,3.+b,0,5.,6.+a,0,1.+a,5.,0,3.+b,6.+a,0,5.,5.,0,6.+a,8.+b,0,3.+b,8.+b,0,6.+a,6.+a,10,1.+a,8.+b,10,3.+b,6.+a,10,5.,8.+b,10,6.+a,3.+b,10,1.+a,5.,10,3.+b,3.+b,10,5.,5.,10,6.+a,3.+b,1.+a,0,5.,3.+b,0,6.+a,1.+a,0,8.+b,3.+b,0,3.+b,5.,0,5.,6.+a,0,6.+a,5.,0,8.+b,6.+a,0,3.+b,8.+b,0,6.+a,8.+b,0,3.+b,1.+a,6.+a,6.+a,1.+a,6.+a,5.,3.+b,6.+a,8.+b,3.+b,6.+a,3.+b,5.,6.+a,6.+a,5.,6.+a,5.,6.+a,6.+a,8.+b,6.+a,6.+a,3.+b,8.+b,6.+a,6.+a,8.+b,6.+a,3.+b,3.+b,5,3.+b,1.+a,3.+b,6.+a,3.+b,5.,6.+a,1.+a,3.+b,5.,3.+b,3.+b,8.+b,3.+b,3.+b,3.+b,6.+a,5.,3.+b,5.,3.+b,6.+a,6.+a,5.,6.+a,5.,3.+b,5.,6.+a,3.+b,8.+b,6.+a,3.+b,3.+b,8.+b,3.+b,6.+a,8.+b,3.+b,3.+b,3.+b,1.+a,6.+a,3.+b,1.+a,3.+b,6.+a,1.+a,6.+a,6.+a,1.+a],111,3)
+ conn=DataArrayInt([30,17,28,32,16,19,29,33,18,83,93,94,58,84,95,96,61,62,85,97,60,30,19,29,33,18,3,12,14,2,84,95,96,61,47,51,50,37,64,86,98,63,30,28,30,34,32,29,31,35,33,87,99,100,93,88,101,102,95,85,89,103,97,30,29,31,35,33,12,13,15,14,88,101,102,95,48,53,52,51,86,90,104,98,30,30,23,22,34,31,21,20,35,91,71,105,99,92,67,106,101,89,72,70,103,30,31,21,20,35,13,5,4,15,92,67,106,101,49,39,54,53,90,68,66,104,30,16,32,24,8,18,33,25,9,94,107,73,57,96,108,75,59,60,97,74,43,30,18,33,25,9,2,14,6,0,96,108,75,59,50,55,40,36,63,98,76,44,30,32,34,26,24,33,35,27,25,100,109,77,107,102,110,79,108,97,103,78,74,30,33,35,27,25,14,15,7,6,102,110,79,108,52,56,41,55,98,104,80,76,30,34,22,10,26,35,20,11,27,105,69,81,109,106,65,82,110,103,70,45,78,30,35,20,11,27,15,4,1,7,106,65,82,110,54,38,42,56,104,66,46,80])
+ connI=DataArrayInt([0,21,42,63,84,105,126,147,168,189,210,231,252])
+ m=MEDCouplingUMesh("mesh",3)
+ m.setConnectivity(conn,connI,True)
+ m.setCoords(coo.deepCopy())# deep copy coo because next line is going to modify it, if it works normaly
+ m.attractSeg3MidPtsAroundNodes(0.1,DataArrayInt([33,35])) # ze call is here !
+ self.assertTrue(not m.getCoords().isEqual(coo,eps)) # some points have had their position changed...
+ ptsExpNotToBeModified=ptsExpToBeModified.buildComplement(len(coo))
+ self.assertTrue(m.getCoords()[ptsExpNotToBeModified].isEqual(coo[ptsExpNotToBeModified],eps))
+ self.assertTrue((m.getCoords()[ptsExpToBeModified]-coo[ptsExpToBeModified]).magnitude().isUniform(4./3.,1e-12))
+ ptsPosExp=DataArrayDouble([6.+a,3.+b,3.+a,6.+a,3.,3.+b,6.+b,3.+b,3.+b,7.,3.+b,3.+b,6.+a,6.+a,3.+a,6.+b,6.+a,3.+b,7.,6.+a,3.+b,6.+a,7.,3.+b,6.+a,3.+b,3.,6.+a,6.+a,3.],10,3)
+ self.assertTrue(m.getCoords()[ptsExpToBeModified].isEqual(ptsPosExp,1e-12))
+ pass
+
pass
if __name__ == '__main__':
