00001 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00015 
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 
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 
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 
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 
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 
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   
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 }