self.assertEqual((1.,0.5,2.),m3.getDXYZ())
# now playing with 3D surf
m4=MEDCouplingIMesh("",3,DataArrayInt([3,1,4]),DataArrayDouble([1.5,2.5,3.5]),DataArrayDouble((0.5,1.,0.25))) ; m4.setAxisUnit("km")
+ self.assertEqual([(1.5,2.5),(2.5,3.5),(3.5,4.25)],m4.getBoundingBox())
self.assertEqual(3,m4.getSpaceDimension())
self.assertEqual(2,m4.getMeshDimension())
self.assertEqual(12,m4.getNumberOfNodes())
self.assertEqual(1,amr.getNumberOfPatches())
self.assertEqual(2,amr.getMaxNumberOfLevelsRelativeToThis())
self.assertEqual(2,amr.getSpaceDimension())
- amr[0].addPatch([(2,3),(1,3)],[2,2])
+ amr[0].addPatch([(2,3),(1,3)],[3,2])
self.assertEqual(amr[0].getBLTRRange(),[(1,2),(0,1)])
self.assertEqual(4,amr.getNumberOfCellsAtCurrentLevel())
- self.assertEqual(28,amr.getNumberOfCellsRecursiveWithOverlap())
- self.assertEqual(25,amr.getNumberOfCellsRecursiveWithoutOverlap())
+ self.assertEqual(32,amr.getNumberOfCellsRecursiveWithOverlap())
+ self.assertEqual(29,amr.getNumberOfCellsRecursiveWithoutOverlap())
self.assertEqual(1,amr.getNumberOfPatches())
self.assertEqual(3,amr.getMaxNumberOfLevelsRelativeToThis())
self.assertEqual(2,amr.getSpaceDimension())
- amr[0].addPatch([(0,2),(3,4)],[3,3])
+ amr[0].addPatch([(0,2),(3,4)],[3,2])
self.assertEqual(16,amr[0].getMesh().getNumberOfCellsAtCurrentLevel())
- self.assertEqual(46,amr.getNumberOfCellsRecursiveWithOverlap())
- self.assertEqual(41,amr.getNumberOfCellsRecursiveWithoutOverlap())
+ self.assertEqual(44,amr.getNumberOfCellsRecursiveWithOverlap())
+ self.assertEqual(39,amr.getNumberOfCellsRecursiveWithoutOverlap())
self.assertEqual(2,amr[0].getMesh().getNumberOfPatches())
self.assertEqual(3,amr.getMaxNumberOfLevelsRelativeToThis())
self.assertEqual(2,amr.getSpaceDimension())
del amr[0][1]
self.assertEqual(amr[0].getBLTRRange(),[(1,2),(0,1)])
self.assertEqual(4,amr.getNumberOfCellsAtCurrentLevel())
- self.assertEqual(28,amr.getNumberOfCellsRecursiveWithOverlap())
- self.assertEqual(25,amr.getNumberOfCellsRecursiveWithoutOverlap())
+ self.assertEqual(32,amr.getNumberOfCellsRecursiveWithOverlap())
+ self.assertEqual(29,amr.getNumberOfCellsRecursiveWithoutOverlap())
self.assertEqual(1,amr.getNumberOfPatches())
self.assertEqual(3,amr.getMaxNumberOfLevelsRelativeToThis())
self.assertEqual(2,amr.getSpaceDimension())
self.assertTrue(m2.getData().isEqual(DataArrayDouble([29,37,37,62]),1e-12))
pass
- def testAMR2(self):
+ def testSwig2AMR2(self):
""" Test condensation of fine IMesh instance into a coarse one, with a factor. See testRemapperAMR1 in MEDCouplingRemapperTest.py file to see how the expected value is obtained."""
coarse=DataArrayDouble(35) ; coarse.iota(0) #X=5,Y=7
fine=DataArrayDouble(3*2*4*4) ; fine.iota(0) #X=3,Y=2 refined by 4
- MEDCouplingIMesh.CondenseFineToCoarse(coarse,[5,7],fine,[(1,4),(2,4)],[4,4])
+ MEDCouplingIMesh.CondenseFineToCoarse([5,7],fine,[(1,4),(2,4)],[4,4],coarse)
self.assertTrue(coarse.isEqual(DataArrayDouble([0,1,2,3,4,5,6,7,8,9,10,312,376,440,14,15,1080,1144,1208,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34]),1e-12))
# 3D
coarse=DataArrayDouble(175) ; coarse.iota(0) #X=5,Y=7,Z=5
fine=DataArrayDouble(3*2*3*4*4*4) ; fine.iota(0) #X=3,Y=2,Z=3 refined by 4
- MEDCouplingIMesh.CondenseFineToCoarse(coarse,[5,7,5],fine,[(1,4),(2,4),(1,4)],[4,4,4])
+ MEDCouplingIMesh.CondenseFineToCoarse([5,7,5],fine,[(1,4),(2,4),(1,4)],[4,4,4],coarse)
self.assertTrue(coarse.isEqual(DataArrayDouble([0.,1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24.,25.,26.,27.,28.,29.,30.,31.,32.,33.,34.,35.,36.,37.,38.,39.,40.,41.,42.,43.,44.,45.,10464.,10720.,10976.,49.,50.,13536.,13792.,14048.,54.,55.,56.,57.,58.,59.,60.,61.,62.,63.,64.,65.,66.,67.,68.,69.,70.,71.,72.,73.,74.,75.,76.,77.,78.,79.,80.,35040.,35296.,35552.,84.,85.,38112.,38368.,38624.,89.,90.,91.,92.,93.,94.,95.,96.,97.,98.,99.,100.,101.,102.,103.,104.,105.,106.,107.,108.,109.,110.,111.,112.,113.,114.,115.,59616.,59872.,60128.,119.,120.,62688.,62944.,63200.,124.,125.,126.,127.,128.,129.,130.,131.,132.,133.,134.,135.,136.,137.,138.,139.,140.,141.,142.,143.,144.,145.,146.,147.,148.,149.,150.,151.,152.,153.,154.,155.,156.,157.,158.,159.,160.,161.,162.,163.,164.,165.,166.,167.,168.,169.,170.,171.,172.,173.,174.]),1e-12))
# 1D
coarse=DataArrayDouble(5) ; coarse.iota(0) #X=5
fine=DataArrayDouble(3*4) ; fine.iota(0) #X=3 refined by 4
- MEDCouplingIMesh.CondenseFineToCoarse(coarse,[5],fine,[(1,4)],[4])
+ MEDCouplingIMesh.CondenseFineToCoarse([5],fine,[(1,4)],[4],coarse)
self.assertTrue(coarse.isEqual(DataArrayDouble([0,6,22,38,4]),1e-12))
pass
- def testAMR3(self):
+ def testSwig2AMR3(self):
""" Test spread of coarse IMesh instance into a fine one, with a factor."""
coarse=DataArrayDouble(35) ; coarse.iota(0) #X=5,Y=7
fine=DataArrayDouble(3*2*4*4) ; fine.iota(0) #X=3,Y=2 refined by 4
self.assertTrue(fine.isEqual(DataArrayDouble([1.,1.,1.,1.,2.,2.,2.,2.,3.,3.,3.,3.]),1e-12))
pass
- def testAMR4(self):
+ def testSwig2AMR4(self):
"""This test focuses on MEDCouplingCartesianAMRMesh.createPatchesFromCriterion method. To test it a field containing 0 everywhere except in the annulus (centered on the center of the mesh) value is 1."""
im=MEDCouplingIMesh("mesh",2,[51,51],[0.,0.],[0.04,0.04])
b=im.getBarycenterAndOwner() ; b-=[1.,1.] ; b=b.magnitude()
#
self.assertEqual(MEDCouplingStructuredMesh.ChangeReferenceToGlobalOfCompactFrmt([(8,32),(4,17)],[(0,24),(2,12)]),[(8,32),(6,16)])
self.assertEqual(MEDCouplingStructuredMesh.ChangeReferenceFromGlobalOfCompactFrmt([(8,32),(4,17)],[(8,32),(6,16)]),[(0,24),(2,12)])
+ self.assertTrue(amr.getImageMesh().isEqual(im,1e-12))
+ m=amr.getImageMesh().asSingleCell().build1SGTUnstructured()
+ self.assertTrue(m.getNodalConnectivity().isEqual(DataArrayInt([1,0,2,3])))
+ self.assertTrue(m.getCoords().isEqualWithoutConsideringStr(DataArrayDouble([(0,0),(2,0),(0,2),(2,2)]),1e-12))
+ pass
+
+ def testSwig2AMR5(self):
+ """ Idem testAMR3, test spread of coarse IMesh instance into a fine one, with a factor, but here ghost is used !"""
+ # 1D
+ coarse=DataArrayDouble(5+2) ; coarse.iota(-1) #X=5 with ghostLev=1
+ fine=DataArrayDouble(3*4+2) ; fine.iota(1000) #X=3 refined by 4 with ghostLev=1
+ MEDCouplingIMesh.SpreadCoarseToFineGhost(coarse,[5],fine,[(1,4)],[4],1)
+ self.assertTrue(fine.isEqual(DataArrayDouble([0,1,1,1,1,2,2,2,2,3,3,3,3,4]),1e-12))
+ coarse.iota(-1000)
+ MEDCouplingIMesh.CondenseFineToCoarseGhost([5],fine,[(1,4)],[4],coarse,1)
+ self.assertTrue(coarse.isEqual(DataArrayDouble([-1000.,-999.,4.,8.,12.,-995.,-994.]),1e-12))
+ # 2D
+ coarse=DataArrayDouble((5+2*1)*(7+2*1)) ; coarse.iota(0) #X=5,Y=7 with ghostLev=1
+ fine=DataArrayDouble((3*4+2*1)*(2*4+2*1)) ; fine.iota(1000) #X=3,Y=2 refined by 4
+ MEDCouplingIMesh.SpreadCoarseToFineGhost(coarse,[5,7],fine,[(1,4),(2,4)],[4,4],1)
+ self.assertTrue(fine.isEqual(DataArrayDouble([15.,16.,16.,16.,16.,17.,17.,17.,17.,18.,18.,18.,18.,19.,22.,23.,23.,23.,23.,24.,24.,24.,24.,25.,25.,25.,25.,26.,22.,23.,23.,23.,23.,24.,24.,24.,24.,25.,25.,25.,25.,26.,22.,23.,23.,23.,23.,24.,24.,24.,24.,25.,25.,25.,25.,26.,22.,23.,23.,23.,23.,24.,24.,24.,24.,25.,25.,25.,25.,26.,29.,30.,30.,30.,30.,31.,31.,31.,31.,32.,32.,32.,32.,33.,29.,30.,30.,30.,30.,31.,31.,31.,31.,32.,32.,32.,32.,33.,29.,30.,30.,30.,30.,31.,31.,31.,31.,32.,32.,32.,32.,33.,29.,30.,30.,30.,30.,31.,31.,31.,31.,32.,32.,32.,32.,33.,36.,37.,37.,37.,37.,38.,38.,38.,38.,39.,39.,39.,39.,40.]),1e-12))
+ f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(MEDCouplingIMesh("",2,DataArrayInt([8,10]),[0.,0.],DataArrayDouble((1.,1.)))) ; f.setArray(coarse) ; f.setName("tutu") ; f.checkCoherency() ; f.writeVTK("coarse.vti")
+ coarse.iota(-1000)
+ fine2=DataArrayDouble.Meld(fine,3*fine) ; coarse2=DataArrayDouble.Meld(coarse,3*coarse)
+ MEDCouplingIMesh.CondenseFineToCoarseGhost([5,7],fine,[(1,4),(2,4)],[4,4],coarse,1)
+ MEDCouplingIMesh.CondenseFineToCoarseGhost([5,7],fine2,[(1,4),(2,4)],[4,4],coarse2,1)
+ f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(MEDCouplingIMesh("",2,DataArrayInt([8,10]),[0.,0.],DataArrayDouble((1.,1.)))) ; f.setArray(coarse) ; f.setName("tutu") ; f.checkCoherency() ; f.writeVTK("coarse.vti")
+ coarseExp=DataArrayDouble([-1000.,-999.,-998.,-997.,-996.,-995.,-994.,-993.,-992.,-991.,-990.,-989.,-988.,-987.,-986.,-985.,-984.,-983.,-982.,-981.,-980.,-979.,-978.,368.,384.,400.,-974.,-973.,-972.,-971.,480.,496.,512.,-967.,-966.,-965.,-964.,-963.,-962.,-961.,-960.,-959.,-958.,-957.,-956.,-955.,-954.,-953.,-952.,-951.,-950.,-949.,-948.,-947.,-946.,-945.,-944.,-943.,-942.,-941.,-940.,-939.,-938.])
+ self.assertTrue(coarse.isEqual(coarseExp,1e-12))
+ self.assertTrue(coarse2[:,0].isEqual(coarseExp,1e-12))
+ self.assertTrue(coarse2[:,1].isEqual(3*coarseExp,1e-12))
+ pass
+
+ def testSwig2AMR6(self):
+ """ Idem testSwig2AMR5, except that only 2D is considered here, and fine to fine is considered here. At the end of the test some checks about typing with AMR structs."""
+ amr=MEDCouplingCartesianAMRMesh("",2,[6,6],[0,0],[1,1])
+ da=DataArrayDouble((5+2)*(5+2)) ; da.iota() ; da+=0.9
+ amr.addPatch([(1,4),(2,4)],[4,4])
+ amr.addPatch([(0,1),(0,1)],[4,4])
+ amr.addPatch([(4,5),(3,4)],[4,4])
+ amr.addPatch([(4,5),(1,3)],[4,4])
+ amr.addPatch([(0,1),(1,4)],[4,4])
+ da0=DataArrayDouble((3*4+2)*(2*4+2)) ; da0.iota() ; da0[:]+=0.2
+ da1=DataArrayDouble((1*4+2)*(1*4+2)) ; da1.iota() ; da1[:]+=0.4
+ da2=DataArrayDouble((1*4+2)*(1*4+2)) ; da2.iota() ; da2[:]+=0.6
+ da3=DataArrayDouble((1*4+2)*(2*4+2)) ; da3.iota() ; da3[:]+=0.7
+ da4=DataArrayDouble((1*4+2)*(3*4+2)) ; da4.iota() ; da4[:]+=0.8
+ self.assertEqual(5,amr.getNumberOfPatches())
+ l=[da0,da1,da2,da3,da4]
+ lCpy=[elt.deepCpy() for elt in l]
+ l2=[DataArrayDouble.Meld(elt,3*elt) for elt in l]
+ amr.fillCellFieldOnPatchGhostAdv(0,da,1,l)
+ amr.fillCellFieldOnPatchGhostAdv(0,DataArrayDouble.Meld(da,3*da),1,l2)
+ amr.fillCellFieldOnPatchOnlyOnGhostZone(0,da,lCpy[0],1)
+ #
+ f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(amr.getImageMesh().buildWithGhost(1)) ; f.setArray(da) ; f.setName("all")
+ f0=MEDCouplingFieldDouble(ON_CELLS) ; f0.setMesh(amr[0].getMesh().getImageMesh().buildWithGhost(1)) ; f0.setArray(da0) ; f0.setName("p0") ; f0.checkCoherency()
+ f1=MEDCouplingFieldDouble(ON_CELLS) ; f1.setMesh(amr[1].getMesh().getImageMesh().buildWithGhost(1)) ; f1.setArray(da1) ; f1.setName("p1") ; f1.checkCoherency()
+ f2=MEDCouplingFieldDouble(ON_CELLS) ; f2.setMesh(amr[2].getMesh().getImageMesh().buildWithGhost(1)) ; f2.setArray(da2) ; f2.setName("p2") ; f2.checkCoherency()
+ f3=MEDCouplingFieldDouble(ON_CELLS) ; f3.setMesh(amr[3].getMesh().getImageMesh().buildWithGhost(1)) ; f3.setArray(da3) ; f3.setName("p3") ; f3.checkCoherency()
+ f4=MEDCouplingFieldDouble(ON_CELLS) ; f4.setMesh(amr[4].getMesh().getImageMesh().buildWithGhost(1)) ; f4.setArray(da4) ; f4.setName("p4") ; f4.checkCoherency()
+ #
+ da0Exp=DataArrayDouble([28.8,16.9,16.9,16.9,16.9,17.9,17.9,17.9,17.9,18.9,18.9,18.9,18.9,25.7,34.8,23.9,23.9,23.9,23.9,24.9,24.9,24.9,24.9,25.9,25.9,25.9,25.9,31.7,40.8,23.9,23.9,23.9,23.9,24.9,24.9,24.9,24.9,25.9,25.9,25.9,25.9,37.7,46.8,23.9,23.9,23.9,23.9,24.9,24.9,24.9,24.9,25.9,25.9,25.9,25.9,43.7,52.8,23.9,23.9,23.9,23.9,24.9,24.9,24.9,24.9,25.9,25.9,25.9,25.9,49.7,58.8,30.9,30.9,30.9,30.9,31.9,31.9,31.9,31.9,32.9,32.9,32.9,32.9,7.6,64.8,30.9,30.9,30.9,30.9,31.9,31.9,31.9,31.9,32.9,32.9,32.9,32.9,13.6,70.8,30.9,30.9,30.9,30.9,31.9,31.9,31.9,31.9,32.9,32.9,32.9,32.9,19.6,76.8,30.9,30.9,30.9,30.9,31.9,31.9,31.9,31.9,32.9,32.9,32.9,32.9,25.6,36.9,37.9,37.9,37.9,37.9,38.9,38.9,38.9,38.9,39.9,39.9,39.9,39.9,40.9])
+ da0Exp2=DataArrayDouble([15.9,16.9,16.9,16.9,16.9,17.9,17.9,17.9,17.9,18.9,18.9,18.9,18.9,19.9,22.9,15.2,16.2,17.2,18.2,19.2,20.2,21.2,22.2,23.2,24.2,25.2,26.2,26.9,22.9,29.2,30.2,31.2,32.2,33.2,34.2,35.2,36.2,37.2,38.2,39.2,40.2,26.9,22.9,43.2,44.2,45.2,46.2,47.2,48.2,49.2,50.2,51.2,52.2,53.2,54.2,26.9,22.9,57.2,58.2,59.2,60.2,61.2,62.2,63.2,64.2,65.2,66.2,67.2,68.2,26.9,29.9,71.2,72.2,73.2,74.2,75.2,76.2,77.2,78.2,79.2,80.2,81.2,82.2,33.9,29.9,85.2,86.2,87.2,88.2,89.2,90.2,91.2,92.2,93.2,94.2,95.2,96.2,33.9,29.9,99.2,100.2,101.2,102.2,103.2,104.2,105.2,106.2,107.2,108.2,109.2,110.2,33.9,29.9,113.2,114.2,115.2,116.2,117.2,118.2,119.2,120.2,121.2,122.2,123.2,124.2,33.9,36.9,37.9,37.9,37.9,37.9,38.9,38.9,38.9,38.9,39.9,39.9,39.9,39.9,40.9])
+ self.assertTrue(da0.isEqual(da0Exp,1e-12))
+ self.assertTrue(l2[0][:,0].isEqual(da0Exp,1e-12))
+ self.assertTrue(l2[0][:,1].isEqual(3*da0Exp,1e-12))
+ self.assertTrue(lCpy[0].isEqual(da0Exp2,1e-12))
+ #
+ g0=amr.retrieveGridsAt(0)
+ self.assertEqual(1,len(g0))
+ self.assertTrue(isinstance(g0[0],MEDCouplingCartesianAMRPatchGF))
+ g1=amr.retrieveGridsAt(1)
+ self.assertEqual(5,len(g1))
+ for i in xrange(5):
+ self.assertTrue(isinstance(g1[i],MEDCouplingCartesianAMRPatch))
+ pass
+ pass
+
+ def testSwig2AMR7(self):
+ """Idem testSwig2AMR6 except that we are in 1D"""
+ amr=MEDCouplingCartesianAMRMesh("",1,[6],[0],[1])
+ da=DataArrayDouble(5+2) ; da.iota() ; da+=0.9
+ amr.addPatch([(1,4)],[4])
+ amr.addPatch([(0,1)],[4])
+ da0=DataArrayDouble(3*4+2) ; da0.iota() ; da0[:]+=0.2
+ da1=DataArrayDouble(1*4+2) ; da1.iota() ; da1[:]+=0.4
+ self.assertEqual(2,amr.getNumberOfPatches())
+ l=[da0,da1]
+ lCpy=[elt.deepCpy() for elt in l]
+ l2=[DataArrayDouble.Meld(elt,3*elt) for elt in l]
+ amr.fillCellFieldOnPatchGhostAdv(0,da,1,l)
+ amr.fillCellFieldOnPatchGhostAdv(0,DataArrayDouble.Meld(da,3*da),1,l2)
+ amr.fillCellFieldOnPatchOnlyOnGhostZone(0,da,lCpy[0],1)
+ #
+ f=MEDCouplingFieldDouble(ON_CELLS) ; f.setMesh(amr.getImageMesh().buildWithGhost(1)) ; f.setArray(da) ; f.setName("all")
+ f0=MEDCouplingFieldDouble(ON_CELLS) ; f0.setMesh(amr[0].getMesh().getImageMesh().buildWithGhost(1)) ; f0.setArray(da0) ; f0.setName("p0") ; f0.checkCoherency()
+ f1=MEDCouplingFieldDouble(ON_CELLS) ; f1.setMesh(amr[1].getMesh().getImageMesh().buildWithGhost(1)) ; f1.setArray(da1) ; f1.setName("p1") ; f1.checkCoherency()
+ #
+ da0Exp=DataArrayDouble([4.4,2.9,2.9,2.9,2.9,3.9,3.9,3.9,3.9,4.9,4.9,4.9,4.9,5.9])
+ da0Exp2=DataArrayDouble([1.9,1.2,2.2,3.2,4.2,5.2,6.2,7.2,8.2,9.2,10.2,11.2,12.2,5.9])
+ self.assertTrue(da0.isEqual(da0Exp,1e-12))
+ self.assertTrue(l2[0][:,0].isEqual(da0Exp,1e-12))
+ self.assertTrue(l2[0][:,1].isEqual(3*da0Exp,1e-12))
+ self.assertTrue(lCpy[0].isEqual(da0Exp2,1e-12))
pass
def setUp(self):
