# -*- 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
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):
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()));
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
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)
#
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
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);
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);
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
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]
#
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
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());
#
eltsV=[1,3];
mesh.convertToPolyTypes(eltsV);
- mesh.checkCoherency();
+ mesh.checkConsistencyLight();
#
desc=DataArrayInt.New();
descIndx=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());
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());
#
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());
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));
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):
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)
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())
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)
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();
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());
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());
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));
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();
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
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();
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);
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);
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();
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
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());
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());
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());
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();
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());
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());
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());
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());
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();
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());
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
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());
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
#
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
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);
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);
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
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,
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();
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
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,
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
#
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
#
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);
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);
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];
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
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.]
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
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
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);
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);
self.assertTrue(mesh1.isEqual(mesh2, 1e-5));
self.assertTrue(not mesh1.isEqual(mesh2, 1e-7));
- self.assertRaises(InterpKernelException, mesh3.checkCoherency1, 1e-12);
- mesh1.checkCoherency1(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]);
-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
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]
#
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]
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
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());
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;
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);
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):
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);
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);
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
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]
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
#
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
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
#
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
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);
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);
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);
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
#
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
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);
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
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
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
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]
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.
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
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]
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));
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
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
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);
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");
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);
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);
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
#
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
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);
#
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