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