X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FMEDCoupling_Swig%2FMEDCouplingBasicsTest1.py;h=4d3692134593b76200104c5e8e4b4b4f04422286;hb=5b299d506fcc63ff652df59a13f5b082fd45dabf;hp=ad43c4491b56e8694ea682057c850e2fc8caa096;hpb=1ff18bfcccb60eb7954f23a427024379f70eb3a9;p=tools%2Fmedcoupling.git diff --git a/src/MEDCoupling_Swig/MEDCouplingBasicsTest1.py b/src/MEDCoupling_Swig/MEDCouplingBasicsTest1.py index ad43c4491..4d3692134 100644 --- a/src/MEDCoupling_Swig/MEDCouplingBasicsTest1.py +++ b/src/MEDCoupling_Swig/MEDCouplingBasicsTest1.py @@ -1,5 +1,5 @@ # -*- coding: utf-8 -*- -# Copyright (C) 2007-2015 CEA/DEN, EDF R&D +# Copyright (C) 2007-2016 CEA/DEN, EDF R&D # # This library is free software; you can redistribute it and/or # modify it under the terms of the GNU Lesser General Public @@ -23,7 +23,7 @@ import unittest from math import pi,e,sqrt,cos,sin from datetime import datetime from MEDCouplingDataForTest import MEDCouplingDataForTest -import rlcompleter,readline # this line has to be here, to ensure a usability of MEDCoupling/MEDLoader. B4 removing it please notify to anthony.geay@cea.fr +import rlcompleter,readline # this line has to be here, to ensure a usability of MEDCoupling/MEDLoader. B4 removing it please notify to anthony.geay@edf.fr class MEDCouplingBasicsTest1(unittest.TestCase): def testArray2(self): @@ -45,11 +45,11 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(7,arr1.getNumberOfTuples()); self.assertEqual(2,arr1.getNumberOfComponents()); self.assertEqual(arr1Ref,list(arr1.getValues())); - arr2=arr1.substr(3); + arr2=arr1.subArray(3); self.assertEqual(4,arr2.getNumberOfTuples()); self.assertEqual(2,arr2.getNumberOfComponents()); self.assertEqual(arr1Ref[6:],list(arr2.getValues())); - arr3=arr1.substr(2,5); + arr3=arr1.subArray(2,5); self.assertEqual(3,arr3.getNumberOfTuples()); self.assertEqual(2,arr3.getNumberOfComponents()); self.assertEqual(arr1Ref[4:10],list(arr3.getValues())); @@ -60,21 +60,21 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(7,arr4.getNumberOfTuples()); self.assertEqual(2,arr4.getNumberOfComponents()); tmp=arr4.getValues() - for i in xrange(14): + for i in range(14): self.assertTrue(abs(arr4Ref[i]-tmp[i])<1e-14); pass - arr5=arr4.substr(3); + arr5=arr4.subArray(3); self.assertEqual(4,arr5.getNumberOfTuples()); self.assertEqual(2,arr5.getNumberOfComponents()); tmp=arr5.getValues() - for i in xrange(8): + for i in range(8): self.assertTrue(abs(arr4Ref[6+i]-tmp[i])<1e-14); pass - arr6=arr4.substr(2,5); + arr6=arr4.subArray(2,5); self.assertEqual(3,arr6.getNumberOfTuples()); self.assertEqual(2,arr6.getNumberOfComponents()); tmp=arr6.getValues() - for i in xrange(6): + for i in range(6): self.assertTrue(abs(arr4Ref[4+i]-tmp[i])<1e-14); pass pass @@ -122,7 +122,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): myCoords.setValues(coords,nbOfNodes,3); self.assertTrue(myCoords.getIJ(3,2)==-0.305) mesh.setCoords(myCoords); - mesh.checkCoherency(); + mesh.checkConsistencyLight(); self.assertTrue(mesh.getAllGeoTypes()==[4]) myFalseConn=DataArrayInt.New() myFalseConn.setValues(tab4,6,4) @@ -130,15 +130,15 @@ class MEDCouplingBasicsTest1(unittest.TestCase): # field=MEDCouplingFieldDouble.New(ON_CELLS) field.setMesh(mesh) - field.setNature(Integral) + field.setNature(ExtensiveMaximum) myCoords=DataArrayDouble.New() sampleTab=[] - for i in range(nbOfCells*9): + for i in range(nbOfCells * 9): sampleTab.append(float(i)) myCoords.setValues(sampleTab,nbOfCells,9); field.setArray(myCoords) self.assertTrue(3==mesh.getSpaceDimension()) - field.checkCoherency() + field.checkConsistencyLight() mesh2=mesh.clone(False) mesh3=mesh.clone(True) mesh3=0 @@ -166,7 +166,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): targetMesh.insertNextCell(NORM_POINT1,1,[7]); targetMesh.insertNextCell(NORM_POINT1,1,[6]); targetMesh.finishInsertingCells(); - self.assertRaises(InterpKernelException,targetMesh.checkCoherency); + self.assertRaises(InterpKernelException,targetMesh.checkConsistencyLight); myCoords=DataArrayDouble.New(); myCoords.setValues(targetCoords,9,3); targetMesh.setCoords(myCoords); @@ -184,7 +184,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertRaises(InterpKernelException,meshM1D.setMeshDimension,-2) self.assertRaises(InterpKernelException,meshM1D.setMeshDimension,-10) meshM1D.setMeshDimension(-1); - meshM1D.checkCoherency(); + meshM1D.checkConsistencyLight(); self.assertEqual(meshM1D.getMeshDimension(),-1); self.assertEqual(meshM1D.getNumberOfCells(),1); self.assertRaises(InterpKernelException,meshM1D.getNumberOfNodes); @@ -196,20 +196,20 @@ class MEDCouplingBasicsTest1(unittest.TestCase): array=DataArrayDouble.New(); array.setValues(6*[7.],1,6); fieldOnCells.setArray(array); - fieldOnCells.checkCoherency(); + fieldOnCells.checkConsistencyLight(); pass def testDeepCopy(self): array=DataArrayDouble.New(); array.setValues(5*3*[7.],5,3); self.assertEqual(array.getIJ(3,2),7.); - array2=array.deepCpy(); + array2=array.deepCopy(); self.assertEqual(array2.getIJ(3,2),7.) # array3=DataArrayInt.New(); array3.setValues(5*3*[17],5,3); self.assertEqual(array3.getIJ(3,2),17); - array4=array3.deepCpy(); + array4=array3.deepCopy(); self.assertEqual(array4.getIJ(3,2),17); pass @@ -228,7 +228,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): mesh=MEDCouplingDataForTest.build2DTargetMesh_1(); elts=[1,3]; mesh.convertToPolyTypes(elts); - mesh.checkCoherency(); + mesh.checkConsistencyLight(); self.assertEqual(5,mesh.getNumberOfCells()); self.assertEqual(23,mesh.getNodalConnectivity().getNumberOfTuples()); expected1=[4, 0, 3, 4, 1, 5, 1, 4, 2, 3, 4, 5, 2, 5, 6, 7, 4, 3, 4, 7, 8, 5, 4] @@ -236,11 +236,11 @@ class MEDCouplingBasicsTest1(unittest.TestCase): # mesh=MEDCouplingDataForTest.build3DTargetMesh_1(); mesh.convertToPolyTypes(elts); - mesh.checkCoherency(); + mesh.checkConsistencyLight(); self.assertEqual(8,mesh.getNumberOfCells()); self.assertEqual(114,mesh.getNodalConnectivity().getNumberOfTuples()); mesh.convertToPolyTypes(elts); - mesh.checkCoherency(); + mesh.checkConsistencyLight(); self.assertEqual(8,mesh.getNumberOfCells()); self.assertEqual(114,mesh.getNodalConnectivity().getNumberOfTuples()); pass @@ -252,7 +252,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): revDesc=DataArrayInt.New(); revDescIndx=DataArrayInt.New(); mesh2=mesh.buildDescendingConnectivity(desc,descIndx,revDesc,revDescIndx); - mesh2.checkCoherency(); + mesh2.checkConsistencyLight(); self.assertEqual(1,mesh2.getMeshDimension()); self.assertEqual(13,mesh2.getNumberOfCells()); self.assertEqual(14,revDescIndx.getNbOfElems()); self.assertEqual(14,revDescIndx.getNumberOfTuples()); @@ -276,7 +276,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): # eltsV=[1,3]; mesh.convertToPolyTypes(eltsV); - mesh.checkCoherency(); + mesh.checkConsistencyLight(); # desc=DataArrayInt.New(); descIndx=DataArrayInt.New(); @@ -284,7 +284,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): revDescIndx=DataArrayInt.New(); # mesh2=mesh.buildDescendingConnectivity(desc,descIndx,revDesc,revDescIndx); - mesh2.checkCoherency(); + mesh2.checkConsistencyLight(); self.assertEqual(1,mesh2.getMeshDimension()); self.assertEqual(13,mesh2.getNumberOfCells()); self.assertEqual(14,revDescIndx.getNbOfElems()); self.assertEqual(14,revDescIndx.getNumberOfTuples()); @@ -309,7 +309,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): revDescIndx=DataArrayInt.New(); # mesh2=mesh.buildDescendingConnectivity(desc,descIndx,revDesc,revDescIndx); - mesh2.checkCoherency(); + mesh2.checkConsistencyLight(); self.assertEqual(2,mesh2.getMeshDimension()); self.assertEqual(36,mesh2.getNumberOfCells()); self.assertEqual(37,revDescIndx.getNbOfElems()); self.assertEqual(37,revDescIndx.getNumberOfTuples()); @@ -339,13 +339,13 @@ class MEDCouplingBasicsTest1(unittest.TestCase): # eltsV=[1,3] mesh.convertToPolyTypes(eltsV); - mesh.checkCoherency(); + mesh.checkConsistencyLight(); desc=DataArrayInt.New(); descIndx=DataArrayInt.New(); revDesc=DataArrayInt.New(); revDescIndx=DataArrayInt.New(); mesh2=mesh.buildDescendingConnectivity(desc,descIndx,revDesc,revDescIndx); - mesh2.checkCoherency(); + mesh2.checkConsistencyLight(); self.assertEqual(2,mesh2.getMeshDimension()); self.assertEqual(36,mesh2.getNumberOfCells()); self.assertEqual(37,revDescIndx.getNbOfElems()); self.assertEqual(37,revDescIndx.getNumberOfTuples()); @@ -634,7 +634,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertTrue(fieldOnCells1.isEqual(fieldOnCells2,1e-12,1e-15)); self.assertTrue(fieldOnCells2.isEqual(fieldOnCells1,1e-12,1e-15)); # - arr2=arr.deepCpy(); + arr2=arr.deepCopy(); fieldOnCells2.setArray(arr2); self.assertTrue(fieldOnCells1.isEqual(fieldOnCells2,1e-12,1e-15)); self.assertTrue(fieldOnCells2.isEqual(fieldOnCells1,1e-12,1e-15)); @@ -662,13 +662,13 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testNatureChecking(self): field=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME); - field.setNature(Integral); - field.setNature(ConservativeVolumic); - field.setNature(IntegralGlobConstraint); + field.setNature(ExtensiveMaximum); + field.setNature(IntensiveMaximum); + field.setNature(ExtensiveConservation); field=MEDCouplingFieldDouble.New(ON_NODES,NO_TIME); - field.setNature(ConservativeVolumic); - self.assertRaises(InterpKernelException,field.setNature,Integral); - self.assertRaises(InterpKernelException,field.setNature,IntegralGlobConstraint); + field.setNature(IntensiveMaximum); + self.assertRaises(InterpKernelException,field.setNature,ExtensiveMaximum); + self.assertRaises(InterpKernelException,field.setNature,ExtensiveConservation); pass def testNatureOperations(self): @@ -678,14 +678,14 @@ class MEDCouplingBasicsTest1(unittest.TestCase): 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()) + f1.setNature(ExtensiveMaximum) + f2.setNature(IntensiveMaximum) + self.assertEqual(ExtensiveMaximum, f1.getNature()) + self.assertEqual(IntensiveMaximum, 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()) + f1.setArray(da); f2.setArray(da.deepCopy()) # All this should complain about nature: self.assertRaises(InterpKernelException, f1.__add__, f2) self.assertRaises(InterpKernelException, f1.__iadd__, f2) @@ -702,19 +702,19 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(NoNature, f3.getNature()) f3 = f1*f2 self.assertEqual(NoNature, f3.getNature()) - f1Tmp = f1.deepCpy(); f1Tmp.setMesh(m); f1Tmp *= f2 + f1Tmp = f1.deepCopy(); 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 + f1Tmp = f1.deepCopy(); 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 + f1Tmp = f1.deepCopy(); f1Tmp.setMesh(m); f1Tmp **= f2 self.assertEqual(NoNature, f1Tmp.getNature()) f3 = MEDCouplingFieldDouble.DotFields(f1,f2) self.assertEqual(NoNature, f3.getNature()) @@ -722,7 +722,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(NoNature, f3.getNature()) da = DataArrayDouble.Meld([da, da, da]) - f1.setArray(da); f2.setArray(da.deepCpy()) + f1.setArray(da); f2.setArray(da.deepCopy()) f3 = MEDCouplingFieldDouble.CrossProductFields(f1,f2) self.assertEqual(NoNature, f3.getNature()) f3 = f1.crossProduct(f2) @@ -757,7 +757,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testExtrudedMesh1(self): mesh3D,mesh2D=MEDCouplingDataForTest.build3DExtrudedUMesh_1(); - ext=MEDCouplingExtrudedMesh.New(mesh3D,mesh2D,1); + ext=MEDCouplingMappedExtrudedMesh.New(mesh3D,mesh2D,1); self.assertEqual(18,ext.getNumberOfCells()); self.assertEqual(60,ext.getNumberOfNodes()); ids3D=ext.getMesh3DIds(); @@ -792,16 +792,16 @@ class MEDCouplingBasicsTest1(unittest.TestCase): m2.rotate(center,vector,-pi/2.); m3=m1.buildExtrudedMesh(m2,0); # - m4=MEDCouplingExtrudedMesh.New(m3,m1,0); + m4=MEDCouplingMappedExtrudedMesh.New(m3,m1,0); self.assertEqual(15,m4.getNumberOfCells()); self.assertEqual(5,m4.getMesh2D().getNumberOfCells()); self.assertEqual(3,m4.getMesh1D().getNumberOfCells()); m3DIds=m4.getMesh3DIds().getValues(); - self.assertEqual(range(15),list(m3DIds)); + self.assertEqual(list(range(15)), list(m3DIds)); #some random in cells to check that extrusion alg find it correctly expected1=[1,3,2,0,6,5,7,10,11,8,12,9,14,13,4] m3.renumberCells(expected1,False); - m4=MEDCouplingExtrudedMesh.New(m3,m1,0); + m4=MEDCouplingMappedExtrudedMesh.New(m3,m1,0); self.assertEqual(15,m4.getNumberOfCells()); self.assertEqual(5,m4.getMesh2D().getNumberOfCells()); self.assertEqual(3,m4.getMesh1D().getNumberOfCells()); @@ -817,7 +817,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(NORM_POLYHED,m3.getTypeOfCell(3)); self.assertEqual(NORM_HEXA8,m3.getTypeOfCell(4)); m3.renumberCells(expected1,False); - m4=MEDCouplingExtrudedMesh.New(m3,m1,0); + m4=MEDCouplingMappedExtrudedMesh.New(m3,m1,0); self.assertEqual(15,m4.getNumberOfCells()); self.assertEqual(5,m4.getMesh2D().getNumberOfCells()); self.assertEqual(3,m4.getMesh1D().getNumberOfCells()); @@ -840,7 +840,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): expected1=[1,3,2,0,6,5,7,10,11,8,12,9,14,13,4] rexpected1=[3, 0, 2, 1, 14, 5, 4, 6, 9, 11, 7, 8, 10, 13, 12] m3.renumberCells(expected1,False); - m4=MEDCouplingExtrudedMesh.New(m3,m1,0); + m4=MEDCouplingMappedExtrudedMesh.New(m3,m1,0); self.assertEqual(NORM_HEXA8,m4.getTypeOfCell(0)); self.assertEqual(NORM_HEXA8,m4.getTypeOfCell(1)); self.assertEqual(NORM_POLYHED,m4.getTypeOfCell(2)); @@ -853,7 +853,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): expected2=[0.075,0.0375,0.0375,0.075,0.075, 0.1125,0.05625,0.05625,0.1125,0.1125, 0.0625,0.03125,0.03125,0.0625,0.0625] - for i in xrange(15): + for i in range(15): self.assertAlmostEqual(expected2[rexpected1[i]],arrPtr[i],16); pass m5=m4.build3DUnstructuredMesh(); @@ -861,10 +861,10 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertTrue(m5.isEqual(m3,1e-12)); f=m5.getMeasureField(True); f.setMesh(m4) - self.assertTrue(isinstance(f.getMesh(),MEDCouplingExtrudedMesh)) + self.assertTrue(isinstance(f.getMesh(),MEDCouplingMappedExtrudedMesh)) arr=f.getArray(); arrPtr=arr.getValues(); - for i in xrange(15): + for i in range(15): self.assertAlmostEqual(expected2[rexpected1[i]],arrPtr[i],15); pass pass @@ -877,7 +877,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): o2n,newNbOfNodes=targetMesh.buildNewNumberingFromCommonNodesFormat(comm,commI); self.assertEqual(27,newNbOfNodes); self.assertEqual(27,o2n.getNumberOfTuples()); - o2nExp1=range(27) + o2nExp1 = list(range(27)) self.assertEqual(o2nExp1,list(o2n.getValues())); # targetMesh=MEDCouplingDataForTest.build3DTargetMeshMergeNode_1(); @@ -981,7 +981,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): m2.translate(vec); m3=m1.mergeMyselfWith(m2); self.assertTrue(isinstance(m3,MEDCouplingUMesh)); - m3.checkCoherency(); + m3.checkConsistencyLight(); m4=MEDCouplingDataForTest.build2DTargetMeshMerged_1(); self.assertTrue(m3.isEqual(m4,1.e-12)); da,isMerged,newNbOfNodes=m3.mergeNodes(1.e-12); @@ -992,14 +992,14 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testMergeMeshOnSameCoords1(self): m1=MEDCouplingDataForTest.build2DTargetMesh_1(); m2=MEDCouplingDataForTest.build2DTargetMesh_1(); - cells=range(5); + cells = list(range(5)); m2.convertToPolyTypes(cells); m1.tryToShareSameCoords(m2,1e-12); m3=MEDCouplingDataForTest.build2DTargetMesh_1(); m3.tryToShareSameCoords(m2,1e-12); meshes=[m1,m2,m3] m4=MEDCouplingUMesh.MergeUMeshesOnSameCoords(meshes); - m4.checkCoherency(); + m4.checkConsistencyLight(); self.assertEqual(15,m4.getNumberOfCells()); cells1=[0,1,2,3,4] m1_1=m4.buildPartOfMySelf(cells1,True); @@ -1023,7 +1023,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f1=m1.getMeasureField(True); f2=m2.getMeasureField(True); f3=MEDCouplingFieldDouble.MergeFields(f1,f2); - f3.checkCoherency(); + f3.checkConsistencyLight(); m4=MEDCouplingDataForTest.build2DTargetMeshMerged_1(); self.assertTrue(f3.getMesh().isEqual(m4,1.e-12)); name=f3.getName(); @@ -1035,7 +1035,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values=[0.25,0.125,0.125,0.25,0.25,0.5,0.5] tmp=f3.getArray().getValues(); self.assertEqual(len(values),len(tmp)) - for i in xrange(7): + for i in range(7): self.assertTrue(abs(values[i]-tmp[i])<1e-12) pass pass @@ -1046,7 +1046,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f1=m.fillFromAnalytic(ON_CELLS,1,"x+y"); self.assertAlmostEqual(3.4,f1.getTime()[0],12) ; self.assertEqual(5,f1.getTime()[1]) ; self.assertEqual(6,f1.getTime()[2]) self.assertEqual("us",f1.getTimeUnit()) - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_CELLS); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(1,f1.getNumberOfComponents()); @@ -1054,12 +1054,12 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values1=[-0.1,0.23333333333333336,0.56666666666666665,0.4,0.9] tmp=f1.getArray().getValues(); self.assertEqual(len(values1),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values1[i])<1.e-12) pass # f1=m.fillFromAnalytic(ON_NODES,1,"x+y"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(1,f1.getNumberOfComponents()); @@ -1067,12 +1067,12 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values2=[-0.6,-0.1,0.4,-0.1,0.4,0.9,0.4,0.9,1.4] tmp=f1.getArray().getValues(); self.assertEqual(len(values2),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values2[i])<1.e-12) pass # f1=m.fillFromAnalytic(ON_NODES,2,"(x+y)*IVec+(2*(x+y))*JVec"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(2,f1.getNumberOfComponents()); @@ -1080,7 +1080,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values3=[-0.6,-1.2,-0.1,-0.2,0.4,0.8,-0.1,-0.2,0.4,0.8,0.9,1.8,0.4,0.8,0.9,1.8,1.4,2.8] tmp=f1.getArray().getValues(); self.assertEqual(len(values3),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values3[i])<1.e-12) pass values4=f1.accumulate(); @@ -1098,7 +1098,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testFillFromAnalytic2(self): m=MEDCouplingDataForTest.build2DTargetMesh_1(); f1=m.fillFromAnalytic(ON_CELLS,1,"y+x"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_CELLS); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(1,f1.getNumberOfComponents()); @@ -1106,12 +1106,12 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values1=[-0.1,0.23333333333333336,0.56666666666666665,0.4,0.9] tmp=f1.getArray().getValues(); self.assertEqual(len(values1),len(tmp)) - for i in xrange(len(values1)): + for i in range(len(values1)): self.assertTrue(abs(values1[i]-tmp[i])<1.e-12); pass # f1=m.fillFromAnalytic(ON_NODES,1,"y+2*x"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(1,f1.getNumberOfComponents()); @@ -1119,11 +1119,11 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values2=[-0.9,0.1,1.1,-0.4,0.6,1.6,0.1,1.1,2.1] tmp=f1.getArray().getValues(); self.assertEqual(len(values2),len(tmp)) - for i in xrange(len(values2)): + for i in range(len(values2)): self.assertTrue(abs(values2[i]-tmp[i])<1.e-12); pass f1=m.fillFromAnalytic(ON_NODES,1,"2.*x+y"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(1,f1.getNumberOfComponents()); @@ -1131,12 +1131,12 @@ class MEDCouplingBasicsTest1(unittest.TestCase): tmp=f1.getArray().getValues(); values2Bis=[-0.9,0.1,1.1,-0.4,0.6,1.6,0.1,1.1,2.1] self.assertEqual(len(values2Bis),len(tmp)) - for i in xrange(len(values2Bis)): + for i in range(len(values2Bis)): self.assertTrue(abs(values2Bis[i]-tmp[i])<1.e-12); pass # f1=m.fillFromAnalytic(ON_NODES,2,"(x+y)*IVec+2*(x+y)*JVec"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(2,f1.getNumberOfComponents()); @@ -1144,7 +1144,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values3=[-0.6,-1.2,-0.1,-0.2,0.4,0.8,-0.1,-0.2,0.4,0.8,0.9,1.8,0.4,0.8,0.9,1.8,1.4,2.8] tmp=f1.getArray().getValues(); self.assertEqual(len(values3),len(tmp)) - for i in xrange(len(values3)): + for i in range(len(values3)): self.assertTrue(abs(values3[i]-tmp[i])<1.e-12); pass values4=f1.accumulate(); @@ -1158,7 +1158,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testApplyFunc(self): m=MEDCouplingDataForTest.build2DTargetMesh_1(); f1=m.fillFromAnalytic(ON_NODES,2,"(x+y)*IVec+(2*(x+y))*JVec"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(2,f1.getNumberOfComponents()); @@ -1171,7 +1171,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values1=[-1.8,-0.3,1.2,-0.3,1.2,2.7,1.2,2.7,4.2] tmp=f1.getArray().getValues(); self.assertEqual(len(values1),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values1[i])<1.e-12) pass pass @@ -1179,7 +1179,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testApplyFunc2(self): m=MEDCouplingDataForTest.build2DTargetMesh_1(); f1=m.fillFromAnalytic(ON_NODES,2,"(x+y)*IVec+2*(x+y)*JVec"); - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); self.assertEqual(2,f1.getNumberOfComponents()); @@ -1201,7 +1201,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): 5.0423700574830965, 17.435300118916864] tmp=f2.getArray().getValues(); self.assertEqual(len(tmp),len(values2)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values2[i])<1.e-12) pass # @@ -1213,7 +1213,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values1=[-1.8,-0.3,1.2,-0.3,1.2,2.7,1.2,2.7,4.2] tmp=f1.getArray().getValues(); self.assertEqual(len(tmp),len(values1)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values1[i])<1.e-12) pass pass @@ -1222,54 +1222,54 @@ class MEDCouplingBasicsTest1(unittest.TestCase): m=MEDCouplingDataForTest.build2DTargetMesh_1(); f1=m.fillFromAnalytic(ON_NODES,1,"x+y"); f2=m.fillFromAnalytic(ON_NODES,1,"x+y"); - f1.checkCoherency(); - f2.checkCoherency(); + f1.checkConsistencyLight(); + f2.checkConsistencyLight(); f3=f1+f2; - f3.checkCoherency(); + f3.checkConsistencyLight(); self.assertEqual(f3.getTypeOfField(),ON_NODES); self.assertEqual(f3.getTimeDiscretization(),ONE_TIME); values1=[-1.2,-0.2,0.8,-0.2,0.8,1.8,0.8,1.8,2.8] tmp=f3.getArray().getValues(); self.assertEqual(len(values1),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values1[i])<1.e-12) pass # f3=f1*f2; - f3.checkCoherency(); + f3.checkConsistencyLight(); self.assertEqual(f3.getTypeOfField(),ON_NODES); self.assertEqual(f3.getTimeDiscretization(),ONE_TIME); values2=[0.36,0.01,0.16,0.01,0.16,0.81,0.16,0.81,1.96] tmp=f3.getArray().getValues(); self.assertEqual(len(values2),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values2[i])<1.e-12) pass # f3=f1+f2; f4=f1-f3; - f4.checkCoherency(); + f4.checkConsistencyLight(); self.assertEqual(f4.getTypeOfField(),ON_NODES); self.assertEqual(f4.getTimeDiscretization(),ONE_TIME); values3=[0.6,0.1,-0.4,0.1,-0.4,-0.9,-0.4,-0.9,-1.4] tmp=f4.getArray().getValues(); self.assertEqual(len(values3),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values3[i])<1.e-12) pass # f3=f1+f2; f4=f3/f2; - f4.checkCoherency(); + f4.checkConsistencyLight(); self.assertEqual(f4.getTypeOfField(),ON_NODES); self.assertEqual(f4.getTimeDiscretization(),ONE_TIME); tmp=f4.getArray().getValues(); - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-2.)<1.e-12) pass # f4=f2.buildNewTimeReprFromThis(NO_TIME,False); - f4.checkCoherency(); + f4.checkConsistencyLight(); self.assertEqual(f4.getTypeOfField(),ON_NODES); self.assertEqual(f4.getTimeDiscretization(),NO_TIME); self.assertRaises(InterpKernelException,f1.__add__,f4); @@ -1280,12 +1280,12 @@ class MEDCouplingBasicsTest1(unittest.TestCase): tmp=f3.getArray().getValues(); values4=[-1.2,-0.2,0.8,-0.2,0.8,1.8,0.8,1.8,2.8] self.assertEqual(len(values3),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values4[i])<1.e-12) pass # f4=f2.buildNewTimeReprFromThis(NO_TIME,True); - f4.checkCoherency(); + f4.checkConsistencyLight(); self.assertEqual(f4.getTypeOfField(),ON_NODES); self.assertEqual(f4.getTimeDiscretization(),NO_TIME); self.assertRaises(InterpKernelException,f1.__add__,f4); @@ -1296,7 +1296,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): tmp=f3.getArray().getValues(); values5=[-1.2,-0.2,0.8,-0.2,0.8,1.8,0.8,1.8,2.8] self.assertEqual(len(values5),len(tmp)) - for i in xrange(len(tmp)): + for i in range(len(tmp)): self.assertTrue(abs(tmp[i]-values5[i])<1.e-12) pass pass @@ -1307,7 +1307,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f1=m.fillFromAnalytic(ON_NODES,1,"x+y+z"); f2=m.fillFromAnalytic(ON_NODES,1,"a*a+b+c*c"); f3=f1/f2; - f3.checkCoherency(); + f3.checkConsistencyLight(); self.assertEqual(f3.getTypeOfField(),ON_NODES); self.assertEqual(f3.getTimeDiscretization(),ONE_TIME); expected1=[-2.4999999999999991, 1.2162162162162162, 0.77868852459016391, @@ -1316,7 +1316,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(1,f3.getNumberOfComponents()); self.assertEqual(9,f3.getNumberOfTuples()); val=f3.getArray().getValues(); - for i in xrange(9): + for i in range(9): self.assertTrue(abs(expected1[i]-val[i])<1.e-12); # f1=m.buildOrthogonalField(); @@ -1326,13 +1326,13 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f3=f1*f2; expected2=[-0.035355339059327376,0.,0.035355339059327376, 0.2592724864350674,0.,-0.2592724864350674, 0.37712361663282529,0.,-0.37712361663282529, -0.035355339059327376,0.,0.035355339059327376, 0.31819805153394637,0.,-0.31819805153394637] val=f3.getArray().getValues(); - for i in xrange(15): + for i in range(15): self.assertTrue(abs(expected2[i]-val[i])<1.e-12); pass # f3=f2*f1; val=f3.getArray().getValues(); - for i in xrange(15): + for i in range(15): self.assertTrue(abs(expected2[i]-val[i])<1.e-12); pass pass @@ -1342,7 +1342,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f1=m.fillFromAnalytic(ON_NODES,1,"x+y+z"); f2=m.fillFromAnalytic(ON_NODES,1,"a*a+b+c*c"); f1/=f2 - f1.checkCoherency(); + f1.checkConsistencyLight(); self.assertEqual(f1.getTypeOfField(),ON_NODES); self.assertEqual(f1.getTimeDiscretization(),ONE_TIME); expected1=[-2.4999999999999991, 1.2162162162162162, 0.77868852459016391, @@ -1351,7 +1351,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(1,f1.getNumberOfComponents()); self.assertEqual(9,f1.getNumberOfTuples()); val=f1.getArray().getValues(); - for i in xrange(9): + for i in range(9): self.assertTrue(abs(expected1[i]-val[i])<1.e-12); pass # @@ -1360,7 +1360,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f1*=f2 expected2=[-0.035355339059327376,0.,0.035355339059327376, 0.2592724864350674,0.,-0.2592724864350674, 0.37712361663282529,0.,-0.37712361663282529, -0.035355339059327376,0.,0.035355339059327376, 0.31819805153394637,0.,-0.31819805153394637] val=f1.getArray().getValues(); - for i in xrange(15): + for i in range(15): self.assertTrue(abs(expected2[i]-val[i])<1.e-12); pass # @@ -1383,7 +1383,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): array.setValues(arr1,nbOfCells,3); f1.setStartTime(2.,0,0); f1.setEndTime(3.,0,0); - f1.checkCoherency(); + f1.checkConsistencyLight(); pos=[0.3,-0.2] res=f1.getValueOn(pos); self.assertTrue(abs(arr1[3]-res[0])<1.e-12); @@ -1401,11 +1401,11 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f2.setArray(f1.getArray()); f2.setStartTime(2.,3,0); f2.setEndTime(4.,13,0); - self.assertRaises(InterpKernelException,f2.checkCoherency) + self.assertRaises(InterpKernelException,f2.checkConsistencyLight) array2=DataArrayDouble.New(); array2.setValues(arr2,nbOfCells,3); f2.setEndArray(array2); - f2.checkCoherency(); + f2.checkConsistencyLight(); # res=None res=f2.getValueOn(pos,3.21); @@ -1546,7 +1546,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(3,len(corr)); expectedVals1=[3,3,2] expectedVals2=[[0,1,2],[3,0,2],[3,0]] - for i in xrange(3): + for i in range(3): arr=corr[i]; self.assertEqual(1,arr.getNumberOfComponents()); nbOfVals=expectedVals1[i]; @@ -1561,7 +1561,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(1,arr2.getNumberOfComponents()); self.assertEqual(4,arr2.getNumberOfTuples()); self.assertEqual(fidExp,list(arr2.getValues())); - for i in xrange(3): + for i in range(3): nbOfVals=expectedVals1[i]; self.assertEqual(list(fidsOfGroups[i]),fidsGrp[i]); pass @@ -1595,7 +1595,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(5,field.getNumberOfTuples()); self.assertEqual(3,field.getNumberOfComponents()); vals=field.getArray().getValues(); - for i in xrange(15): + for i in range(15): self.assertTrue(abs(expected[i%3]-vals[i])<1e-12); # testing targetCoords=[0.,0.,0.,0.5,0.,0.5,1.,0.,1.,0.,1.,0.] @@ -1699,7 +1699,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): fieldOnCells.setMesh(targetMesh); array=DataArrayDouble.New(); tmp=2*nbOfCells*[None] - for i in xrange(nbOfCells): + for i in range(nbOfCells): tmp[2*i]=7.+float(i); tmp[2*i+1]=17.+float(i) pass @@ -1719,7 +1719,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): fieldOnNodes.setMesh(targetMesh); array=DataArrayDouble.New(); tmp=2*nbOfNodes*[None] - for i in xrange(nbOfNodes): + for i in range(nbOfNodes): tmp[2*i]=17.+float(i); tmp[2*i+1]=27.+float(i) pass @@ -1765,7 +1765,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): 7., 4., 6., 7., 9., 5., 7., 8., 10., 7., 9., 10., 12.]; val=fieldOnNodes.getArray().getValues(); - for i in xrange(64): + for i in range(64): self.assertAlmostEqual(expected1[i], val[i], 12) res=fieldOnNodes.getValueOnPos(1, 3, 2); self.assertAlmostEqual(7., res[0], 12); @@ -1776,12 +1776,12 @@ class MEDCouplingBasicsTest1(unittest.TestCase): val=fieldOnCells.getArray().getValues(); expected2=[0, 1.5, 3, 1.5, 3, 4.5, 3, 4.5, 6, 1.5, 3, 4.5, 3, 4.5, 6, 4.5, 6, 7.5, 3, 4.5, 6, 4.5, 6, 7.5, 6, 7.5, 9]; - for i in xrange(27): + for i in range(27): self.assertAlmostEqual(expected2[i], val[i], 12); #res=fieldOnCells.getValueOnPos(1,2,1); #self.assertAlmostEqual(6.,res,12); # - meshDeepCopy=mesh.deepCpy(); + meshDeepCopy=mesh.deepCopy(); meshClone=mesh.clone(False); meshEmpty.copyTinyStringsFrom(mesh); @@ -1845,8 +1845,8 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertTrue(mesh1.isEqual(mesh2, 1e-5)); self.assertTrue(not mesh1.isEqual(mesh2, 1e-7)); - self.assertRaises(InterpKernelException, mesh3.checkCoherency1, 1e-12); - mesh1.checkCoherency2(1e-12); + self.assertRaises(InterpKernelException, mesh3.checkConsistency, 1e-12); + mesh1.checkConsistency(1e-12); self.assertEqual(NORM_HEXA8, mesh1.getTypeOfCell(1)); self.assertEqual(NORM_HEXA8, mesh1.getAllGeoTypes()[0]); @@ -1962,7 +1962,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): -0.05,0.45, 0.2,0.45, 0.45,0.45] val=mesh.getCoords().getValues(); self.assertEqual(18,len(val)) - for i in xrange(18): + for i in range(18): self.assertTrue(abs(expected1[i]-val[i])<1e-12); pass pass @@ -2002,7 +2002,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(9,len(n)); m3dSurf=mesh.buildFacePartOfMySelfNode(n,True); self.assertTrue(isinstance(m3dSurf,MEDCouplingUMesh)) - me=MEDCouplingExtrudedMesh.New(mesh,m3dSurf,0); + me=MEDCouplingMappedExtrudedMesh.New(mesh,m3dSurf,0); da=me.getMesh3DIds(); self.assertEqual(8,me.getNumberOfCells()); expected=[0,1,2,3,4,5,6,7] @@ -2011,7 +2011,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): # m3dSurf=mesh.buildFacePartOfMySelfNode(n,True); self.assertTrue(isinstance(m3dSurf,MEDCouplingUMesh)) - me=MEDCouplingExtrudedMesh.New(mesh,m3dSurf,0); + me=MEDCouplingMappedExtrudedMesh.New(mesh,m3dSurf,0); da=me.getMesh3DIds(); self.assertEqual(8,me.getNumberOfCells()); expected=[0,1,2,3,4,5,6,7] @@ -2049,7 +2049,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(3,m1.getSpaceDimension()); expected=[-0.3,-0.3,0., 0.2,-0.3,0., 0.7,-0.3,0., -0.3,0.2,0., 0.2,0.2,0., 0.7,0.2,0., -0.3,0.7,0., 0.2,0.7,0., 0.7,0.7,0.] val=m1.getCoords().getValues(); - for i in xrange(27): + for i in range(27): self.assertTrue(abs(expected[i]-val[i])<1e-14); pass pass @@ -2084,19 +2084,19 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(2,f.getNbOfGaussLocalization()); array=DataArrayDouble.New(); ptr=18*2*[None] - for i in xrange(18*2): + for i in range(18 * 2): ptr[i]=float(i+1) array.setValues(ptr,18,2); ptr=array.getPointer(); f.setArray(array); f.setName("MyFirstFieldOnGaussPoint"); - f.checkCoherency(); + f.checkConsistencyLight(); self.assertAlmostEqual(27.,f.getIJK(2,5,0),14); self.assertAlmostEqual(16.,f.getIJK(1,5,1),14); # f.clearGaussLocalizations(); self.assertEqual(0,f.getNbOfGaussLocalization()); - self.assertRaises(InterpKernelException,f.checkCoherency); + self.assertRaises(InterpKernelException,f.checkConsistencyLight); ids1=[0,1,3,4] self.assertRaises(InterpKernelException,f.setGaussLocalizationOnCells,ids1,_refCoo2,_gsCoo1,_wg1); self.assertEqual(0,f.getNbOfGaussLocalization()); @@ -2111,7 +2111,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(0,f.getGaussLocalizationIdOfOneCell(0)); self.assertEqual(1,f.getGaussLocalizationIdOfOneCell(1)); self.assertEqual(1,f.getGaussLocalizationIdOfOneCell(2)); - self.assertRaises(InterpKernelException,f.checkCoherency);#<- cell 3 has no localization + self.assertRaises(InterpKernelException,f.checkConsistencyLight);#<- cell 3 has no localization ids4=[3] _gsCoo2=_gsCoo1; _wg2=_wg1; @@ -2121,10 +2121,10 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(3,f.getNbOfGaussLocalization()); tmpIds=f.getCellIdsHavingGaussLocalization(0); self.assertEqual(ids2,list(tmpIds.getValues())); - self.assertRaises(InterpKernelException,f.checkCoherency);#<- it's always not ok because undelying array not with the good size. - array2=f.getArray().substr(0,10); + self.assertRaises(InterpKernelException,f.checkConsistencyLight);#<- it's always not ok because undelying array not with the good size. + array2=f.getArray().subArray(0,10); f.setArray(array2); - f.checkCoherency();#<- here it is OK + f.checkConsistencyLight();#<- here it is OK f2=f.clone(True); self.assertTrue(f.isEqual(f2,1e-14,1e-14)); gl1=f2.getGaussLocalization(0); @@ -2134,7 +2134,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertTrue(not f.isEqual(f2,1e-14,1e-14)); gl1.setGaussCoord(1,1,tmp); self.assertTrue(f.isEqual(f2,1e-14,1e-14)); - f2.checkCoherency(); + f2.checkConsistencyLight(); pass def testGaussPointNEField1(self): @@ -2146,14 +2146,14 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f.setDescription("MyDescriptionNE"); array=DataArrayDouble.New(); tmp=18*2*[None] - for i in xrange(18*2): + for i in range(18 * 2): tmp[i]=float(i+7) pass array.setValues(tmp,18,2); ptr=array.getPointer(); f.setArray(array); # - f.checkCoherency(); + f.checkConsistencyLight(); f2=f.clone(True); self.assertTrue(f.isEqual(f2,1e-14,1e-14)); self.assertAlmostEqual(21.,f.getIJK(2,0,0),14); @@ -2196,7 +2196,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): m2.orientCorrectlyPolyhedrons(); res1=m2.arePolyhedronsNotCorrectlyOriented(); self.assertTrue(len(res1)==0); - m2.checkCoherency(); + m2.checkConsistencyLight(); self.assertEqual(18,m2.getNumberOfCells()); cellIds2=[0,6,12] m2.convertToPolyTypes(cellIds2); @@ -2228,15 +2228,15 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(1,f3.getNumberOfComponents()); f3Ptr=f3.getArray().getValues(); expected1=[0.075,0.0375,0.0375,0.075,0.075, 0.1125,0.05625,0.05625,0.1125,0.1125, 0.0625,0.03125,0.03125,0.0625,0.0625]; - for i in xrange(15): + for i in range(15): self.assertTrue(abs(expected1[i]-f3Ptr[i])<1e-12); pass - f4=m5.getBarycenterAndOwner(); + f4=m5.computeCellCenterOfMass(); self.assertEqual(15,f4.getNumberOfTuples()); self.assertEqual(3,f4.getNumberOfComponents()); f4Ptr=f4.getValues(); expected2=[-0.05,-0.05,0.15, 0.3666666666666667,-0.13333333333333333,0.15, 0.53333333333333333,0.033333333333333333,0.15, -0.05,0.45,0.15, 0.45,0.45,0.15,-0.05,-0.05,0.525, 0.3666666666666667,-0.13333333333333333,0.525, 0.53333333333333333,0.033333333333333333,0.525, -0.05,0.45,0.525, 0.45,0.45,0.525,-0.05,-0.05,0.875, 0.3666666666666667,-0.13333333333333333,0.875, 0.53333333333333333,0.033333333333333333,0.875, -0.05,0.45,0.875, 0.45,0.45,0.875]; - for i in xrange(45): + for i in range(45): self.assertTrue(abs(expected2[i]-f4Ptr[i])<1e-12); pass pass @@ -2244,7 +2244,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testCellOrientation3(self): from cmath import rect - c = [rect(1.0, i*pi/4.0) for i in range(8)] + c = [rect(1.0, i * pi / 4.0) for i in range(8)] coords = [c[-1].real,c[-1].imag, c[3].real,c[3].imag, c[5].real,c[5].imag, c[1].real,c[1].imag] connec = [0,1,2,3] @@ -2274,35 +2274,35 @@ class MEDCouplingBasicsTest1(unittest.TestCase): myCoords=DataArrayDouble.New(); myCoords.setValues(coords,9,3); meshN.setCoords(myCoords); - meshN.checkCoherency(); + meshN.checkConsistencyLight(); # res1=meshN.arePolyhedronsNotCorrectlyOriented(); meshN.orientCorrectlyPolyhedrons(); self.assertTrue(len(res1)==0); - da=meshN.getBarycenterAndOwner(); + da=meshN.computeCellCenterOfMass(); self.assertEqual(1,da.getNumberOfTuples()); self.assertEqual(3,da.getNumberOfComponents()); daPtr=da.getValues(); ref=meshN.getCoords().getValues()[24:]; - for i in xrange(3): + for i in range(3): self.assertTrue(abs(ref[i]-daPtr[i])<1e-12); pass # center=[0.,0.,0.] vec=[0.,2.78,0.] - da=meshN.getBarycenterAndOwner(); + da=meshN.computeCellCenterOfMass(); daPtr=da.getValues(); ref=meshN.getCoords().getValues()[24:]; - for i in xrange(3): + for i in range(3): self.assertTrue(abs(ref[i]-daPtr[i])<1e-12); pass # meshN.rotate(center,vec,pi/7.); meshN.translate(vec); - da=meshN.getBarycenterAndOwner(); + da=meshN.computeCellCenterOfMass(); daPtr=da.getValues(); ref=meshN.getCoords().getValues()[24:]; - for i in xrange(3): + for i in range(3): self.assertTrue(abs(ref[i]-daPtr[i])<1e-12); pass # @@ -2310,10 +2310,10 @@ class MEDCouplingBasicsTest1(unittest.TestCase): vec2=[4.5,9.3,2.8] meshN.rotate(center2,vec2,e); meshN.translate(vec2); - da=meshN.getBarycenterAndOwner(); + da=meshN.computeCellCenterOfMass(); daPtr=da.getValues(); ref=meshN.getCoords().getValues()[24:]; - for i in xrange(3): + for i in range(3): self.assertTrue(abs(ref[i]-daPtr[i])<1e-10); pass pass @@ -2327,12 +2327,12 @@ class MEDCouplingBasicsTest1(unittest.TestCase): array.setValues(arr,m1.getNumberOfCells(),3); f1.setArray(array); # - f3=m1.getBarycenterAndOwner(); + f3=m1.computeCellCenterOfMass(); self.assertEqual(4,f3.getNumberOfTuples()); self.assertEqual(1,f3.getNumberOfComponents()); expected9=[0.75,5.105,0.8,5.155] ptr=f3.getValues(); - for i in xrange(4): + for i in range(4): self.assertTrue(abs(expected9[i]-ptr[i])<1e-12); pass # @@ -2341,13 +2341,13 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(1,f2.getNumberOfComponents()); expected1=[0.5,0.21,-0.6,-0.31] ptr=f2.getArray().getValues(); - for i in xrange(4): + for i in range(4): self.assertTrue(abs(expected1[i]-ptr[i])<1e-12); pass expected2=[0.5,0.21,0.6,0.31] f2=m1.getMeasureField(True); ptr=f2.getArray().getValues(); - for i in xrange(4): + for i in range(4): self.assertTrue(abs(expected2[i]-ptr[i])<1e-12); pass #integral @@ -2355,7 +2355,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): res=f1.integral(False); self.assertTrue(3,len(res)) expected3=[0.9866,-0.3615,0.4217] - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected3[i]-res[i])<1e-12); pass self.assertTrue(abs(expected3[0]-f1.integral(0,False))<1e-12); @@ -2363,14 +2363,14 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertTrue(abs(expected3[2]-f1.integral(2,False))<1e-12); res=f1.integral(True); expected4=[-3.4152,8.7639,-14.6879] - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected4[i]-res[i])<1e-12); pass #normL1 res=f1.normL1(); self.assertTrue(3,len(res)) expected5=[6.979506172839505, 16.89018518518518, 27.02969135802469] - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected5[i]-res[i])<1e-12); pass self.assertTrue(abs(expected5[0]-f1.normL1(0))<1e-12); @@ -2380,7 +2380,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): res=f1.normL2(); self.assertTrue(3,len(res)) expected7=[7.090910979452395, 16.9275542960123, 27.053271464160858] - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected7[i]-res[i])<1e-9); pass self.assertTrue(abs(expected7[0]-f1.normL2(0))<1e-9); @@ -2397,23 +2397,23 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(4,f2.getArray().getNumberOfTuples()); self.assertEqual(1,f2.getNumberOfComponents()); ptr=f2.getArray().getValues(); - for i in xrange(4): + for i in range(4): self.assertTrue(abs(sqrt(2.)*expected2[i]-ptr[i])<1e-12); pass f2=m1.getMeasureField(True); self.assertEqual(4,f2.getArray().getNumberOfTuples()); self.assertEqual(1,f2.getNumberOfComponents()); ptr=f2.getArray().getValues(); - for i in xrange(4): + for i in range(4): self.assertTrue(abs(expected2[i]*sqrt(2.)-ptr[i])<1e-12); pass #bary - f3=m1.getBarycenterAndOwner(); + f3=m1.computeCellCenterOfMass(); self.assertEqual(4,f3.getNumberOfTuples()); self.assertEqual(2,f3.getNumberOfComponents()); expected10=[0.75,0.75,5.105,5.105,0.8,0.8,5.155,5.155] ptr=f3.getValues(); - for i in xrange(8): + for i in range(8): self.assertTrue(abs(expected10[i]-ptr[i])<1e-12); pass # @@ -2423,19 +2423,19 @@ class MEDCouplingBasicsTest1(unittest.TestCase): array.setValues(arr,m1.getNumberOfCells(),3); f1.setArray(array); res=f1.integral(False); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(sqrt(2.)*expected4[i]-res[i])<1e-12); pass res=f1.integral(True); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(sqrt(2.)*expected4[i]-res[i])<1e-12); pass res=f1.normL1(); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected5[i]-res[i])<1e-12); pass res=f1.normL2(); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected7[i]-res[i])<1e-12); pass pass @@ -2447,20 +2447,20 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(1,f1.getNumberOfComponents()); expected1=[-0.5,-1,-1.5,-0.5,-1, 0.5,1,1.5,0.5,1] ptr=f1.getArray().getValues(); - for i in xrange(10): + for i in range(10): self.assertTrue(abs(expected1[i]-ptr[i])<1e-12); pass f1=m1.getMeasureField(True); ptr=f1.getArray().getValues(); - for i in xrange(10): + for i in range(10): self.assertTrue(abs(abs(expected1[i])-ptr[i])<1e-12); pass - f2=m1.getBarycenterAndOwner(); + f2=m1.computeCellCenterOfMass(); self.assertEqual(10,f2.getNumberOfTuples()); self.assertEqual(2,f2.getNumberOfComponents()); expected2=[0.5,0.3333333333333333,0.5,0.5,0.5,0.77777777777777777,0.5,0.3333333333333333,0.5,0.5,0.5,0.3333333333333333,0.5,0.5,0.5,0.77777777777777777,0.5,0.3333333333333333,0.5,0.5] ptr=f2.getValues(); - for i in xrange(20): + for i in range(20): self.assertTrue(abs(expected2[i]-ptr[i])<1e-12); pass m1.changeSpaceDimension(3); @@ -2468,15 +2468,15 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(10,f1.getArray().getNumberOfTuples()); self.assertEqual(1,f1.getNumberOfComponents()); ptr=f1.getArray().getValues(); - for i in xrange(10): + for i in range(10): self.assertTrue(abs(abs(expected1[i])-ptr[i])<1e-12); pass - f2=m1.getBarycenterAndOwner(); + f2=m1.computeCellCenterOfMass(); self.assertEqual(10,f2.getNumberOfTuples()); self.assertEqual(3,f2.getNumberOfComponents()); ptr=f2.getValues(); expected3=[0.5,0.3333333333333333,0.,0.5,0.5,0.,0.5,0.77777777777777777,0.,0.5,0.3333333333333333,0.,0.5,0.5,0., 0.5,0.3333333333333333,0.,0.5,0.5,0.,0.5,0.77777777777777777,0.,0.5,0.3333333333333333,0.,0.5,0.5,0.] - for i in xrange(30): + for i in range(30): self.assertTrue(abs(expected3[i]-ptr[i])<1e-12); pass pass @@ -2553,17 +2553,17 @@ class MEDCouplingBasicsTest1(unittest.TestCase): myCoords=DataArrayDouble.New(); myCoords.setValues(coords,69,3); meshN.setCoords(myCoords); - meshN.checkCoherency(); + meshN.checkConsistencyLight(); res1=meshN.arePolyhedronsNotCorrectlyOriented(); meshN.orientCorrectlyPolyhedrons(); res1=meshN.arePolyhedronsNotCorrectlyOriented(); self.assertTrue(len(res1)==0); # - da=meshN.getBarycenterAndOwner(); + da=meshN.computeCellCenterOfMass(); self.assertEqual(4,da.getNumberOfTuples()); self.assertEqual(3,da.getNumberOfComponents()); daPtr=da.getValues(); - for i in xrange(12): + for i in range(12): self.assertTrue(abs(barys[i]-daPtr[i])<1e-12); pass pass @@ -2579,22 +2579,22 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f.setArray(arr); renumber1=[3,1,0,4,2] loc=[-0.05,-0.05, 0.55,-0.25, 0.55,0.15, -0.05,0.45, 0.45,0.45] - for j in xrange(5): + for j in range(5): res=f.getValueOn(loc[2*j:2*j+2]); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(values1[i+3*j]-res[i])<1e-12); pass pass f.renumberCells(renumber1,False); ptr=f.getArray().getValues(); expected1=[9.,109.,10009.,8.,108.,10008.,11.,111.,10011.,7.,107.,10007.,10.,110.,10010.] - for i in xrange(15): + for i in range(15): self.assertTrue(abs(expected1[i]-ptr[i])<1e-12); pass #check that fields remains the same geometrically - for j in xrange(5): + for j in range(5): res=f.getValueOn(loc[2*j:2*(j+1)]); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(values1[i+3*j]-res[i])<1e-12); pass pass @@ -2621,14 +2621,14 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values2=[1.,1001.,2.,1002., 11.,1011.,12.,1012.,13.,1013.,14.,1014.,15.,1015.,16.,1016., 21.,1021.,22.,1022.,23.,1023.,24.,1024.,25.,1025.,26.,1026., 31.,1031.,32.,1032., 41.,1041.,42.,1042.] arr.setValues(values2,18,2); f.setArray(arr); - f.checkCoherency(); + f.checkConsistencyLight(); fCpy=f.clone(True); self.assertTrue(f.isEqual(fCpy,1e-12,1e-12)); f.renumberCells(renumber1,False); self.assertTrue(not f.isEqual(fCpy,1e-12,1e-12)); expected2=[21.,1021.,22.,1022.,23.,1023.,24.,1024.,25.,1025.,26.,1026., 11.,1011.,12.,1012.,13.,1013.,14.,1014.,15.,1015.,16.,1016., 41.,1041.,42.,1042., 1.,1001.,2.,1002., 31.,1031.,32.,1032.] ptr=f.getArray().getValues(); - for i in xrange(36): + for i in range(36): self.assertTrue(abs(expected2[i]-ptr[i])<1e-12); pass renumber2=[2,1,4,0,3] @@ -2641,14 +2641,14 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values3=[1.,1001.,2.,1002.,3.,1003.,4.,1004., 11.,1011.,12.,1012.,13.,1013., 21.,1021.,22.,1022.,23.,1023., 31.,1031.,32.,1032.,33.,1033.,34.,1034., 41.,1041.,42.,1042.,43.,1043.,44.,1044.] arr.setValues(values3,18,2); f.setArray(arr); - f.checkCoherency(); + f.checkConsistencyLight(); fCpy=f.clone(True); self.assertTrue(f.isEqual(fCpy,1e-12,1e-12)); f.renumberCells(renumber1,False); self.assertTrue(not f.isEqual(fCpy,1e-12,1e-12)); expected3=[21.,1021.,22.,1022.,23.,1023.,11.,1011.,12.,1012.,13.,1013.,41.,1041.,42.,1042.,43.,1043.,44.,1044.,1.,1001.,2.,1002.,3.,1003.,4.,1004.,31.,1031.,32.,1032.,33.,1033.,34.,1034.] ptr=f.getArray().getValues(); - for i in xrange(36): + for i in range(36): self.assertTrue(abs(expected3[i]-ptr[i])<1e-12); pass f.renumberCells(renumber2,False);#perform reverse operation of renumbering to check that the resulting field is equal. @@ -2666,13 +2666,13 @@ class MEDCouplingBasicsTest1(unittest.TestCase): values1=[7.,107.,10007.,8.,108.,10008.,9.,109.,10009.,10.,110.,10010.,11.,111.,10011.,12.,112.,10012.,13.,113.,10013.,14.,114.,10014.,15.,115.,10015.] arr.setValues(values1,nbOfNodes,3); f.setArray(arr); - f.checkCoherency(); + f.checkConsistencyLight(); renumber1=[0,4,1,3,5,2,6,7,8] loc=[0.5432,-0.2432, 0.5478,0.1528] expected1=[9.0272, 109.0272, 10009.0272, 11.4124,111.4124,10011.4124] - for j in xrange(2): + for j in range(2): res=f.getValueOn(loc[2*j:2*j+2]); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected1[i+3*j]-res[i])<1e-12); pass pass @@ -2680,14 +2680,14 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertTrue(f.isEqual(fCpy,1e-12,1e-12)); f.renumberNodes(renumber1); self.assertTrue(not f.isEqual(fCpy,1e-12,1e-12)); - for j in xrange(2): + for j in range(2): res=f.getValueOn(loc[2*j:2*j+2]); - for i in xrange(3): + for i in range(3): self.assertTrue(abs(expected1[i+3*j]-res[i])<1e-12); pass pass expected2=[7.,107.,10007.,9.,109.,10009.,12.,112.,10012.,10.,110.,10010.,8.,108.,10008.,11.,111.,10011.,13.,113.,10013.,14.,114.,10014.,15.,115.,10015.] - for i in xrange(27): + for i in range(27): self.assertTrue(abs(expected2[i]-f.getArray().getValues()[i])<1e-12); pass renumber2=[0,2,5,3,1,4,6,7,8] @@ -2697,24 +2697,24 @@ class MEDCouplingBasicsTest1(unittest.TestCase): def testConvertQuadraticCellsToLinear(self): mesh=MEDCouplingDataForTest.build2DTargetMesh_3(); - mesh.checkCoherency(); + mesh.checkConsistencyLight(); types=mesh.getAllGeoTypes(); types.sort() self.assertEqual(5,len(types)); - expected1=[NORM_POLYGON, NORM_TRI3, NORM_QUAD4, NORM_TRI6, NORM_QUAD8] + expected1 = [NORM_POLYGON, NORM_TRI3, NORM_QUAD4, NORM_TRI6, NORM_QUAD8] expected1.sort() self.assertEqual(expected1,types); self.assertTrue(mesh.isPresenceOfQuadratic()); - self.assertEqual(62,mesh.getMeshLength()); + self.assertEqual(62,mesh.getNodalConnectivityArrayLen()); f1=mesh.getMeasureField(False); # mesh.convertQuadraticCellsToLinear(); self.assertTrue(not mesh.isPresenceOfQuadratic()); # - mesh.checkCoherency(); + mesh.checkConsistencyLight(); f2=mesh.getMeasureField(False); self.assertTrue(f1.getArray().isEqual(f2.getArray(),1e-12)); - self.assertEqual(48,mesh.getMeshLength()); + self.assertEqual(48,mesh.getNodalConnectivityArrayLen()); types2=mesh.getAllGeoTypes(); types2.sort() self.assertEqual(3,len(types2)); @@ -2824,7 +2824,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): self.assertEqual(None,cellCor); self.assertNotEqual(None,nodeCor); expected1=[0, 1, 3, 4, 5, 6, 7, 8, 9] - for i in xrange(9): + for i in range(9): self.assertEqual(expected1[i],nodeCor.getIJ(i,0)); pass pass @@ -2846,7 +2846,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): m = MEDCouplingUMesh("tst", 2) m.setCoords(DataArrayDouble([], 0,2)) m.setConnectivity(DataArrayInt([]), DataArrayInt([0])) - m2 = m.deepCpy() + m2 = m.deepCopy() m.checkDeepEquivalWith(m2, 0, eps) # Should not raise! pass @@ -2863,7 +2863,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): a1.fillWithZero(); a1.setInfoOnComponent(0,"c"); a1.setInfoOnComponent(1,"d"); - a2=a1.deepCpy(); + a2=a1.deepCopy(); a2.setInfoOnComponent(0,"e"); a2.setInfoOnComponent(1,"f"); f.setArray(a1); @@ -2874,7 +2874,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): m.getCoords().setInfoOnComponent(1,"i"); m.getCoords().setInfoOnComponent(2,"j"); # - f.checkCoherency(); + f.checkConsistencyLight(); f2=f.clone(True); self.assertTrue(f2.isEqual(f,1e-12,1e-12)); f2.setName("smth"); @@ -2935,15 +2935,15 @@ class MEDCouplingBasicsTest1(unittest.TestCase): myCoords=DataArrayDouble.New(); myCoords.setValues(targetCoords,4,2); m2.setCoords(myCoords); - m2.checkCoherency(); - m1.checkCoherency(); + m2.checkConsistencyLight(); + m1.checkConsistencyLight(); # expected1=[0.25,0.125,0.125,0.25,0.25] f1=m1.getMeasureField(False); f2=m2.getMeasureField(False); self.assertEqual(5,f1.getArray().getNumberOfTuples()); self.assertEqual(1,f2.getArray().getNumberOfTuples()); - for i in xrange(5): + for i in range(5): self.assertAlmostEqual(expected1[i],f1.getIJ(i,0),12); pass self.assertAlmostEqual(expected1[0],f2.getIJ(0,0),12); @@ -2954,7 +2954,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f2=m2.getMeasureField(False); self.assertEqual(5,f1.getArray().getNumberOfTuples()); self.assertEqual(1,f2.getArray().getNumberOfTuples()); - for i in xrange(5): + for i in range(5): self.assertAlmostEqual(expected1[i],f1.getIJ(i,0),12); pass self.assertAlmostEqual(expected1[0],f2.getIJ(0,0),12); @@ -2978,7 +2978,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f1.changeUnderlyingMesh(mesh2,10,1e-12); #self.assertTrue(f1.getMesh()==mesh2); expected1=[7.,107.,9.,109.,8.,108.,10.,110.,11.,111.,12.,112.,13.,113.,15.,115.,14.,114.,16.,116.] - for i in xrange(20): + for i in range(20): self.assertAlmostEqual(expected1[i],f1.getArray().getIJ(0,i),12); pass # @@ -2995,7 +2995,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f1.changeUnderlyingMesh(mesh2,10,1e-12); #self.assertTrue(f1.getMesh()==mesh2); expected2=[7.,107.,17.,117.,8.,108.,10.,110.,11.,111.,12.,112.,13.,113.,15.,115.,14.,114.,16.,116.,9.,109.] - for i in xrange(22): + for i in range(22): self.assertAlmostEqual(expected2[i],f1.getArray().getIJ(0,i),12); pass pass @@ -3014,7 +3014,7 @@ class MEDCouplingBasicsTest1(unittest.TestCase): f.setArray(a1); f.setEndArray(a2); f.setEndTime(3.,3,4); - f.checkCoherency(); + f.checkConsistencyLight(); # self.assertAlmostEqual(8.,f.getMaxValue(),14); self.assertAlmostEqual(0.,f.getMinValue(),14); @@ -3076,13 +3076,13 @@ class MEDCouplingBasicsTest1(unittest.TestCase): # f3=f1.dot(f2); expected1=[842.,1820.,2816.,3830.,4862.,5912.,6980.,8066.,9170.,10292.] - for i in xrange(10): + for i in range(10): self.assertAlmostEqual(expected1[i],f3.getIJ(i,0),9); pass # f4=f1.crossProduct(f2); expected2=[-93., 186., -93., -392., 784., -392., -691., 1382., -691., -990., 1980., -990., -1289., 2578., -1289., -1588., 3176., -1588., -1887., 3774., -1887., -2186., 4372., -2186., -2485., 4970., -2485., -2784., 5568., -2784.] - for i in xrange(30): + for i in range(30): self.assertAlmostEqual(expected2[i],f4.getIJ(0,i),9); pass pass