self.assertRaises(InterpKernelException,field.setNature,Integral);
self.assertRaises(InterpKernelException,field.setNature,IntegralGlobConstraint);
pass
+
+ def testNatureOperations(self):
+ """ Check nature constraints on field operations """
+ m = MEDCouplingCMesh()
+ m.setCoordsAt(0, DataArrayDouble([1.0,2.0,3.0]))
+ m.setCoordsAt(1, DataArrayDouble([1.0,2.0,3.0]))
+ m = m.buildUnstructured()
+ f1, f2 = MEDCouplingFieldDouble.New(ON_CELLS, NO_TIME), MEDCouplingFieldDouble.New(ON_CELLS, NO_TIME)
+ f1.setNature(Integral)
+ f2.setNature(ConservativeVolumic)
+ self.assertEqual(Integral, f1.getNature())
+ self.assertEqual(ConservativeVolumic, f2.getNature())
+
+ da = DataArrayDouble([1.0,2.0,3.0,4.0])
+ f1.setMesh(m); f2.setMesh(m)
+ f1.setArray(da); f2.setArray(da.deepCpy())
+ # All this should complain about nature:
+ self.assertRaises(InterpKernelException, f1.__add__, f2)
+ self.assertRaises(InterpKernelException, f1.__iadd__, f2)
+ self.assertRaises(InterpKernelException, f1.__sub__, f2)
+ self.assertRaises(InterpKernelException, f1.__isub__, f2)
+ self.assertRaises(InterpKernelException, f1.__radd__, f2)
+ self.assertRaises(InterpKernelException, f1.__rsub__, f2)
+ self.assertRaises(InterpKernelException, MEDCouplingFieldDouble.AddFields, f1, f2)
+ self.assertRaises(InterpKernelException, MEDCouplingFieldDouble.SubstractFields, f1, f2)
+ self.assertRaises(InterpKernelException, MEDCouplingFieldDouble.MaxFields, f1, f2)
+ self.assertRaises(InterpKernelException, MEDCouplingFieldDouble.MinFields, f1, f2)
+ # Not those ones:
+ f3 = MEDCouplingFieldDouble.MultiplyFields(f1,f2)
+ self.assertEqual(NoNature, f3.getNature())
+ f3 = f1*f2
+ self.assertEqual(NoNature, f3.getNature())
+ f1Tmp = f1.deepCpy(); f1Tmp.setMesh(m); f1Tmp *= f2
+ self.assertEqual(NoNature, f1Tmp.getNature())
+ f3 = MEDCouplingFieldDouble.DivideFields(f1,f2)
+ self.assertEqual(NoNature, f3.getNature())
+ f3 = f1/f2
+ self.assertEqual(NoNature, f3.getNature())
+ f1Tmp = f1.deepCpy(); f1Tmp.setMesh(m); f1Tmp /= f2
+ self.assertEqual(NoNature, f1Tmp.getNature())
+# f3 = MEDCouplingFieldDouble.PowFields(f1,f2)
+# self.assertEqual(NoNature, f3.getNature())
+ f3 = f1**f2
+ self.assertEqual(NoNature, f3.getNature())
+ f1Tmp = f1.deepCpy(); f1Tmp.setMesh(m); f1Tmp **= f2
+ self.assertEqual(NoNature, f1Tmp.getNature())
+ f3 = MEDCouplingFieldDouble.DotFields(f1,f2)
+ self.assertEqual(NoNature, f3.getNature())
+ f3 = f1.dot(f2)
+ self.assertEqual(NoNature, f3.getNature())
+
+ da = DataArrayDouble.Meld([da, da, da])
+ f1.setArray(da); f2.setArray(da.deepCpy())
+ f3 = MEDCouplingFieldDouble.CrossProductFields(f1,f2)
+ self.assertEqual(NoNature, f3.getNature())
+ f3 = f1.crossProduct(f2)
def testBuildSubMeshData(self):
targetMesh=MEDCouplingDataForTest.build2DTargetMesh_1()
f.setArray(arr)
f.checkCoherency()
self.assertTrue(DataArrayDouble(f.integral(False)).isEqual(DataArrayDouble([-211.66121638700983,-4863.9563007698835]),1e-11))
- self.assertTrue(DataArrayDouble(f.getWeightedAverageValue()).isEqual(DataArrayDouble([45.496085813113595,1045.496085813114]),1e-11))
+ self.assertTrue(DataArrayDouble(f.getWeightedAverageValue()).isEqual(DataArrayDouble([45.4960858131136,1045.496085813114]),1e-11))
self.assertTrue(DataArrayDouble(f.normL1()).isEqual(DataArrayDouble([45.49608581311362,1045.496085813114]),1e-11))
self.assertTrue(DataArrayDouble(f.normL2()).isEqual(DataArrayDouble([58.16846378340898,1046.1241521947334]),1e-11))
pass
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.getNodalConnectivity().isEqual(DataArrayInt([NORM_POLYGON,0,2,3,4])))
self.assertTrue(m.getNodalConnectivityIndex().isEqual(DataArrayInt([0,5])))
pass
+ def testSwig2Colinearize2D3(self):
+ """ colinearize was too agressive, potentially producing cells with one edge """
+ # Flat polygon with 3 edges - nothing should happen (min number of edges for a linear polyg)
+ coo = DataArrayDouble([0.0,0.0, 2.0,0.0, 1.5,0.0, 1.0,0.0, 0.5,0.0], 5,2)
+ m = MEDCouplingUMesh("m", 2)
+ c, cI = [DataArrayInt(l) for l in [[NORM_POLYGON, 0,1,2], [0,4]] ]
+ m.setCoords(coo); m.setConnectivity(c, cI)
+ m.colinearize2D(1e-10)
+ m.checkCoherency2()
+ self.assertEqual(c.getValues(), m.getNodalConnectivity().getValues())
+ self.assertEqual(cI.getValues(), m.getNodalConnectivityIndex().getValues())
+
+ # Flat quad polygon, 2 edges - nothing should happen (min number of edges for a quad polyg)
+ m = MEDCouplingUMesh("m", 2)
+ c, cI = [DataArrayInt(l) for l in [[NORM_QPOLYG, 0,1, 2,3], [0,5]] ]
+ m.setCoords(coo); m.setConnectivity(c, cI)
+ m.colinearize2D(1e-10)
+ m.checkCoherency2()
+ self.assertEqual(c.getValues(), m.getNodalConnectivity().getValues())
+ self.assertEqual(cI.getValues(), m.getNodalConnectivityIndex().getValues())
+
+ # Flat polygon, 4 edges - one reduction should happen
+ m = MEDCouplingUMesh("m", 2)
+ c, cI = [DataArrayInt(l) for l in [[NORM_POLYGON, 0,1,2,3], [0,5]] ]
+ m.setCoords(coo); m.setConnectivity(c, cI)
+ m.colinearize2D(1e-10)
+ m.checkCoherency2()
+ self.assertEqual([NORM_POLYGON, 3,1,2], m.getNodalConnectivity().getValues())
+ self.assertEqual([0,4], m.getNodalConnectivityIndex().getValues())
+
+ # Flat quad polygon, 3 edges - one reduction expected
+ m = MEDCouplingUMesh("m", 2)
+ c, cI = [DataArrayInt(l) for l in [[NORM_QPOLYG, 0,1,3, 3,2,4], [0,7]] ]
+ m.setCoords(coo); m.setConnectivity(c, cI)
+ m.colinearize2D(1e-10)
+ m.checkCoherency2()
+ self.assertEqual([NORM_QPOLYG, 3,1, 5,2], m.getNodalConnectivity().getValues())
+ self.assertTrue( m.getCoords()[5].isEqual( DataArrayDouble([(1.5,0.0)]), 1.0e-12 ) )
+ self.assertEqual([0,5], m.getNodalConnectivityIndex().getValues())
+
+ # Now an actual (neutronic) case: circle made of 4 SEG3. Should be reduced to 2 SEG3
+ m = MEDCouplingDataForTest.buildCircle2(0.0, 0.0, 1.0)
+ c, cI = [DataArrayInt(l) for l in [[NORM_QPOLYG, 7,5,3,1, 6,4,2,0], [0,9]] ]
+ m.colinearize2D(1e-10)
+ m.checkCoherency2()
+ self.assertEqual([NORM_QPOLYG, 3,5, 8,4], m.getNodalConnectivity().getValues())
+ self.assertTrue( m.getCoords()[8].isEqual( DataArrayDouble([(1.0,0.0)]), 1.0e-12 ) )
+ self.assertEqual([0,5], m.getNodalConnectivityIndex().getValues())
+
def testSwig2CheckAndPreparePermutation2(self):
a=DataArrayInt([10003,9999999,5,67])
self.assertTrue(DataArrayInt.CheckAndPreparePermutation(a).isEqual(DataArrayInt([2,3,0,1])))
m2.zipCoords()
self.assertTrue(m2.areAllNodesFetched())
pass
+
+ def testMEDCouplingPointSetComputeDiameterField1(self):
+ arrX=DataArrayDouble([0.,1.1,1.7,2.1])
+ arrY=DataArrayDouble([0.,0.7,0.8,1.9])
+ arrZ=DataArrayDouble([0.,1.3,2.1,2.4])
+ m=MEDCouplingCMesh() ; m.setCoords(arrX,arrY,arrZ) ; m=m.buildUnstructured()
+ f=m.computeDiameterField()
+ f.checkCoherency()
+ exp=DataArrayDouble([1.8411952639521971,1.5937377450509227,1.5297058540778357,1.705872210923198,1.4352700094407325,1.3638181696985856,2.0273134932713295,1.8055470085267789,1.7492855684535902,1.5297058540778357,1.2206555615733703,1.1357816691600546,1.3638181696985856,1.004987562112089,0.9,1.7492855684535902,1.4866068747318506,1.4177446878757824,1.3379088160259651,0.9695359714832656,0.8602325267042626,1.1445523142259597,0.6782329983125266,0.5099019513592785,1.5842979517754858,1.2884098726725124,1.208304597359457])
+ self.assertTrue(exp.isEqual(f.getArray(),1e-12))
+ m1=m[::2]
+ m2=m[1::2]
+ m2.simplexize(PLANAR_FACE_5)
+ m3=MEDCouplingUMesh.MergeUMeshesOnSameCoords(m1,m2)
+ f=m3.computeDiameterField()
+ f.checkCoherency()
+ exp2=DataArrayDouble([1.8411952639521971,1.5297058540778357,1.4352700094407325,2.0273134932713295,1.7492855684535902,1.2206555615733703,1.3638181696985856,0.9,1.4866068747318506,1.3379088160259651,0.8602325267042626,0.6782329983125266,1.5842979517754858,1.208304597359457,1.47648230602334,1.47648230602334,1.47648230602334,1.47648230602334,1.47648230602334,1.7029386365926402,1.7029386365926402,1.7029386365926402,1.7029386365926402,1.7029386365926402,1.3601470508735445,1.3601470508735445,1.3601470508735445,1.3601470508735445,1.3601470508735445,1.70293863659264,1.70293863659264,1.70293863659264,1.70293863659264,1.70293863659264,1.3601470508735445,1.3601470508735445,1.3601470508735445,1.3601470508735445,1.3601470508735445,1.063014581273465,1.063014581273465,1.063014581273465,1.063014581273465,1.063014581273465,1.0,1.0,1.0,1.0,1.0,1.5556349186104046,1.5556349186104046,1.5556349186104046,1.5556349186104046,1.5556349186104046,1.3601470508735443,1.3601470508735443,1.3601470508735443,1.3601470508735443,1.3601470508735443,0.9219544457292886,0.9219544457292886,0.9219544457292886,0.9219544457292886,0.9219544457292886,1.140175425099138,1.140175425099138,1.140175425099138,1.140175425099138,1.140175425099138,0.5,0.5,0.5,0.5,0.5,1.2529964086141667,1.2529964086141667,1.2529964086141667,1.2529964086141667,1.2529964086141667])
+ self.assertTrue(exp2.isEqual(f.getArray(),1e-12))
+ # TRI3 - spacedim = 2
+ coo=DataArrayDouble([(1,1),(5,1.9),(2.1,3)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_TRI3) ; m.setCoords(coo)
+ for c in [[0,1,2],[0,2,1],[2,1,0]]:
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],4.1,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],4.1,12)
+ m3=m.buildUnstructured() ; m3.convertLinearCellsToQuadratic(1)
+ self.assertAlmostEqual(m3.computeDiameterField().getArray()[0],4.1,12)
+ # TRI3 - spacedim = 3
+ coo=DataArrayDouble([(1.3198537928820775,1.0991902391274959,-0.028645697595823361),(5.2486835106806335,2.2234012799688281,0.30368935050077939),(2.2973688139447361,3.1572023778066649,0.10937756365410012)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_TRI3) ; m.setCoords(coo)
+ for c in [[0,1,2],[0,2,1],[2,1,0]]:
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],4.1,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],4.1,12)
+ m3=m.buildUnstructured() ; m3.convertLinearCellsToQuadratic(1)
+ self.assertAlmostEqual(m3.computeDiameterField().getArray()[0],4.1,12)
+ # QUAD4 - spacedim = 2
+ coo=DataArrayDouble([(0,2),(2,0),(6,4),(4,9)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_QUAD4) ; m.setCoords(coo)
+ exp3=sqrt(85.)
+ for delta in xrange(4):
+ c=[(elt+delta)%4 for elt in xrange(4)]
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp3,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp3,12)
+ m3=m.buildUnstructured() ; m3.convertLinearCellsToQuadratic(1)
+ self.assertAlmostEqual(m3.computeDiameterField().getArray()[0],exp3,12)
+ c.reverse()
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp3,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp3,12)
+ m3=m.buildUnstructured() ; m3.convertLinearCellsToQuadratic(1)
+ self.assertAlmostEqual(m3.computeDiameterField().getArray()[0],exp3,12)
+ # QUAD4 - spacedim = 3
+ coo=DataArrayDouble([(0.26570992384234871,2.0405889913271817,-0.079134238105786903),(2.3739976619218064,0.15779148692781009,0.021842842914139737),(6.1207841448393197,4.3755532938679655,0.43666375769970678),(3.8363255342943359,9.2521096041694229,0.41551170895942313)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_QUAD4) ; m.setCoords(coo)
+ for delta in xrange(4):
+ c=[(elt+delta)%4 for elt in xrange(4)]
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp3,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp3,12)
+ m3=m.buildUnstructured() ; m3.convertLinearCellsToQuadratic(1)
+ self.assertAlmostEqual(m3.computeDiameterField().getArray()[0],exp3,12)
+ c.reverse()
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp3,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp3,12)
+ m3=m.buildUnstructured() ; m3.convertLinearCellsToQuadratic(1)
+ self.assertAlmostEqual(m3.computeDiameterField().getArray()[0],exp3,12)
+ # PENTA6
+ # noise of coo=DataArrayDouble([(0,0,0),(1,0,0),(0,1,0),(0,0,2),(1,0,2),(0,1,2)]) + rotation([0.7,-1.2,0.6],[-4,-1,10],0.3)
+ coo=DataArrayDouble([(-0.28594726851554486,-0.23715005500928255,-0.10268080010083136),(0.6167364988633947,-0.008923258436324799,-0.08574087516687756),(-0.6132873463333834,0.6943403970881654,-0.2806118260037991),(-0.40705974936532896,-0.05868487929989308,1.7724055544436323),(0.5505955507861958,0.19145393798144705,1.8788156352163994),(-0.6092686217773406,0.812502961290914,1.685712743757831)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_PENTA6) ; m.setCoords(coo)
+ exp4=2.5041256256889888
+ self.assertAlmostEqual(exp4,coo.buildEuclidianDistanceDenseMatrix().getMaxValue()[0],12)# <- the definition of diameter
+ for delta in xrange(3):
+ c=[(elt+delta)%3 for elt in xrange(3)]
+ c+=[elt+3 for elt in c]
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp4,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp4,12)
+ c.reverse()
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp4,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp4,12)
+ # HEXA8
+ # noise of coo=DataArrayDouble([(0,0,0),(1,0,0),(1,1,0),(0,1,0),(0,0,2),(1,0,2),(1,1,2),(0,1,2)]) + rotation([0.7,-1.2,0.6],[-4,-1,10],0.3)
+ coo=DataArrayDouble([(-0.21266406388867243,-0.3049569460042527,-0.11012394815006032),(0.7641037943272584,-0.06990814759929553,-0.0909613877456491),(0.47406560768559974,0.8681310650341907,-0.2577311403703061),(-0.5136830410871793,0.644390554940524,-0.21319015989794698),(-0.4080167737381202,-0.12853761670628505,1.7869166291979348),(0.5650318811550441,0.20476257733110748,1.8140158890821603),(0.3230844436386215,1.1660778242678538,1.7175073141333406),(-0.6656588358432984,0.918357550969698,1.7566470691880265)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_HEXA8) ; m.setCoords(coo)
+ exp5=2.5366409441884215
+ self.assertAlmostEqual(exp5,coo.buildEuclidianDistanceDenseMatrix().getMaxValue()[0],12)# <- the definition of diameter
+ for delta in xrange(4):
+ c=[(elt+delta)%4 for elt in xrange(4)]
+ c+=[elt+4 for elt in c]
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp5,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp5,12)
+ c.reverse()
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp5,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp5,12)
+ # PYRA5 (1) 5th node is further
+ # noise of coo=DataArrayDouble([(0,0,0),(1,0,0),(1,1,0),(0,1,0),(0.5,0.5,2)]) + rotation([0.7,-1.2,0.6],[-4,-1,10],0.3)
+ coo=DataArrayDouble([(-0.31638393672228626,-0.3157865246451914,-0.12555467233075002),(0.7281379795666488,0.03836511217237115,-0.08431662762197323),(0.4757967840735147,0.8798897996143908,-0.2680890320119049),(-0.5386339871809047,0.5933159894201252,-0.2975311238319419),(0.012042592988768974,0.534282135495012,1.7859521682027926)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_PYRA5) ; m.setCoords(coo)
+ exp6=2.1558368027391386
+ self.assertAlmostEqual(exp6,coo.buildEuclidianDistanceDenseMatrix().getMaxValue()[0],12)# <- the definition of diameter
+ for delta in xrange(4):
+ c=[(elt+delta)%4 for elt in xrange(4)]
+ c+=[4]
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp6,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp6,12)
+ pass
+ # PYRA5 (2) 5th node is closer
+ # noise of coo=DataArrayDouble([(0,0,0),(1,0,0),(1,1,0),(0,1,0),(0.5,0.5,0.1)]) + rotation([0.7,-1.2,0.6],[-4,-1,10],0.3)
+ coo=DataArrayDouble([(-0.31638393672228626,-0.3157865246451914,-0.12555467233075002),(0.7281379795666488,0.03836511217237115,-0.08431662762197323),(0.4757967840735147,0.8798897996143908,-0.2680890320119049),(-0.5386339871809047,0.5933159894201252,-0.2975311238319419),(0.092964408350795,0.33389670321297005,-0.10171764888060142)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_PYRA5) ; m.setCoords(coo)
+ exp7=1.4413563787228953
+ self.assertAlmostEqual(exp7,coo.buildEuclidianDistanceDenseMatrix().getMaxValue()[0],12)# <- the definition of diameter
+ for delta in xrange(4):
+ c=[(elt+delta)%4 for elt in xrange(4)]
+ c+=[4]
+ m.setNodalConnectivity(DataArrayInt(c))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp7,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp7,12)
+ pass
+ # TETRA4
+ # noise of coo=DataArrayDouble([(0,0,0),(1,0,0),(0,1,0),(1,1,1)]) + rotation([0.7,-1.2,0.6],[-4,-1,10],0.3)
+ coo=DataArrayDouble([(-0.2256894071281369,-0.27631691290428106,-0.20266086543995965),(0.655458695100186,-0.08173323565551605,-0.19254662462061933),(-0.49893490718947264,0.5848097154568599,-0.3039928255382145),(0.2988102920828487,1.0582266398878504,0.7347375047372364)])
+ m=MEDCoupling1SGTUMesh("mesh",NORM_TETRA4) ; m.setCoords(coo)
+ exp8=1.7131322579364157
+ self.assertAlmostEqual(exp8,coo.buildEuclidianDistanceDenseMatrix().getMaxValue()[0],12)# <- the definition of diameter
+ for c in [[0,1,2,3],[0,3,2,1],[0,1,3,2],[0,2,3,1],[0,3,1,2],[0,2,1,3]]:
+ for i in xrange(4):
+ m.setNodalConnectivity(DataArrayInt([(elt+i)%4 for elt in c]))
+ self.assertAlmostEqual(m.computeDiameterField().getArray()[0],exp8,12)
+ m2=m.buildUnstructured() ; m2.convertLinearCellsToQuadratic(0)
+ self.assertAlmostEqual(m2.computeDiameterField().getArray()[0],exp8,12)
+ pass
+ pass
+ pass
+
+ def testMEDCouplingSkyLineArray(self):
+ index = DataArrayInt([ 0, 3, 5, 6, 6 ])
+ value = DataArrayInt([ 1, 2, 3, 2, 3, 3 ])
+
+ sla0 = MEDCouplingSkyLineArray()
+ self.assertEqual( -1, sla0.getNumberOf() )
+ self.assertEqual( 0, sla0.getLength() )
+ sla0.set( index, value )
+ self.assertTrue( index.isEqual( sla0.getIndexArray() ))
+ self.assertTrue( value.isEqual( sla0.getValueArray() ))
+ self.assertEqual( 4, sla0.getNumberOf() )
+ self.assertEqual( 6, sla0.getLength() )
+
+ sla1 = MEDCouplingSkyLineArray( index, value )
+ self.assertTrue( index.isEqual( sla1.getIndexArray() ))
+ self.assertTrue( value.isEqual( sla1.getValueArray() ))
+ self.assertEqual( 4, sla1.getNumberOf() )
+ self.assertEqual( 6, sla1.getLength() )
+
+ sla2 = MEDCouplingSkyLineArray( sla1 )
+ self.assertTrue( index.isEqual( sla2.getIndexArray() ))
+ self.assertTrue( value.isEqual( sla2.getValueArray() ))
+ self.assertEqual( 4, sla2.getNumberOf() )
+ self.assertEqual( 6, sla2.getLength() )
+
+ indexVec = ivec(); indexVec.reserve( len( index ))
+ for i in index: indexVec.push_back( i[0] )
+ valueVec = ivec(); valueVec.reserve( len( value ))
+ for i in value: valueVec.push_back( i[0] )
+ sla3 = MEDCouplingSkyLineArray( indexVec, valueVec )
+ self.assertTrue( index.isEqual( sla3.getIndexArray() ))
+ self.assertTrue( value.isEqual( sla3.getValueArray() ))
+ self.assertEqual( 4, sla3.getNumberOf() )
+ self.assertEqual( 6, sla3.getLength() )
+
+ pass
+
+ def testMEDCouplingUMeshgenerateGraph(self):
+ # cartesian mesh 3x3
+ arr=DataArrayDouble(4) ; arr.iota()
+ c=MEDCouplingCMesh() ; c.setCoords(arr,arr)
+ m=c.buildUnstructured()
+ graph = m.generateGraph()
+ # 0 1 2
+ # 3 4 5
+ # 6 7 8
+ valRef=[ 0,1,3,
+ 0,1,2,4,
+ 1,2,5,
+ 0,3,4,6,
+ 1,3,4,5,7,
+ 2,4,5,8,
+ 3,6,7,
+ 4,6,7,8,
+ 5,7,8]
+ self.assertEqual(valRef,list(graph.getValueArray().getValues()));
+
+ indRef=[0, 3, 7, 10, 14, 19, 23, 26, 30, 33]
+ self.assertEqual(indRef,list(graph.getIndexArray().getValues()));
+ pass
+
+ def testSwig2MEDCouplingCurveLinearReprQuick1(self):
+ """Non regression test. Error in m.__str__ when m is a MEDCouplingCurveLinear with spaceDim != meshDim."""
+ arr=DataArrayDouble(12) ; arr.iota() ; arr.rearrange(2)
+ m=MEDCouplingCurveLinearMesh()
+ m.setCoords(arr)
+ m.setNodeGridStructure([3,2])
+ m.checkCoherency()
+ self.assertEqual(m.getMeshDimension(),2)
+ self.assertEqual(m.getSpaceDimension(),2)
+ self.assertTrue(not "mismatch" in m.__str__())
+ self.assertTrue(not "mismatch" in m.__repr__())
+ #
+ arr=DataArrayDouble(18) ; arr.iota() ; arr.rearrange(3)
+ m.setCoords(arr)
+ self.assertEqual(m.getMeshDimension(),2)
+ self.assertEqual(m.getSpaceDimension(),3)
+ self.assertTrue(not "mismatch" in m.__str__())
+ self.assertTrue(not "mismatch" in m.__repr__())# bug was here !
+ pass
+
pass
if __name__ == '__main__':
