MED fichier
test10_f32_i32.c
Aller à la documentation de ce fichier.
1/* This file is part of MED.
2 *
3 * COPYRIGHT (C) 1999 - 2020 EDF R&D, CEA/DEN
4 * MED is free software: you can redistribute it and/or modify
5 * it under the terms of the GNU Lesser General Public License as published by
6 * the Free Software Foundation, either version 3 of the License, or
7 * (at your option) any later version.
8 *
9 * MED 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
12 * GNU Lesser General Public License for more details.
13 *
14 * You should have received a copy of the GNU Lesser General Public License
15 * along with MED. If not, see <http://www.gnu.org/licenses/>.
16 */
17#define _a 0.446948490915965
18#define _b 0.091576213509771
19#define _p1 0.11169079483905
20#define _p2 0.0549758718227661
21
22/******************************************************************************
23 * - Nom du fichier : test10.c
24 *
25 * - Description : ecriture de champs de resultats MED
26 *
27 *****************************************************************************/
28
29#include <med.h>
30#define MESGERR 1
31#include "med_utils.h"
32#include <string.h>
33
34#ifdef DEF_LECT_ECR
35#define MODE_ACCES MED_ACC_RDWR
36#elif DEF_LECT_AJOUT
37#define MODE_ACCES MED_ACC_RDEXT
38#else
39#define MODE_ACCES MED_ACC_CREAT
40#endif
41
42#ifndef USER_INTERLACE
43#define USER_INTERLACE MED_FULL_INTERLACE
44#endif
45
46#define USER_MODE MED_COMPACT_STMODE
47
48#ifdef DEF_MED_FLOAT32
49#define FTYPECHA MED_FLOAT32
50#define ftypecha med_float32
51#define filename_ "test10_f32"
52#else
53#define FTYPECHA MED_DOUBLE
54#define ftypecha med_double
55#define filename_ "test10"
56#endif
57
58#ifdef DEF_MED_INT32
59#define ITYPECHA MED_INT32
60#define itypecha med_int32
61#define filename__ filename_ "_i32"
62#elif DEF_MED_INT64
63#define ITYPECHA MED_INT64
64#define itypecha med_int64
65#define filename__ filename_ "_i64"
66#else
67#define ITYPECHA MED_INT
68#define itypecha med_int
69#define filename__ filename_
70#endif
71
72
73#define filename filename__ ".med"
74
75int main (int argc, char **argv)
76
77
78{
79 med_err ret=0;
80 med_idt fid;
81
82
83
84 /* Maillage support aux champs*/
85 /* Ces maillages sont vides*/
86 char maa1[MED_NAME_SIZE+1]= "maa1";
87 char maa2[MED_NAME_SIZE+1]= "maa2";
88 char * lien_maa2 = "./testfoo.med";
89 char maa3[MED_NAME_SIZE+1]= "maa3";
90
91
92 /* Caractéristiques du champ n° 1 sur TRIA6 */
93 char nomcha1[MED_NAME_SIZE+1] = "champ reel";
94 char comp1[2*MED_SNAME_SIZE+1] = "comp1 comp2 ";
95 /*12345678901234561234567890123456*/
96 char unit1[2*MED_SNAME_SIZE+1] = "unit1 unit2 ";
97 med_int ncomp1 = 2;
98 /* Caractéristiques du model n° 1 de localisation des points de gauss pour le champ n°1*/
99 med_int ngauss1_1 = 6;
100 char gauss1_1[MED_NAME_SIZE+1] = "Model n1";
101 med_float refcoo1[12] = { -1.0,1.0, -1.0,-1.0, 1.0,-1.0, -1.0,0.0, 0.0,-1.0, 0.0,0.0 };
102
103 /* Constantes */
104
105 med_float gscoo1_1[12] = { 2*_b-1, 1-4*_b, 2*_b-1, 2*_b-1, 1-4*_b,
106 2*_b-1, 1-4*_a, 2*_a-1, 2*_a-1, 1-4*_a, 2*_a-1, 2*_a-1 };
107 med_float wg1_1[6] = { 4*_p2, 4*_p2, 4*_p2, 4*_p1, 4*_p1, 4*_p1 };
108
109 med_int nval1_1= 1*6; /*1 valeurs et 6 points de gauss par valeur */
110 med_int _nent1_1= 1; /*1 valeurs et 6 points de gauss par valeur */
111 ftypecha valr1_1[1*6*2] = {0.0,1.0, 2.0,3.0, 10.0,11.0, 12.0,13.0, 20.0,21.0, 22.0,23.0}; /* 2 composantes*/
112 /* Caractéristiques du model n° 2 de localisation des points de gauss pour le champ n°1*/
113 med_int ngauss1_2 = 3;
114 char gauss1_2[MED_NAME_SIZE+1] = "Model n2";
115 med_float gscoo1_2[6] = { -2.0/3,1.0/3, -2.0/3,-2.0/3, 1.0/3,-2.0/3 };
116 med_float wg1_2[3] = { 2.0/3, 2.0/3, 2.0/3 };
117 med_int nval1_2= 2*3; /*2 valeurs et 3 points de gauss par valeur */
118 med_int _nent1_2= 2; /*2 valeurs et 3 points de gauss par valeur */
119 ftypecha valr1_2[2*3*2] = {0.0,1.0, 2.0,3.0, 10.0,11.0, 12.0,13.0, 20.0,21.0, 22.0,23.0}; /* 2 composantes*/
120 ftypecha valr1_2p[2*3*2] = { 12.0,13.0, 20.0,21.0, 22.0,23.0}; /* 2 composantes*/
121 /* Caractéristiques du model n° 3 sans points de gauss pour le champ n°1*/
122 med_int nval1_3= 6; /*6 valeurs et pas de points de gauss */
123 med_int _nent1_3= 6; /*6 valeurs et pas de points de gauss */
124 ftypecha valr1_3[2*3*2] = {0.0,1.0, 2.0,3.0, 10.0,11.0, 12.0,13.0, 20.0,21.0, 22.0,23.0}; /* 2 composantes*/
125 ftypecha valr1_3p[2*2*2] = { 2.0,3.0, 10.0,11.0 }; /* 2 composantes profil1 */
126
127 /* Caractéristiques du champ n° 2 */
128 char nomcha2[MED_NAME_SIZE+1] = "champ entier";
129 char comp2[3*MED_SNAME_SIZE+1] = "comp1 comp2 comp3 ";
130 /*123456789012345612345678901234561234567890123456*/
131 char unit2[3*MED_SNAME_SIZE+1] = "unit1 unit2 unit3 ";
132 med_int ncomp2 = 3;
133 med_int nval2 = 5; /*5 valeurs */
134 itypecha valr2[5*3 ] = {0,1,2, 10,11,12, 20,21,22, 30,31,32, 40,41,42}; /* 3 composantes*/
135 /* u32b : 4294967296 s32b : 2147483647 -2147483647 */
136 itypecha valr2_[5*3 ] = {2147483647, -2147483648, 4294967296, 10,11,12, 20,21,22, 30,31,32, 40,41,42};
137 itypecha valr2p[3*3 ] = {0,1,2, 20,21,22, 40,41,42}; /* 3 composantes*/
138
139 /* Profils utilisés */
140 char nomprofil1[MED_NAME_SIZE+1] = "PROFIL(champ(1))";
141 char nomprofil1b[MED_NAME_SIZE+1] = "PROFIL(champ(1b))";
142 char nomprofil2[MED_NAME_SIZE+1] = "PROFIL(champ2)";
143 med_int profil1[2] = { 2, 3 };
144 med_int profil2[3] = { 1, 3, 5 };
145
146
147 /* Caractéristiques du champ n° 3 */
148 char nomcha3[MED_NAME_SIZE+1] = "champ entier 3";
149 char comp3[3*MED_SNAME_SIZE+1] = "comp1 comp2 ";
150 /*123456789012345612345678901234561234567890123456*/
151 char unit3[3*MED_SNAME_SIZE+1] = "unit1 unit2 ";
152 char dtunit[MED_SNAME_SIZE+1] = "s";
153 med_int ncomp3 = 2;
154 med_int nval3 = 5*4; /*5 valeurs et 4 noeuds par element*/
155 med_int _nent3 = 5; /*5 valeurs et 4 noeuds par element*/
156 itypecha valr3[5*4*2] = {0,1, 10,11, 20,21, 30,31,
157 40,41, 50,51, 60,61, 70,71,
158 80,81, 90,91, 100,101, 110,111,
159 120,121, 130,131, 140,141, 150,151,
160 160,161, 170,171, 180,181, 190,191}; /* 2 composantes*/
161 itypecha valr3p[3*4*2] = {0,1, 10,11, 20,21, 30,31,
162 80,81, 90,91, 100,101, 110,111,
163 160,161, 170,171, 180,181, 190,191}; /* 2 composantes*/
164
165
166 char nomcoo[3*MED_SNAME_SIZE+1] = "x y z ";
167 char unicoo[3*MED_SNAME_SIZE+1] = "cm cm cm ";
168
169
170
171
172 /* ouverture du fichier */
174 MESSAGE("Erreur à l'ouverture du fichier : ");
175 return -1;
176 }
177
178 /* creation de maa1 de dimension 3*/
179 if (MEDmeshCr( fid, maa1, 3, 3, MED_UNSTRUCTURED_MESH,
180 "Maillage vide","s", MED_SORT_DTIT,
181 MED_CARTESIAN, nomcoo, unicoo) < 0) {
182 MESSAGE("Erreur a la creation du maillage : "); SSCRUTE(maa1);
183 ret = -1;
184 }
185
186
187 /* creation de maa3 de dimension 3*/
188 if (MEDmeshCr( fid, maa3, 3, 3, MED_UNSTRUCTURED_MESH,
189 "Maillage vide","s", MED_SORT_DTIT,
190 MED_CARTESIAN, nomcoo, unicoo) < 0) {
191 MESSAGE("Erreur a la creation du maillage : "); SSCRUTE(maa3);
192 ret = -1;
193 }
194
195
196 /* creation du champ réel n°1 */
197 if ( MEDfieldCr(fid,nomcha1,FTYPECHA,ncomp1,comp1,unit1,dtunit,maa1 ) < 0) {
198 MESSAGE("Erreur à la création du champ : ");SSCRUTE(nomcha1);
199 ret = -1;
200 };
201
202 /* creation du champ entier n°2 */
203 if ( MEDfieldCr(fid,nomcha2,ITYPECHA,ncomp2,comp2,unit2,dtunit,maa2 ) < 0) {
204 MESSAGE("Erreur à la création du champ : ");SSCRUTE(nomcha2);
205 ret = -1;
206 };
207
208 /* creation du lien au fichier distant contenant maa2 */
209 if (MEDlinkWr(fid,maa2,lien_maa2) < 0) {
210 MESSAGE("Erreur à la création du lien : ");SSCRUTE(lien_maa2);
211 ret = -1;
212 };
213
214 /* creation de la localisation des points de Gauss modèle n°1 */
215 if (MEDlocalizationWr(fid, gauss1_1, MED_TRIA6, MED_TRIA6/100, refcoo1, USER_INTERLACE,
216 ngauss1_1, gscoo1_1, wg1_1,
218 MESSAGE("Erreur à la création du modèle n°1 : ");SSCRUTE(gauss1_1);
219 ret = -1;
220 };
221
222 /* creation de la localisation des points de Gauss modèle n°2 */
223 if (MEDlocalizationWr(fid, gauss1_2, MED_TRIA6, MED_TRIA6/100, refcoo1, USER_INTERLACE,
224 ngauss1_2, gscoo1_2, wg1_2,
226 MESSAGE("Erreur à la création du modèle n°1 : ");SSCRUTE(gauss1_2);
227 ret = -1;
228 };
229
230 /* ecriture du champ n°1*/
231 /* enregistre uniquement les composantes n°2 de valr1_1, et n'utilise ni pas de temps ni n° d'ordre*/
232
234 gauss1_1,USER_INTERLACE, 2, _nent1_1, (unsigned char*)valr1_1 ) < 0) {
235 MESSAGE("Erreur à l'écriture du champ : ");
237 SSCRUTE(maa1);
238 ret = -1;
239 };
240
241
242
243 /* enregistre uniquement les composantes n°1 de valr1_1, et n'utilise ni pas de temps ni n° d'ordre */
244
246 gauss1_1,USER_INTERLACE, 1, _nent1_1, (unsigned char*)valr1_1 ) < 0) {
247 MESSAGE("Erreur à l'écriture du champ : ");
249 SSCRUTE(maa1);
250 ret = -1;
251 };
252
253 /* enregistre uniquement les composantes n°1 de valr1_2, au pas de temps n°1(5.5), n'utilise pas de n°d'ordre*/
254 /* ce champ repose sur le maillage maa2 qui est distant */
255
257 gauss1_2,USER_INTERLACE, 1, _nent1_2, (unsigned char*)valr1_2 ) < 0) {
258 MESSAGE("Erreur à l'écriture du champ : ");
260 SSCRUTE(maa1);
261 ret = -1;
262 };
263
264/*Ce test utilise un deuxième maillage pour un même champ, ceci n'existe plus en 3.0*/
265 /* enregistre uniquement les composantes n°2 de valr1_2, au pas de temps n°1(5.5), n'utilise pas de n°d'ordre*/
266 /* ce champ repose sur le maillage maa1 qui est local */
267
268 /*Ce test utilise un deuxième maillage pour un même champ, ceci n'existe plus en 3.0*/
270 gauss1_1,USER_INTERLACE, 2, _nent1_1, (unsigned char*)valr1_1 ) < 0) {
271 MESSAGE("Erreur à l'écriture du champ : ");
273 SSCRUTE(maa1);
274 ret = -1;
275 };
276
277
278 /* enregistre uniquement les composantes n°1 de valr1_1, au pas de temps n°1(5.5), et n°d'itération n°2*/
279 /* ce champ repose sur le maillage maa3 qui est local */
280
282 gauss1_2,USER_INTERLACE, 1, _nent1_2, (unsigned char*)valr1_2 ) < 0) {
283 MESSAGE("Erreur à l'écriture du champ : ");
285 SSCRUTE(maa1);
286 ret = -1;
287 };
288
289 /* Creation d'un profil (selection du deuxieme élément de valr1_1) */
290 /* On n'utilise que la première valeur (2) du profil */
291 if ( MEDprofileWr(fid,nomprofil1,1,profil1) < 0) {
292 MESSAGE("Erreur à l'écriture du profil : ");
293 SSCRUTE(nomprofil1);
294 ret = -1;
295 };
296
297
298 if ( MEDprofileWr(fid,nomprofil1b,1,profil1) < 0) {
299 MESSAGE("Erreur à l'écriture du profil : ");
300 SSCRUTE(nomprofil1b);
301 ret = -1;
302 };
303
304
305
306
307 /* enregistre toutes les composantes du deuxième élément de valr1_1 (premier élément en stockage compact de valr1p),
308 au pas de temps n°2(5.6), et n°d'itération n°2*/
309 if ( MEDfieldValueWithProfileWr(fid, nomcha1,2,2,5.6,MED_CELL,MED_TRIA6,USER_MODE,nomprofil1,
310 MED_NO_LOCALIZATION,USER_INTERLACE, MED_ALL_CONSTITUENT, nval1_3, (unsigned char*)valr1_3p ) < 0) {
311 MESSAGE("Erreur à l'écriture du champ : ");
313 SSCRUTE(maa1);
314 ret = -1;
315 };
316
317
318 /* enregistre toutes les composantes du deuxième élément de valr1_1 (premier élément en stockage compact de valr1p),
319 au pas de temps n°2(5.6), et n°d'itération n°2 */
320
321 if ( MEDfieldValueWithProfileWr(fid, nomcha1,2,2,5.6,MED_CELL,MED_TRIA6,USER_MODE,nomprofil1b,
322 gauss1_2,USER_INTERLACE, MED_ALL_CONSTITUENT, _nent1_2, (unsigned char*)valr1_2p ) < 0) {
323 MESSAGE("Erreur à l'écriture du champ : ");
325 SSCRUTE(maa1);
326 ret = -1;
327 };
328
329
330 if ( MEDfieldValueWithProfileWr(fid, nomcha1,3,2,5.7,MED_CELL,MED_TRIA6,USER_MODE,nomprofil1,
331 MED_NO_LOCALIZATION,USER_INTERLACE, 2, _nent1_3, (unsigned char*)valr1_3p ) < 0) {
332 MESSAGE("Erreur à l'écriture du champ : ");
334 SSCRUTE(maa1);
335 ret = -1;
336 };
337
338
339 /* Ecriture du champ n° 2 */
340
341
342 if ( MEDfieldValueWr(fid, nomcha2,MED_NO_DT,MED_NO_IT,0,
344 USER_INTERLACE, 1, nval2, (unsigned char*)valr2 ) < 0) {
345 MESSAGE("Erreur à l'écriture du champ : ");
347 SSCRUTE(maa1);
348 ret = -1;
349 };
350
351
353 USER_INTERLACE, 2, nval2, (unsigned char*)valr2 ) < 0) {
354 MESSAGE("Erreur à l'écriture du champ : ");
356 SSCRUTE(maa1);
357 ret = -1;
358 };
359
360
362 USER_INTERLACE, 3, nval2, (unsigned char*)valr2 ) < 0) {
363 MESSAGE("Erreur à l'écriture du champ : ");
365 SSCRUTE(maa1);
366 ret = -1;
367 };
368
369 /* Creation d'un profil (selection des éléments 1,3,5 de valr2) */
370 /* On utilise les trois valeurs du profil */
371 if ( MEDprofileWr(fid,nomprofil2,3,profil2) < 0) {
372 MESSAGE("Erreur à l'écriture du profil : ");
373 SSCRUTE(nomprofil2);
374 ret = -1;
375 };
376
377
379 MED_NO_LOCALIZATION,USER_INTERLACE, 3, nval2, (unsigned char*)valr2p ) < 0) {
380 MESSAGE("Erreur à l'écriture du champ : ");
382 SSCRUTE(maa1);
383 ret = -1;
384 };
385
386 /* creation du champ entier n°3 */
387 if ( MEDfieldCr(fid,nomcha3,ITYPECHA,ncomp3,comp3,unit3,dtunit,maa1) < 0) {
388 MESSAGE("Erreur à la création du champ : ");SSCRUTE(nomcha3);
389 ret = -1;
390 };
391
392 /* Ecriture du champ n° 3 */
393
395 USER_INTERLACE, 1, nval3, (unsigned char*)valr3 ) < 0) {
396 MESSAGE("Erreur à l'écriture du champ : ");
398 SSCRUTE(maa1);
399 ret = -1;
400 };
401
403 USER_INTERLACE, MED_ALL_CONSTITUENT, _nent3, (unsigned char*)valr3 ) < 0) {
404 MESSAGE("Erreur à l'écriture du champ : ");
406 SSCRUTE(maa1);
407 ret = -1;
408 };
409
411 MED_NO_LOCALIZATION,USER_INTERLACE, MED_ALL_CONSTITUENT, _nent3, (unsigned char*)valr3p ) < 0) {
412 MESSAGE("Erreur à l'écriture du champ : ");
414 SSCRUTE(maa1);
415 ret = -1;
416 };
417
418
419 /* fermeture du fichier */
420 if ( MEDfileClose(fid) < 0 ) ret=-1;
421
422 return ret;
423}
424
425
426
427
MEDC_EXPORT med_err MEDfieldCr(const med_idt fid, const char *const fieldname, const med_field_type fieldtype, const med_int ncomponent, const char *const componentname, const char *const componentunit, const char *const dtunit, const char *const meshname)
Cette fonction crée un champ dans un fichier.
Definition MEDfieldCr.c:44
MEDC_EXPORT med_err MEDfieldValueWithProfileWr(const med_idt fid, const char *const fieldname, const med_int numdt, const med_int numit, const med_float dt, const med_entity_type entitype, const med_geometry_type geotype, const med_storage_mode storagemode, const char *const profilename, const char *const localizationname, const med_switch_mode switchmode, const med_int componentselect, const med_int nentity, const unsigned char *const value)
Cette fonction permet d'écrire les valeurs d'un champ définies sur des entités d'un maillage pour une...
MEDC_EXPORT med_err MEDfieldValueWr(const med_idt fid, const char *const fieldname, const med_int numdt, const med_int numit, const med_float dt, const med_entity_type entitype, const med_geometry_type geotype, const med_switch_mode switchmode, const med_int componentselect, const med_int nentity, const unsigned char *const value)
Cette fonction permet d'écrire les valeurs d'un champ définies sur des entités d'un maillage pour une...
MEDC_EXPORT med_err MEDfileClose(med_idt fid)
Fermeture d'un fichier MED.
MEDC_EXPORT med_idt MEDfileVersionOpen(const char *const filename, const med_access_mode accessmode, const med_int major, const med_int minor, const med_int release)
Ouverture d'un fichier MED en indiquant la version du modèle à utiliser en cas de création d'un nouve...
MEDC_EXPORT med_err MEDlocalizationWr(const med_idt fid, const char *const localizationname, const med_geometry_type geotype, const med_int spacedimension, const med_float *const elementcoordinate, const med_switch_mode switchmode, const med_int nipoint, const med_float *const ipointcoordinate, const med_float *const weight, const char *const geointerpname, const char *const ipointstructmeshname)
Cette routine permet l'écriture d'une localisation localizationname de points d'intégration dans/auto...
MEDC_EXPORT med_err MEDmeshCr(const med_idt fid, const char *const meshname, const med_int spacedim, const med_int meshdim, const med_mesh_type meshtype, const char *const description, const char *const dtunit, const med_sorting_type sortingtype, const med_axis_type axistype, const char *const axisname, const char *const axisunit)
Cette routine permet de créer un maillage dans un fichier.
Definition MEDmeshCr.c:45
MEDC_EXPORT med_err MEDprofileWr(const med_idt fid, const char *const profilename, const med_int profilesize, const med_int *const profilearray)
Cette routine permet d'écrire un profil dans un fichier MED.
#define MED_NAME_SIZE
Definition med.h:81
#define MED_SNAME_SIZE
Definition med.h:82
#define MED_QUAD4
Definition med.h:204
#define MED_MAJOR_NUM
Definition med.h:57
#define MED_SEG2
Definition med.h:200
#define MED_NO_LOCALIZATION
Definition med.h:273
#define MED_ALLENTITIES_PROFILE
Definition med.h:286
#define MED_NO_INTERPOLATION
Definition med.h:275
@ MED_CARTESIAN
Definition med.h:258
@ MED_SORT_DTIT
Definition med.h:300
#define MED_MINOR_NUM
Definition med.h:58
#define MED_ALL_CONSTITUENT
Definition med.h:293
@ MED_UNSTRUCTURED_MESH
Definition med.h:131
#define MED_TRIA6
Definition med.h:205
int med_int
Definition med.h:333
#define MED_NO_DT
Definition med.h:311
#define MED_NO_IT
Definition med.h:312
#define MED_NONE
Definition med.h:231
#define MED_RELEASE_NUM
Definition med.h:59
#define MED_NO_PROFILE
Definition med.h:279
@ MED_NODE
Definition med.h:143
@ MED_CELL
Definition med.h:143
@ MED_NODE_ELEMENT
Definition med.h:144
@ MED_DESCENDING_FACE
Definition med.h:143
@ MED_DESCENDING_EDGE
Definition med.h:143
double med_float
Definition med.h:327
herr_t med_err
Definition med.h:323
#define MED_NO_MESH_SUPPORT
Definition med.h:271
hid_t med_idt
Definition med.h:322
#define SSCRUTE(chaine)
Definition med_utils.h:323
#define MESSAGE(chaine)
Definition med_utils.h:324
#define ISCRUTE_int(entier)
Definition med_utils.h:314
#define FTYPECHA
#define ftypecha
#define MODE_ACCES
int main(int argc, char **argv)
#define _b
#define _p2
#define USER_MODE
#define _p1
#define _a
#define itypecha
#define filename
#define USER_INTERLACE
#define ITYPECHA