MEDinterpBaseFunctionWr.c

Aller à la documentation de ce fichier.
00001 /*  This file is part of MED.
00002  *
00003  *  COPYRIGHT (C) 1999 - 2015  EDF R&D, CEA/DEN
00004  *  MED is free software: you can redistribute it and/or modify
00005  *  it under the terms of the GNU Lesser General Public License as published by
00006  *  the Free Software Foundation, either version 3 of the License, or
00007  *  (at your option) any later version.
00008  *
00009  *  MED is distributed in the hope that it will be useful,
00010  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
00011  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00012  *  GNU Lesser General Public License for more details.
00013  *
00014  *  You should have received a copy of the GNU Lesser General Public License
00015  *  along with MED.  If not, see <http://www.gnu.org/licenses/>.
00016  */
00017 
00018 
00019 #include <med.h>
00020 #include <med_config.h>
00021 #include <med_outils.h>
00022 #include <hdf5.h>
00023 #include <string.h>
00024 
00040 med_err
00041 MEDinterpBaseFunctionWr( const med_idt          fid,
00042                          const char*      const interpname,
00043                          const med_int          basisfuncit,
00044                          const med_int          ncoef,
00045                          const med_int*   const power,
00046                          const med_float* const coefficient)
00047 {
00048  med_access_mode   _MED_ACCESS_MODE;
00049  med_idt           _root=0,_interpid=0,_bid=0;
00050   med_int           _ret=-1;
00051   char              _path[MED_TAILLE_INTERP+MED_NAME_SIZE+1]=MED_INTERP;
00052   char              _basisfuncname[MED_MAX_PARA+1]="";
00053   med_int           _nvariable  =0;
00054   med_filter        _filter        = MED_FILTER_INIT;
00055 
00056   /*
00057    * On inhibe le gestionnaire d'erreur
00058    */
00059   _MEDmodeErreurVerrouiller();
00060  if (_MEDcheckVersion30(fid) < 0) goto ERROR;
00061 
00062   if ( (_MED_ACCESS_MODE = _MEDmodeAcces(fid) ) == MED_ACC_UNDEF ) {
00063     MED_ERR_(_ret,MED_ERR_UNRECOGNIZED,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
00064     goto ERROR;
00065   }
00066 
00067   if ( _MED_ACCESS_MODE == MED_ACC_RDONLY) {
00068     MED_ERR_(_ret,MED_ERR_INVALID,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
00069     ISCRUTE_int(_MED_ACCESS_MODE);
00070     goto ERROR;
00071   }
00072 
00073   if (basisfuncit < 1 ) {
00074     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_PARAMETER,"");
00075     ISCRUTE(basisfuncit);
00076     goto ERROR;
00077   }
00078 
00079   NOFINALBLANK(interpname,ERROR);
00080   strcat(_path,interpname);
00081   /*
00082    * Ouverture du  DataGroup /INTERP/<interpname>
00083    */
00084   if ((_interpid = _MEDdatagroupOpen(fid,_path)) < 0) {
00085     MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,interpname);
00086     SSCRUTE(_path);goto ERROR;
00087   }
00088 
00089   /*
00090    * Lecture de l'attribut MED_NOM_NBV
00091    */
00092   if ( _MEDattrEntierLire(_interpid,MED_NOM_NBV,&_nvariable) < 0 ) {
00093     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,_path);
00094     SSCRUTE(MED_NOM_NBV);ISCRUTE(_nvariable);
00095     goto ERROR;
00096   }
00097 
00098 
00099   sprintf(_basisfuncname,"%0*li",MED_MAX_PARA,(long ) basisfuncit);
00100   /*
00101    * Si le DataGroup /INTERP/<interpname>/<basisfuncit> n'existe pas, on le cree
00102    */
00103   if ((_bid = _MEDdatagroupOpen(_interpid,_basisfuncname)) < 0)
00104     if ((_bid = _MEDdatagroupCreer(_interpid,_basisfuncname)) < 0) {
00105       MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_DATAGROUP,_path);SSCRUTE(_basisfuncname);
00106       goto ERROR;
00107     }
00108 
00109   /*
00110    * Creation/Ecriture de l'attribut MED_NOM_MCF
00111    */
00112   if ( _MEDattributeIntWr(_bid,MED_NOM_MCF,&ncoef) < 0 ) {
00113     MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_ATTRIBUTE,_path);SSCRUTE(_basisfuncname);
00114     SSCRUTE(MED_NOM_MCF);ISCRUTE(ncoef);
00115     goto ERROR;
00116   }
00117 
00118   /*TODO : Créer un _MEDdatasetSimpleWr */
00119   if ( MEDfilterEntityCr(fid, ncoef, 1, 1, MED_ALL_CONSTITUENT,
00120                          MED_NO_INTERLACE,MED_UNDEF_PFLMODE,
00121                          MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
00122     MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_FILTER,MED_ERR_INTERNAL_MSG);
00123     goto ERROR;
00124   }
00125 
00126   if ( _MEDdatasetWr(_bid,MED_NOM_COE,MED_INTERNAL_FLOAT64,&_filter, coefficient) < 0) {
00127     MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_DATASET,MED_NOM_COE);
00128     SSCRUTE(_path);SSCRUTE(_basisfuncname);
00129     goto ERROR;
00130   }
00131 
00132   if ( MEDfilterClose(&_filter) < 0 ) {
00133     MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_FILTER,MED_ERR_STRUCT_ELEMENT_MSG);
00134     SSCRUTE(_path);
00135     goto ERROR;
00136   }
00137 
00138 
00139   if ( MEDfilterEntityCr(fid, ncoef*_nvariable, 1, 1, MED_ALL_CONSTITUENT,
00140                          MED_NO_INTERLACE,MED_UNDEF_PFLMODE,
00141                          MED_NO_PROFILE, MED_UNDEF_SIZE, NULL, &_filter) < 0 ) {
00142     MED_ERR_(_ret,MED_ERR_CREATE,MED_ERR_FILTER,MED_ERR_INTERNAL_MSG);
00143     goto ERROR;
00144   }
00145 
00146   if ( _MEDdatasetWr(_bid,MED_NOM_POW,MED_INTERNAL_INT,&_filter, power) < 0) {
00147     MED_ERR_(_ret,MED_ERR_WRITE,MED_ERR_DATASET,MED_NOM_POW);
00148     SSCRUTE(_path);SSCRUTE(_basisfuncname);
00149     goto ERROR;
00150   }
00151 
00152   if ( MEDfilterClose(&_filter) < 0 ) {
00153     MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_FILTER,MED_ERR_STRUCT_ELEMENT_MSG);
00154     SSCRUTE(_path);
00155     goto ERROR;
00156   }
00157 
00158 
00159   _ret=0;
00160  ERROR:
00161 
00162   if (_interpid>0)            if (_MEDdatagroupFermer(_interpid) < 0) {
00163     MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_path);
00164     ISCRUTE_id(_interpid);
00165   }
00166 
00167   if (_bid>0)            if (_MEDdatagroupFermer(_bid) < 0) {
00168     MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_basisfuncname);
00169     ISCRUTE_id(_bid);SSCRUTE(_path);
00170   }
00171 
00172   return _ret;
00173 }

Généré le Thu Oct 8 14:26:16 2015 pour MED fichier par  doxygen 1.6.1