src/MEDCoupling_Swig/MEDCouplingRemapperTest.py

   1 #  -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2007-2012  CEA/DEN, EDF R&D
   3 #
   4 # This library is free software; you can redistribute it and/or
   5 # modify it under the terms of the GNU Lesser General Public
   6 # License as published by the Free Software Foundation; either
   7 # version 2.1 of the License.
   8 #
   9 # This library is distributed in the hope that it will be useful,
  10 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 # Lesser General Public License for more details.
  13 #
  14 # You should have received a copy of the GNU Lesser General Public
  15 # License along with this library; if not, write to the Free Software
  16 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 #
  18 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  19 #
  20
  21 from MEDCouplingRemapper import *
  22 from math import *
  23 import unittest
  24
  25 class MEDCouplingBasicsTest(unittest.TestCase):
  26     def testRemapper1(self):
  27         sourceMesh=self.build2DSourceMesh_1();
  28         targetMesh=self.build2DTargetMesh_1();
  29         remapper=MEDCouplingRemapper()
  30         remapper.setPrecision(1e-12);
  31         remapper.setIntersectionType(Triangulation);
  32         self.failUnless(remapper.prepare(sourceMesh,targetMesh,"P0P0")==1);
  33         srcField=MEDCouplingFieldDouble.New(ON_CELLS);
  34         srcField.setNature(ConservativeVolumic);
  35         srcField.setMesh(sourceMesh);
  36         array=DataArrayDouble.New();
  37         ptr=sourceMesh.getNumberOfCells()*[None]
  38         for i in xrange(sourceMesh.getNumberOfCells()):
  39             ptr[i]=float(i+7)
  40             pass
  41         array.setValues(ptr,sourceMesh.getNumberOfCells(),1);
  42         srcField.setArray(array);
  43         srcField.setName("abc") ; srcField.setDescription("def")
  44         srcField.setTime(7.7,9,10)
  45         trgfield=remapper.transferField(srcField,4.57);
  46         self.assertEqual("abc",trgfield.getName())
  47         self.assertEqual("def",trgfield.getDescription())
  48         a,b,c=trgfield.getTime()
  49         self.assertAlmostEqual(7.7,a,14)
  50         self.assertEqual(b,9)
  51         self.assertEqual(c,10)
  52         values=trgfield.getArray().getValues();
  53         valuesExpected=[7.5 ,7. ,7.,8.,7.5];
  54         for i in xrange(targetMesh.getNumberOfCells()):
  55             self.failUnless(abs(values[i]-valuesExpected[i])<1e-12);
  56             pass
  57         self.failUnless(1==trgfield.getArray().getNumberOfComponents());
  58         pass
  59
  60     def testPrepareEx1(self):
  61         sourceMesh=self.build2DSourceMesh_1();
  62         targetMesh=self.build2DTargetMesh_3();
  63         #
  64         remapper=MEDCouplingRemapper();
  65         remapper.setPrecision(1e-12);
  66         remapper.setIntersectionType(Triangulation);
  67         srcFt=MEDCouplingFieldTemplate.New(ON_CELLS);
  68         trgFt=MEDCouplingFieldTemplate.New(ON_CELLS);
  69         srcFt.setMesh(sourceMesh);
  70         trgFt.setMesh(targetMesh);
  71         self.assertEqual(1,remapper.prepareEx(srcFt,trgFt));
  72         srcField=MEDCouplingFieldDouble.New(ON_CELLS);
  73         srcField.setNature(ConservativeVolumic);
  74         srcField.setMesh(sourceMesh);
  75         array=DataArrayDouble.New();
  76         ptr=sourceMesh.getNumberOfCells()*[None]
  77         for i in xrange(sourceMesh.getNumberOfCells()):
  78             ptr[i]=float(i+7);
  79             pass
  80         array.setValues(ptr,sourceMesh.getNumberOfCells(),1);
  81         srcField.setArray(array);
  82         trgfield=remapper.transferField(srcField,4.220173);
  83         values=trgfield.getArray().getValues();
  84         valuesExpected=[7.75, 7.0625, 4.220173,8.0]
  85         self.assertEqual(4,trgfield.getArray().getNumberOfTuples());
  86         self.assertEqual(1,trgfield.getArray().getNumberOfComponents());
  87         for i0 in xrange(4):
  88             self.assertAlmostEqual(valuesExpected[i0],values[i0],12);
  89             pass
  90         pass
  91
  92     def testPartialTransfer1(self):
  93         sourceMesh=self.build2DSourceMesh_1();
  94         targetMesh=self.build2DTargetMesh_3();
  95         #
  96         remapper=MEDCouplingRemapper();
  97         remapper.setPrecision(1e-12);
  98         remapper.setIntersectionType(Triangulation);
  99         srcFt=MEDCouplingFieldTemplate.New(ON_CELLS);
 100         trgFt=MEDCouplingFieldTemplate.New(ON_CELLS);
 101         srcFt.setMesh(sourceMesh);
 102         trgFt.setMesh(targetMesh);
 103         self.assertEqual(1,remapper.prepareEx(srcFt,trgFt));
 104         srcField=MEDCouplingFieldDouble.New(ON_CELLS);
 105         srcField.setNature(ConservativeVolumic);
 106         srcField.setMesh(sourceMesh);
 107         array=DataArrayDouble.New();
 108         ptr=sourceMesh.getNumberOfCells()*[None]
 109         for i in xrange(sourceMesh.getNumberOfCells()):
 110             ptr[i]=float(i+7);
 111             pass
 112         array.setValues(ptr,sourceMesh.getNumberOfCells(),1);
 113         srcField.setArray(array);
 114         trgfield=MEDCouplingFieldDouble.New(ON_CELLS);
 115         trgfield.setNature(ConservativeVolumic);
 116         trgfield.setMesh(targetMesh);
 117         array=DataArrayDouble.New();
 118         ptr=targetMesh.getNumberOfCells()*[None]
 119         for i in xrange(targetMesh.getNumberOfCells()):
 120             ptr[i]=4.220173;
 121             pass
 122         array.setValues(ptr,targetMesh.getNumberOfCells(),1);
 123         trgfield.setArray(array);
 124         remapper.partialTransfer(srcField,trgfield);
 125         values=trgfield.getArray().getValues();
 126         valuesExpected=[7.75, 7.0625, 4.220173,8.0]
 127         self.assertEqual(4,trgfield.getArray().getNumberOfTuples());
 128         self.assertEqual(1,trgfield.getArray().getNumberOfComponents());
 129         for i0 in xrange(4):
 130             self.assertAlmostEqual(valuesExpected[i0],values[i0],12);
 131             pass
 132         pass
 133
 134     def testPrepareUC(self):
 135         # 1D
 136         coords=DataArrayDouble([0.,0.5,0.7])
 137         src=MEDCouplingUMesh("",1) ; src.setCoords(coords)
 138         src.allocateCells(2) ; src.insertNextCell(NORM_SEG2,[0,1]) ; src.insertNextCell(NORM_SEG2,[1,2]) ; src.finishInsertingCells()
 139         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 140         trg.setCoordsAt(0,arr)
 141         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
 142         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 143         rem=MEDCouplingRemapper()
 144         rem.prepare(src,trg,"P0P0")
 145         trgField=rem.transferField(fieldSrc,-7.)
 146         expected1=[-7.,4.,36.,-7.,-7.]
 147         self.assertEqual(5,trgField.getArray().getNumberOfTuples())
 148         self.assertEqual(5,len(expected1))
 149         for i,val in enumerate(expected1):
 150             self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
 151             pass
 152         # 2D
 153         coords=DataArrayDouble([0.,0.,0.,1.,1.,1.,1.,0.,0.5,-0.2],5,2)
 154         src=MEDCouplingUMesh("",2) ; src.setCoords(coords)
 155         src.allocateCells(2) ; src.insertNextCell(NORM_TRI3,[0,1,2]) ; src.insertNextCell(NORM_TRI3,[3,4,0]) ; src.finishInsertingCells()
 156         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 157         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr)
 158         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
 159         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 160         rem=MEDCouplingRemapper()
 161         rem.prepare(src,trg,"P0P0")
 162         trgField=rem.transferField(fieldSrc,-7.)
 163         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.]
 164         self.assertEqual(25,trgField.getArray().getNumberOfTuples())
 165         self.assertEqual(25,len(expected2))
 166         for i,val in enumerate(expected2):
 167             self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
 168             pass
 169         # 3D
 170         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)
 171         src=MEDCouplingUMesh("",3) ; src.setCoords(coords)
 172         src.allocateCells(2) ; src.insertNextCell(NORM_TETRA4,[0,1,2,5]) ; src.insertNextCell(NORM_TETRA4,[3,4,0,6]) ; src.finishInsertingCells()
 173         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.])
 174         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr) ; trg.setCoordsAt(2,arr2)
 175         src.checkCoherency2(1e-10)
 176         trg.checkCoherency()
 177         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
 178         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 179         rem=MEDCouplingRemapper()
 180         rem.prepare(src,trg,"P0P0")
 181         trgField=rem.transferField(fieldSrc,-7.)
 182         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.]
 183         self.assertEqual(100,trgField.getArray().getNumberOfTuples())
 184         self.assertEqual(100,len(expected3))
 185         for i,val in enumerate(expected3):
 186             self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
 187             pass
 188         pass
 189
 190     def testPrepareCU(self):
 191         # 1D
 192         coords=DataArrayDouble([0.,0.5,0.7])
 193         trg=MEDCouplingUMesh("",1) ; trg.setCoords(coords)
 194         trg.allocateCells(2) ; trg.insertNextCell(NORM_SEG2,[0,1]) ; trg.insertNextCell(NORM_SEG2,[1,2]) ; trg.finishInsertingCells()
 195         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 196         src.setCoordsAt(0,arr)
 197         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrTrg=DataArrayDouble([10.,30.,40.,70.,80.])
 198         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrTrg)
 199         rem=MEDCouplingRemapper()
 200         rem.prepare(src,trg,"P0P0")
 201         trgField=rem.transferField(fieldSrc,-7.)
 202         expected1=[44.,16.]
 203         self.assertEqual(2.,trgField.getArray().getNumberOfTuples())
 204         self.assertEqual(2,len(expected1))
 205         for i,val in enumerate(expected1):
 206             self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
 207             pass
 208         # 2D
 209         coords=DataArrayDouble([0.,0.,0.,1.,1.,1.,1.,0.,0.5,-0.2],5,2)
 210         trg=MEDCouplingUMesh("",2) ; trg.setCoords(coords)
 211         trg.allocateCells(2) ; trg.insertNextCell(NORM_TRI3,[0,1,2]) ; trg.insertNextCell(NORM_TRI3,[3,4,0]) ; trg.finishInsertingCells()
 212         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 213         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr)
 214         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.])
 215         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 216         rem=MEDCouplingRemapper()
 217         rem.prepare(src,trg,"P0P0")
 218         trgField=rem.transferField(fieldSrc,-7.)
 219         expected2=[441.3050624589086,68.69529914529915]
 220         self.assertEqual(2,trgField.getArray().getNumberOfTuples())
 221         self.assertEqual(2,len(expected2))
 222         for i,val in enumerate(expected2):
 223             self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
 224             pass
 225         # 3D
 226         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)
 227         trg=MEDCouplingUMesh("",3) ; trg.setCoords(coords)
 228         trg.allocateCells(2) ; trg.insertNextCell(NORM_TETRA4,[0,1,2,5]) ; trg.insertNextCell(NORM_TETRA4,[3,4,0,6]) ; trg.finishInsertingCells()
 229         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.])
 230         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr) ; src.setCoordsAt(2,arr2)
 231         trg.checkCoherency2(1e-10)
 232         src.checkCoherency()
 233         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble(100) ; arrSrc.iota(7.7)
 234         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 235         rem=MEDCouplingRemapper()
 236         rem.prepare(src,trg,"P0P0")
 237         trgField=rem.transferField(fieldSrc,-7.)
 238         expected3=[39.635196634558845,12.13422356758468]
 239         self.assertEqual(2,trgField.getArray().getNumberOfTuples())
 240         self.assertEqual(2,len(expected3))
 241         for i,val in enumerate(expected3):
 242             self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
 243             pass
 244         pass
 245
 246     def testPrepareCC(self):
 247         # 1D
 248         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 249         src.setCoordsAt(0,arr)
 250         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
 251         trg.setCoordsAt(0,arr)
 252         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrTrg=DataArrayDouble([10.,30.,40.,70.,80.])
 253         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrTrg)
 254         rem=MEDCouplingRemapper()
 255         rem.prepare(src,trg,"P0P0")
 256         trgField=rem.transferField(fieldSrc,-7.)
 257         expected1=[10.,25.]
 258         self.assertEqual(2.,trgField.getArray().getNumberOfTuples())
 259         self.assertEqual(2,len(expected1))
 260         for i,val in enumerate(expected1):
 261             self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
 262             pass
 263         # 2D
 264         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 265         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr)
 266         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
 267         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr)
 268         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.])
 269         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 270         rem=MEDCouplingRemapper()
 271         rem.prepare(src,trg,"P0P0")
 272         trgField=rem.transferField(fieldSrc,-7.)
 273         expected2=[10.,25.,91.66666666666666,90.27777777777777]
 274         self.assertEqual(4,trgField.getArray().getNumberOfTuples())
 275         self.assertEqual(4,len(expected2))
 276         for i,val in enumerate(expected2):
 277             self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
 278             pass
 279         # 3D
 280         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 281         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr) ; src.setCoordsAt(2,arr)
 282         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
 283         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr) ; trg.setCoordsAt(2,arr)
 284         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble(125) ; arrSrc.iota(7.7)
 285         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc) ; fieldSrc.checkCoherency()
 286         rem=MEDCouplingRemapper()
 287         rem.prepare(src,trg,"P0P0")
 288         trgField=rem.transferField(fieldSrc,-7.)
 289         expected3=[7.7, 7.249999999999999, 10.583333333333332, 9.513888888888886, 27.25, 23.40277777777777, 26.180555555555546, 22.39583333333333]
 290         self.assertEqual(8,trgField.getArray().getNumberOfTuples())
 291         self.assertEqual(8,len(expected3))
 292         for i,val in enumerate(expected3):
 293             self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
 294             pass
 295         pass
 296
 297     # Bug when source mesh is not homogeneously oriented in source mesh
 298     def testNonRegressionNonHomegenousOrriented3DCells(self):
 299         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)
 300         src1=MEDCouplingUMesh("src",3) ; src1.allocateCells(0) ; src1.insertNextCell(NORM_TETRA4,[0,1,4,3]) ; src1.insertNextCell(NORM_TETRA4,[2,0,4,3])
 301         src2=MEDCouplingUMesh("src",3) ; src2.allocateCells(0) ; src2.insertNextCell(NORM_TETRA4,[0,4,1,3]) ; src2.insertNextCell(NORM_TETRA4,[2,0,4,3])
 302         src1.setCoords(csrc) ; src2.setCoords(csrc)
 303         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)
 304         trg=MEDCouplingUMesh("trg",3) ; trg.allocateCells(0) ; trg.insertNextCell(NORM_HEXA8,[0,1,3,2,4,5,7,6])
 305         trg.setCoords(ctrg)
 306         rem1=MEDCouplingRemapper() ; rem1.setSplittingPolicy(PLANAR_FACE_5) ; rem1.prepare(src1,trg,"P0P0")
 307         rem2=MEDCouplingRemapper() ; rem2.setSplittingPolicy(PLANAR_FACE_5) ; rem2.prepare(src1,trg,"P0P0")
 308         mat1=rem1.getCrudeMatrix() ; mat2=rem2.getCrudeMatrix()
 309         self.assertEqual(1,len(mat1)) ; self.assertEqual(1,len(mat2))
 310         self.assertEqual(mat1[0].keys(),mat2[0].keys()) ; self.assertEqual([0,1],mat1[0].keys())
 311         self.assertAlmostEqual(1.25884108122e-06,mat1[0][0],16) ; self.assertAlmostEqual(1.25884108122e-06,mat2[0][0],16)
 312         self.assertAlmostEqual(1.25884086663e-06,mat1[0][1],16) ; self.assertAlmostEqual(1.25884086663e-06,mat2[0][1],16)
 313         #
 314         d=DataArrayDouble([13.45,27.67],2,1)
 315         f1=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f1.setMesh(src1) ; f1.setArray(d) ; f1.setNature(RevIntegral)
 316         f2=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f2.setMesh(src2) ; f2.setArray(d) ; f2.setNature(RevIntegral)
 317         f11=rem1.transferField(f1,1e300) ; f22=rem2.transferField(f2,1e300)
 318         expected1=DataArrayDouble([0.012480539537637884])
 319         self.assertTrue(f11.getArray().isEqual(expected1,1e-15))
 320         self.assertTrue(f22.getArray().isEqual(expected1,1e-15))
 321         #
 322         f1.setNature(Integral) ; f2.setNature(Integral)
 323         f11=rem1.transferField(f1,1e300) ; f22=rem2.transferField(f2,1e300)
 324         #
 325         expected2=DataArrayDouble([41.12])
 326         self.assertTrue(f11.getArray().isEqual(expected2,1e-13))
 327         self.assertTrue(f22.getArray().isEqual(expected2,1e-13))
 328         pass
 329
 330     def testCellToNodeReverse3D(self):
 331         c=DataArrayDouble([0.,1.,2.5])
 332         cc=MEDCouplingCMesh()
 333         cc.setCoords(c,c,c)
 334         um=cc.buildUnstructured()
 335         f=um.getMeasureField(ON_CELLS)
 336         #
 337         n2o=um.simplexize(PLANAR_FACE_5)
 338         f.setArray(f.getArray()[n2o])
 339         f.checkCoherency()
 340         f.setNature(ConservativeVolumic)
 341         f.setTime(5.6,7,8)
 342         f.setName("toto") ; f.setDescription("aDescription")
 343         p=MEDCouplingRemapper()
 344         p.setP1P0BaryMethod(True)
 345         p.prepare(um,um,"P1P0")
 346         fNode=p.reverseTransferField(f,1e300)
 347         self.assertEqual("toto",fNode.getName())
 348         self.assertEqual("aDescription",fNode.getDescription())
 349         a,b,c=fNode.getTime()
 350         self.assertAlmostEqual(5.6,a,14)
 351         self.assertEqual(7,b) ; self.assertEqual(8,c)
 352         #
 353         integExpected=34.328125
 354         self.assertAlmostEqual(fNode.integral(False)[0],integExpected,14)
 355         self.assertAlmostEqual(f.integral(False)[0],integExpected,14)
 356         pass
 357
 358     def build2DSourceMesh_1(self):
 359         sourceCoords=[-0.3,-0.3, 0.7,-0.3, -0.3,0.7, 0.7,0.7]
 360         sourceConn=[0,3,1,0,2,3]
 361         sourceMesh=MEDCouplingUMesh.New("my name of mesh 2D",2)
 362         sourceMesh.allocateCells(2);
 363         sourceMesh.insertNextCell(NORM_TRI3,3,sourceConn[0:3]);
 364         sourceMesh.insertNextCell(NORM_TRI3,3,sourceConn[3:6]);
 365         sourceMesh.finishInsertingCells();
 366         myCoords=DataArrayDouble.New();
 367         myCoords.setValues(sourceCoords,4,2);
 368         sourceMesh.setCoords(myCoords);
 369         return sourceMesh;
 370
 371     def build2DTargetMesh_1(self):
 372         targetCoords=[-0.3,-0.3, 0.2,-0.3, 0.7,-0.3, -0.3,0.2, 0.2,0.2, 0.7,0.2, -0.3,0.7, 0.2,0.7, 0.7,0.7 ]
 373         targetConn=[0,3,4,1, 1,4,2, 4,5,2, 6,7,4,3, 7,8,5,4]
 374         targetMesh=MEDCouplingUMesh.New();
 375         targetMesh.setMeshDimension(2);
 376         targetMesh.allocateCells(5);
 377         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[0:4]);
 378         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[4:7]);
 379         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[7:10]);
 380         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[10:14]);
 381         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[14:18]);
 382         targetMesh.finishInsertingCells();
 383         myCoords=DataArrayDouble.New();
 384         myCoords.setValues(targetCoords,9,2);
 385         targetMesh.setCoords(myCoords);
 386         return targetMesh;
 387
 388     def build2DTargetMesh_3(self):
 389         targetCoords=[-0.6,-0.4, -0.1,-0.4, 1.1,-0.4, 2.1,-0.4, -0.6,0.1,  -0.1,0.1,  1.1,0.1,  2.1,0.1, -0.6,1.1,  -0.1,1.1]
 390         targetConn=[0,4,5,1, 1,5,6,2, 2,6,7,3, 4,8,9,5]
 391         targetMesh=MEDCouplingUMesh.New();
 392         targetMesh.setMeshDimension(2);
 393         targetMesh.allocateCells(4);
 394         for i in xrange(4):
 395             targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[4*i:4*(i+1)])
 396             pass
 397         targetMesh.finishInsertingCells();
 398         myCoords=DataArrayDouble.New();
 399         myCoords.setValues(targetCoords,10,2);
 400         targetMesh.setCoords(myCoords);
 401         return targetMesh;
 402         pass
 403
 404     def setUp(self):
 405         pass
 406     pass
 407
 408 unittest.main()