src/MEDCouplingCorba_Swig/MEDCouplingCorbaSwigTest.py

   1 # Copyright (C) 2007-2013  CEA/DEN, EDF R&D
   2 #
   3 # This library is free software; you can redistribute it and/or
   4 # modify it under the terms of the GNU Lesser General Public
   5 # License as published by the Free Software Foundation; either
   6 # version 2.1 of the License.
   7 #
   8 # This library is distributed in the hope that it will be useful,
   9 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 # Lesser General Public License for more details.
  12 #
  13 # You should have received a copy of the GNU Lesser General Public
  14 # License along with this library; if not, write to the Free Software
  15 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 #
  17 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 #
  19 # Author : Anthony Geay (CEA/DEN)
  20
  21 from MEDCoupling import *
  22 import math
  23 import os
  24
  25 class MEDCouplingCorbaServBasicsTest:
  26     def build1DMesh(self):
  27         coords=[ 0.0, 0.3, 0.75, 1.0 ]
  28         conn=[ 0,1, 1,2, 2,3 ]
  29         mesh=MEDCouplingUMesh.New("1DMeshForCorba",1);
  30         mesh.setDescription("build1DMesh");
  31         mesh.allocateCells(3);
  32         mesh.setTime(5.6,7,8);
  33         mesh.setTimeUnit("ms");
  34         mesh.insertNextCell(NORM_SEG2,2,conn[0:2]);
  35         mesh.insertNextCell(NORM_SEG2,2,conn[2:4]);
  36         mesh.insertNextCell(NORM_SEG2,2,conn[4:6]);
  37         mesh.finishInsertingCells();
  38         myCoords=DataArrayDouble.New();
  39         myCoords.setValues(coords,4,1);
  40         mesh.setCoords(myCoords);
  41         mesh.changeSpaceDimension(3);
  42         myCoords=mesh.getCoords();
  43         myCoords.setInfoOnComponent(0,"X1D [m]");
  44         myCoords.setInfoOnComponent(1,"Y1D [dm]");
  45         myCoords.setInfoOnComponent(2,"Z1D [pm]");
  46         center=[0.,0.,0.]
  47         vector=[0.,1.,0.]
  48         mesh.rotate(center,vector,-math.pi/2.);
  49         return mesh
  50
  51     def build2DMesh(self):
  52         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 ];
  53         targetConn=[0,3,4,1, 1,4,2, 4,5,2, 6,7,4,3, 7,8,5,4]
  54         targetMesh=MEDCouplingUMesh.New();
  55         targetMesh.setMeshDimension(2);
  56         targetMesh.setName("MyMesh2D");
  57         targetMesh.setDescription("build2DMesh");
  58         targetMesh.allocateCells(5);
  59         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[0:4]);
  60         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[4:7]);
  61         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[7:10]);
  62         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[10:14]);
  63         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[14:18]);
  64         targetMesh.finishInsertingCells();
  65         myCoords=DataArrayDouble.New();
  66         myCoords.setValues(targetCoords,9,2);
  67         targetMesh.setCoords(myCoords);
  68         return targetMesh;
  69
  70     def build3DMesh(self):
  71         targetCoords=[ 0., 0., 0., 50., 0., 0. , 200., 0., 0.  , 0., 50., 0., 50., 50., 0. , 200., 50., 0.,   0., 200., 0., 50., 200., 0. , 200., 200., 0. ,
  72                        0., 0., 50., 50., 0., 50. , 200., 0., 50.  , 0., 50., 50., 50., 50., 50. , 200., 50., 50.,   0., 200., 50., 50., 200., 50. , 200., 200., 50. ,
  73                        0., 0., 200., 50., 0., 200. , 200., 0., 200.  , 0., 50., 200., 50., 50., 200. , 200., 50., 200.,   0., 200., 200., 50., 200., 200. , 200., 200., 200. ];
  74         targetConn=[0,1,4,3,9,10,13,12, 1,2,5,4,10,11,14,13, 3,4,7,6,12,13,16,15, 4,5,8,7,13,14,17,16,
  75                     9,10,13,12,18,19,22,21, 10,11,14,13,19,20,23,22, 12,13,16,15,21,22,25,24, 13,14,17,16,22,23,26,25];
  76         targetMesh=MEDCouplingUMesh.New();
  77         targetMesh.setMeshDimension(3);
  78         targetMesh.setName("MyMesh3D");
  79         targetMesh.setDescription("build3DMesh");
  80         targetMesh.allocateCells(12);
  81         for i in xrange(8):
  82             targetMesh.insertNextCell(NORM_HEXA8,8,targetConn[8*i:8*(i+1)]);
  83             pass
  84         targetMesh.finishInsertingCells();
  85         myCoords=DataArrayDouble.New();
  86         myCoords.setValues(targetCoords,27,3);
  87         targetMesh.setCoords(myCoords)
  88         myCoords.setName("check in case")
  89         return targetMesh;
  90
  91     def build3DSurfMesh(self):
  92         targetCoords=[-0.3,-0.3,0.5, 0.2,-0.3,1., 0.7,-0.3,1.5, -0.3,0.2,0.5, 0.2,0.2,1., 0.7,0.2,1.5, -0.3,0.7,0.5, 0.2,0.7,1., 0.7,0.7,1.5];
  93         targetConn=[0,3,4,1, 1,4,2, 4,5,2, 6,7,4,3, 7,8,5,4];
  94         targetMesh=MEDCouplingUMesh.New();
  95         targetMesh.setMeshDimension(2);
  96         targetMesh.setName("MyMesh3DSurf");
  97         targetMesh.setDescription("build3DSurfMesh");
  98         targetMesh.allocateCells(5);
  99         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[0:4]);
 100         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[4:7]);
 101         targetMesh.insertNextCell(NORM_TRI3,3,targetConn[7:10]);
 102         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[10:14]);
 103         targetMesh.insertNextCell(NORM_QUAD4,4,targetConn[14:18]);
 104         targetMesh.finishInsertingCells();
 105         myCoords=DataArrayDouble.New();
 106         myCoords.setValues(targetCoords,9,3);
 107         targetMesh.setCoords(myCoords);
 108         myCoords.setInfoOnComponent(0,"X [m]");
 109         myCoords.setInfoOnComponent(1,"X [dm]");
 110         myCoords.setInfoOnComponent(2,"X [m]");
 111         return targetMesh;
 112
 113     def build0DMesh(self):
 114         targetCoords=[-0.3,-0.3,0.5, 0.2,-0.3,1., 0.7,-0.3,1.5, -0.3,0.2,0.5, 0.2,0.2,1., 0.7,0.2,1.5, -0.3,0.7,0.5, 0.2,0.7,1., 0.7,0.7,1.5];
 115         targetMesh=MEDCouplingUMesh.New();
 116         targetMesh.setMeshDimension(0);
 117         targetMesh.allocateCells(8);
 118         targetMesh.setName("Wonderfull 0D mesh");
 119         targetMesh.setDescription("build0DMesh");
 120         targetConn=[]
 121         targetMesh.insertNextCell(NORM_POINT1,1,[0]);
 122         targetMesh.insertNextCell(NORM_POINT1,1,[1]);
 123         targetMesh.insertNextCell(NORM_POINT1,1,[2]);
 124         targetMesh.insertNextCell(NORM_POINT1,1,[3]);
 125         targetMesh.insertNextCell(NORM_POINT1,1,[4]);
 126         targetMesh.insertNextCell(NORM_POINT1,1,[5]);
 127         targetMesh.insertNextCell(NORM_POINT1,1,[7]);
 128         targetMesh.insertNextCell(NORM_POINT1,1,[6]);
 129         targetMesh.finishInsertingCells();
 130         myCoords=DataArrayDouble.New();
 131         myCoords.setValues(targetCoords,9,3);
 132         targetMesh.setCoords(myCoords);
 133         myCoords.setInfoOnComponent(0,"X [m]");
 134         myCoords.setInfoOnComponent(1,"YY [Pm]");
 135         myCoords.setInfoOnComponent(2,"ZZZ [m]");
 136         targetMesh.checkCoherency();
 137         return targetMesh;
 138
 139     def buildM1DMesh(self):
 140         meshM1D=MEDCouplingUMesh.New("wonderfull -1 D mesh",-1);
 141         meshM1D.setDescription("buildM1DMesh");
 142         meshM1D.checkCoherency();
 143         return meshM1D;
 144
 145     def buildExtrudedMesh(self):
 146         m2D=self.build2DMesh();
 147         m2D.changeSpaceDimension(3);
 148         m1D=self.build1DMesh();
 149         retu=m2D.buildExtrudedMesh(m1D,0);
 150         ret=MEDCouplingExtrudedMesh.New(retu,m2D,2);
 151         ret.setName("ExtrudedTestForCorbaTest");
 152         ret.setDescription("buildExtrudedMesh");
 153         return ret;
 154
 155     def buildCMesh(self):
 156         targetMesh=MEDCouplingCMesh.New();
 157         targetMesh.setTime(2.3,4,5);
 158         targetMesh.setTimeUnit("us");
 159         targetMesh.setName("Example of CMesh");
 160         targetMesh.setDescription("buildCMesh");
 161         a1=DataArrayDouble.New();
 162         a1Data=[3.,4.,5.,6.,7.]
 163         a1.setValues(a1Data,5,1);
 164         a1.setInfoOnComponent(0,"SmthX");
 165         a2=DataArrayDouble.New();
 166         a2Data=[2.78,3.,4.,5.,6.,7.]
 167         a2.setValues(a2Data,6,1);
 168         a2.setInfoOnComponent(0,"SmthZ");
 169         #
 170         targetMesh.setCoordsAt(0,a1);
 171         targetMesh.setCoordsAt(2,a2);
 172         #
 173         #
 174         targetMesh.checkCoherency();
 175         #
 176         return targetMesh;
 177
 178     def buildCLMesh(self):
 179         targetMesh=MEDCouplingCurveLinearMesh();
 180         targetMesh.setTime(2.3,4,5);
 181         targetMesh.setTimeUnit("us");
 182         targetMesh.setName("Example of Cuve linear mesh");
 183         targetMesh.setDescription("buildCLMesh");
 184         a1=DataArrayDouble(3*20,1);
 185         a1.iota(7.) ; a1.rearrange(3);
 186         targetMesh.setCoords(a1);
 187         targetMesh.setNodeGridStructure([4,5]);
 188         #
 189         targetMesh.checkCoherency();
 190         #
 191         return targetMesh;
 192
 193     def buildFieldScalarOn2DNT(self):
 194         mesh=self.build2DMesh();
 195         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
 196         fieldOnCells.setTimeUnit("ms");
 197         fieldOnCells.setName("toto");
 198         fieldOnCells.setMesh(mesh);
 199         array=DataArrayDouble.New();
 200         tmp=mesh.getNumberOfCells()*6*[7.]
 201         array.setValues(tmp,mesh.getNumberOfCells(),6);
 202         fieldOnCells.setArray(array)
 203         fieldOnCells.checkCoherency();
 204         return fieldOnCells;
 205
 206     def buildFieldNodeScalarOn2DNT(self):
 207         mesh=self.build2DMesh();
 208         fieldOnNodes=MEDCouplingFieldDouble.New(ON_NODES,NO_TIME);
 209         fieldOnNodes.setName("toto2");
 210         fieldOnNodes.setTimeUnit("s");
 211         fieldOnNodes.setDescription("my wonderful field toto2");
 212         fieldOnNodes.setMesh(mesh);
 213         array=DataArrayDouble.New();
 214         tmp=mesh.getNumberOfNodes()*5*[7.1234]
 215         array.setValues(tmp,mesh.getNumberOfNodes(),5);
 216         fieldOnNodes.setArray(array);
 217         fieldOnNodes.checkCoherency();
 218         return fieldOnNodes;
 219
 220     def buildFieldScalarOn3DNT(self):
 221         mesh=self.build3DMesh();
 222         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
 223         fieldOnCells.setNature(ConservativeVolumic);
 224         fieldOnCells.setName("toto");
 225         fieldOnCells.setDescription("my wonderful 3D field toto2");
 226         fieldOnCells.setMesh(mesh);
 227         array=DataArrayDouble.New();
 228         tmp=mesh.getNumberOfCells()*6*[7.]
 229         array.setValues(tmp,mesh.getNumberOfCells(),6);
 230         fieldOnCells.setArray(array);
 231         fieldOnCells.checkCoherency();
 232         return fieldOnCells;
 233
 234     def buildFieldScalarOn3DSurfWT(self):
 235         mesh=self.build3DSurfMesh();
 236         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 237         fieldOnCells.setName("toto25");
 238         fieldOnCells.setDescription("my wonderful 3D surf field toto25");
 239         fieldOnCells.setTimeUnit("us");
 240         fieldOnCells.setMesh(mesh);
 241         array=DataArrayDouble.New();
 242         tmp=mesh.getNumberOfCells()*3*[7.]
 243         array.setValues(tmp,mesh.getNumberOfCells(),3);
 244         array.setInfoOnComponents(["aaa","bbbb","ccccc"])
 245         fieldOnCells.setArray(array);
 246         fieldOnCells.setTime(6.7,1,4);
 247         fieldOnCells.checkCoherency();
 248         return fieldOnCells;
 249
 250     def buildFieldScalarOn3DSurfCOTI(self):
 251         mesh=self.build3DSurfMesh();
 252         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,CONST_ON_TIME_INTERVAL);
 253         fieldOnCells.setName("toto26");
 254         fieldOnCells.setDescription("my wonderful 3D surf field toto26");
 255         fieldOnCells.setTimeUnit("us");
 256         fieldOnCells.setMesh(mesh);
 257         array=DataArrayDouble.New();
 258         tmp=mesh.getNumberOfCells()*3*[7.]
 259         array.setValues(tmp,mesh.getNumberOfCells(),3);
 260         fieldOnCells.setArray(array);
 261         fieldOnCells.setStartTime(6.7,1,4);
 262         fieldOnCells.setEndTime(7.2,2,8);
 263         fieldOnCells.checkCoherency();
 264         return fieldOnCells;
 265
 266     def buildFieldScalarOn2DLT(self):
 267         mesh=self.build2DMesh()
 268         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,LINEAR_TIME)
 269         fieldOnCells.setName("toto27");
 270         fieldOnCells.setDescription("my wonderful 2D field toto27");
 271         fieldOnCells.setTimeUnit("ms");
 272         fieldOnCells.setMesh(mesh);
 273         array=DataArrayDouble.New();
 274         arr1=[1.2,1.02,1.002,1.0002, 3.4,3.04,3.004,3.0004, 5.6,5.06,5.006,5.0006, 7.8,7.08,7.008,7.0008, 9.1,9.01,9.001,9.0001]
 275         array.setValues(arr1,mesh.getNumberOfCells(),4);
 276         fieldOnCells.setArray(array);
 277         array=DataArrayDouble.New();
 278         arr2=[71.2,71.02,71.002,71.0002, 73.4,73.04,73.004,73.0004, 75.6,75.06,75.006,75.0006, 77.8,77.08,77.008,77.0008, 79.1,79.01,79.001,79.0001]
 279         array.setValues(arr2,mesh.getNumberOfCells(),4);
 280         fieldOnCells.setEndArray(array)
 281         fieldOnCells.setStartTime(6.7,25,26);
 282         fieldOnCells.setEndTime(17.2,125,126);
 283         fieldOnCells.checkCoherency();
 284         return fieldOnCells;
 285
 286     def buildFieldGaussPt2DWT(self):
 287         _a=0.446948490915965;
 288         _b=0.091576213509771;
 289         _p1=0.11169079483905;
 290         _p2=0.0549758718227661;
 291         refCoo1=[ 0.,0., 1.,0., 0.,1. ]
 292         gsCoo1=[ 2*_b-1, 1-4*_b, 2*_b-1, 2.07*_b-1, 1-4*_b,
 293                  2*_b-1, 1-4*_a, 2*_a-1, 2*_a-1, 1-4*_a, 2*_a-1, 2*_a-1]
 294         wg1=[ 4*_p2, 4*_p2, 4*_p2, 4*_p1, 4*_p1, 4*_p1 ]
 295         _refCoo1=refCoo1
 296         _gsCoo1=gsCoo1
 297         _wg1=wg1
 298         m=self.build2DMesh();
 299         f=MEDCouplingFieldDouble.New(ON_GAUSS_PT,ONE_TIME);
 300         f.setTime(6.7,1,4);
 301         f.setMesh(m);
 302         f.setGaussLocalizationOnType(NORM_TRI3,_refCoo1,_gsCoo1,_wg1);
 303         refCoo2=[ 0.,0., 1.,0., 1.,1., 0.,1. ]
 304         _refCoo2=refCoo2
 305         _gsCoo1=_gsCoo1[0:4]
 306         _wg1=_wg1[0:2]
 307         f.setGaussLocalizationOnType(NORM_QUAD4,_refCoo2,_gsCoo1,_wg1);
 308         array=DataArrayDouble.New();
 309         ptr=18*2*[None]
 310         for i in xrange(18*2):
 311             ptr[i]=float(i+1);
 312             pass
 313         array.setValues(ptr,18,2);
 314         array.setInfoOnComponent(0,"Power [MW]");
 315         array.setInfoOnComponent(1,"Density [kg/m^3]");
 316         f.setArray(array);
 317         f.setName("MyFirstFieldOnGaussPoint");
 318         f.setTimeUnit("ms");
 319         f.setDescription("mmmmmmmmmmmm");
 320         f.checkCoherency();
 321         return f;
 322
 323     def buildFieldGaussPtNE2DWT(self):
 324         m=self.build2DMesh();
 325         f=MEDCouplingFieldDouble.New(ON_GAUSS_NE,ONE_TIME);
 326         f.setTime(6.8,11,8);
 327         f.setMesh(m);
 328         f.setTimeUnit("ms");
 329         f.setName("MyFirstFieldOnNE");
 330         f.setDescription("MyDescriptionNE");
 331         array=DataArrayDouble.New();
 332         ptr=18*2*[None]
 333         for i in xrange(18*2):
 334             ptr[i]=float(i+7)
 335         array.setValues(ptr,18,2);
 336         array.setInfoOnComponent(0,"Power [MW]");
 337         array.setInfoOnComponent(1,"Density [kg/m^3]");
 338         f.setArray(array);
 339         #
 340         f.checkCoherency();
 341         return f;
 342
 343     def buildFieldVectorOnExtrudedWT(self):
 344         ext=self.buildExtrudedMesh();
 345         #
 346         f=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 347         f.setTime(6.8,11,8);
 348         f.setMesh(ext);
 349         f.setName("MyFieldOnExtruM");
 350         f.setDescription("desc of MyFiOnExtruM");
 351         array=DataArrayDouble.New();
 352         nbOfCells=ext.getNumberOfCells();
 353         ptr=2*nbOfCells*[None]
 354         for i in range(nbOfCells*2):
 355             ptr[i]=float(i/2+7)+float((i%2)*1000);
 356             pass
 357         array.setValues(ptr,nbOfCells,2);
 358         array.setInfoOnComponent(0,"Power [MW]");
 359         array.setInfoOnComponent(1,"Density [kg/m^3]");
 360         f.setArray(array);
 361         #
 362         f.checkCoherency();
 363         return f
 364
 365     def buildFieldVectorOnCMeshWT(self):
 366         ext=self.buildCMesh();
 367         #
 368         f=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 369         f.setTime(6.8,11,8);
 370         f.setMesh(ext);
 371         f.setName("MyFieldOnCMesh");
 372         f.setDescription("desc of MyFiOnCMesh");
 373         array=DataArrayDouble.New();
 374         nbOfCells=ext.getNumberOfCells();
 375         ptr=2*nbOfCells*[None]
 376         for i in range(nbOfCells*2):
 377             ptr[i]=float(i/2+7)+float((i%2)*1000);
 378             pass
 379         array.setValues(ptr,nbOfCells,2);
 380         array.setInfoOnComponent(0,"Power [GW]");
 381         array.setInfoOnComponent(1,"Density [kg/m^3]");
 382         f.setArray(array);
 383         #
 384         f.checkCoherency();
 385         return f
 386
 387     def buildFieldTemplateCellOn2D(self):
 388         f1=self.buildFieldScalarOn2DNT();
 389         f2=MEDCouplingFieldTemplate.New(f1);
 390         f2.setNature(NoNature);
 391         return f2
 392
 393     def buildFieldTemplateNodeOn2D(self):
 394         f1=self.buildFieldNodeScalarOn2DNT();
 395         f2=MEDCouplingFieldTemplate.New(f1);
 396         f2.setNature(ConservativeVolumic);
 397         return f2
 398
 399     def buildFieldTemplateGaussPtOn2D(self):
 400         f1=self.buildFieldGaussPt2DWT();
 401         f2=MEDCouplingFieldTemplate.New(f1);
 402         f2.setNature(Integral);
 403         return f2
 404
 405     def buildFieldTemplateGaussNEOn2D(self):
 406         f1=self.buildFieldGaussPtNE2DWT();
 407         f2=MEDCouplingFieldTemplate.New(f1);
 408         f2.setNature(IntegralGlobConstraint);
 409         return f2
 410
 411     def buildMultiFields1(self):
 412         m1=self.build2DMesh();
 413         m1.setName("m1");
 414         m2=self.build2DMesh();
 415         m2.setName("m2");
 416         vals0=[-0.7,-1.,-2.,-3.,-4.];
 417         vals1=[0.,1.,2.,3.,4.,0.1,0.2,0.3,0.4];
 418         vals1_1=[170.,171.,172.,173.,174.,170.1,170.2,170.3,170.4];
 419         vals2=[5.,6.,7.,8.,9.];
 420         vals4=[15.,16.,17.,18.,19.];
 421         d0=DataArrayDouble.New();
 422         d0.setValues(vals0,5,1);
 423         d1=DataArrayDouble.New();
 424         d1.setValues(vals1,9,1);
 425         d1_1=DataArrayDouble.New();
 426         d1_1.setValues(vals1_1,9,1);
 427         d2=DataArrayDouble.New();
 428         d2.setValues(vals2,5,1);
 429         d4=DataArrayDouble.New();
 430         d4.setValues(vals4,5,1);
 431         d0.setName("d0"); d1.setName("d1"); d1_1.setName("d1_1"); d2.setName("d2"); d4.setName("d4");
 432         d0.setInfoOnComponent(0,"c1");
 433         d1.setInfoOnComponent(0,"c6");
 434         d1_1.setInfoOnComponent(0,"c9");
 435         d2.setInfoOnComponent(0,"c5");
 436         d4.setInfoOnComponent(0,"c7");
 437         f0=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 438         f0.setMesh(m1);
 439         f0.setArray(d0);
 440         f0.setTime(0.2,5,6);
 441         f0.setName("f0");
 442         f1=MEDCouplingFieldDouble.New(ON_NODES,LINEAR_TIME);
 443         f1.setMesh(m1);
 444         f1.setArrays([d1,d1_1]);
 445         f1.setStartTime(0.7,7,8);
 446         f1.setEndTime(1.2,9,10);
 447         f1.setName("f1");
 448         f2=MEDCouplingFieldDouble.New(ON_CELLS,CONST_ON_TIME_INTERVAL);
 449         f2.setMesh(m2);
 450         f2.setArray(d2);
 451         f2.setTime(1.2,11,12);
 452         f2.setEndTime(1.5,13,14);
 453         f2.setName("f2");
 454         f3=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 455         f3.setMesh(m1);
 456         f3.setArray(d2);
 457         f3.setTime(1.7,15,16);
 458         f3.setName("f3");
 459         f4=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
 460         f4.setMesh(m2);
 461         f4.setArray(d4);
 462         f4.setName("f4");
 463         ret=MEDCouplingMultiFields.New([f0,f1,f2,f3,f4]);
 464         return ret;
 465
 466     def buildArrayDouble1(self):
 467         ret=DataArrayDouble.New();
 468         vals=[2.4,3.2,5.6,9.6,47.6,20.4,24.6,278.1,2.01,3.3,2.4,9.4];
 469         ret.setValues(vals,4,3);
 470         ret.setName("toto");
 471         ret.setInfoOnComponent(0,"sss");
 472         ret.setInfoOnComponent(1,"ppp");
 473         ret.setInfoOnComponent(2,"ttt");
 474         return ret;
 475
 476     def buildArrayDouble2(self):
 477         ret=DataArrayDouble.New();
 478         ret.setName("titi");
 479         return ret;
 480
 481     def buildArrayDouble3(self):
 482         ret=DataArrayDouble.New();
 483         ret.setName("titi");
 484         ret.alloc(0,3);
 485         ret.setInfoOnComponent(0,"sss");
 486         ret.setInfoOnComponent(1,"ppp");
 487         ret.setInfoOnComponent(2,"ttt");
 488         return ret;
 489
 490     def buildArrayInt1(self):
 491         ret=DataArrayInt.New();
 492         vals=[2,3,5,9,47,20,24,278,2,3,2,9];
 493         ret.setValues(vals,4,3);
 494         ret.setName("toto");
 495         ret.setInfoOnComponent(0,"sss");
 496         ret.setInfoOnComponent(1,"ppp");
 497         ret.setInfoOnComponent(2,"ttt");
 498         return ret;
 499
 500     def buildArrayInt2(self):
 501         ret=DataArrayInt.New();
 502         ret.setName("titi");
 503         return ret;
 504
 505     def buildArrayInt3(self):
 506         ret=DataArrayInt.New();
 507         ret.setName("titi");
 508         ret.alloc(0,3);
 509         ret.setInfoOnComponent(0,"sss");
 510         ret.setInfoOnComponent(1,"ppp");
 511         ret.setInfoOnComponent(2,"ttt");
 512         return ret;
 513
 514     def buildMultiFields2(self):
 515         m1=self.build2DMesh();
 516         m1.setName("m1");
 517         m2=self.build2DMesh();
 518         m2.setName("m2");
 519         vals0=[-0.7,-1.,-2.,-3.,-4.];
 520         vals1=[0.,1.,2.,3.,4.,0.1,0.2,0.3,0.4];
 521         vals1_1=[170.,171.,172.,173.,174.];
 522         vals2=[5.,6.,7.,8.,9.];
 523         vals4=[15.,16.,17.,18.,19.];
 524         d0=DataArrayDouble.New();
 525         d0.setValues(vals0,5,1);
 526         d1=DataArrayDouble.New();
 527         d1.setValues(vals1[:5],5,1);
 528         d1_1=DataArrayDouble.New();
 529         d1_1.setValues(vals1_1,5,1);
 530         d2=DataArrayDouble.New();
 531         d2.setValues(vals2,5,1);
 532         d4=DataArrayDouble.New();
 533         d4.setValues(vals4,5,1);
 534         d0.setName("d0"); d1.setName("d1"); d1_1.setName("d1_1"); d2.setName("d2"); d4.setName("d4");
 535         d0.setInfoOnComponent(0,"c1");
 536         d1.setInfoOnComponent(0,"c6");
 537         d1_1.setInfoOnComponent(0,"c9");
 538         d2.setInfoOnComponent(0,"c5");
 539         d4.setInfoOnComponent(0,"c7");
 540         f0=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 541         f0.setMesh(m1);
 542         f0.setArray(d0);
 543         f0.setTime(0.2,5,6);
 544         f0.setName("f0");
 545         f1=MEDCouplingFieldDouble.New(ON_CELLS,LINEAR_TIME);
 546         f1.setMesh(m1);
 547         f1.setArrays([d1,d1_1]);
 548         f1.setStartTime(0.7,7,8);
 549         f1.setEndTime(1.2,9,10);
 550         f1.setName("f1");
 551         f2=MEDCouplingFieldDouble.New(ON_CELLS,CONST_ON_TIME_INTERVAL);
 552         f2.setMesh(m2);
 553         f2.setArray(d2);
 554         f2.setTime(1.2,11,12);
 555         f2.setEndTime(1.5,13,14);
 556         f2.setName("f2");
 557         f3=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 558         f3.setMesh(m1);
 559         f3.setArray(d2);
 560         f3.setTime(1.7,15,16);
 561         f3.setName("f3");
 562         f4=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
 563         f4.setMesh(m2);
 564         f4.setArray(d4);
 565         f4.setName("f4");
 566         f4.setTime(2.7,25,26);
 567         ret=MEDCouplingFieldOverTime.New([f0,f1,f2,f3,f4]);
 568         ret.checkCoherency();
 569         return ret;
 570
 571     def buildFileNameForIOR(self):
 572         ret=os.getenv("TMP");
 573         ret+="/entryPointMEDCouplingCorba.ior";
 574         return ret;
 575     pass
 576
 577 def testMesh():
 578     tab4=[1, 2, 8, 7, 2, 3, 9, 8, 3,
 579           4, 10, 9, 4, 5, 11, 10, 5,
 580           0, 6, 11, 0, 1, 7, 6 ]
 581     nbOfNodes=12
 582     nbOfCells=6
 583     coords=[ 0.024155, 0.04183768725682622, -0.305, 0.04831000000000001, -1.015761910347357e-17,
 584              -0.305, 0.09662000000000001, -1.832979297858306e-18, -0.305, 0.120775, 0.04183768725682623,
 585              -0.305, 0.09662000000000001, 0.08367537451365245, -0.305, 0.04831000000000001, 0.08367537451365246,
 586              -0.305, 0.024155, 0.04183768725682622, -0.2863, 0.04831000000000001, -1.015761910347357e-17, -0.2863,
 587              0.09662000000000001, -1.832979297858306e-18, -0.2863, 0.120775, 0.04183768725682623, -0.2863, 0.09662000000000001,
 588              0.08367537451365245, -0.2863, 0.04831000000000001, 0.08367537451365246, -0.2863 ]
 589     mesh=MEDCouplingUMesh.New()
 590     mesh.setMeshDimension(2)
 591     mesh.allocateCells(8);
 592     mesh.setName("mesh1")
 593     for i in range(nbOfCells):
 594         mesh.insertNextCell(NORM_QUAD4,4,tab4[4*i:4*(i+1)]);
 595         pass
 596     mesh.finishInsertingCells()
 597     myCoords=DataArrayDouble.New()
 598     myCoords.setValues(coords,nbOfNodes,3);
 599     mesh.setCoords(myCoords);
 600     mesh.checkCoherency();
 601     myFalseConn=DataArrayInt.New()
 602     myFalseConn.setValues(tab4,6,4)
 603     return mesh
 604
 605 def testField():
 606     #
 607     mesh=testMesh()
 608     nbOfCells=mesh.getNumberOfCells()
 609     field=MEDCouplingFieldDouble.New(ON_CELLS)
 610     field.setMesh(mesh)
 611     field.setNature(Integral)
 612     myCoords=DataArrayDouble.New()
 613     sampleTab=[]
 614     for i in range(nbOfCells*9):
 615         sampleTab.append(float(i))
 616         myCoords.setValues(sampleTab,nbOfCells,9);
 617         pass
 618     field.setArray(myCoords)
 619     return field
 620