MED fichier
MEDinterpBaseFunctionWr.c
Aller à la documentation de ce fichier.
1 /* This file is part of MED.
2  *
3  * COPYRIGHT (C) 1999 - 2019 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 #include <hdf5.h>
23 #include <string.h>
24 
40 med_err
42  const char* const interpname,
43  const med_int basisfuncit,
44  const med_int ncoef,
45  const med_int* const power,
46  const med_float* const coefficient)
47 {
48  med_access_mode _MED_ACCESS_MODE;
49  med_idt _root=0,_interpid=0,_bid=0;
50  med_int _ret=-1;
52  char _basisfuncname[MED_MAX_PARA+1]="";
53  med_int _nvariable =0;
54  med_filter _filter = MED_FILTER_INIT;
55 
56  /*
57  * On inhibe le gestionnaire d'erreur
58  */
60  if (_MEDcheckVersion30(fid) < 0) goto ERROR;
61 
62  if ( (_MED_ACCESS_MODE = _MEDmodeAcces(fid) ) == MED_ACC_UNDEF ) {
64  goto ERROR;
65  }
66 
67  if ( _MED_ACCESS_MODE == MED_ACC_RDONLY) {
69  ISCRUTE_int(_MED_ACCESS_MODE);
70  goto ERROR;
71  }
72 
73  if (basisfuncit < 1 ) {
75  ISCRUTE(basisfuncit);
76  goto ERROR;
77  }
78 
79  NOFINALBLANK(interpname,ERROR);
80  strcat(_path,interpname);
81  /*
82  * Ouverture du DataGroup /INTERP/<interpname>
83  */
84  if ((_interpid = _MEDdatagroupOpen(fid,_path)) < 0) {
85  MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,interpname);
86  SSCRUTE(_path);goto ERROR;
87  }
88 
89  /*
90  * Lecture de l'attribut MED_NOM_NBV
91  */
92  if ( _MEDattrEntierLire(_interpid,MED_NOM_NBV,&_nvariable) < 0 ) {
94  SSCRUTE(MED_NOM_NBV);ISCRUTE(_nvariable);
95  goto ERROR;
96  }
97 
98 
99  sprintf(_basisfuncname,"%0*li",MED_MAX_PARA,(long ) basisfuncit);
100  /*
101  * Si le DataGroup /INTERP/<interpname>/<basisfuncit> n'existe pas, on le cree
102  */
103  if ((_bid = _MEDdatagroupOpen(_interpid,_basisfuncname)) < 0)
104  if ((_bid = _MEDdatagroupCreer(_interpid,_basisfuncname)) < 0) {
105  MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_path);SSCRUTE(_basisfuncname);
106  goto ERROR;
107  }
108 
109  /*
110  * Creation/Ecriture de l'attribut MED_NOM_MCF
111  */
112  if ( _MEDattributeIntWr(_bid,MED_NOM_MCF,&ncoef) < 0 ) {
113  MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,_path);SSCRUTE(_basisfuncname);
114  SSCRUTE(MED_NOM_MCF);ISCRUTE(ncoef);
115  goto ERROR;
116  }
117 
118  /*TODO : Créer un _MEDdatasetSimpleWr */
119  if ( MEDfilterEntityCr(fid, ncoef, 1, 1, MED_ALL_CONSTITUENT,
121  MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
123  goto ERROR;
124  }
125 
126  if ( _MEDdatasetWr(_bid,MED_NOM_COE,MED_INTERNAL_FLOAT64,&_filter, coefficient) < 0) {
128  SSCRUTE(_path);SSCRUTE(_basisfuncname);
129  goto ERROR;
130  }
131 
132  if ( MEDfilterClose(&_filter) < 0 ) {
134  SSCRUTE(_path);
135  goto ERROR;
136  }
137 
138 
139  if ( MEDfilterEntityCr(fid, ncoef*_nvariable, 1, 1, MED_ALL_CONSTITUENT,
141  MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
143  goto ERROR;
144  }
145 
146  if ( _MEDdatasetWr(_bid,MED_NOM_POW,MED_INTERNAL_INT,&_filter, power) < 0) {
148  SSCRUTE(_path);SSCRUTE(_basisfuncname);
149  goto ERROR;
150  }
151 
152  if ( MEDfilterClose(&_filter) < 0 ) {
154  SSCRUTE(_path);
155  goto ERROR;
156  }
157 
158 
159  _ret=0;
160  ERROR:
161 
162  if (_interpid>0) if (_MEDdatagroupFermer(_interpid) < 0) {
164  ISCRUTE_id(_interpid);
165  }
166 
167  if (_bid>0) if (_MEDdatagroupFermer(_bid) < 0) {
168  MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_basisfuncname);
169  ISCRUTE_id(_bid);SSCRUTE(_path);
170  }
171 
172  return _ret;
173 }
Filtre de sélection.
Definition: med.h:350
#define MED_ERR_ATTRIBUTE
Definition: med_err.h:99
herr_t med_err
Definition: med.h:327
#define NOFINALBLANK(chaine, label)
Definition: med_utils.h:304
#define MED_NOM_POW
Definition: med_outils.h:150
#define MED_ERR_(rt, r1, r2, r3)
Definition: med_utils.h:158
#define MED_ERR_FILE_MSG
Definition: med_err.h:129
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)
#define MED_ERR_INVALID
Definition: med_err.h:37
#define MED_ERR_READ
Definition: med_err.h:25
MEDC_EXPORT med_err _MEDcheckVersion30(med_idt fid)
MEDC_EXPORT med_access_mode _MEDmodeAcces(med_idt oid)
#define MED_INTERPOLATION_GRP
Definition: med_outils.h:176
med_access_mode
Definition: med.h:118
#define MED_ERR_PARAMETER
Definition: med_err.h:102
#define MED_NOM_COE
Definition: med_outils.h:149
double med_float
Definition: med.h:331
#define MED_FILTER_INIT
Definition: med.h:369
#define MED_ALL_CONSTITUENT
Definition: med.h:296
#define ISCRUTE_int(entier)
Definition: med_utils.h:307
#define MED_ERR_ACCESSMODE
Definition: med_err.h:103
#define MED_ERR_UNRECOGNIZED
Definition: med_err.h:36
#define MED_ERR_WRITE
Definition: med_err.h:26
#define MED_MAX_PARA
Definition: med.h:74
int med_int
Definition: med.h:337
#define MED_ERR_CREATE
Definition: med_err.h:24
#define ISCRUTE(entier)
Definition: med_utils.h:306
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)
#define MED_ERR_CLOSE
Definition: med_err.h:27
#define MED_NOM_NBV
Definition: med_outils.h:145
#define _MEDattributeIntWr(w, x, y)
Definition: med_hdfi.h:45
med_err MEDinterpBaseFunctionWr(const med_idt fid, const char *const interpname, const med_int basisfuncit, const med_int ncoef, const med_int *const power, const med_float *const coefficient)
Cette routine permet l'écriture d'une fonction de base/forme de l'interpolation interpname.
#define MED_INTERPOLATION_GRP_SIZE
Definition: med_outils.h:177
#define MED_ERR_RANGE
Definition: med_err.h:31
#define MED_NO_PROFILE
Definition: med.h:278
MEDC_EXPORT med_err MEDfilterClose(med_filter *const filter)
Désalloue les ressources hdf détenues par un filtre.
#define MED_ERR_INTERNAL_MSG
Definition: med_err.h:191
#define MED_ERR_DATAGROUP
Definition: med_err.h:98
#define _MEDattrEntierLire(x, y, z)
Definition: med_hdfi.h:64
hid_t med_idt
Definition: med.h:326
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)
#define MED_ERR_STRUCT_ELEMENT_MSG
Definition: med_err.h:133
#define MED_NAME_SIZE
Definition: med.h:79
#define MED_NOM_MCF
Definition: med_outils.h:147
#define ISCRUTE_id(entier)
Definition: med_utils.h:311
MEDC_EXPORT med_idt _MEDdatagroupCreer(med_idt pid, const char *const nom)
#define MED_UNDEF_SIZE
Definition: med.h:303
#define MED_ERR_DATASET
Definition: med_err.h:97
#define SSCRUTE(chaine)
Definition: med_utils.h:315
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
#define MED_ERR_FILTER
Definition: med_err.h:92
MEDC_EXPORT med_idt _MEDdatagroupOpen(const med_idt pid, const char *const name)