src/MEDCoupling_Swig/MEDCouplingRemapperTest.py

   1 #  -*- coding: iso-8859-1 -*-
   2 # Copyright (C) 2007-2015  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, or (at your option) any later version.
   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 MEDCouplingDataForTest import MEDCouplingDataForTest
  22 from MEDCouplingRemapper import *
  23 from math import *
  24 import unittest
  25
  26 class MEDCouplingBasicsTest(unittest.TestCase):
  27     def testRemapper1(self):
  28         sourceMesh=self.build2DSourceMesh_1();
  29         targetMesh=self.build2DTargetMesh_1();
  30         remapper=MEDCouplingRemapper()
  31         remapper.setPrecision(1e-12);
  32         remapper.setIntersectionType(Triangulation);
  33         self.failUnless(remapper.prepare(sourceMesh,targetMesh,"P0P0")==1);
  34         srcField=MEDCouplingFieldDouble.New(ON_CELLS);
  35         srcField.setNature(ConservativeVolumic);
  36         srcField.setMesh(sourceMesh);
  37         array=DataArrayDouble.New();
  38         ptr=sourceMesh.getNumberOfCells()*[None]
  39         for i in xrange(sourceMesh.getNumberOfCells()):
  40             ptr[i]=float(i+7)
  41             pass
  42         array.setValues(ptr,sourceMesh.getNumberOfCells(),1);
  43         srcField.setArray(array);
  44         srcField.setName("abc") ; srcField.setDescription("def")
  45         srcField.setTime(7.7,9,10)
  46         trgfield=remapper.transferField(srcField,4.57);
  47         self.assertEqual("abc",trgfield.getName())
  48         self.assertEqual("def",trgfield.getDescription())
  49         a,b,c=trgfield.getTime()
  50         self.assertAlmostEqual(7.7,a,14)
  51         self.assertEqual(b,9)
  52         self.assertEqual(c,10)
  53         values=trgfield.getArray().getValues();
  54         valuesExpected=[7.5 ,7. ,7.,8.,7.5];
  55         for i in xrange(targetMesh.getNumberOfCells()):
  56             self.failUnless(abs(values[i]-valuesExpected[i])<1e-12);
  57             pass
  58         self.failUnless(1==trgfield.getArray().getNumberOfComponents());
  59         pass
  60
  61     def testPrepareEx1(self):
  62         sourceMesh=self.build2DSourceMesh_1();
  63         targetMesh=self.build2DTargetMesh_3();
  64         #
  65         remapper=MEDCouplingRemapper();
  66         remapper.setPrecision(1e-12);
  67         remapper.setIntersectionType(Triangulation);
  68         srcFt=MEDCouplingFieldTemplate.New(ON_CELLS);
  69         trgFt=MEDCouplingFieldTemplate.New(ON_CELLS);
  70         srcFt.setMesh(sourceMesh);
  71         trgFt.setMesh(targetMesh);
  72         self.assertEqual(1,remapper.prepareEx(srcFt,trgFt));
  73         srcField=MEDCouplingFieldDouble.New(ON_CELLS);
  74         srcField.setNature(ConservativeVolumic);
  75         srcField.setMesh(sourceMesh);
  76         array=DataArrayDouble.New();
  77         ptr=sourceMesh.getNumberOfCells()*[None]
  78         for i in xrange(sourceMesh.getNumberOfCells()):
  79             ptr[i]=float(i+7);
  80             pass
  81         array.setValues(ptr,sourceMesh.getNumberOfCells(),1);
  82         srcField.setArray(array);
  83         trgfield=remapper.transferField(srcField,4.220173);
  84         values=trgfield.getArray().getValues();
  85         valuesExpected=[7.75, 7.0625, 4.220173,8.0]
  86         self.assertEqual(4,trgfield.getArray().getNumberOfTuples());
  87         self.assertEqual(1,trgfield.getArray().getNumberOfComponents());
  88         for i0 in xrange(4):
  89             self.assertAlmostEqual(valuesExpected[i0],values[i0],12);
  90             pass
  91         pass
  92
  93     def testPartialTransfer1(self):
  94         sourceMesh=self.build2DSourceMesh_1();
  95         targetMesh=self.build2DTargetMesh_3();
  96         #
  97         remapper=MEDCouplingRemapper();
  98         remapper.setPrecision(1e-12);
  99         remapper.setIntersectionType(Triangulation);
 100         srcFt=MEDCouplingFieldTemplate.New(ON_CELLS);
 101         trgFt=MEDCouplingFieldTemplate.New(ON_CELLS);
 102         srcFt.setMesh(sourceMesh);
 103         trgFt.setMesh(targetMesh);
 104         self.assertEqual(1,remapper.prepareEx(srcFt,trgFt));
 105         srcField=MEDCouplingFieldDouble.New(ON_CELLS);
 106         srcField.setNature(ConservativeVolumic);
 107         srcField.setMesh(sourceMesh);
 108         array=DataArrayDouble.New();
 109         ptr=sourceMesh.getNumberOfCells()*[None]
 110         for i in xrange(sourceMesh.getNumberOfCells()):
 111             ptr[i]=float(i+7);
 112             pass
 113         array.setValues(ptr,sourceMesh.getNumberOfCells(),1);
 114         srcField.setArray(array);
 115         trgfield=MEDCouplingFieldDouble.New(ON_CELLS);
 116         trgfield.setNature(ConservativeVolumic);
 117         trgfield.setMesh(targetMesh);
 118         array=DataArrayDouble.New();
 119         ptr=targetMesh.getNumberOfCells()*[None]
 120         for i in xrange(targetMesh.getNumberOfCells()):
 121             ptr[i]=4.220173;
 122             pass
 123         array.setValues(ptr,targetMesh.getNumberOfCells(),1);
 124         trgfield.setArray(array);
 125         remapper.partialTransfer(srcField,trgfield);
 126         values=trgfield.getArray().getValues();
 127         valuesExpected=[7.75, 7.0625, 4.220173,8.0]
 128         self.assertEqual(4,trgfield.getArray().getNumberOfTuples());
 129         self.assertEqual(1,trgfield.getArray().getNumberOfComponents());
 130         for i0 in xrange(4):
 131             self.assertAlmostEqual(valuesExpected[i0],values[i0],12);
 132             pass
 133         pass
 134
 135     def testPrepareUC(self):
 136         # 1D
 137         coords=DataArrayDouble([0.,0.5,0.7])
 138         src=MEDCouplingUMesh("",1) ; src.setCoords(coords)
 139         src.allocateCells(2) ; src.insertNextCell(NORM_SEG2,[0,1]) ; src.insertNextCell(NORM_SEG2,[1,2]) ; src.finishInsertingCells()
 140         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 141         trg.setCoordsAt(0,arr)
 142         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
 143         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 144         rem=MEDCouplingRemapper()
 145         rem.prepare(src,trg,"P0P0")
 146         trgField=rem.transferField(fieldSrc,-7.)
 147         expected1=[-7.,4.,36.,-7.,-7.]
 148         self.assertEqual(5,trgField.getArray().getNumberOfTuples())
 149         self.assertEqual(5,len(expected1))
 150         for i,val in enumerate(expected1):
 151             self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
 152             pass
 153         # 2D
 154         coords=DataArrayDouble([0.,0.,0.,1.,1.,1.,1.,0.,0.5,-0.2],5,2)
 155         src=MEDCouplingUMesh("",2) ; src.setCoords(coords)
 156         src.allocateCells(2) ; src.insertNextCell(NORM_TRI3,[0,1,2]) ; src.insertNextCell(NORM_TRI3,[3,4,0]) ; src.finishInsertingCells()
 157         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 158         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr)
 159         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
 160         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 161         rem=MEDCouplingRemapper()
 162         rem.prepare(src,trg,"P0P0")
 163         trgField=rem.transferField(fieldSrc,-7.)
 164         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.]
 165         self.assertEqual(25,trgField.getArray().getNumberOfTuples())
 166         self.assertEqual(25,len(expected2))
 167         for i,val in enumerate(expected2):
 168             self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
 169             pass
 170         # 3D
 171         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)
 172         src=MEDCouplingUMesh("",3) ; src.setCoords(coords)
 173         src.allocateCells(2) ; src.insertNextCell(NORM_TETRA4,[0,1,2,5]) ; src.insertNextCell(NORM_TETRA4,[3,4,0,6]) ; src.finishInsertingCells()
 174         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.])
 175         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr) ; trg.setCoordsAt(2,arr2)
 176         src.checkCoherency2(1e-10)
 177         trg.checkCoherency()
 178         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble([10.,30.])
 179         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 180         rem=MEDCouplingRemapper()
 181         rem.prepare(src,trg,"P0P0")
 182         trgField=rem.transferField(fieldSrc,-7.)
 183         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.]
 184         self.assertEqual(100,trgField.getArray().getNumberOfTuples())
 185         self.assertEqual(100,len(expected3))
 186         for i,val in enumerate(expected3):
 187             self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
 188             pass
 189         pass
 190
 191     def testPrepareCU(self):
 192         # 1D
 193         coords=DataArrayDouble([0.,0.5,0.7])
 194         trg=MEDCouplingUMesh("",1) ; trg.setCoords(coords)
 195         trg.allocateCells(2) ; trg.insertNextCell(NORM_SEG2,[0,1]) ; trg.insertNextCell(NORM_SEG2,[1,2]) ; trg.finishInsertingCells()
 196         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 197         src.setCoordsAt(0,arr)
 198         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrTrg=DataArrayDouble([10.,30.,40.,70.,80.])
 199         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrTrg)
 200         rem=MEDCouplingRemapper()
 201         rem.prepare(src,trg,"P0P0")
 202         trgField=rem.transferField(fieldSrc,-7.)
 203         expected1=[44.,16.]
 204         self.assertEqual(2.,trgField.getArray().getNumberOfTuples())
 205         self.assertEqual(2,len(expected1))
 206         for i,val in enumerate(expected1):
 207             self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
 208             pass
 209         # 2D
 210         coords=DataArrayDouble([0.,0.,0.,1.,1.,1.,1.,0.,0.5,-0.2],5,2)
 211         trg=MEDCouplingUMesh("",2) ; trg.setCoords(coords)
 212         trg.allocateCells(2) ; trg.insertNextCell(NORM_TRI3,[0,1,2]) ; trg.insertNextCell(NORM_TRI3,[3,4,0]) ; trg.finishInsertingCells()
 213         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 214         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr)
 215         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.])
 216         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 217         rem=MEDCouplingRemapper()
 218         rem.prepare(src,trg,"P0P0")
 219         trgField=rem.transferField(fieldSrc,-7.)
 220         expected2=[441.3050624589086,68.69529914529915]
 221         self.assertEqual(2,trgField.getArray().getNumberOfTuples())
 222         self.assertEqual(2,len(expected2))
 223         for i,val in enumerate(expected2):
 224             self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
 225             pass
 226         # 3D
 227         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)
 228         trg=MEDCouplingUMesh("",3) ; trg.setCoords(coords)
 229         trg.allocateCells(2) ; trg.insertNextCell(NORM_TETRA4,[0,1,2,5]) ; trg.insertNextCell(NORM_TETRA4,[3,4,0,6]) ; trg.finishInsertingCells()
 230         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.])
 231         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr) ; src.setCoordsAt(2,arr2)
 232         trg.checkCoherency2(1e-10)
 233         src.checkCoherency()
 234         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble(100) ; arrSrc.iota(7.7)
 235         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 236         rem=MEDCouplingRemapper()
 237         rem.prepare(src,trg,"P0P0")
 238         trgField=rem.transferField(fieldSrc,-7.)
 239         expected3=[39.635196634558845,12.13422356758468]
 240         self.assertEqual(2,trgField.getArray().getNumberOfTuples())
 241         self.assertEqual(2,len(expected3))
 242         for i,val in enumerate(expected3):
 243             self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
 244             pass
 245         pass
 246
 247     def testPrepareCC(self):
 248         # 1D
 249         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 250         src.setCoordsAt(0,arr)
 251         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
 252         trg.setCoordsAt(0,arr)
 253         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrTrg=DataArrayDouble([10.,30.,40.,70.,80.])
 254         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrTrg)
 255         rem=MEDCouplingRemapper()
 256         rem.prepare(src,trg,"P0P0")
 257         trgField=rem.transferField(fieldSrc,-7.)
 258         expected1=[10.,25.]
 259         self.assertEqual(2.,trgField.getArray().getNumberOfTuples())
 260         self.assertEqual(2,len(expected1))
 261         for i,val in enumerate(expected1):
 262             self.assertAlmostEqual(expected1[i],trgField.getArray().getIJ(i,0),12);
 263             pass
 264         # 2D
 265         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 266         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr)
 267         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
 268         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr)
 269         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.])
 270         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc)
 271         rem=MEDCouplingRemapper()
 272         rem.prepare(src,trg,"P0P0")
 273         trgField=rem.transferField(fieldSrc,-7.)
 274         expected2=[10.,25.,91.66666666666666,90.27777777777777]
 275         self.assertEqual(4,trgField.getArray().getNumberOfTuples())
 276         self.assertEqual(4,len(expected2))
 277         for i,val in enumerate(expected2):
 278             self.assertAlmostEqual(expected2[i],trgField.getArray().getIJ(i,0),12);
 279             pass
 280         # 3D
 281         src=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.7,-0.1,0.2,0.7,2.,2.3])
 282         src.setCoordsAt(0,arr) ; src.setCoordsAt(1,arr) ; src.setCoordsAt(2,arr)
 283         trg=MEDCouplingCMesh() ; arr=DataArrayDouble([-0.9,-0.1,0.15])
 284         trg.setCoordsAt(0,arr) ; trg.setCoordsAt(1,arr) ; trg.setCoordsAt(2,arr)
 285         fieldSrc=MEDCouplingFieldDouble(ON_CELLS,NO_TIME) ; fieldSrc.setMesh(src) ; arrSrc=DataArrayDouble(125) ; arrSrc.iota(7.7)
 286         fieldSrc.setNature(Integral) ;  fieldSrc.setArray(arrSrc) ; fieldSrc.checkCoherency()
 287         rem=MEDCouplingRemapper()
 288         rem.prepare(src,trg,"P0P0")
 289         trgField=rem.transferField(fieldSrc,-7.)
 290         expected3=[7.7, 7.249999999999999, 10.583333333333332, 9.513888888888886, 27.25, 23.40277777777777, 26.180555555555546, 22.39583333333333]
 291         self.assertEqual(8,trgField.getArray().getNumberOfTuples())
 292         self.assertEqual(8,len(expected3))
 293         for i,val in enumerate(expected3):
 294             self.assertAlmostEqual(expected3[i],trgField.getArray().getIJ(i,0),12);
 295             pass
 296         pass
 297
 298     # Bug when source mesh is not homogeneously oriented in source mesh
 299     def testNonRegressionNonHomegenousOrriented3DCells(self):
 300         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)
 301         src1=MEDCouplingUMesh("src",3) ; src1.allocateCells(0) ; src1.insertNextCell(NORM_TETRA4,[0,1,4,3]) ; src1.insertNextCell(NORM_TETRA4,[2,0,4,3])
 302         src2=MEDCouplingUMesh("src",3) ; src2.allocateCells(0) ; src2.insertNextCell(NORM_TETRA4,[0,4,1,3]) ; src2.insertNextCell(NORM_TETRA4,[2,0,4,3])
 303         src1.setCoords(csrc) ; src2.setCoords(csrc)
 304         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)
 305         trg=MEDCouplingUMesh("trg",3) ; trg.allocateCells(0) ; trg.insertNextCell(NORM_HEXA8,[0,1,3,2,4,5,7,6])
 306         trg.setCoords(ctrg)
 307         rem1=MEDCouplingRemapper() ; rem1.setSplittingPolicy(PLANAR_FACE_5) ; rem1.prepare(src1,trg,"P0P0")
 308         rem2=MEDCouplingRemapper() ; rem2.setSplittingPolicy(PLANAR_FACE_5) ; rem2.prepare(src1,trg,"P0P0")
 309         mat1=rem1.getCrudeMatrix() ; mat2=rem2.getCrudeMatrix()
 310         self.assertEqual(1,len(mat1)) ; self.assertEqual(1,len(mat2))
 311         self.assertEqual(mat1[0].keys(),mat2[0].keys()) ; self.assertEqual([0,1],mat1[0].keys())
 312         self.assertAlmostEqual(1.25884108122e-06,mat1[0][0],16) ; self.assertAlmostEqual(1.25884108122e-06,mat2[0][0],16)
 313         self.assertAlmostEqual(1.25884086663e-06,mat1[0][1],16) ; self.assertAlmostEqual(1.25884086663e-06,mat2[0][1],16)
 314         #
 315         d=DataArrayDouble([13.45,27.67],2,1)
 316         f1=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f1.setMesh(src1) ; f1.setArray(d) ; f1.setNature(RevIntegral)
 317         f2=MEDCouplingFieldDouble(ON_CELLS,ONE_TIME) ; f2.setMesh(src2) ; f2.setArray(d) ; f2.setNature(RevIntegral)
 318         f11=rem1.transferField(f1,1e300) ; f22=rem2.transferField(f2,1e300)
 319         expected1=DataArrayDouble([0.012480539537637884])
 320         self.assertTrue(f11.getArray().isEqual(expected1,1e-15))
 321         self.assertTrue(f22.getArray().isEqual(expected1,1e-15))
 322         #
 323         f1.setNature(Integral) ; f2.setNature(Integral)
 324         f11=rem1.transferField(f1,1e300) ; f22=rem2.transferField(f2,1e300)
 325         #
 326         expected2=DataArrayDouble([41.12])
 327         self.assertTrue(f11.getArray().isEqual(expected2,1e-13))
 328         self.assertTrue(f22.getArray().isEqual(expected2,1e-13))
 329         pass
 330
 331     def testCellToNodeReverse3D(self):
 332         c=DataArrayDouble([0.,1.,2.5])
 333         cc=MEDCouplingCMesh()
 334         cc.setCoords(c,c,c)
 335         um=cc.buildUnstructured()
 336         f=um.getMeasureField(ON_CELLS)
 337         #
 338         n2o=um.simplexize(PLANAR_FACE_5)
 339         f.setArray(f.getArray()[n2o])
 340         f.checkCoherency()
 341         f.setNature(ConservativeVolumic)
 342         f.setTime(5.6,7,8)
 343         f.setName("toto") ; f.setDescription("aDescription")
 344         p=MEDCouplingRemapper()
 345         p.setIntersectionType(Barycentric)
 346         p.prepare(um,um,"P1P0")
 347         fNode=p.reverseTransferField(f,1e300)
 348         self.assertEqual("toto",fNode.getName())
 349         self.assertEqual("aDescription",fNode.getDescription())
 350         a,b,c=fNode.getTime()
 351         self.assertAlmostEqual(5.6,a,14)
 352         self.assertEqual(7,b) ; self.assertEqual(8,c)
 353         #
 354         integExpected=34.328125
 355         self.assertAlmostEqual(fNode.integral(False)[0],integExpected,14)
 356         self.assertAlmostEqual(f.integral(False)[0],integExpected,14)
 357         pass
 358
 359     def testGauss2Gauss2DValidated(self):
 360         srcFt=MEDCouplingDataForTest.buildFieldOnGauss_1()
 361         trgFt=MEDCouplingDataForTest.buildFieldOnGauss_2()
 362         src=MEDCouplingFieldDouble(srcFt)
 363         self.assertEqual(srcFt.getMesh().getHiddenCppPointer(),src.getMesh().getHiddenCppPointer())
 364         self.assertEqual(srcFt.getDiscretization().getHiddenCppPointer(),src.getDiscretization().getHiddenCppPointer())
 365         #values given by ASTER usecase
 366         src.setArray(DataArrayDouble([1.,1.,0.,0.,1.,1.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.]))
 367         src.getArray().setInfoOnComponents(["DOMA"])
 368         rem=MEDCouplingRemapper()
 369         rem.setIntersectionType(PointLocator)
 370         rem.prepareEx(srcFt,trgFt)
 371         trg=rem.transferField(src,1e300)
 372         self.assertEqual(trg.getMesh().getHiddenCppPointer(),trgFt.getMesh().getHiddenCppPointer())
 373         self.assertEqual(trg.getDiscretization().getHiddenCppPointer(),trgFt.getDiscretization().getHiddenCppPointer())
 374         #values given after interpolation in ASTER
 375         arrExpected=DataArrayDouble([1.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.]) ; arrExpected.setInfoOnComponents(["DOMA"])
 376         self.assertTrue(trg.getArray().isEqual(arrExpected,1e-12))
 377         #
 378         # second part of the test : reverse source and target
 379         #
 380         rem.prepareEx(trgFt,srcFt)# sorry trgFt is in the place of source and srcFt in the place of target it is not a bug
 381         trg=MEDCouplingFieldDouble(trgFt)
 382         #values given after interpolation in ASTER
 383         trg.setArray(DataArrayDouble([1.,1.,0.,0.,1.,0.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,1.,1.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,0.,1.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,0.,0.,0.,1.,0.,0.,1.,1.,0.,0.,1.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.]))
 384         trg.getArray().setInfoOnComponents(["DOMA"])
 385         src=rem.transferField(trg,1e300)
 386         #values given after interpolation in ASTER
 387         arrExpected2=DataArrayDouble([1.,1.,0.,0.,1.,1.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,1.,1.,1., 1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,1.,0.,1.,1.,1.,1.,0.,0.,0.,0.,1.,0.,0.,0.,1.,1.,1.,0.,1.,1.,1.,1.,1.,1.,0.,0.,1.,1.,0.,1.,1.,1.,0.,1.,0.,0.,0.,1.,0.,0.,0.,1.,1.,1.,0.,1.,1.,1.,1.,1.]) ; arrExpected2.setInfoOnComponents(["DOMA"])
 388         # modification of values in ASTER due to modification of algorithm
 389         # target PG 82 in target cell 32(C)/36 PG 1(C)/9 is in source cell 58(C)/120 source Gauss point 113 (1(C)/4). Values must be 1. and not 0.
 390         arrExpected2.setIJ(82,0,1.)
 391         self.assertTrue(src.getArray().isEqual(arrExpected2,1e-12))
 392         pass
 393
 394     def testGauss2Gauss3DValidated(self):
 395         srcFt=MEDCouplingDataForTest.buildFieldOnGauss_3()
 396         trgFt=MEDCouplingDataForTest.buildFieldOnGauss_4()
 397         src=MEDCouplingFieldDouble(srcFt)
 398         self.assertEqual(srcFt.getMesh().getHiddenCppPointer(),src.getMesh().getHiddenCppPointer())
 399         self.assertEqual(srcFt.getDiscretization().getHiddenCppPointer(),src.getDiscretization().getHiddenCppPointer())
 400         #values given by ASTER usecase
 401         src.setArray(DataArrayDouble([0.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,1.,1.,0.,0.,1.,1.,1.,1.,0.,0.,1.,1.,0.,0.]))
 402         src.getArray().setInfoOnComponents(["DOMA"])
 403         rem=MEDCouplingRemapper()
 404         rem.setIntersectionType(PointLocator)
 405         rem.prepareEx(srcFt,trgFt)
 406         trg=rem.transferField(src,1e300)
 407         self.assertEqual(trg.getMesh().getHiddenCppPointer(),trgFt.getMesh().getHiddenCppPointer())
 408         self.assertEqual(trg.getDiscretization().getHiddenCppPointer(),trgFt.getDiscretization().getHiddenCppPointer())
 409         #values given after interpolation in ASTER
 410         arrExpected=DataArrayDouble([0.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,0.,0.,1.,1.,0.,0.,1.,1.,1.,1.,0.,1.,1.,1.,0.,1.]) ; arrExpected.setInfoOnComponents(["DOMA"])
 411         self.assertTrue(trg.getArray().isEqual(arrExpected,1e-12))
 412         #
 413         # second part of the test : reverse source and target
 414         #
 415         rem.prepareEx(trgFt,srcFt)# sorry trgFt is in the place of source and srcFt in the place of target it is not a bug
 416         trg=MEDCouplingFieldDouble(trgFt)
 417         #values given after interpolation in ASTER
 418         trg.setArray(DataArrayDouble([0.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,0.,1.,1.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.]))
 419         trg.getArray().setInfoOnComponents(["DOMA"])
 420         src=rem.transferField(trg,1e300)
 421         #values given after interpolation in ASTER
 422         arrExpected2=DataArrayDouble([0.,1.,1.,1.,1.,1.,0.,0.,0.,0.,0.,0.,0.,0.,0.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,0.,1.,0.,1.,0.,1.,1.,1.,0.,1.,1.,1.,1.,0.,1.,1.,1.,0.,1.]) ; arrExpected2.setInfoOnComponents(["DOMA"])
 423         self.assertTrue(src.getArray().isEqual(arrExpected2,1e-12))
 424         pass
 425
 426     def testSwig2MixOfUMesh(self):
 427         arr0=DataArrayDouble([0,1,1.5]) ; arr1=DataArrayDouble([0,1])
 428         sc=MEDCouplingCMesh() ; sc.setCoords(arr0,arr1,arr1)
 429         tc=sc.deepCpy() ; tc.translate([0.4,0.3,0.3])
 430         # umesh-umesh
 431         # 90 (umesh-1sgtumesh)
 432         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 433         s=sc.buildUnstructured() ; t=tc.build1SGTUnstructured()
 434         self.assertTrue(isinstance(s,MEDCouplingUMesh))
 435         self.assertTrue(isinstance(t,MEDCoupling1SGTUMesh))
 436         rem.prepare(s,t,"P0P0")
 437         mat=rem.getCrudeMatrix()
 438         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 439         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 440         del s,t
 441         # 91 (umesh-1dgtumesh)
 442         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 443         s=sc.buildUnstructured() ; t=tc.buildUnstructured() ; t.convertAllToPoly() ; t=MEDCoupling1DGTUMesh(t)
 444         self.assertTrue(isinstance(s,MEDCouplingUMesh))
 445         self.assertTrue(isinstance(t,MEDCoupling1DGTUMesh))
 446         rem.prepare(s,t,"P0P0")
 447         mat=rem.getCrudeMatrix()
 448         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 449         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 450         del s,t
 451         # 165 (1sgtumesh-umesh)
 452         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 453         s=sc.build1SGTUnstructured() ; t=tc.buildUnstructured()
 454         self.assertTrue(isinstance(s,MEDCoupling1SGTUMesh))
 455         self.assertTrue(isinstance(t,MEDCouplingUMesh))
 456         rem.prepare(s,t,"P0P0")
 457         mat=rem.getCrudeMatrix()
 458         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 459         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 460         del s,t
 461         # 181 (1dgtumesh-umesh
 462         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 463         s=sc.buildUnstructured() ; s.convertAllToPoly() ; s=MEDCoupling1DGTUMesh(s) ; t=tc.buildUnstructured()
 464         self.assertTrue(isinstance(s,MEDCoupling1DGTUMesh))
 465         self.assertTrue(isinstance(t,MEDCouplingUMesh))
 466         rem.prepare(s,t,"P0P0")
 467         mat=rem.getCrudeMatrix()
 468         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 469         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 470         del s,t
 471         # 170 (1sgtumesh-1sgtumesh)
 472         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 473         s=sc.build1SGTUnstructured() ; t=tc.build1SGTUnstructured()
 474         self.assertTrue(isinstance(s,MEDCoupling1SGTUMesh))
 475         self.assertTrue(isinstance(t,MEDCoupling1SGTUMesh))
 476         rem.prepare(s,t,"P0P0")
 477         mat=rem.getCrudeMatrix()
 478         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 479         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 480         del s,t
 481         # 171 (1sgtumesh-1dgtumesh)
 482         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 483         s=sc.build1SGTUnstructured() ; t=tc.buildUnstructured() ; t.convertAllToPoly() ; t=MEDCoupling1DGTUMesh(t)
 484         self.assertTrue(isinstance(s,MEDCoupling1SGTUMesh))
 485         self.assertTrue(isinstance(t,MEDCoupling1DGTUMesh))
 486         rem.prepare(s,t,"P0P0")
 487         mat=rem.getCrudeMatrix()
 488         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 489         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 490         del s,t
 491         # 186 (1dgtumesh-1sgtumesh)
 492         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 493         s=sc.buildUnstructured() ; s.convertAllToPoly() ; s=MEDCoupling1DGTUMesh(s) ; t=tc.build1SGTUnstructured()
 494         self.assertTrue(isinstance(s,MEDCoupling1DGTUMesh))
 495         self.assertTrue(isinstance(t,MEDCoupling1SGTUMesh))
 496         rem.prepare(s,t,"P0P0")
 497         mat=rem.getCrudeMatrix()
 498         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 499         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 500         del s,t
 501         # 187 (1dgtumesh-1dgtumesh)
 502         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 503         s=sc.buildUnstructured() ; s.convertAllToPoly() ; s=MEDCoupling1DGTUMesh(s) ; t=tc.buildUnstructured() ; t.convertAllToPoly() ; t=MEDCoupling1DGTUMesh(t)
 504         self.assertTrue(isinstance(s,MEDCoupling1DGTUMesh))
 505         self.assertTrue(isinstance(t,MEDCoupling1DGTUMesh))
 506         rem.prepare(s,t,"P0P0")
 507         mat=rem.getCrudeMatrix()
 508         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 509         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 510         del s,t
 511         # (umesh-cmesh)
 512         # 167 (1sgtumesh-cmesh)
 513         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 514         s=sc.build1SGTUnstructured() ; t=tc.deepCpy()
 515         self.assertTrue(isinstance(s,MEDCoupling1SGTUMesh))
 516         self.assertTrue(isinstance(t,MEDCouplingCMesh))
 517         rem.prepare(s,t,"P0P0")
 518         mat=rem.getCrudeMatrix()
 519         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 520         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 521         del s,t
 522         # 183 (1dgtumesh-cmesh)
 523         #rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 524         #s=sc.buildUnstructured() ; s.convertAllToPoly() ; s=MEDCoupling1DGTUMesh(s) ; t=tc.deepCpy()
 525         #self.assertTrue(isinstance(s,MEDCoupling1DGTUMesh))
 526         #self.assertTrue(isinstance(t,MEDCouplingCMesh))
 527         #rem.prepare(s,t,"P0P0")
 528         #mat=rem.getCrudeMatrix()
 529         #self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 530         #self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 531         #del s,t
 532         # (cmesh-umesh)
 533         # 122 (cmesh-1sgtumesh)
 534         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 535         s=sc.deepCpy() ; t=tc.build1SGTUnstructured()
 536         self.assertTrue(isinstance(s,MEDCouplingCMesh))
 537         self.assertTrue(isinstance(t,MEDCoupling1SGTUMesh))
 538         rem.prepare(s,t,"P0P0")
 539         mat=rem.getCrudeMatrix()
 540         self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 541         self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 542         del s,t
 543         # 123 (cmesh-1dgtumesh)
 544         #rem=MEDCouplingRemapper() ; rem.setIntersectionType(Triangulation)
 545         #s=sc.deepCpy() ; t=tc.buildUnstructured() ; t.convertAllToPoly() ; t=MEDCoupling1DGTUMesh(t)
 546         #self.assertTrue(isinstance(s,MEDCouplingCMesh))
 547         #self.assertTrue(isinstance(t,MEDCoupling1DGTUMesh))
 548         #rem.prepare(s,t,"P0P0")
 549         #mat=rem.getCrudeMatrix()
 550         #self.assertEqual(2,len(mat)) ; self.assertEqual(2,len(mat[0])) ; self.assertEqual(1,len(mat[1]))
 551         #self.assertAlmostEqual(0.294,mat[0][0],14) ; self.assertAlmostEqual(0.196,mat[0][1],14) ; self.assertAlmostEqual(0.049,mat[1][1],14)
 552         #del s,t
 553         pass
 554
 555     def testSwig2BarycentricP1P13D_1(self):
 556         sCoo=DataArrayDouble([0.313,0.00218,6.90489,0.313,0.10692667,6.90489,0.313,0.10692667,6.96790167,0.313,0.00218,6.9773125,0.313,0.21167333,6.90489,0.313,0.21167333,6.95849083,0.313,0.31642,6.90489,0.313,0.31642,6.94908,0.313,0.09383333,7.04891667,0.313,0.00218,7.049735,0.313,0.18548667,7.04809833,0.313,0.27714,7.04728,0.313,0.05782667,7.133205,0.313,0.00218,7.1221575,0.313,0.11347333,7.1442525,0.313,0.16912,7.1553,0.313,0.02509333,7.19458,0.313,0.00218,7.19458,0.313,0.04800667,7.19458,0.313,0.07092,7.19458,0.31005609,0.00218,6.90460005,0.31005609,0.10692667,6.90460005,0.29776312,0.10692667,6.96640097,0.29592716,0.00218,6.97563097,0.31005609,0.21167333,6.90460005,0.29959908,0.21167333,6.95717096,0.31005609,0.31642,6.90460005,0.30143505,0.31642,6.94794095,0.28195788,0.09383333,7.04585928,0.28179823,0.00218,7.04666189,0.28211753,0.18548667,7.04505668,0.28227718,0.27714,7.04425407,0.26551404,0.05782667,7.12852804,0.2676693,0.00218,7.11769282,0.26335878,0.11347333,7.13936327,0.26120352,0.16912,7.15019849,0.25354037,0.02509333,7.18872374,0.25354037,0.00218,7.18872374,0.25354037,0.04800667,7.18872374,0.25354037,0.07092,7.18872374,0.30722531,0.00218,6.90374134,0.30722531,0.10692667,6.90374134,0.28311179,0.10692667,6.96195653,0.27951042,0.00218,6.97065101,0.30722531,0.21167333,6.90374134,0.28671316,0.21167333,6.95326205,0.30722531,0.31642,6.90374134,0.29031453,0.31642,6.94456758,0.25210869,0.09383333,7.03680463,0.25179553,0.00218,7.03756067,0.25242185,0.18548667,7.03604859,0.25273501,0.27714,7.03529255,0.21985294,0.05782667,7.1146769,0.22408063,0.00218,7.10447034,0.21562524,0.11347333,7.12488346,0.21139755,0.16912,7.13509002,0.19636574,0.02509333,7.17138,0.19636574,0.00218,7.17138,0.19636574,0.04800667,7.17138,0.19636574,0.07092,7.17138,0.30461645,0.00218,6.90234688,0.30461645,0.10692667,6.90234688,0.26960904,0.10692667,6.95473916,0.26438066,0.00218,6.96256398,0.30461645,0.21167333,6.90234688,0.27483742,0.21167333,6.94691434,0.30461645,0.31642,6.90234688,0.2800658,0.31642,6.93908952,0.22459952,0.09383333,7.02210067,0.22414487,0.00218,7.02278109,0.22505416,0.18548667,7.02142025,0.2255088,0.27714,7.02073983,0.17777143,0.05782667,7.09218386,0.18390909,0.00218,7.0829982,0.17163377,0.11347333,7.10136952,0.1654961,0.16912,7.11055518,0.1436733,0.02509333,7.14321531,0.1436733,0.00218,7.14321531,0.1436733,0.04800667,7.14321531,0.1436733,0.07092,7.14321531,0.30232976,0.00218,6.90047024,0.30232976,0.10692667,6.90047024,0.25777378,0.10692667,6.94502622,0.25111932,0.00218,6.95168068,0.30232976,0.21167333,6.90047024,0.26442825,0.21167333,6.93837175,0.30232976,0.31642,6.90047024,0.27108271,0.31642,6.93171729,0.20048753,0.09383333,7.00231247,0.19990888,0.00218,7.00289112,0.20106618,0.18548667,7.00173382,0.20164482,0.27714,7.00115518,0.14088667,0.05782667,7.06191333,0.14869844,0.00218,7.05410156,0.13307491,0.11347333,7.06972509,0.12526315,0.16912,7.07753685,0.097488,0.02509333,7.105312,0.097488,0.00218,7.105312,0.097488,0.04800667,7.105312,0.097488,0.07092,7.105312,0.30045312,0.00218,6.89818355,0.30045312,0.10692667,6.89818355,0.24806084,0.10692667,6.93319096,0.24023602,0.00218,6.93841934,0.30045312,0.21167333,6.89818355,0.25588566,0.21167333,6.92796258,0.30045312,0.31642,6.89818355,0.26371048,0.31642,6.9227342,0.18069933,0.09383333,6.97820048,0.18001891,0.00218,6.97865513,0.18137975,0.18548667,6.97774584,0.18206017,0.27714,6.9772912,0.11061614,0.05782667,7.02502857,0.1198018,0.00218,7.01889091,0.10143048,0.11347333,7.03116623,0.09224482,0.16912,7.0373039,0.05958469,0.02509333,7.0591267,0.05958469,0.00218,7.0591267,0.05958469,0.04800667,7.0591267,0.05958469,0.07092,7.0591267,0.29905866,0.00218,6.89557469,0.29905866,0.10692667,6.89557469,0.24084347,0.10692667,6.91968821,0.23214899,0.00218,6.92328958,0.29905866,0.21167333,6.89557469,0.24953795,0.21167333,6.91608684,0.29905866,0.31642,6.89557469,0.25823242,0.31642,6.91248547,0.16599537,0.09383333,6.95069131,0.16523933,0.00218,6.95100447,0.16675141,0.18548667,6.95037815,0.16750745,0.27714,6.95006499,0.0881231,0.05782667,6.98294706,0.09832966,0.00218,6.97871937,0.07791654,0.11347333,6.98717476,0.06770998,0.16912,6.99140245,0.03142,0.02509333,7.00643426,0.03142,0.00218,7.00643426,0.03142,0.04800667,7.00643426,0.03142,0.07092,7.00643426,0.29819995,0.00218,6.89274391,0.29819995,0.10692667,6.89274391,0.23639903,0.10692667,6.90503688,0.22716903,0.00218,6.90687284,0.29819995,0.21167333,6.89274391,0.24562904,0.21167333,6.90320092,0.29819995,0.31642,6.89274391,0.25485905,0.31642,6.90136495,0.15694072,0.09383333,6.92084212,0.15613811,0.00218,6.92100177,0.15774332,0.18548667,6.92068247,0.15854593,0.27714,6.92052282,0.07427196,0.05782667,6.93728596,0.08510718,0.00218,6.9351307,0.06343673,0.11347333,6.93944122,0.05260151,0.16912,6.94159648,0.01407626,0.02509333,6.94925963,0.01407626,0.00218,6.94925963,0.01407626,0.04800667,6.94925963,0.01407626,0.07092,6.94925963,0.29792818,0.00218,6.89054043,0.29792818,0.10692667,6.89054043,0.23499241,0.10692667,6.89363227,0.22559291,0.00218,6.89409403,0.29792818,0.21167333,6.89054043,0.24439191,0.21167333,6.8931705,0.29792818,0.31642,6.89054043,0.25379141,0.31642,6.89270873,0.154075,0.09383333,6.89760748,0.15325765,0.00218,6.89764764,0.15489234,0.18548667,6.89756733,0.15570969,0.27714,6.89752718,0.06988819,0.05782667,6.90174332,0.08092238,0.00218,6.90120124,0.058854,0.11347333,6.90228539,0.04781981,0.16912,6.90282747,0.00858712,0.02509333,6.90475485,0.00858712,0.00218,6.90475485,0.00858712,0.04800667,6.90475485,0.00858712,0.07092,6.90475485,0.29791,0.00218,6.820902,0.29791,0.10692667,6.820902,0.23489833,0.10692667,6.820902,0.2254875,0.00218,6.820902,0.29791,0.21167333,6.820902,0.24430917,0.21167333,6.820902,0.29791,0.31642,6.820902,0.25372,0.31642,6.820902,0.15388333,0.09383333,6.820902,0.153065,0.00218,6.820902,0.15470167,0.18548667,6.820902,0.15552,0.27714,6.820902,0.069595,0.05782667,6.820902,0.0806425,0.00218,6.820902,0.0585475,0.11347333,6.820902,0.0475,0.16912,6.820902,0.00822,0.02509333,6.820902,0.00822,0.00218,6.820902,0.00822,0.04800667,6.820902,0.00822,0.07092,6.820902],200,3)
 557         sConn=DataArrayInt([0,1,2,3,20,21,22,23,1,4,5,2,21,24,25,22,4,6,7,5,24,26,27,25,3,2,8,9,23,22,28,29,2,5,10,8,22,25,30,28,5,7,11,10,25,27,31,30,9,8,12,13,29,28,32,33,8,10,14,12,28,30,34,32,10,11,15,14,30,31,35,34,13,12,16,17,33,32,36,37,12,14,18,16,32,34,38,36,14,15,19,18,34,35,39,38,20,21,22,23,40,41,42,43,21,24,25,22,41,44,45,42,24,26,27,25,44,46,47,45,23,22,28,29,43,42,48,49,22,25,30,28,42,45,50,48,25,27,31,30,45,47,51,50,29,28,32,33,49,48,52,53,28,30,34,32,48,50,54,52,30,31,35,34,50,51,55,54,33,32,36,37,53,52,56,57,32,34,38,36,52,54,58,56,34,35,39,38,54,55,59,58,40,41,42,43,60,61,62,63,41,44,45,42,61,64,65,62,44,46,47,45,64,66,67,65,43,42,48,49,63,62,68,69,42,45,50,48,62,65,70,68,45,47,51,50,65,67,71,70,49,48,52,53,69,68,72,73,48,50,54,52,68,70,74,72,50,51,55,54,70,71,75,74,53,52,56,57,73,72,76,77,52,54,58,56,72,74,78,76,54,55,59,58,74,75,79,78,60,61,62,63,80,81,82,83,61,64,65,62,81,84,85,82,64,66,67,65,84,86,87,85,63,62,68,69,83,82,88,89,62,65,70,68,82,85,90,88,65,67,71,70,85,87,91,90,69,68,72,73,89,88,92,93,68,70,74,72,88,90,94,92,70,71,75,74,90,91,95,94,73,72,76,77,93,92,96,97,72,74,78,76,92,94,98,96,74,75,79,78,94,95,99,98,80,81,82,83,100,101,102,103,81,84,85,82,101,104,105,102,84,86,87,85,104,106,107,105,83,82,88,89,103,102,108,109,82,85,90,88,102,105,110,108,85,87,91,90,105,107,111,110,89,88,92,93,109,108,112,113,88,90,94,92,108,110,114,112,90,91,95,94,110,111,115,114,93,92,96,97,113,112,116,117,92,94,98,96,112,114,118,116,94,95,99,98,114,115,119,118,100,101,102,103,120,121,122,123,101,104,105,102,121,124,125,122,104,106,107,105,124,126,127,125,103,102,108,109,123,122,128,129,102,105,110,108,122,125,130,128,105,107,111,110,125,127,131,130,109,108,112,113,129,128,132,133,108,110,114,112,128,130,134,132,110,111,115,114,130,131,135,134,113,112,116,117,133,132,136,137,112,114,118,116,132,134,138,136,114,115,119,118,134,135,139,138,120,121,122,123,140,141,142,143,121,124,125,122,141,144,145,142,124,126,127,125,144,146,147,145,123,122,128,129,143,142,148,149,122,125,130,128,142,145,150,148,125,127,131,130,145,147,151,150,129,128,132,133,149,148,152,153,128,130,134,132,148,150,154,152,130,131,135,134,150,151,155,154,133,132,136,137,153,152,156,157,132,134,138,136,152,154,158,156,134,135,139,138,154,155,159,158,140,141,142,143,160,161,162,163,141,144,145,142,161,164,165,162,144,146,147,145,164,166,167,165,143,142,148,149,163,162,168,169,142,145,150,148,162,165,170,168,145,147,151,150,165,167,171,170,149,148,152,153,169,168,172,173,148,150,154,152,168,170,174,172,150,151,155,154,170,171,175,174,153,152,156,157,173,172,176,177,152,154,158,156,172,174,178,176,154,155,159,158,174,175,179,178,160,161,162,163,180,181,182,183,161,164,165,162,181,184,185,182,164,166,167,165,184,186,187,185,163,162,168,169,183,182,188,189,162,165,170,168,182,185,190,188,165,167,171,170,185,187,191,190,169,168,172,173,189,188,192,193,168,170,174,172,188,190,194,192,170,171,175,174,190,191,195,194,173,172,176,177,193,192,196,197,172,174,178,176,192,194,198,196,174,175,179,178,194,195,199,198])
 558         s=MEDCoupling1SGTUMesh("target",NORM_HEXA8) ; s.setCoords(sCoo)
 559         s.setNodalConnectivity(sConn)
 560         #
 561         tCoo=DataArrayDouble([0.328,0.012,6.8598,0.328,0.168320184237353,6.8598,0.328,0.324640368474706,6.8598,0.328,0.0,6.8598,0.298,0.012,6.8598,0.1565,0.012,6.8598,0.180205346493166,0.144794653506834,6.8598,0.298,0.168320184237353,6.8598,0.0,0.012,6.8598,0.0916755774886107,0.233324422511389,6.8598,0.298,0.324640368474706,6.8598,0.298,0.0,6.8598,0.1565,0.0,6.8598,0.0,0.0,6.8598,0.328,0.012,7.2298,0.328,0.168320184237353,7.2298,0.328,0.324640368474706,7.2298,0.328,0.0,7.2298,0.298,0.012,7.2298,0.1565,0.012,7.2298,0.180205346493166,0.144794653506834,7.2298,0.298,0.168320184237353,7.2298,0.0,0.012,7.2298,0.0916755774886107,0.233324422511389,7.2298,0.298,0.324640368474706,7.2298,0.298,0.0,7.2298,0.1565,0.0,7.2298,0.0,0.0,7.2298],28,3)
 562         tConn=DataArrayInt([4,5,6,7,18,19,20,21,5,8,9,6,19,22,23,20,6,9,10,7,20,23,24,21,11,12,5,4,25,26,19,18,12,13,8,5,26,27,22,19,3,11,4,0,17,25,18,14,0,4,7,1,14,18,21,15,1,7,10,2,15,21,24,16])
 563         t=MEDCoupling1SGTUMesh("target",NORM_HEXA8) ; t.setCoords(tCoo)
 564         t.setNodalConnectivity(tConn)
 565         #
 566         s.simplexize(PLANAR_FACE_5)
 567         aRemapper=MEDCouplingRemapper()
 568         aRemapper.setPrecision(1e-12)
 569         aRemapper.setIntersectionType(Barycentric)
 570         self.assertEqual(aRemapper.prepare(s,t,'P1P1'),1)
 571         m=aRemapper.getCrudeMatrix()
 572         self.assertEqual(len(m),28)
 573         for i in xrange(28):
 574             if i not in [5,6]:
 575                 self.assertEqual(len(m[i]),0)
 576                 pass
 577             pass
 578         self.assertEqual(len(m[5]),4)
 579         self.assertEqual(len(m[6]),4)
 580         self.assertAlmostEqual(0.10714286103952797,m[5][168],12)
 581         self.assertAlmostEqual(0.35691534416938014,m[5][169],12)
 582         self.assertAlmostEqual(0.04492099619713096,m[5][163],12)
 583         self.assertAlmostEqual(0.49102079859396097,m[5][189],12)
 584         self.assertAlmostEqual(0.14039089397104254,m[6][185],12)
 585         self.assertAlmostEqual(0.16362822318261033,m[6][162],12)
 586         self.assertAlmostEqual(0.3438363717836785 ,m[6][188],12)
 587         self.assertAlmostEqual(0.3521445110626687 ,m[6][170],12)
 588         pass
 589
 590     @unittest.skipUnless(MEDCouplingHasNumPyBindings() and MEDCouplingHasSciPyBindings(),"requires numpy AND scipy")
 591     def testGetCrudeCSRMatrix1(self):
 592         """ testing CSR matrix output using numpy/scipy.
 593         """
 594         from scipy.sparse import spdiags #diags
 595         import scipy
 596         from numpy import array
 597         arr=DataArrayDouble(3) ; arr.iota()
 598         m=MEDCouplingCMesh() ; m.setCoords(arr,arr)
 599         src=m.buildUnstructured()
 600         trg=src.deepCpy() ; trg=trg[[0,1,3]]
 601         trg.getCoords()[:]*=0.5 ; trg.getCoords()[:]+=[0.3,0.25]
 602         # Let's interpolate.
 603         rem=MEDCouplingRemapper()
 604         rem.prepare(src,trg,"P0P0")
 605         # Internal crude sparse matrix computed. Let's manipulate it using CSR matrix in scipy.
 606         for i in xrange(10):
 607             m=rem.getCrudeCSRMatrix()
 608             pass
 609         m2=rem.getCrudeCSRMatrix()
 610         diff=m-m2
 611         self.assertTrue(isinstance(m,scipy.sparse.csr.csr_matrix))
 612         self.assertEqual(m.getnnz(),7)
 613         self.assertAlmostEqual(m[0,0],0.25,12)
 614         self.assertAlmostEqual(m[1,0],0.1,12)
 615         self.assertAlmostEqual(m[1,1],0.15,12)
 616         self.assertAlmostEqual(m[2,0],0.05,12)
 617         self.assertAlmostEqual(m[2,1],0.075,12)
 618         self.assertAlmostEqual(m[2,2],0.05,12)
 619         self.assertAlmostEqual(m[2,3],0.075,12)
 620         self.assertEqual(diff.getnnz(),0)
 621         # IntegralGlobConstraint (division by sum of cols)
 622         colSum=m.sum(axis=0)
 623         # version 0.12.0 # m_0=m*diags(array(1/colSum),[0])
 624         m_0=m*spdiags(array(1/colSum),[0],colSum.shape[1],colSum.shape[1])
 625         del colSum
 626         self.assertAlmostEqual(m_0[0,0],0.625,12)
 627         self.assertAlmostEqual(m_0[1,0],0.25,12)
 628         self.assertAlmostEqual(m_0[1,1],0.6666666666666667,12)
 629         self.assertAlmostEqual(m_0[2,0],0.125,12)
 630         self.assertAlmostEqual(m_0[2,1],0.3333333333333333,12)
 631         self.assertAlmostEqual(m_0[2,2],1.,12)
 632         self.assertAlmostEqual(m_0[2,3],1.,12)
 633         self.assertEqual(m_0.getnnz(),7)
 634         # ConservativeVolumic (division by sum of rows)
 635         rowSum=m.sum(axis=1)
 636         # version 0.12.0 # m_1=diags(array(1/rowSum.transpose()),[0])*m
 637         m_1=spdiags(array(1/rowSum.transpose()),[0],rowSum.shape[0],rowSum.shape[0])*m
 638         del rowSum
 639         self.assertAlmostEqual(m_1[0,0],1.,12)
 640         self.assertAlmostEqual(m_1[1,0],0.4,12)
 641         self.assertAlmostEqual(m_1[1,1],0.6,12)
 642         self.assertAlmostEqual(m_1[2,0],0.2,12)
 643         self.assertAlmostEqual(m_1[2,1],0.3,12)
 644         self.assertAlmostEqual(m_1[2,2],0.2,12)
 645         self.assertAlmostEqual(m_1[2,3],0.3,12)
 646         self.assertEqual(m_1.getnnz(),7)
 647         pass
 648
 649     @unittest.skipUnless(MEDCouplingHasNumPyBindings() and MEDCouplingHasSciPyBindings(),"requires numpy AND scipy")
 650     def testP0P1Bary_1(self):
 651         a=MEDCouplingUMesh("a",2)
 652         a.allocateCells()
 653         conna=[0,1,3,2,1,4,5,3,4,6,7,5,6,8,9,7,8,10,11,9,10,12,13,11,12,14,15,13,14,16,17,15,16,18,19,17,18,20,21,19,20,22,23,21,22,24,25,23,24,26,27,25]
 654         a.setCoords(DataArrayDouble([1.54,0,-0.01,1.54,0.02,-0.01,1.54,0,0.01,1.54,0.02,0.01,1.54,0.04,-0.01,1.54,0.04,0.01,1.54,0.06,-0.01,1.54,0.06,0.01,1.54,0.08,-0.01,1.54,0.08,0.01,1.54,0.1,-0.01,1.54,0.1,0.01,1.54,0.12,-0.01,1.54,0.12,0.01,1.54,0.14,-0.01,1.54,0.14,0.01,1.54,0.16,-0.01,1.54,0.16,0.01,1.54,0.18,-0.01,1.54,0.18,0.01,1.54,0.2,-0.01,1.54,0.2,0.01,1.54,0.22,-0.01,1.54,0.22,0.01,1.54,0.24,-0.01,1.54,0.24,0.01,1.54,0.26,-0.01,1.54,0.26,0.01],28,3))
 655         for i in xrange(13):
 656             a.insertNextCell(NORM_QUAD4,conna[4*i:4*(i+1)])
 657             pass
 658         a.finishInsertingCells() ; a.simplexize(0)
 659         #
 660         connb=[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,0,2,39,3,5,40,6,8,41,9,11,42,12,14,43,15,17,44,18,20,45,21,23,46,24,26,47,27,29,48,30,32,49,33,35,50,36,38,51,52,2,39,53,5,40,54,8,41,55,11,42,56,14,43,57,17,44,58,20,45,59,23,46,60,26,47,61,29,48,62,32,49,63,35,50,64,38,51,52,2,65,53,5,66,54,8,67,55,11,68,56,14,69,57,17,70,58,20,71,59,23,72,60,26,73,61,29,74,62,32,75,63,35,76,64,38,77,53,2,65,54,5,66,55,8,67,56,11,68,57,14,69,58,17,70,59,20,71,60,23,72,61,26,73,62,29,74,63,32,75,64,35,76,78,38,77,53,2,40,54,5,41,55,8,42,56,11,43,57,14,44,58,17,45,59,20,46,60,23,47,61,26,48,62,29,49,63,32,50,64,35,51,78,38,79,3,2,40,6,5,41,9,8,42,12,11,43,15,14,44,18,17,45,21,20,46,24,23,47,27,26,48,30,29,49,33,32,50,36,35,51,80,38,79,3,2,1,6,5,4,9,8,7,12,11,10,15,14,13,18,17,16,21,20,19,24,23,22,27,26,25,30,29,28,33,32,31,36,35,34,80,38,37]
 661         b=MEDCouplingUMesh("b",2)
 662         b.allocateCells()
 663         for i in xrange(104):
 664             b.insertNextCell(NORM_TRI3,connb[3*i:3*(i+1)])
 665             pass
 666         b.setCoords(DataArrayDouble([1.54,0,-0.01,1.54,0.01,-0.01,1.54,0.01,0,1.54,0.02,-0.01,1.54,0.03,-0.01,1.54,0.03,0,1.54,0.04,-0.01,1.54,0.05,-0.01,1.54,0.05,0,1.54,0.06,-0.01,1.54,0.07,-0.01,1.54,0.07,0,1.54,0.08,-0.01,1.54,0.09,-0.01,1.54,0.09,0,1.54,0.1,-0.01,1.54,0.11,-0.01,1.54,0.11,0,1.54,0.12,-0.01,1.54,0.13,-0.01,1.54,0.13,0,1.54,0.14,-0.01,1.54,0.15,-0.01,1.54,0.15,0,1.54,0.16,-0.01,1.54,0.17,-0.01,1.54,0.17,0,1.54,0.18,-0.01,1.54,0.19,-0.01,1.54,0.19,0,1.54,0.2,-0.01,1.54,0.21,-0.01,1.54,0.21,0,1.54,0.22,-0.01,1.54,0.23,-0.01,1.54,0.23,0,1.54,0.24,-0.01,1.54,0.25,-0.01,1.54,0.25,0,1.54,0,0,1.54,0.02,0,1.54,0.04,0,1.54,0.06,0,1.54,0.08,0,1.54,0.1,0,1.54,0.12,0,1.54,0.14,0,1.54,0.16,0,1.54,0.18,0,1.54,0.2,0,1.54,0.22,0,1.54,0.24,0,1.54,0,0.01,1.54,0.02,0.01,1.54,0.04,0.01,1.54,0.06,0.01,1.54,0.08,0.01,1.54,0.1,0.01,1.54,0.12,0.01,1.54,0.14,0.01,1.54,0.16,0.01,1.54,0.18,0.01,1.54,0.2,0.01,1.54,0.22,0.01,1.54,0.24,0.01,1.54,0.01,0.01,1.54,0.03,0.01,1.54,0.05,0.01,1.54,0.07,0.01,1.54,0.09,0.01,1.54,0.11,0.01,1.54,0.13,0.01,1.54,0.15,0.01,1.54,0.17,0.01,1.54,0.19,0.01,1.54,0.21,0.01,1.54,0.23,0.01,1.54,0.25,0.01,1.54,0.26,0.01,1.54,0.26,0,1.54,0.26,-0.01],81,3))
 667         #
 668         rem=MEDCouplingRemapper() ; rem.setIntersectionType(Barycentric)
 669         rem.prepare(a,b,"P1P0")
 670         m0=rem.getCrudeCSRMatrix()
 671         self.assertEqual(m0.nnz,312)
 672         #
 673         ids=4*[None] ; vs=4*[None]
 674         ids[0]=DataArrayInt([0,3,6,9,12,15,18,21,24,27,30,33,36,39,42,45,48,51,54,57,60,63,66,69,72,75,158,161,164,167,170,173,176,179,182,185,188,191,194,197,200,203,206,209,212,215,218,221,224,227,230,233]) ; vs[0]=10./3.
 675         ids[1]=DataArrayInt([1,2,4,5,7,8,10,11,13,14,16,17,19,20,22,23,25,26,28,29,31,32,34,35,37,38,40,41,43,44,46,47,49,50,52,53,55,56,58,59,61,62,64,65,67,68,70,71,73,74,76,77,80,83,86,89,92,95,98,101,104,107,110,113,116,117,120,123,126,129,132,135,138,141,144,147,150,153,156,157,159,160,162,163,165,166,168,169,171,172,174,175,177,178,180,181,183,184,186,187,189,190,192,193,195,196,198,199,201,202,204,205,207,208,210,211,213,214,216,217,219,220,222,223,225,226,228,229,231,232,234,237,240,243,246,249,252,255,258,261,264,267,270,275,278,281,284,287,290,293,296,299,302,305,308,311]) ; vs[1]=5./6.
 676         ids[2]=DataArrayInt([78,81,84,87,90,93,96,99,102,105,108,111,114,119,122,125,128,131,134,137,140,143,146,149,152,155,236,239,242,245,248,251,254,257,260,263,266,269,272,273,276,279,282,285,288,291,294,297,300,303,306,309]) ; vs[2]=5./3.
 677         ids[3]=DataArrayInt([79,82,85,88,91,94,97,100,103,106,109,112,115,118,121,124,127,130,133,136,139,142,145,148,151,154,235,238,241,244,247,250,253,256,259,262,265,268,271,274,277,280,283,286,289,292,295,298,301,304,307,310]) ; vs[3]=2.5
 678         vals=DataArrayDouble(312,1)
 679         for idd,v in zip(ids,vs):
 680             vals[idd]=v
 681             pass
 682         vals*=1e-5
 683         eps0=DataArrayDouble(m0.data)-vals ; eps0.abs()
 684         self.assertTrue(eps0.getIdsInRange(1e-17,1e300).empty())
 685         self.assertTrue(DataArrayInt(m0.indices).isEqual(DataArrayInt([0,1,3,1,4,5,4,6,7,6,8,9,8,10,11,10,12,13,12,14,15,14,16,17,16,18,19,18,20,21,20,22,23,22,24,25,24,26,27,0,2,3,1,3,5,4,5,7,6,7,9,8,9,11,10,11,13,12,13,15,14,15,17,16,17,19,18,19,21,20,21,23,22,23,25,24,25,27,0,2,3,1,3,5,4,5,7,6,7,9,8,9,11,10,11,13,12,13,15,14,15,17,16,17,19,18,19,21,20,21,23,22,23,25,24,25,27,0,2,3,1,3,5,4,5,7,6,7,9,8,9,11,10,11,13,12,13,15,14,15,17,16,17,19,18,19,21,20,21,23,22,23,25,24,25,27,0,2,3,1,3,5,4,5,7,6,7,9,8,9,11,10,11,13,12,13,15,14,15,17,16,17,19,18,19,21,20,21,23,22,23,25,24,25,27,0,1,3,1,4,5,4,6,7,6,8,9,8,10,11,10,12,13,12,14,15,14,16,17,16,18,19,18,20,21,20,22,23,22,24,25,24,26,27,0,1,3,1,4,5,4,6,7,6,8,9,8,10,11,10,12,13,12,14,15,14,16,17,16,18,19,18,20,21,20,22,23,22,24,25,24,26,27,0,1,3,1,4,5,4,6,7,6,8,9,8,10,11,10,12,13,12,14,15,14,16,17,16,18,19,18,20,21,20,22,23,22,24,25,24,26,27])))
 686         self.assertTrue(DataArrayInt(m0.indptr).isEqual(DataArrayInt([0,3,6,9,12,15,18,21,24,27,30,33,36,39,42,45,48,51,54,57,60,63,66,69,72,75,78,81,84,87,90,93,96,99,102,105,108,111,114,117,120,123,126,129,132,135,138,141,144,147,150,153,156,159,162,165,168,171,174,177,180,183,186,189,192,195,198,201,204,207,210,213,216,219,222,225,228,231,234,237,240,243,246,249,252,255,258,261,264,267,270,273,276,279,282,285,288,291,294,297,300,303,306,309,312])))
 687         #
 688         rem2=MEDCouplingRemapper() ; rem2.setIntersectionType(Barycentric)
 689         rem2.prepare(b,a,"P0P1")
 690         m1=rem2.getCrudeCSRMatrix()
 691         self.assertEqual(m1.nnz,312)
 692         #
 693         m1=rem2.getCrudeCSRMatrix()
 694         m1t=m1.transpose()
 695         delta=m0-m1t
 696         self.assertTrue(DataArrayDouble(delta.data).isUniform(0.,1e-17))
 697         pass
 698
 699     @unittest.skipUnless(MEDCouplingHasNumPyBindings() and MEDCouplingHasSciPyBindings(),"requires numpy AND scipy")
 700     def testNonConformWithRemapper_1(self):
 701         coo=DataArrayDouble([-0.396700000780411,-0.134843245350081,-0.0361311386958691,-0.407550009429364,-0.13484324535008,-0.0361311386958923,-0.396700000780411,-0.132191446077668,-0.0448729493559049,-0.407550009429364,-0.132191446077666,-0.0448729493559254,-0.396700000780411,-0.128973582738749,-0.0534226071577727,-0.407550009429364,-0.128973582738747,-0.0534226071577904,-0.396700000780411,-0.128348829636458,-0.0346583696473619,-0.407550009429364,-0.128348829636457,-0.0346583696473822,-0.396700000780411,-0.125874740261886,-0.0430683597970123,-0.407550009429364,-0.125874740261885,-0.0430683597970302,-0.396700000780411,-0.122905344829122,-0.051310216195766,-0.407550009429364,-0.12290534482912,-0.0513102161957814],12,3)
 702         conn=DataArrayInt([2,9,3,11,2,3,5,11,2,8,9,11,2,10,8,11,2,5,4,11,2,4,10,11,3,0,1,6,3,1,7,6,3,2,0,6,3,8,2,6,3,7,9,6,3,9,8,6])
 703         m=MEDCoupling1SGTUMesh("mesh",NORM_TETRA4)
 704         m.setNodalConnectivity(conn)
 705         m.setCoords(coo)
 706         # m is ready
 707         m1,d,di,rd,rdi=m.buildUnstructured().buildDescendingConnectivity()
 708         rdi2=rdi.deltaShiftIndex()
 709         cellIds=rdi2.getIdsEqual(1)
 710         skinAndNonConformCells=m1[cellIds]
 711         skinAndNonConformCells.zipCoords() # at this point skinAndNonConformCells contains non conform cells and skin cells. Now trying to split them in two parts.
 712         #
 713         rem=MEDCouplingRemapper()
 714         rem.setMaxDistance3DSurfIntersect(1e-12)
 715         rem.setMinDotBtwPlane3DSurfIntersect(0.99)# this line is important it is to tell to remapper to select only cells with very close orientation
 716         rem.prepare(skinAndNonConformCells,skinAndNonConformCells,"P0P0")
 717         mat=rem.getCrudeCSRMatrix()
 718         indptr=DataArrayInt(mat.indptr)
 719         indptr2=indptr.deltaShiftIndex()
 720         cellIdsOfNonConformCells=indptr2.getIdsNotEqual(1)
 721         cellIdsOfSkin=indptr2.getIdsEqual(1)
 722         self.assertTrue(cellIdsOfSkin.isEqual(DataArrayInt([1,2,3,5,6,7,8,9,10,11,12,13,14,15,16,17,19,20,21,23])))
 723         self.assertTrue(cellIdsOfNonConformCells.isEqual(DataArrayInt([0,4,18,22])))
 724         pass
 725
 726     def test3D1DOnP1P0_1(self):
 727         """ This test focused on P1P0 interpolation with a source with meshDim=1 spaceDim=3 and a target with meshDim=3.
 728         This test has revealed a bug in remapper. A reverse matrix is computed so a reverse method should be given in input.
 729         """
 730         target=MEDCouplingCMesh()
 731         arrX=DataArrayDouble([0,1]) ; arrY=DataArrayDouble([0,1]) ; arrZ=DataArrayDouble(11) ; arrZ.iota()
 732         target.setCoords(arrX,arrY,arrZ)
 733         target=target.buildUnstructured() ; target.setName("TargetSecondaire")
 734         #
 735         sourceCoo=DataArrayDouble([(0.5,0.5,0.1),(0.5,0.5,1.2),(0.5,0.5,1.6),(0.5,0.5,1.8),(0.5,0.5,2.43),(0.5,0.5,2.55),(0.5,0.5,4.1),(0.5,0.5,4.4),(0.5,0.5,4.9),(0.5,0.5,5.1),(0.5,0.5,7.6),(0.5,0.5,7.7),(0.5,0.5,8.2),(0.5,0.5,8.4),(0.5,0.5,8.6),(0.5,0.5,8.8),(0.5,0.5,9.2),(0.5,0.5,9.6),(0.5,0.5,11.5)])
 736         source=MEDCoupling1SGTUMesh("SourcePrimaire",NORM_SEG2)
 737         source.setCoords(sourceCoo)
 738         source.allocateCells()
 739         for i in xrange(len(sourceCoo)-1):
 740             source.insertNextCell([i,i+1])
 741             pass
 742         source=source.buildUnstructured()
 743         fsource=MEDCouplingFieldDouble(ON_NODES) ; fsource.setName("field")
 744         fsource.setMesh(source)
 745         arr=DataArrayDouble(len(sourceCoo)) ; arr.iota(0.7) ; arr*=arr
 746         fsource.setArray(arr)
 747         fsource.setNature(ConservativeVolumic)
 748         #
 749         rem=MEDCouplingRemapper()
 750         rem.setIntersectionType(PointLocator)
 751         rem.prepare(source,target,"P1P0")
 752         f2Test=rem.transferField(fsource,-27)
 753         self.assertEqual(f2Test.getName(),fsource.getName())
 754         self.assertEqual(f2Test.getMesh().getHiddenCppPointer(),target.getHiddenCppPointer())
 755         expArr=DataArrayDouble([0.49,7.956666666666667,27.29,-27,59.95666666666667,94.09,-27,125.69,202.89,296.09])
 756         self.assertTrue(f2Test.getArray().isEqual(expArr,1e-12))
 757         f2Test=rem.reverseTransferField(f2Test,-36)
 758         self.assertEqual(f2Test.getName(),fsource.getName())
 759         self.assertEqual(f2Test.getMesh().getHiddenCppPointer(),source.getHiddenCppPointer())
 760         expArr2=DataArrayDouble([0.49,7.956666666666667,7.956666666666667,7.956666666666667,27.29,27.29,59.95666666666667,59.95666666666667,59.95666666666667,94.09,125.69,125.69,202.89,202.89,202.89,202.89,296.09,296.09,-36.])
 761         self.assertTrue(f2Test.getArray().isEqual(expArr2,1e-12))
 762         pass
 763
 764     def testRemapperAMR1(self):
 765         """ This test is the origin of the ref values for MEDCouplingBasicsTest.testAMR2"""
 766         coarse=DataArrayDouble(35) ; coarse.iota(0) #X=5,Y=7
 767         fine=DataArrayDouble(3*2*4*4) ; fine.iota(0) #X=3,Y=2 refined by 4
 768         MEDCouplingIMesh.CondenseFineToCoarse([5,7],fine,[(1,4),(2,4)],[4,4],coarse)
 769         #
 770         m=MEDCouplingCartesianAMRMesh("mesh",2,[6,8],[0.,0.],[1.,1.])
 771         trgMesh=m.buildUnstructured()
 772         m.addPatch([(1,4),(2,4)],[4,4])
 773         srcMesh=m[0].getMesh().buildUnstructured()
 774         srcField=MEDCouplingFieldDouble(ON_CELLS)
 775         fine2=DataArrayDouble(3*2*4*4) ; fine2.iota(0) ; srcField.setArray(fine2)
 776         srcField.setMesh(srcMesh) ; srcField.setNature(Integral)
 777         #
 778         trgField=MEDCouplingFieldDouble(ON_CELLS)
 779         coarse2=DataArrayDouble(35) ; coarse2.iota(0) ; trgField.setArray(coarse2)
 780         trgField.setMesh(trgMesh) ; trgField.setNature(Integral)
 781         #
 782         rem=MEDCouplingRemapper()
 783         rem.prepare(srcMesh,trgMesh,"P0P0")
 784         rem.partialTransfer(srcField,trgField)
 785         #
 786         self.assertTrue(coarse.isEqual(trgField.getArray(),1e-12))
 787         pass
 788
 789     @unittest.skipUnless(MEDCouplingHasNumPyBindings() and MEDCouplingHasSciPyBindings(),"requires numpy AND scipy")
 790     def test1DPointLocator1(self):
 791         """This test focuses on PointLocator for P1P1 in 1D and 2DCurve."""
 792         from numpy import array
 793         from scipy.sparse import diags,csr_matrix,identity
 794         ## basic case 1D
 795         arrS=DataArrayInt.Range(0,11,1).convertToDblArr()
 796         arrT=DataArrayDouble([0.1,1.7,5.5,9.6])
 797         mS=MEDCouplingCMesh() ; mS.setCoords(arrS)
 798         mT=MEDCouplingCMesh() ; mT.setCoords(arrT)
 799         rem=MEDCouplingRemapper()
 800         rem.setIntersectionType(PointLocator)
 801         self.assertEqual(rem.prepare(mS.buildUnstructured(),mT.buildUnstructured(),"P1P1"),1)
 802         m=rem.getCrudeCSRMatrix()
 803         rowSum=m.sum(axis=1)
 804         m=diags(array(1/rowSum.transpose()),[0])*m
 805         # expected matrix
 806         row=array([0,0,1,1,2,2,3,3])
 807         col=array([0,1,1,2,5,6,9,10])
 808         data=array([0.9,0.1,0.3,0.7,0.5,0.5,0.4,0.6])
 809         mExp0=csr_matrix((data,(row,col)),shape=(4,11))
 810         # compute diff and check
 811         diff=abs(m-mExp0)
 812         self.assertAlmostEqual(diff.max(),0.,14)
 813         ## full specific case 1D where target=source
 814         rem=MEDCouplingRemapper()
 815         rem.setIntersectionType(PointLocator)
 816         self.assertEqual(rem.prepare(mS.buildUnstructured(),mS.buildUnstructured(),"P1P1"),1)
 817         m=rem.getCrudeCSRMatrix()
 818         rowSum=m.sum(axis=1)
 819         m=diags(array(1/rowSum.transpose()),[0])*m
 820         # expected matrix
 821         mExp1=identity(11)
 822         diff=abs(m-mExp1)
 823         self.assertAlmostEqual(diff.max(),0.,14)
 824         ## case where some points in target are not in source
 825         arrT=DataArrayDouble([-0.2,0.1,1.7,5.5,10.3])
 826         mT=MEDCouplingCMesh() ; mT.setCoords(arrT)
 827         mT=mT.buildUnstructured()
 828         rem=MEDCouplingRemapper()
 829         rem.setIntersectionType(PointLocator)
 830         self.assertEqual(rem.prepare(mS.buildUnstructured(),mT,"P1P1"),1)
 831         m=rem.getCrudeCSRMatrix()
 832         row=array([1,1,2,2,3,3])
 833         col=array([0,1,1,2,5,6])
 834         data=array([0.9,0.1,0.3,0.7,0.5,0.5])
 835         mExp2=csr_matrix((data,(row,col)),shape=(5,11))
 836         diff=abs(m-mExp2)
 837         self.assertAlmostEqual(diff.max(),0.,14)
 838         ## basic case 2D Curve
 839         arrS=DataArrayInt.Range(0,11,1).convertToDblArr()
 840         arrT=DataArrayDouble([0.1,1.7,5.5,9.6])
 841         mS=MEDCouplingCMesh() ; mS.setCoords(arrS)
 842         mT=MEDCouplingCMesh() ; mT.setCoords(arrT)
 843         mS=mS.buildUnstructured() ; mS.changeSpaceDimension(2)
 844         mT=mT.buildUnstructured() ; mT.changeSpaceDimension(2)
 845         mS.rotate([-1.,-1.],1.2)
 846         mT.rotate([-1.,-1.],1.2)
 847         rem=MEDCouplingRemapper()
 848         rem.setIntersectionType(PointLocator)
 849         self.assertEqual(rem.prepare(mS,mT,"P1P1"),1)
 850         m=rem.getCrudeCSRMatrix()
 851         rowSum=m.sum(axis=1)
 852         m=diags(array(1/rowSum.transpose()),[0])*m
 853         diff=abs(m-mExp0)
 854         self.assertAlmostEqual(diff.max(),0.,14)
 855         pass
 856
 857     def build2DSourceMesh_1(self):
 858         sourceCoords=[-0.3,-0.3, 0.7,-0.3, -0.3,0.7, 0.7,0.7]
 859         sourceConn=[0,3,1,0,2,3]
 860         sourceMesh=MEDCouplingUMesh.New("my name of mesh 2D",2)
 861         sourceMesh.allocateCells(2);
 862         sourceMesh.insertNextCell(NORM_TRI3,3,sourceConn[0:3]);
 863         sourceMesh.insertNextCell(NORM_TRI3,3,sourceConn[3:6]);
 864         sourceMesh.finishInsertingCells();
 865         myCoords=DataArrayDouble.New();
 866         myCoords.setValues(sourceCoords,4,2);
 867         sourceMesh.setCoords(myCoords);
 868         return sourceMesh;
 869
 870     def build2DTargetMesh_1(self):
 871         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 ]
 872         targetConn=[0,3,4,1, 1,4,2, 4,5,2, 6,7,4,3, 7,8,5,4]
 873         targetMesh=MEDCouplingUMesh.New();
 874         targetMesh.setMeshDimension(2);
 875         targetMesh.allocateCells(5);
 876         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[0:4]);
 877         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[4:7]);
 878         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[7:10]);
 879         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[10:14]);
 880         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[14:18]);
 881         targetMesh.finishInsertingCells();
 882         myCoords=DataArrayDouble.New();
 883         myCoords.setValues(targetCoords,9,2);
 884         targetMesh.setCoords(myCoords);
 885         return targetMesh;
 886
 887     def build2DTargetMesh_3(self):
 888         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]
 889         targetConn=[0,4,5,1, 1,5,6,2, 2,6,7,3, 4,8,9,5]
 890         targetMesh=MEDCouplingUMesh.New();
 891         targetMesh.setMeshDimension(2);
 892         targetMesh.allocateCells(4);
 893         for i in xrange(4):
 894             targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[4*i:4*(i+1)])
 895             pass
 896         targetMesh.finishInsertingCells();
 897         myCoords=DataArrayDouble.New();
 898         myCoords.setValues(targetCoords,10,2);
 899         targetMesh.setCoords(myCoords);
 900         return targetMesh;
 901         pass
 902
 903     def setUp(self):
 904         pass
 905     pass
 906
 907 unittest.main()