Salome HOME
Merge branch 'abn/fix_intersec' into V7_main
[modules/med.git] / src / MEDCouplingCorba_Swig / MEDCouplingCorbaSwigTest.py
1 # Copyright (C) 2007-2014  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, or (at your option) any later version.
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(1,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 build1SGTUMesh(self):
194         targetMesh=MEDCoupling1SGTUMesh("Mesh1SGT",NORM_QUAD4)
195         targetMesh.setTime(2.3,44,-55)
196         targetMesh.setTimeUnit("us")
197         targetMesh.setDescription("My Description of 1SGTU");
198         a1=DataArrayDouble([1.,1.,0.,2.,1.,0.,3.,1.,0.,1.,0.,0.,2.,0.,0.,0.,0.,0.,0.,1.,0.,3.,0.,0.,4.,0.,0.,4.,1.,0.],10,3)
199         a1.setInfoOnComponents(["X1 [m]","YY2 [km]","ZZZ3 [km]"])
200         targetMesh.setCoords(a1)
201         a2=DataArrayInt([6,0,3,5,3,0,1,4,1,2,7,4,8,7,2,9],4*4,1)
202         targetMesh.setNodalConnectivity(a2)
203         #
204         targetMesh.checkCoherency();
205         #
206         return targetMesh;
207
208     def build1DGTUMesh(self):
209         targetMesh=MEDCoupling1DGTUMesh("Mesh1DGT",NORM_POLYGON);
210         targetMesh.setTime(2.3,55,-66)
211         targetMesh.setTimeUnit("us")
212         targetMesh.setDescription("My Description of 1DGTU");
213         a1=DataArrayDouble([1.,1.,0.,2.,1.,0.,3.,1.,0.,1.,0.,0.,2.,0.,0.,0.,0.,0.,0.,1.,0.,3.,0.,0.,4.,0.,0.,4.,1.,0.],10,3)
214         a1.setInfoOnComponents(["X1 [m]","YY2 [km]","ZZZ3 [km]"])
215         targetMesh.setCoords(a1)
216         a2=DataArrayInt([6,0,3,5,3,0,1,4,1,2,7,4,8,7,2],15,1)
217         a3=DataArrayInt([0,4,8,12,15],5,1)
218         targetMesh.setNodalConnectivity(a2,a3)
219         #
220         targetMesh.checkCoherency();
221         #
222         return targetMesh;
223         
224     def buildFieldScalarOn2DNT(self):
225         mesh=self.build2DMesh();
226         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
227         fieldOnCells.setTimeUnit("ms");
228         fieldOnCells.setName("toto");
229         fieldOnCells.setMesh(mesh);
230         array=DataArrayDouble.New();
231         tmp=mesh.getNumberOfCells()*6*[7.]
232         array.setValues(tmp,mesh.getNumberOfCells(),6);
233         fieldOnCells.setArray(array)
234         fieldOnCells.checkCoherency();
235         return fieldOnCells;
236
237     def buildFieldNodeScalarOn2DNT(self):
238         mesh=self.build2DMesh();
239         fieldOnNodes=MEDCouplingFieldDouble.New(ON_NODES,NO_TIME);
240         fieldOnNodes.setName("toto2");
241         fieldOnNodes.setTimeUnit("s");
242         fieldOnNodes.setDescription("my wonderful field toto2");
243         fieldOnNodes.setMesh(mesh);
244         array=DataArrayDouble.New();
245         tmp=mesh.getNumberOfNodes()*5*[7.1234]
246         array.setValues(tmp,mesh.getNumberOfNodes(),5);
247         fieldOnNodes.setArray(array);
248         fieldOnNodes.checkCoherency();
249         return fieldOnNodes;
250
251     def buildFieldScalarOn3DNT(self):
252         mesh=self.build3DMesh();
253         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
254         fieldOnCells.setNature(ConservativeVolumic);
255         fieldOnCells.setName("toto");
256         fieldOnCells.setDescription("my wonderful 3D field toto2");
257         fieldOnCells.setMesh(mesh);
258         array=DataArrayDouble.New();
259         tmp=mesh.getNumberOfCells()*6*[7.]
260         array.setValues(tmp,mesh.getNumberOfCells(),6);
261         fieldOnCells.setArray(array);
262         fieldOnCells.checkCoherency();
263         return fieldOnCells;
264
265     def buildFieldScalarOn3DSurfWT(self):
266         mesh=self.build3DSurfMesh();
267         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
268         fieldOnCells.setName("toto25");
269         fieldOnCells.setDescription("my wonderful 3D surf field toto25");
270         fieldOnCells.setTimeUnit("us");
271         fieldOnCells.setMesh(mesh);
272         array=DataArrayDouble.New();
273         tmp=mesh.getNumberOfCells()*3*[7.]
274         array.setValues(tmp,mesh.getNumberOfCells(),3);
275         array.setInfoOnComponents(["aaa","bbbb","ccccc"])
276         fieldOnCells.setArray(array);
277         fieldOnCells.setTime(6.7,1,4);
278         fieldOnCells.checkCoherency();
279         return fieldOnCells;
280
281     def buildFieldScalarOn3DSurfCOTI(self):
282         mesh=self.build3DSurfMesh();
283         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,CONST_ON_TIME_INTERVAL);
284         fieldOnCells.setName("toto26");
285         fieldOnCells.setDescription("my wonderful 3D surf field toto26");
286         fieldOnCells.setTimeUnit("us");
287         fieldOnCells.setMesh(mesh);
288         array=DataArrayDouble.New();
289         tmp=mesh.getNumberOfCells()*3*[7.]
290         array.setValues(tmp,mesh.getNumberOfCells(),3);
291         fieldOnCells.setArray(array);
292         fieldOnCells.setStartTime(6.7,1,4);
293         fieldOnCells.setEndTime(7.2,2,8);
294         fieldOnCells.checkCoherency();
295         return fieldOnCells;
296
297     def buildFieldScalarOn2DLT(self):
298         mesh=self.build2DMesh()
299         fieldOnCells=MEDCouplingFieldDouble.New(ON_CELLS,LINEAR_TIME)
300         fieldOnCells.setName("toto27");
301         fieldOnCells.setDescription("my wonderful 2D field toto27");
302         fieldOnCells.setTimeUnit("ms");
303         fieldOnCells.setMesh(mesh);
304         array=DataArrayDouble.New();
305         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]
306         array.setValues(arr1,mesh.getNumberOfCells(),4);
307         fieldOnCells.setArray(array);
308         array=DataArrayDouble.New();
309         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]
310         array.setValues(arr2,mesh.getNumberOfCells(),4);
311         fieldOnCells.setEndArray(array)
312         fieldOnCells.setStartTime(6.7,25,26);
313         fieldOnCells.setEndTime(17.2,125,126);
314         fieldOnCells.checkCoherency();
315         return fieldOnCells;
316
317     def buildFieldGaussPt2DWT(self):
318         _a=0.446948490915965;
319         _b=0.091576213509771;
320         _p1=0.11169079483905;
321         _p2=0.0549758718227661;
322         refCoo1=[ 0.,0., 1.,0., 0.,1. ]
323         gsCoo1=[ 2*_b-1, 1-4*_b, 2*_b-1, 2.07*_b-1, 1-4*_b,
324                  2*_b-1, 1-4*_a, 2*_a-1, 2*_a-1, 1-4*_a, 2*_a-1, 2*_a-1]
325         wg1=[ 4*_p2, 4*_p2, 4*_p2, 4*_p1, 4*_p1, 4*_p1 ]
326         _refCoo1=refCoo1
327         _gsCoo1=gsCoo1
328         _wg1=wg1
329         m=self.build2DMesh();
330         f=MEDCouplingFieldDouble.New(ON_GAUSS_PT,ONE_TIME);
331         f.setTime(6.7,1,4);
332         f.setMesh(m);
333         f.setGaussLocalizationOnType(NORM_TRI3,_refCoo1,_gsCoo1,_wg1);
334         refCoo2=[ 0.,0., 1.,0., 1.,1., 0.,1. ]
335         _refCoo2=refCoo2
336         _gsCoo1=_gsCoo1[0:4]
337         _wg1=_wg1[0:2]
338         f.setGaussLocalizationOnType(NORM_QUAD4,_refCoo2,_gsCoo1,_wg1);
339         array=DataArrayDouble.New();
340         ptr=18*2*[None]
341         for i in xrange(18*2):
342             ptr[i]=float(i+1);
343             pass
344         array.setValues(ptr,18,2);
345         array.setInfoOnComponent(0,"Power [MW]");
346         array.setInfoOnComponent(1,"Density [kg/m^3]");
347         f.setArray(array);
348         f.setName("MyFirstFieldOnGaussPoint");
349         f.setTimeUnit("ms");
350         f.setDescription("mmmmmmmmmmmm");
351         f.checkCoherency();
352         return f;
353
354     def buildFieldGaussPtNE2DWT(self):
355         m=self.build2DMesh();
356         f=MEDCouplingFieldDouble.New(ON_GAUSS_NE,ONE_TIME);
357         f.setTime(6.8,11,8);
358         f.setMesh(m);
359         f.setTimeUnit("ms");
360         f.setName("MyFirstFieldOnNE");
361         f.setDescription("MyDescriptionNE");
362         array=DataArrayDouble.New();
363         ptr=18*2*[None]
364         for i in xrange(18*2):
365             ptr[i]=float(i+7)
366         array.setValues(ptr,18,2);
367         array.setInfoOnComponent(0,"Power [MW]");
368         array.setInfoOnComponent(1,"Density [kg/m^3]");
369         f.setArray(array);
370         #
371         f.checkCoherency();
372         return f;
373
374     def buildFieldVectorOnExtrudedWT(self):
375         ext=self.buildExtrudedMesh();
376         #
377         f=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
378         f.setTime(6.8,11,8);
379         f.setMesh(ext);
380         f.setName("MyFieldOnExtruM");
381         f.setDescription("desc of MyFiOnExtruM");
382         array=DataArrayDouble.New();
383         nbOfCells=ext.getNumberOfCells();
384         ptr=2*nbOfCells*[None]
385         for i in range(nbOfCells*2):
386             ptr[i]=float(i/2+7)+float((i%2)*1000);
387             pass
388         array.setValues(ptr,nbOfCells,2);
389         array.setInfoOnComponent(0,"Power [MW]");
390         array.setInfoOnComponent(1,"Density [kg/m^3]");
391         f.setArray(array);
392         #
393         f.checkCoherency();
394         return f
395
396     def buildFieldVectorOnCMeshWT(self):
397         ext=self.buildCMesh();
398         #
399         f=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
400         f.setTime(6.8,11,8);
401         f.setMesh(ext);
402         f.setName("MyFieldOnCMesh");
403         f.setDescription("desc of MyFiOnCMesh");
404         array=DataArrayDouble.New();
405         nbOfCells=ext.getNumberOfCells();
406         ptr=2*nbOfCells*[None]
407         for i in range(nbOfCells*2):
408             ptr[i]=float(i/2+7)+float((i%2)*1000);
409             pass
410         array.setValues(ptr,nbOfCells,2);
411         array.setInfoOnComponent(0,"Power [GW]");
412         array.setInfoOnComponent(1,"Density [kg/m^3]");
413         f.setArray(array);
414         #
415         f.checkCoherency();
416         return f
417
418     def buildFieldTemplateCellOn2D(self):
419         f1=self.buildFieldScalarOn2DNT();
420         f2=MEDCouplingFieldTemplate.New(f1);
421         f2.setNature(NoNature);
422         return f2
423     
424     def buildFieldTemplateNodeOn2D(self):
425         f1=self.buildFieldNodeScalarOn2DNT();
426         f2=MEDCouplingFieldTemplate.New(f1);
427         f2.setNature(ConservativeVolumic);
428         return f2
429             
430     def buildFieldTemplateGaussPtOn2D(self):
431         f1=self.buildFieldGaussPt2DWT();
432         f2=MEDCouplingFieldTemplate.New(f1);
433         f2.setNature(Integral);
434         return f2
435
436     def buildFieldTemplateGaussNEOn2D(self):
437         f1=self.buildFieldGaussPtNE2DWT();
438         f2=MEDCouplingFieldTemplate.New(f1);
439         f2.setNature(IntegralGlobConstraint);
440         return f2
441
442     def buildMultiFields1(self):
443         m1=self.build2DMesh();
444         m1.setName("m1");
445         m2=self.build2DMesh();
446         m2.setName("m2");
447         vals0=[-0.7,-1.,-2.,-3.,-4.];
448         vals1=[0.,1.,2.,3.,4.,0.1,0.2,0.3,0.4];
449         vals1_1=[170.,171.,172.,173.,174.,170.1,170.2,170.3,170.4];
450         vals2=[5.,6.,7.,8.,9.];
451         vals4=[15.,16.,17.,18.,19.];
452         d0=DataArrayDouble.New();
453         d0.setValues(vals0,5,1);
454         d1=DataArrayDouble.New();
455         d1.setValues(vals1,9,1);
456         d1_1=DataArrayDouble.New();
457         d1_1.setValues(vals1_1,9,1);
458         d2=DataArrayDouble.New();
459         d2.setValues(vals2,5,1);
460         d4=DataArrayDouble.New();
461         d4.setValues(vals4,5,1);
462         d0.setName("d0"); d1.setName("d1"); d1_1.setName("d1_1"); d2.setName("d2"); d4.setName("d4");
463         d0.setInfoOnComponent(0,"c1");
464         d1.setInfoOnComponent(0,"c6");
465         d1_1.setInfoOnComponent(0,"c9");
466         d2.setInfoOnComponent(0,"c5");
467         d4.setInfoOnComponent(0,"c7");
468         f0=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
469         f0.setMesh(m1);
470         f0.setArray(d0);
471         f0.setTime(0.2,5,6);
472         f0.setName("f0");
473         f1=MEDCouplingFieldDouble.New(ON_NODES,LINEAR_TIME);
474         f1.setMesh(m1);
475         f1.setArrays([d1,d1_1]);
476         f1.setStartTime(0.7,7,8);
477         f1.setEndTime(1.2,9,10);
478         f1.setName("f1");
479         f2=MEDCouplingFieldDouble.New(ON_CELLS,CONST_ON_TIME_INTERVAL);
480         f2.setMesh(m2);
481         f2.setArray(d2);
482         f2.setTime(1.2,11,12);
483         f2.setEndTime(1.5,13,14);
484         f2.setName("f2");
485         f3=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
486         f3.setMesh(m1);
487         f3.setArray(d2);
488         f3.setTime(1.7,15,16);
489         f3.setName("f3");
490         f4=MEDCouplingFieldDouble.New(ON_CELLS,NO_TIME);
491         f4.setMesh(m2);
492         f4.setArray(d4);
493         f4.setName("f4");
494         ret=MEDCouplingMultiFields.New([f0,f1,f2,f3,f4]);
495         return ret;
496     
497     def buildArrayDouble1(self):
498         ret=DataArrayDouble.New();
499         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];
500         ret.setValues(vals,4,3);
501         ret.setName("toto");
502         ret.setInfoOnComponent(0,"sss");
503         ret.setInfoOnComponent(1,"ppp");
504         ret.setInfoOnComponent(2,"ttt");
505         return ret;
506
507     def buildArrayDouble2(self):
508         ret=DataArrayDouble.New();
509         ret.setName("titi");
510         return ret;
511
512     def buildArrayDouble3(self):
513         ret=DataArrayDouble.New();
514         ret.setName("titi");
515         ret.alloc(0,3);
516         ret.setInfoOnComponent(0,"sss");
517         ret.setInfoOnComponent(1,"ppp");
518         ret.setInfoOnComponent(2,"ttt");
519         return ret;
520
521     def buildArrayInt1(self):
522         ret=DataArrayInt.New();
523         vals=[2,3,5,9,47,20,24,278,2,3,2,9];
524         ret.setValues(vals,4,3);
525         ret.setName("toto");
526         ret.setInfoOnComponent(0,"sss");
527         ret.setInfoOnComponent(1,"ppp");
528         ret.setInfoOnComponent(2,"ttt");
529         return ret;
530
531     def buildArrayInt2(self):
532         ret=DataArrayInt.New();
533         ret.setName("titi");
534         return ret;
535
536     def buildArrayInt3(self):
537         ret=DataArrayInt.New();
538         ret.setName("titi");
539         ret.alloc(0,3);
540         ret.setInfoOnComponent(0,"sss");
541         ret.setInfoOnComponent(1,"ppp");
542         ret.setInfoOnComponent(2,"ttt");
543         return ret;
544
545     def buildMultiFields2(self):
546         m1=self.build2DMesh();
547         m1.setName("m1");
548         m2=self.build2DMesh();
549         m2.setName("m2");
550         vals0=[-0.7,-1.,-2.,-3.,-4.];
551         vals1=[0.,1.,2.,3.,4.,0.1,0.2,0.3,0.4];
552         vals1_1=[170.,171.,172.,173.,174.];
553         vals2=[5.,6.,7.,8.,9.];
554         vals4=[15.,16.,17.,18.,19.];
555         d0=DataArrayDouble.New();
556         d0.setValues(vals0,5,1);
557         d1=DataArrayDouble.New();
558         d1.setValues(vals1[:5],5,1);
559         d1_1=DataArrayDouble.New();
560         d1_1.setValues(vals1_1,5,1);
561         d2=DataArrayDouble.New();
562         d2.setValues(vals2,5,1);
563         d4=DataArrayDouble.New();
564         d4.setValues(vals4,5,1);
565         d0.setName("d0"); d1.setName("d1"); d1_1.setName("d1_1"); d2.setName("d2"); d4.setName("d4");
566         d0.setInfoOnComponent(0,"c1");
567         d1.setInfoOnComponent(0,"c6");
568         d1_1.setInfoOnComponent(0,"c9");
569         d2.setInfoOnComponent(0,"c5");
570         d4.setInfoOnComponent(0,"c7");
571         f0=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
572         f0.setMesh(m1);
573         f0.setArray(d0);
574         f0.setTime(0.2,5,6);
575         f0.setName("f0");
576         f1=MEDCouplingFieldDouble.New(ON_CELLS,LINEAR_TIME);
577         f1.setMesh(m1);
578         f1.setArrays([d1,d1_1]);
579         f1.setStartTime(0.7,7,8);
580         f1.setEndTime(1.2,9,10);
581         f1.setName("f1");
582         f2=MEDCouplingFieldDouble.New(ON_CELLS,CONST_ON_TIME_INTERVAL);
583         f2.setMesh(m2);
584         f2.setArray(d2);
585         f2.setTime(1.2,11,12);
586         f2.setEndTime(1.5,13,14);
587         f2.setName("f2");
588         f3=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
589         f3.setMesh(m1);
590         f3.setArray(d2);
591         f3.setTime(1.7,15,16);
592         f3.setName("f3");
593         f4=MEDCouplingFieldDouble.New(ON_CELLS,ONE_TIME);
594         f4.setMesh(m2);
595         f4.setArray(d4);
596         f4.setName("f4");
597         f4.setTime(2.7,25,26);
598         ret=MEDCouplingFieldOverTime.New([f0,f1,f2,f3,f4]);
599         ret.checkCoherency();
600         return ret;
601
602     def buildFileNameForIOR(self):
603         tmpdir=os.getenv("TMP", "/tmp");
604         username="";
605         if os.getenv("USERNAME"):
606             username = os.getenv("USERNAME")+"_";
607         elif os.getenv("USER"):
608             username = os.getenv("USER")+"_";
609         ret=tmpdir+"/"+username+"entryPointMEDCouplingCorba.ior";
610         return ret;
611     pass
612
613 def testMesh():
614     tab4=[1, 2, 8, 7, 2, 3, 9, 8, 3,
615           4, 10, 9, 4, 5, 11, 10, 5,
616           0, 6, 11, 0, 1, 7, 6 ]
617     nbOfNodes=12
618     nbOfCells=6
619     coords=[ 0.024155, 0.04183768725682622, -0.305, 0.04831000000000001, -1.015761910347357e-17,
620              -0.305, 0.09662000000000001, -1.832979297858306e-18, -0.305, 0.120775, 0.04183768725682623,
621              -0.305, 0.09662000000000001, 0.08367537451365245, -0.305, 0.04831000000000001, 0.08367537451365246,
622              -0.305, 0.024155, 0.04183768725682622, -0.2863, 0.04831000000000001, -1.015761910347357e-17, -0.2863, 
623              0.09662000000000001, -1.832979297858306e-18, -0.2863, 0.120775, 0.04183768725682623, -0.2863, 0.09662000000000001,
624              0.08367537451365245, -0.2863, 0.04831000000000001, 0.08367537451365246, -0.2863 ]
625     mesh=MEDCouplingUMesh.New()
626     mesh.setMeshDimension(2)
627     mesh.allocateCells(8);
628     mesh.setName("mesh1")
629     for i in range(nbOfCells):
630         mesh.insertNextCell(NORM_QUAD4,4,tab4[4*i:4*(i+1)]);
631         pass
632     mesh.finishInsertingCells()
633     myCoords=DataArrayDouble.New()
634     myCoords.setValues(coords,nbOfNodes,3);
635     mesh.setCoords(myCoords);
636     mesh.checkCoherency();
637     myFalseConn=DataArrayInt.New()
638     myFalseConn.setValues(tab4,6,4)
639     return mesh
640
641 def testField():
642     #
643     mesh=testMesh()
644     nbOfCells=mesh.getNumberOfCells()
645     field=MEDCouplingFieldDouble.New(ON_CELLS)
646     field.setMesh(mesh)
647     field.setNature(Integral)
648     myCoords=DataArrayDouble.New()
649     sampleTab=[]
650     for i in range(nbOfCells*9):
651         sampleTab.append(float(i))
652         myCoords.setValues(sampleTab,nbOfCells,9);
653         pass
654     field.setArray(myCoords)
655     return field
656