MED fichier
MEDinterpBaseFunctionRd.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
18
19#include <med.h>
20#include <med_config.h>
21#include <med_outils.h>
22#include <hdf5.h>
23#include <string.h>
24
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;
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) {
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}
MEDC_EXPORT med_err MEDfilterClose(med_filter *const filter)
Désalloue les ressources hdf détenues par un filtre.
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)
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.
#define MED_NAME_SIZE
Definition med.h:81
@ MED_NO_INTERLACE
Definition med.h:98
#define MED_FILTER_INIT
Definition med.h:365
@ MED_UNDEF_STMODE
Definition med.h:108
#define MED_MAX_PARA
Definition med.h:76
#define MED_ALL_CONSTITUENT
Definition med.h:293
int med_int
Definition med.h:333
#define MED_UNDEF_SIZE
Definition med.h:297
@ MED_INTERNAL_INT
Definition med.h:158
@ MED_INTERNAL_FLOAT64
Definition med.h:156
#define MED_NO_PROFILE
Definition med.h:279
double med_float
Definition med.h:327
herr_t med_err
Definition med.h:323
hid_t med_idt
Definition med.h:322
#define MED_ERR_OPEN
Definition med_err.h:37
#define MED_ERR_FILTER
Definition med_err.h:93
#define MED_ERR_DATAGROUP
Definition med_err.h:99
#define MED_ERR_CLOSE
Definition med_err.h:30
#define MED_ERR_STRUCT_ELEMENT_MSG
Definition med_err.h:136
#define MED_ERR_ATTRIBUTE
Definition med_err.h:100
#define MED_ERR_CREATE
Definition med_err.h:27
#define MED_ERR_INTERNAL_MSG
Definition med_err.h:196
#define MED_ERR_PARAMETER
Definition med_err.h:103
#define MED_ERR_READ
Definition med_err.h:28
#define MED_ERR_RANGE
Definition med_err.h:34
#define MED_ERR_WRITE
Definition med_err.h:29
#define MED_ERR_DATASET
Definition med_err.h:98
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)
#define _MEDattrEntierLire(x, y, z)
Definition med_hdfi.h:68
MEDC_EXPORT med_idt _MEDdatagroupOpen(const med_idt pid, const char *const name)
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)
#define MED_NOM_COE
Definition med_outils.h:155
#define MED_NOM_NBV
Definition med_outils.h:151
#define MED_INTERPOLATION_GRP_SIZE
Definition med_outils.h:180
#define MED_NOM_MCF
Definition med_outils.h:153
#define MED_INTERPOLATION_GRP
Definition med_outils.h:179
#define MED_NOM_POW
Definition med_outils.h:156
#define MED_ERR_(rt, r1, r2, r3)
Definition med_utils.h:160
#define SSCRUTE(chaine)
Definition med_utils.h:323
#define ISCRUTE(entier)
Definition med_utils.h:313
#define ISCRUTE_id(entier)
Definition med_utils.h:319
#define NOFINALBLANK(chaine, label)
Definition med_utils.h:311
#define ISCRUTE_int(entier)
Definition med_utils.h:314
Filtre de sélection.
Definition med.h:346