MED fichier

MEDinterpBaseFunctionWr.c
Aller à la documentation de ce fichier.
1 /* This file is part of MED.
2 *
3 * COPYRIGHT (C) 1999 - 2021 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
41 MEDinterpBaseFunctionWr( const med_idt fid,
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;
51 char _path[MED_INTERPOLATION_GRP_SIZE+MED_NAME_SIZE+1]=MED_INTERPOLATION_GRP;
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 */
59 _MEDmodeErreurVerrouiller();
60 if (_MEDcheckVersion30(fid) < 0) goto ERROR;
61
62 if ( (_MED_ACCESS_MODE = _MEDmodeAcces(fid) ) == MED_ACC_UNDEF ) {
63 MED_ERR_(_ret,MED_ERR_UNRECOGNIZED,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
64 goto ERROR;
65 }
66
67 if ( _MED_ACCESS_MODE == MED_ACC_RDONLY) {
68 MED_ERR_(_ret,MED_ERR_INVALID,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
69 ISCRUTE_int(_MED_ACCESS_MODE);
70 goto ERROR;
71 }
72
73 if (basisfuncit < 1 ) {
74 MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_PARAMETER,"");
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 ) {
93 MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,_path);
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,
120 MED_NO_INTERLACE,MED_UNDEF_STMODE,
121 MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
122 MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_FILTER,MED_ERR_INTERNAL_MSG);
123 goto ERROR;
124 }
125
126 if ( _MEDdatasetWr(_bid,MED_NOM_COE,MED_INTERNAL_FLOAT64,&_filter, coefficient) < 0) {
127 MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_DATASET,MED_NOM_COE);
128 SSCRUTE(_path);SSCRUTE(_basisfuncname);
129 goto ERROR;
130 }
131
132 if ( MEDfilterClose(&_filter) < 0 ) {
133 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_FILTER,MED_ERR_STRUCT_ELEMENT_MSG);
134 SSCRUTE(_path);
135 goto ERROR;
136 }
137
138
139 if ( MEDfilterEntityCr(fid, ncoef*_nvariable, 1, 1, MED_ALL_CONSTITUENT,
140 MED_NO_INTERLACE,MED_UNDEF_STMODE,
141 MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
142 MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_FILTER,MED_ERR_INTERNAL_MSG);
143 goto ERROR;
144 }
145
146 if ( _MEDdatasetWr(_bid,MED_NOM_POW,MED_INTERNAL_INT,&_filter, power) < 0) {
147 MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_DATASET,MED_NOM_POW);
148 SSCRUTE(_path);SSCRUTE(_basisfuncname);
149 goto ERROR;
150 }
151
152 if ( MEDfilterClose(&_filter) < 0 ) {
153 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_FILTER,MED_ERR_STRUCT_ELEMENT_MSG);
154 SSCRUTE(_path);
155 goto ERROR;
156 }
157
158
159 _ret=0;
160 ERROR:
161
162 if (_interpid>0) if (_MEDdatagroupFermer(_interpid) < 0) {
163 MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_path);
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 }
MED_ERR_CREATE
#define MED_ERR_CREATE
Definition: med_err.h:27
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
_MEDdatagroupFermer
MEDC_EXPORT med_err _MEDdatagroupFermer(med_idt id)
med_filter
Filtre de sélection.
Definition: med.h:355
MED_FILTER_INIT
#define MED_FILTER_INIT
Definition: med.h:374
ISCRUTE_id
#define ISCRUTE_id(entier)
Definition: med_utils.h:319
MED_UNDEF_SIZE
#define MED_UNDEF_SIZE
Definition: med.h:306
_MEDmodeAcces
MEDC_EXPORT med_access_mode _MEDmodeAcces(med_idt oid)
ISCRUTE_int
#define ISCRUTE_int(entier)
Definition: med_utils.h:314
MEDinterpBaseFunctionWr
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.
Definition: MEDinterpBaseFunctionWr.c:41
med_idt
hid_t med_idt
Definition: med.h:331
_MEDdatagroupOpen
MEDC_EXPORT med_idt _MEDdatagroupOpen(const med_idt pid, const char *const name)
MED_ERR_
#define MED_ERR_(rt, r1, r2, r3)
Definition: med_utils.h:160
MED_NOM_POW
#define MED_NOM_POW
Definition: med_outils.h:156
_MEDattributeIntWr
#define _MEDattributeIntWr(w, x, y)
Definition: med_hdfi.h:45
MED_ACC_RDONLY
Definition: med.h:120
MED_ERR_INVALID
#define MED_ERR_INVALID
Definition: med_err.h:40
MED_UNDEF_STMODE
Definition: med.h:108
MED_ERR_ATTRIBUTE
#define MED_ERR_ATTRIBUTE
Definition: med_err.h:100
med_err
herr_t med_err
Definition: med.h:332
_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_MAX_PARA
#define MED_MAX_PARA
Definition: med.h:76
MED_ERR_CLOSE
#define MED_ERR_CLOSE
Definition: med_err.h:30
MED_INTERNAL_FLOAT64
Definition: med.h:156
MED_ERR_FILTER
#define MED_ERR_FILTER
Definition: med_err.h:93
med_int
int med_int
Definition: med.h:342
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.h
ISCRUTE
#define ISCRUTE(entier)
Definition: med_utils.h:313
med_float
double med_float
Definition: med.h:336
MED_ERR_READ
#define MED_ERR_READ
Definition: med_err.h:28
med_access_mode
med_access_mode
Definition: med.h:120
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_INTERNAL_INT
Definition: med.h:158
MED_ERR_INTERNAL_MSG
#define MED_ERR_INTERNAL_MSG
Definition: med_err.h:196
MED_ERR_ACCESSMODE
#define MED_ERR_ACCESSMODE
Definition: med_err.h:104
_MEDdatagroupCreer
MEDC_EXPORT med_idt _MEDdatagroupCreer(med_idt pid, const char *const nom)
MED_INTERPOLATION_GRP
#define MED_INTERPOLATION_GRP
Definition: med_outils.h:179
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:323
MED_ERR_DATASET
#define MED_ERR_DATASET
Definition: med_err.h:98
MED_ERR_DATAGROUP
#define MED_ERR_DATAGROUP
Definition: med_err.h:99
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:81
_MEDcheckVersion30
MEDC_EXPORT med_err _MEDcheckVersion30(med_idt fid)
med_outils.h
med_config.h
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:299
MED_NO_PROFILE
#define MED_NO_PROFILE
Definition: med.h:281
MED_NO_INTERLACE
Definition: med.h:98
_MEDattrEntierLire
#define _MEDattrEntierLire(x, y, z)
Definition: med_hdfi.h:68
MED_ERR_RANGE
#define MED_ERR_RANGE
Definition: med_err.h:34
NOFINALBLANK
#define NOFINALBLANK(chaine, label)
Definition: med_utils.h:311
MED_ERR_UNRECOGNIZED
#define MED_ERR_UNRECOGNIZED
Definition: med_err.h:39
MED_ERR_PARAMETER
#define MED_ERR_PARAMETER
Definition: med_err.h:103
MED_NOM_MCF
#define MED_NOM_MCF
Definition: med_outils.h:153
MED_INTERPOLATION_GRP_SIZE
#define MED_INTERPOLATION_GRP_SIZE
Definition: med_outils.h:180
MED_ERR_FILE_MSG
#define MED_ERR_FILE_MSG
Definition: med_err.h:132
MED_NOM_COE
#define MED_NOM_COE
Definition: med_outils.h:155
_MEDmodeErreurVerrouiller
MEDC_EXPORT void _MEDmodeErreurVerrouiller(void)
MED_ACC_UNDEF
Definition: med.h:124