MED fichier
MEDsubdomainCorrespondenceWr.c
Aller à la documentation de ce fichier.
1 /* This file is part of MED.
2  *
3  * COPYRIGHT (C) 1999 - 2023 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 
18 
19 #include <med.h>
20 #include <med_config.h>
21 #include <med_outils.h>
22 
23 #include <string.h>
24 #include <stdlib.h>
25 
46 med_err
48  const char * const localmeshname,
49  const char * const jointname,
50  const med_int numdt,
51  const med_int numit,
52  const med_entity_type localentitype,
53  const med_geometry_type localgeotype,
54  const med_entity_type remoteentitype,
55  const med_geometry_type remotegeotype,
56  const med_int nentity,
57  const med_int * const correspondence) {
58 
59  med_access_mode _MED_ACCESS_MODE;
60  med_idt _root=0,_jntid=0,_meshid=0,_cstpid=0,_datagroup2=0;
61  med_idt _dataset=0;
62  med_err _ret=-1;
64  char _computationstepname[2*MED_MAX_PARA+1]="";
65  char _datagroupname2[4*MED_TAILLE_NOM_ENTITE+4]="";
66  char _localgeotypename [MED_TAILLE_NOM_ENTITE+1]="";
67  char _remotegeotypename [MED_TAILLE_NOM_ENTITE+1]="";
68  med_sorting_type _sortingtype=0;
69  med_filter _filter = MED_FILTER_INIT;
70  med_int _localgeotype = MED_NONE;
71  med_int _localentitype = localentitype;
72  med_int _remotegeotype = MED_NONE;
73  med_int _remoteentitype = remoteentitype;
74 
75  if ( localentitype == MED_NODE_ELEMENT ) _localgeotype=MED_NODE ;
76 
77 
78  if ( localentitype != MED_NODE ) _localgeotype = localgeotype ;
79  if ( remoteentitype != MED_NODE ) _remotegeotype = remotegeotype ;
80 
81  /*
82  * On inhibe le gestionnaire d'erreur
83  */
85  if (_MEDcheckVersion30(fid) < 0) goto ERROR;
86 
87  if ( (_MED_ACCESS_MODE = _MEDmodeAcces(fid) ) == MED_ACC_UNDEF ) {
89  goto ERROR;
90  }
91 
92  if ( _MED_ACCESS_MODE == MED_ACC_RDONLY) {
94  ISCRUTE_int(_MED_ACCESS_MODE);
95  goto ERROR;
96  }
97 
98  /*
99  * Ouverture du dataGroup /JNT/
100  */
101  if ((_root = _MEDdatagroupOuvrir(fid,_path)) < 0) {
103  goto ERROR;
104  }
105 
106  /*
107  * Ouverture du dataGroup <meshname>
108  */
109  if ((_meshid = _MEDdatagroupOuvrir(_root,localmeshname)) < 0) {
110  MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,localmeshname);
111  SSCRUTE(_path);goto ERROR;
112  }
113 
114  strcat(_path,localmeshname);
115 
116  /*
117  * Ouverture du data Group "/JNT/<localmeshname>/<jointname>"
118  */
119  if ((_jntid = _MEDdatagroupOuvrir(_meshid,jointname)) < 0) {
120  MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,jointname);
121  SSCRUTE(_path);goto ERROR;
122  }
123 
124  strcat(_path,"/");
125  strcat(_path,jointname);
126 
127  if ( MEDmeshSortingTypeRd(fid,localmeshname,&_sortingtype) < 0 ) {
128  MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API," MEDmeshSortingTypeRd");
129  SSCRUTE(localmeshname);ISCRUTE_int(_sortingtype);goto ERROR;
130  }
131 
132  _MEDgetComputationStepName(_sortingtype,numdt,numit,_computationstepname);
133 
134  if ((_cstpid = _MEDdatagroupOuvrir(_jntid,_computationstepname)) < 0)
135  if ((_cstpid = _MEDdatagroupCreer(_jntid,_computationstepname)) < 0 ) {
136  MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_computationstepname);
137  SSCRUTE(_path);goto ERROR;
138  }
139 
140  /*Cree ou ouvre l'attribut MED_NOM_NDT pour �criture */
141  if ( _MEDattributeIntWr(_cstpid,MED_NOM_NDT,&numdt) < 0) {
143  SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);
144  goto ERROR;
145  }
146 
147 
148  /*Cree ou ouvre l'attribut MED_NOM_NOR pour �criture */
149  if ( _MEDattributeIntWr(_cstpid,MED_NOM_NOR,&numit) < 0) {
151  SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);
152  goto ERROR;
153  }
154 
155 
156  /*
157  * Creation/Ouverture du datagroup de niveau
158  * <localentitype>[.<localgeotype>].<remoteentitype>[.<remotegeotype>]
159  */
160  if (_MEDgetEntityTypeName(_datagroupname2,localentitype) < 0) {
162  ISCRUTE_int(localentitype);SSCRUTE(jointname);goto ERROR;
163  }
164  if ( localentitype != MED_NODE ) {
165  if ( _MEDgetInternalGeometryTypeName(fid,_localgeotypename,localgeotype) < 0) {
167  ISCRUTE_int(localgeotype);SSCRUTE(jointname);goto ERROR;
168  }
169  strcat(_datagroupname2,".");
170  strcat(_datagroupname2,_localgeotypename);
171  }
172 
173  strcat(_datagroupname2,".");
174 
175  if (_MEDgetEntityTypeName(&_datagroupname2[strlen(_datagroupname2)],remoteentitype) < 0) {
177  ISCRUTE_int(remoteentitype);SSCRUTE(jointname);goto ERROR;
178  }
179  if ( remoteentitype != MED_NODE ) {
180  if ( _MEDgetInternalGeometryTypeName(fid,_remotegeotypename,remotegeotype) < 0) {
182  ISCRUTE_int(remotegeotype);SSCRUTE(jointname);goto ERROR;
183  }
184  strcat(_datagroupname2,".");
185  strcat(_datagroupname2,_remotegeotypename);
186  }
187 
188 
189  if ( (_datagroup2 = _MEDdatagroupOuvrir(_cstpid,_datagroupname2)) < 0)
190  if ((_datagroup2 = _MEDdatagroupCreer(_cstpid,_datagroupname2)) < 0) {
191  MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_datagroupname2);
192  SSCRUTE(_path);SSCRUTE(jointname);goto ERROR;
193  }
194 
195  /*Cree ou ouvre l'attribut MED_NOM_ENT pour �criture */
196  if (_MEDattributeIntWr(_datagroup2,MED_NOM_ENT,&_localentitype) < 0) {
198  ISCRUTE(_localentitype);goto ERROR;
199  }
200 
201  /*Cree ou ouvre l'attribut MED_NOM_GEO pour �criture */
202  if (_MEDattributeIntWr(_datagroup2,MED_NOM_GEO,&_localgeotype) < 0) {
204  ISCRUTE(_localgeotype);goto ERROR;
205  }
206 
207  /*Cree ou ouvre l'attribut MED_NOM_END pour �criture */
208  if (_MEDattributeIntWr(_datagroup2,MED_NOM_END,&_remoteentitype) < 0) {
210  ISCRUTE(_remoteentitype);goto ERROR;
211  }
212 
213  /*Cree ou ouvre l'attribut MED_NOM_GED pour �criture */
214  if (_MEDattributeIntWr(_datagroup2,MED_NOM_GED,&_remotegeotype) < 0) {
216  ISCRUTE(_remotegeotype);goto ERROR;
217  }
218 
219 
220  if ( MEDfilterEntityCr(fid, nentity, 1, 2, MED_ALL_CONSTITUENT,
222  MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
224  goto ERROR;
225  }
226 
227  if ( _MEDdatasetWr(_datagroup2,MED_NOM_COR,MED_INTERNAL_INT,&_filter, correspondence) < 0) {
229  SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);
230  goto ERROR;
231  }
232 
233  if ( MEDfilterClose(&_filter) < 0 ) {
235  SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);
236  goto ERROR;
237  }
238 
239  if ((_dataset = _MEDdatasetOuvrir(_datagroup2,MED_NOM_COR)) < 0) {
241  SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);SSCRUTE(_datagroupname2);
242  goto ERROR;
243  }
244 
245  if ( _MEDattributeIntWr(_dataset,MED_NOM_NBR,&nentity) < 0) {
247  SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);SSCRUTE(_datagroupname2);
248  SSCRUTE(MED_NOM_NBR);ISCRUTE(nentity);goto ERROR;
249  }
250 
251 
252 
253  _ret=0;
254  ERROR:
255 
256  if (_dataset>0) if (_MEDdatasetFermer(_dataset) < 0) {
258  ISCRUTE_id(_dataset);
259  }
260 
261  if (_datagroup2>0) if (_MEDdatagroupFermer(_datagroup2) < 0) {
262  MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_datagroupname2);
263  ISCRUTE_id(_datagroup2);SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);
264  }
265 
266  if (_cstpid>0) if (_MEDdatagroupFermer(_cstpid) < 0) {
267  MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_computationstepname);
268  ISCRUTE_id(_cstpid);SSCRUTE(_path);ISCRUTE(numdt);ISCRUTE(numit);
269  }
270 
271  if (_jntid>0) if (_MEDdatagroupFermer(_jntid) < 0) {
272  MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,jointname);
273  ISCRUTE_id(_jntid);SSCRUTE(_path);
274  }
275 
276  if (_meshid>0) if (_MEDdatagroupFermer(_meshid) < 0) {
277  MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,localmeshname);
278  ISCRUTE_id(_meshid);
279  }
280 
281  if (_root>0) if (_MEDdatagroupFermer(_root) < 0) {
283  ISCRUTE_id(_root);
284  }
285 
286  return _ret;
287 }
MEDfilterClose
MEDC_EXPORT med_err MEDfilterClose(med_filter *const filter)
Désalloue les ressources hdf détenues par un filtre.
Definition: MEDfilterClose.c:35
MED_ACC_RDONLY
Definition: med.h:122
MED_ERR_API
#define MED_ERR_API
Definition: med_err.h:111
MED_ERR_CREATE
#define MED_ERR_CREATE
Definition: med_err.h:27
ISCRUTE_id
#define ISCRUTE_id(entier)
Definition: med_utils.h:319
MED_UNDEF_STMODE
Definition: med.h:110
MED_FILTER_INIT
#define MED_FILTER_INIT
Definition: med.h:376
MED_UNDEF_SIZE
#define MED_UNDEF_SIZE
Definition: med.h:308
med_geometry_type
int med_geometry_type
Definition: med.h:196
_MEDattributeIntWr
#define _MEDattributeIntWr(w, x, y)
Definition: med_hdfi.h:45
MEDfilterEntityCr
MEDC_EXPORT med_err MEDfilterEntityCr(const med_idt fid, const med_int nentity, const med_int nvaluesperentity, const med_int nconstituentpervalue, const med_int constituentselect, const med_switch_mode switchmode, const med_storage_mode storagemode, const char *const profilename, const med_int filterarraysize, const med_int *const filterarray, med_filter *const filter)
Definition: MEDfilterEntityCr.c:55
ISCRUTE_int
#define ISCRUTE_int(entier)
Definition: med_utils.h:314
_MEDmodeAcces
MEDC_EXPORT med_access_mode _MEDmodeAcces(med_idt oid)
MED_NO_INTERLACE
Definition: med.h:100
med_idt
hid_t med_idt
Definition: med.h:333
MED_ERR_
#define MED_ERR_(rt, r1, r2, r3)
Definition: med_utils.h:160
MED_ERR_SUBDOMAINJOINT_MSG
#define MED_ERR_SUBDOMAINJOINT_MSG
Definition: med_err.h:172
MED_ERR_ATTRIBUTE
#define MED_ERR_ATTRIBUTE
Definition: med_err.h:100
_MEDgetEntityTypeName
MEDC_EXPORT med_err _MEDgetEntityTypeName(char *const entitytypename, const med_entity_type entitytype)
MED_ERR_INVALID
#define MED_ERR_INVALID
Definition: med_err.h:40
_MEDdatasetWr
MEDC_EXPORT med_err _MEDdatasetWr(const med_idt id, const char *const datasetname, const med_internal_type datatype, const med_filter *const filter, const void *const value)
MED_NOM_END
#define MED_NOM_END
Definition: med_outils.h:137
med_err
herr_t med_err
Definition: med.h:334
med_sorting_type
med_sorting_type
Definition: med.h:311
MEDsubdomainCorrespondenceWr
med_err MEDsubdomainCorrespondenceWr(const med_idt fid, const char *const localmeshname, const char *const jointname, const med_int numdt, const med_int numit, const med_entity_type localentitype, const med_geometry_type localgeotype, const med_entity_type remoteentitype, const med_geometry_type remotegeotype, const med_int nentity, const med_int *const correspondence)
Cette routine permet l'écriture d'une correspondance dans un joint pour un type de couple d'entité en...
Definition: MEDsubdomainCorrespondenceWr.c:47
MEDmeshSortingTypeRd
MEDC_EXPORT med_err MEDmeshSortingTypeRd(const med_idt fid, const char *const meshname, med_sorting_type *const sortingtype)
Cette routine lit l'ordre de tri des étapes évolutives du maillage.
Definition: MEDmeshSortingTypeRd.c:35
med_outils.h
med_entity_type
med_entity_type
Definition: med.h:145
MED_MAX_PARA
#define MED_MAX_PARA
Definition: med.h:78
MED_ERR_CLOSE
#define MED_ERR_CLOSE
Definition: med_err.h:30
MED_ERR_FILTER
#define MED_ERR_FILTER
Definition: med_err.h:93
MED_NOM_GEO
#define MED_NOM_GEO
Definition: med_outils.h:133
_MEDgetInternalGeometryTypeName
MEDC_EXPORT med_err _MEDgetInternalGeometryTypeName(const med_idt fid, char *const geotypename, med_geometry_type geotype)
MED_NOM_COR
#define MED_NOM_COR
Definition: med_outils.h:55
MED_ERR_OPEN
#define MED_ERR_OPEN
Definition: med_err.h:37
med_int
int med_int
Definition: med.h:344
MED_ERR_VALUE_MSG
#define MED_ERR_VALUE_MSG
Definition: med_err.h:184
_MEDgetComputationStepName
MEDC_EXPORT med_err _MEDgetComputationStepName(const med_sorting_type sortingtype, const med_int numdt, const med_int numit, char *const datagroupname)
med_filter
Filtre de sélection.
Definition: med.h:357
ISCRUTE
#define ISCRUTE(entier)
Definition: med_utils.h:313
MED_ERR_WRITE
#define MED_ERR_WRITE
Definition: med_err.h:29
med_config.h
_MEDdatasetFermer
MEDC_EXPORT med_err _MEDdatasetFermer(med_idt id)
MED_ERR_CALL
#define MED_ERR_CALL
Definition: med_err.h:48
MED_ERR_ENTITY
#define MED_ERR_ENTITY
Definition: med_err.h:85
MED_JOINT_GRP
#define MED_JOINT_GRP
Definition: med_outils.h:252
MED_ACC_UNDEF
Definition: med.h:126
MED_NOM_NDT
#define MED_NOM_NDT
Definition: med_outils.h:142
med_access_mode
med_access_mode
Definition: med.h:122
MED_NONE
#define MED_NONE
Definition: med.h:233
MED_TAILLE_NOM_ENTITE
#define MED_TAILLE_NOM_ENTITE
Definition: med_outils.h:38
MED_ERR_INTERNAL_MSG
#define MED_ERR_INTERNAL_MSG
Definition: med_err.h:196
MED_ERR_DATAGROUP
#define MED_ERR_DATAGROUP
Definition: med_err.h:99
_MEDdatagroupCreer
MEDC_EXPORT med_idt _MEDdatagroupCreer(med_idt pid, const char *const nom)
MED_ERR_DATASET
#define MED_ERR_DATASET
Definition: med_err.h:98
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:323
MED_ERR_ACCESSMODE
#define MED_ERR_ACCESSMODE
Definition: med_err.h:104
MED_NOM_NOR
#define MED_NOM_NOR
Definition: med_outils.h:145
_MEDdatasetOuvrir
MEDC_EXPORT med_idt _MEDdatasetOuvrir(med_idt pid, char *nom)
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:83
_MEDcheckVersion30
MEDC_EXPORT med_err _MEDcheckVersion30(med_idt fid)
MED_NODE_ELEMENT
Definition: med.h:146
MED_ALL_CONSTITUENT
#define MED_ALL_CONSTITUENT
Definition: med.h:301
MED_NODE
Definition: med.h:145
MED_NO_PROFILE
#define MED_NO_PROFILE
Definition: med.h:283
MED_NOM_NBR
#define MED_NOM_NBR
Definition: med_outils.h:40
med.h
MED_ERR_UNRECOGNIZED
#define MED_ERR_UNRECOGNIZED
Definition: med_err.h:39
_MEDdatagroupOuvrir
MEDC_EXPORT med_idt _MEDdatagroupOuvrir(med_idt pid, const char *const nom)
MED_INTERNAL_INT
Definition: med.h:160
MED_NOM_ENT
#define MED_NOM_ENT
Definition: med_outils.h:135
_MEDmodeErreurVerrouiller
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
MED_JOINT_GRP_SIZE
#define MED_JOINT_GRP_SIZE
Definition: med_outils.h:253
MED_ERR_GEOMETRY
#define MED_ERR_GEOMETRY
Definition: med_err.h:107
MED_ERR_FILE_MSG
#define MED_ERR_FILE_MSG
Definition: med_err.h:132
MED_NOM_GED
#define MED_NOM_GED
Definition: med_outils.h:136
_MEDdatagroupFermer
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)