MED fichier
MEDinterpBaseFunctionRd.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 #include <hdf5.h>
23 #include <string.h>
24 
40 med_err
42  const char* const interpname,
43  const int basisfuncit,
44  med_int* const ncoef,
45  med_int* const power,
46  med_float* const coefficient)
47 {
48  med_idt _interpid=0,_bid=0;
49  med_err _ret=-1;
51  char _basisfuncname[MED_MAX_PARA+1]="";
52  med_int _nvariable =0;
53  med_filter _filter = MED_FILTER_INIT;
54 
55  /*
56  * On inhibe le gestionnaire d'erreur
57  */
59 
60  if (basisfuncit < 1 ) {
62  ISCRUTE_int(basisfuncit);
63  goto ERROR;
64  }
65 
66  NOFINALBLANK(interpname,ERROR);
67  strcat(_path,interpname);
68  /*
69  * Ouverture du DataGroup /INTERP/<interpname>
70  */
71  if ((_interpid = _MEDdatagroupOpen(fid,_path)) < 0) {
72  MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,interpname);
73  SSCRUTE(_path);goto ERROR;
74  }
75 
76  /*
77  * Lecture de l'attribut MED_NOM_NBV
78  */
79  if ( _MEDattrEntierLire(_interpid,MED_NOM_NBV,&_nvariable) < 0 ) {
81  SSCRUTE(MED_NOM_NBV);ISCRUTE(_nvariable);
82  goto ERROR;
83  }
84 
85 
86  sprintf(_basisfuncname,"%0*li",MED_MAX_PARA,(long ) basisfuncit);
87  /*
88  * Ouverture du DataGroup /INTERP/<interpname>/<basisfuncit>
89  */
90  if ((_bid = _MEDdatagroupOpen(_interpid,_basisfuncname)) < 0) {
91  MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,_path);SSCRUTE(_basisfuncname);
92  goto ERROR;
93  }
94 
95  /*
96  * Lecture de l'attribut MED_NOM_MCF
97  */
98  if ( _MEDattrEntierLire(_bid,MED_NOM_MCF,ncoef) < 0 ) {
99  MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,_path);SSCRUTE(_basisfuncname);
100  SSCRUTE(MED_NOM_MCF);ISCRUTE(*ncoef);
101  goto ERROR;
102  }
103 
104  /*TODO : Créer un _MEDdatasetSimpleWr */
105  if ( MEDfilterEntityCr(fid, *ncoef, 1, 1, MED_ALL_CONSTITUENT,
107  MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
109  goto ERROR;
110  }
111 
112  if ( _MEDdatasetRd(_bid,MED_NOM_COE,MED_INTERNAL_FLOAT64,&_filter, (unsigned char * const) coefficient) < 0) {
114  SSCRUTE(_path);SSCRUTE(_basisfuncname);
115  goto ERROR;
116  }
117 
118  if ( MEDfilterClose(&_filter) < 0 ) {
120  SSCRUTE(_path);
121  goto ERROR;
122  }
123 
124 
125  if ( MEDfilterEntityCr(fid, (*ncoef)*_nvariable, 1, 1, MED_ALL_CONSTITUENT,
127  MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
129  goto ERROR;
130  }
131 
132  if ( _MEDdatasetRd(_bid,MED_NOM_POW,MED_INTERNAL_INT,&_filter, (unsigned char * const) power) < 0) {
134  SSCRUTE(_path);SSCRUTE(_basisfuncname);
135  goto ERROR;
136  }
137 
138  if ( MEDfilterClose(&_filter) < 0 ) {
140  SSCRUTE(_path);
141  goto ERROR;
142  }
143 
144 
145  _ret=0;
146  ERROR:
147 
148  if (_interpid>0) if (_MEDdatagroupFermer(_interpid) < 0) {
150  ISCRUTE_id(_interpid);
151  }
152 
153  if (_bid>0) if (_MEDdatagroupFermer(_bid) < 0) {
154  MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_basisfuncname);
155  ISCRUTE_id(_bid);SSCRUTE(_path);
156  }
157 
158  return _ret;
159 }
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_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_NOM_POW
#define MED_NOM_POW
Definition: med_outils.h:156
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
MED_NO_INTERLACE
Definition: med.h:100
med_idt
hid_t med_idt
Definition: med.h:333
_MEDdatasetRd
MEDC_EXPORT med_err _MEDdatasetRd(const med_idt id, const char *const datasetname, const med_internal_type datatype, const med_filter *const filter, unsigned char *const value)
MED_ERR_
#define MED_ERR_(rt, r1, r2, r3)
Definition: med_utils.h:160
_MEDdatagroupOpen
MEDC_EXPORT med_idt _MEDdatagroupOpen(const med_idt pid, const char *const name)
MED_ERR_ATTRIBUTE
#define MED_ERR_ATTRIBUTE
Definition: med_err.h:100
med_err
herr_t med_err
Definition: med.h:334
med_outils.h
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_ERR_OPEN
#define MED_ERR_OPEN
Definition: med_err.h:37
med_int
int med_int
Definition: med.h:344
med_filter
Filtre de sélection.
Definition: med.h:357
ISCRUTE
#define ISCRUTE(entier)
Definition: med_utils.h:313
MED_NOM_NBV
#define MED_NOM_NBV
Definition: med_outils.h:151
MED_ERR_WRITE
#define MED_ERR_WRITE
Definition: med_err.h:29
med_config.h
MED_ERR_READ
#define MED_ERR_READ
Definition: med_err.h:28
med_float
double med_float
Definition: med.h:338
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
MED_INTERPOLATION_GRP
#define MED_INTERPOLATION_GRP
Definition: med_outils.h:179
MED_ERR_DATASET
#define MED_ERR_DATASET
Definition: med_err.h:98
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:323
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:83
_MEDattrEntierLire
#define _MEDattrEntierLire(x, y, z)
Definition: med_hdfi.h:68
MED_ERR_STRUCT_ELEMENT_MSG
#define MED_ERR_STRUCT_ELEMENT_MSG
Definition: med_err.h:136
MED_ALL_CONSTITUENT
#define MED_ALL_CONSTITUENT
Definition: med.h:301
MED_NO_PROFILE
#define MED_NO_PROFILE
Definition: med.h:283
MED_INTERNAL_FLOAT64
Definition: med.h:158
med.h
NOFINALBLANK
#define NOFINALBLANK(chaine, label)
Definition: med_utils.h:311
MED_ERR_RANGE
#define MED_ERR_RANGE
Definition: med_err.h:34
MED_NOM_MCF
#define MED_NOM_MCF
Definition: med_outils.h:153
MED_INTERNAL_INT
Definition: med.h:160
_MEDmodeErreurVerrouiller
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
MED_NOM_COE
#define MED_NOM_COE
Definition: med_outils.h:155
MED_ERR_PARAMETER
#define MED_ERR_PARAMETER
Definition: med_err.h:103
MEDinterpBaseFunctionRd
med_err MEDinterpBaseFunctionRd(const med_idt fid, const char *const interpname, const int basisfuncit, med_int *const ncoef, med_int *const power, med_float *const coefficient)
Cette routine permet la lecture d'une fonction de base/forme de l'interpolation interpname.
Definition: MEDinterpBaseFunctionRd.c:41
MED_INTERPOLATION_GRP_SIZE
#define MED_INTERPOLATION_GRP_SIZE
Definition: med_outils.h:180
_MEDdatagroupFermer
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)