pass
array.setValues(ptr,sourceMesh.getNumberOfCells(),1);
srcField.setArray(array);
+ srcField.setName("abc") ; srcField.setDescription("def")
+ srcField.setTime(7.7,9,10)
trgfield=remapper.transferField(srcField,4.57);
+ self.assertEqual("abc",trgfield.getName())
+ self.assertEqual("def",trgfield.getDescription())
+ a,b,c=trgfield.getTime()
+ self.assertAlmostEqual(7.7,a,14)
+ self.assertEqual(b,9)
+ self.assertEqual(c,10)
values=trgfield.getArray().getValues();
valuesExpected=[7.5 ,7. ,7.,8.,7.5];
for i in xrange(targetMesh.getNumberOfCells()):
self.assertAlmostEqual(valuesExpected[i0],values[i0],12);
pass
pass
+
+ def testPrepareUC(self):
+ # 1D
+ coords=DataArrayDouble([0.,0.5,0.7])
+ src=MEDCouplingUMesh("",1) ; src.setCoords(coords)
+ src.allocateCells(2) ; src.insertNextCell(NORM_SEG2,[0,1]) ; src.insertNextCell(NORM_SEG2,[1,2]) ; src.finishInsertingCells()
+ trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
+ trg.setCoordsAt(0,arr)
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrSrc)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected1=[-7.,4.,36.,-7.,-7.]
+ self.assertEqual(5,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(5,len(expected1))
+ for i,val in enumerate(expected1):
+ self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ # 2D
+ coords=DataArrayDouble([0.,0.,0.,1.,1.,1.,1.,0.,0.5,-0.2],5,2)
+ src=MEDCouplingUMesh("",2) ; src.setCoords(coords)
+ src.allocateCells(2) ; src.insertNextCell(NORM_TRI3,[0,1,2]) ; src.insertNextCell(NORM_TRI3,[3,4,0]) ; src.finishInsertingCells()
+ trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
+ trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr)
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrSrc)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected2=[-7.,-7.,7.35,0.15,-7.,-7.,2.8,14.85,5.25,-7.,-7.,2.,2.5,-7.,-7.,-7.,1.2,3.,0.9,-7.,-7.,-7.,-7.,-7.,-7.]
+ self.assertEqual(25,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(25,len(expected2))
+ for i,val in enumerate(expected2):
+ self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ # 3D
+ coords=DataArrayDouble([0.,0.,0.,0.,1.,0.,1.,1.,0.,1.,0.,0.,0.5,-0.2,0.,0.1,0.8,1.,0.5,0.,1.],7,3)
+ src=MEDCouplingUMesh("",3) ; src.setCoords(coords)
+ src.allocateCells(2) ; src.insertNextCell(NORM_TETRA4,[0,1,2,5]) ; src.insertNextCell(NORM_TETRA4,[3,4,0,6]) ; src.finishInsertingCells()
+ trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3]) ; arr2=DataArrayDouble([-0.7,0.2,0.6,1.2,2.])
+ trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr) ; trg.setCoordsAt(2,arr2)
+ src.checkCoherency2(1e-10)
+ trg.checkCoherency()
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrSrc)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected3=[-7.,-7.,2.925,0.015,-7.,-7.,0.9392,8.595,2.265,-7.,-7.,1.1008,1.1192,-7.,-7.,-7.,0.6392,1.6408,0.2808,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,0.81,-7.,-7.,-7.,0.1208,11.55,0.96,-7.,-7.,1.1752,0.6592,-7.,-7.,-7.,0.8512,1.7744,0.0192,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,1.92,-7.,-7.,-7.,0.12578571428571422,0.007314285714285673,-7.,-7.,-7.,0.3189253968253971,0.1879746031746033,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.,-7.]
+ self.assertEqual(100,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(100,len(expected3))
+ for i,val in enumerate(expected3):
+ self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ pass
+
+ def testPrepareCU(self):
+ # 1D
+ coords=DataArrayDouble([0.,0.5,0.7])
+ trg=MEDCouplingUMesh("",1) ; trg.setCoords(coords)
+ trg.allocateCells(2) ; trg.insertNextCell(NORM_SEG2,[0,1]) ; trg.insertNextCell(NORM_SEG2,[1,2]) ; trg.finishInsertingCells()
+ src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
+ src.setCoordsAt(0,arr)
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrTrg=DataArrayDouble([10.,30.,40.,70.,80.])
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrTrg)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected1=[44.,16.]
+ self.assertEqual(2.,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(2,len(expected1))
+ for i,val in enumerate(expected1):
+ self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ # 2D
+ coords=DataArrayDouble([0.,0.,0.,1.,1.,1.,1.,0.,0.5,-0.2],5,2)
+ trg=MEDCouplingUMesh("",2) ; trg.setCoords(coords)
+ trg.allocateCells(2) ; trg.insertNextCell(NORM_TRI3,[0,1,2]) ; trg.insertNextCell(NORM_TRI3,[3,4,0]) ; trg.finishInsertingCells()
+ src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
+ src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr)
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.,40.,70.,80.,110.,130.,140.,170.,180.,210.,230.,240.,270.,280.,310.,330.,340.,370.,380.,410.,430.,440.,470.,480.])
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrSrc)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected2=[441.3050624589086,68.69529914529915]
+ self.assertEqual(2,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(2,len(expected2))
+ for i,val in enumerate(expected2):
+ self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ # 3D
+ coords=DataArrayDouble([0.,0.,0.,0.,1.,0.,1.,1.,0.,1.,0.,0.,0.5,-0.2,0.,0.1,0.8,1.,0.5,0.,1.],7,3)
+ trg=MEDCouplingUMesh("",3) ; trg.setCoords(coords)
+ trg.allocateCells(2) ; trg.insertNextCell(NORM_TETRA4,[0,1,2,5]) ; trg.insertNextCell(NORM_TETRA4,[3,4,0,6]) ; trg.finishInsertingCells()
+ src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3]) ; arr2=DataArrayDouble([-0.7,0.2,0.6,1.2,2.])
+ src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr) ; src.setCoordsAt(2,arr2)
+ trg.checkCoherency2(1e-10)
+ src.checkCoherency()
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble(100) ; arrSrc.iota(7.7)
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrSrc)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected3=[39.635196634558845,12.13422356758468]
+ self.assertEqual(2,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(2,len(expected3))
+ for i,val in enumerate(expected3):
+ self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ pass
+
+ def testPrepareCC(self):
+ # 1D
+ src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
+ src.setCoordsAt(0,arr)
+ trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
+ trg.setCoordsAt(0,arr)
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrTrg=DataArrayDouble([10.,30.,40.,70.,80.])
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrTrg)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected1=[10.,25.]
+ self.assertEqual(2.,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(2,len(expected1))
+ for i,val in enumerate(expected1):
+ self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ # 2D
+ src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
+ src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr)
+ trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
+ trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr)
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.,40.,70.,80.,110.,130.,140.,170.,180.,210.,230.,240.,270.,280.,310.,330.,340.,370.,380.,410.,430.,440.,470.,480.])
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrSrc)
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected2=[10.,25.,91.66666666666666,90.27777777777777]
+ self.assertEqual(4,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(4,len(expected2))
+ for i,val in enumerate(expected2):
+ self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ # 3D
+ src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
+ src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr) ; src.setCoordsAt(2,arr)
+ trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
+ trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr) ; trg.setCoordsAt(2,arr)
+ fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble(125) ; arrSrc.iota(7.7)
+ fieldSrc.setNature(Integral) ; fieldSrc.setArray(arrSrc) ; fieldSrc.checkCoherency()
+ rem=MEDCouplingRemapper()
+ rem.prepare(src,trg,"P0P0")
+ trgField=rem.transferField(fieldSrc,-7.)
+ expected3=[7.7, 7.249999999999999, 10.583333333333332, 9.513888888888886, 27.25, 23.40277777777777, 26.180555555555546, 22.39583333333333]
+ self.assertEqual(8,trgField.getArray().getNumberOfTuples())
+ self.assertEqual(8,len(expected3))
+ for i,val in enumerate(expected3):
+ self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
+ pass
+ pass
+
+ # Bug when source mesh is not homogeneously oriented in source mesh
+ def testNonRegressionNonHomegenousOrriented3DCells(self):
+ csrc=DataArrayDouble([-0.15240000188350677,0,0,-0.1086929515004158,0,0,-0.15240000188350677,0.018142856657505035,0,-0.13054648041725159,0.0090714283287525177,0.019050000235438347,-0.13054648041725159,0.0090714283287525177,0],5,3)
+ src1=MEDCouplingUMesh("src",3) ; src1.allocateCells(0) ; src1.insertNextCell(NORM_TETRA4,[0,1,4,3]) ; src1.insertNextCell(NORM_TETRA4,[2,0,4,3])
+ src2=MEDCouplingUMesh("src",3) ; src2.allocateCells(0) ; src2.insertNextCell(NORM_TETRA4,[0,4,1,3]) ; src2.insertNextCell(NORM_TETRA4,[2,0,4,3])
+ src1.setCoords(csrc) ; src2.setCoords(csrc)
+ ctrg=DataArrayDouble([-0.15240000188350677,-0.038100000470876694,0,0.32379999756813049,-0.038100000470876694,0,-0.15240000188350677,0.076200000941753387,0,0.32379999756813049,0.076200000941753387,0,-0.15240000188350677,-0.038100000470876694,0.076200000941753387,0.32379999756813049,-0.038100000470876694,0.076200000941753387,-0.15240000188350677,0.076200000941753387,0.076200000941753387,0.32379999756813049,0.076200000941753387,0.076200000941753387],8,3)
+ trg=MEDCouplingUMesh("trg",3) ; trg.allocateCells(0) ; trg.insertNextCell(NORM_HEXA8,[0,1,3,2,4,5,7,6])
+ trg.setCoords(ctrg)
+ rem1=MEDCouplingRemapper() ; rem1.setSplittingPolicy(PLANAR_FACE_5) ; rem1.prepare(src1,trg,"P0P0")
+ rem2=MEDCouplingRemapper() ; rem2.setSplittingPolicy(PLANAR_FACE_5) ; rem2.prepare(src1,trg,"P0P0")
+ mat1=rem1.getCrudeMatrix() ; mat2=rem2.getCrudeMatrix()
+ self.assertEqual(1,len(mat1)) ; self.assertEqual(1,len(mat2))
+ self.assertEqual(mat1[0].keys(),mat2[0].keys()) ; self.assertEqual([0,1],mat1[0].keys())
+ self.assertAlmostEqual(1.25884108122e-06,mat1[0][0],16) ; self.assertAlmostEqual(1.25884108122e-06,mat2[0][0],16)
+ self.assertAlmostEqual(1.25884086663e-06,mat1[0][1],16) ; self.assertAlmostEqual(1.25884086663e-06,mat2[0][1],16)
+ #
+ d=DataArrayDouble([13.45,27.67],2,1)
+ f1=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f1.setMesh(src1) ; f1.setArray(d) ; f1.setNature(RevIntegral)
+ f2=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f2.setMesh(src2) ; f2.setArray(d) ; f2.setNature(RevIntegral)
+ f11=rem1.transferField(f1,1e300) ; f22=rem2.transferField(f2,1e300)
+ expected1=DataArrayDouble([0.012480539537637884])
+ self.assertTrue(f11.getArray().isEqual(expected1,1e-15))
+ self.assertTrue(f22.getArray().isEqual(expected1,1e-15))
+ #
+ f1.setNature(Integral) ; f2.setNature(Integral)
+ f11=rem1.transferField(f1,1e300) ; f22=rem2.transferField(f2,1e300)
+ #
+ expected2=DataArrayDouble([41.12])
+ self.assertTrue(f11.getArray().isEqual(expected2,1e-13))
+ self.assertTrue(f22.getArray().isEqual(expected2,1e-13))
+ pass
+
+ def testCellToNodeReverse3D(self):
+ c=DataArrayDouble([0.,1.,2.5])
+ cc=MEDCouplingCMesh()
+ cc.setCoords(c,c,c)
+ um=cc.buildUnstructured()
+ f=um.getMeasureField(ON_CELLS)
+ #
+ n2o=um.simplexize(PLANAR_FACE_5)
+ f.setArray(f.getArray()[n2o])
+ f.checkCoherency()
+ f.setNature(ConservativeVolumic)
+ f.setTime(5.6,7,8)
+ f.setName("toto") ; f.setDescription("aDescription")
+ p=MEDCouplingRemapper()
+ p.setP1P0BaryMethod(True)
+ p.prepare(um,um,"P1P0")
+ fNode=p.reverseTransferField(f,1e300)
+ self.assertEqual("toto",fNode.getName())
+ self.assertEqual("aDescription",fNode.getDescription())
+ a,b,c=fNode.getTime()
+ self.assertAlmostEqual(5.6,a,14)
+ self.assertEqual(7,b) ; self.assertEqual(8,c)
+ #
+ integExpected=34.328125
+ self.assertAlmostEqual(fNode.integral(False)[0],integExpected,14)
+ self.assertAlmostEqual(f.integral(False)[0],integExpected,14)
+ pass
def build2DSourceMesh_1(self):
sourceCoords=[-0.3,-0.3, 0.7,-0.3, -0.3,0.7, 0.7,0.7]
