Test_MEDinterpBaseFunctionRd.c
Aller à la documentation de ce fichier.
/* This file is part of MED.
*
* COPYRIGHT (C) 1999 - 2021 EDF R&D, CEA/DEN
* MED is free software: you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* MED is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with MED. If not, see <http://www.gnu.org/licenses/>.
*/
#include <med.h>
#define MESGERR 1
#include <med_utils.h>
#include <string.h>
#ifdef DEF_LECT_ECR
#define MODE_ACCES MED_ACC_RDWR
#elif DEF_LECT_AJOUT
#define MODE_ACCES MED_ACC_RDEXT
#else
#define MODE_ACCES MED_ACC_CREAT
#endif
int main (int argc, char **argv)
{
med_err _ret=-1;
med_idt _fid=0;
med_int _ninterp=0;
int _interpit =0;
char _interpname[MED_NAME_SIZE+1]="";
med_geometry_type _geotype =MED_NONE;
med_bool _cellnodes =MED_FALSE;
med_int _nbasisfunc =0;
med_int _nvariable =0;
med_int _maxdegree =0;
med_int _nmaxcoefficient =0;
int _basisfuncit =0;
int _powerit =0;
med_int _ncoefficient =0;
med_int* _power =NULL;
med_float* _coefficient =NULL;
int _coefficientit =0;
/* Ouverture en mode creation du fichier "current.med" */
_fid = MEDfileOpen("current.med",MED_ACC_RDONLY);
if (_fid < 0) {
MESSAGE("Erreur a la creation du fichier current.med");
return -1;
}
if ( (_ninterp = MEDnInterp(_fid)) <0) _ret=_ninterp;
ISCRUTE(_ninterp);
for ( _interpit=1; _interpit<= _ninterp; ++_interpit) {
if ( (_ret = MEDinterpInfo(_fid,
_interpit,
_interpname,
&_geotype,
&_cellnodes,
&_nbasisfunc,
&_nvariable,
&_maxdegree,
&_nmaxcoefficient
) <0) ) {
MESSAGE("Erreur à la demande d'information de la fonction d'interpolation n°");
ISCRUTE(_interpit);
goto ERROR;
}
fprintf(stdout,"Fonction d'interpolation n° %d |%s| sur le type géométrique n° %d\n",
_interpit,_interpname, _geotype);
if ( _cellnodes )
if ( _nbasisfunc == (_geotype % 100) )
fprintf(stdout,"\t Les noeuds de construction sont les noeuds de la maille de référence.\n");
else {
MESSAGE("Erreur : le nombre de noeuds de construction "\
"est différent du nombre de noeuds de la maille de référence.\n");
ISCRUTE(_nbasisfunc); ISCRUTE(_geotype % 100);
}
if ( _nvariable != (_geotype / 100) ) {
MESSAGE("Erreur : le nombre de variables "\
"est différent de la dimension de l'espace de la maille de référence.\n");
ISCRUTE(_nvariable); ISCRUTE (_geotype / 100);
} else
fprintf(stdout,"\t Il y a %d fonctions de base avec %d variables\n ",_nbasisfunc,_nvariable);
fprintf(stdout,"\t Le degré maximum des fonctions de base est %d et possèdent au maximum %d coefficients\n"
,_maxdegree,_nmaxcoefficient);
_coefficient = (med_float*) calloc(sizeof(med_float),_nmaxcoefficient);
_power = (med_int*) calloc(sizeof(med_int),_nvariable*_nmaxcoefficient);
for ( _basisfuncit=1; _basisfuncit<= _nbasisfunc; ++_basisfuncit) {
if ( (_ret = MEDinterpBaseFunctionRd( _fid,
_interpname,
_basisfuncit,
&_ncoefficient,
_power,
_coefficient
) <0) ) {
MESSAGE("Erreur à la lecture de la fonction de base n°");ISCRUTE(_basisfuncit);goto ERROR;
} else {
fprintf(stdout,"\n\t Tableau de coefficients de la fonctions de base n° %d :\n\t",_basisfuncit);
for ( _coefficientit=0; _coefficientit< _ncoefficient; ++_coefficientit)
fprintf(stdout," %4f ",_coefficient[_coefficientit]);
fprintf(stdout,"\n\t Tableau de puissances de la fonctions de base n° %d :\n\t",_basisfuncit);
for ( _powerit=0; _powerit< _nvariable*_ncoefficient; ++_powerit)
fprintf(stdout," %4d ",_power[_powerit]);
}
}
fprintf(stdout,"\n");
free(_coefficient);
free(_power);
}
_ret=0;
ERROR:
if (MEDfileClose(_fid) < 0) {
MESSAGE("ERROR : file closing");
return -1;
}
return _ret;
}
MEDnInterp
MEDC_EXPORT med_int MEDnInterp(const med_idt fid)
Cette routine renvoie le nombre d'interpolations disponibles dans le fichier.
Definition: MEDnInterp.c:34
MEDfileOpen
MEDC_EXPORT med_idt MEDfileOpen(const char *const filename, const med_access_mode accessmode)
Ouverture d'un fichier MED.
Definition: MEDfileOpen.c:42
med_geometry_type
int med_geometry_type
Definition: med.h:194
med_idt
hid_t med_idt
Definition: med.h:331
MED_ACC_RDONLY
Definition: med.h:120
med_err
herr_t med_err
Definition: med.h:332
MED_FALSE
Definition: med.h:260
MESSAGE
#define MESSAGE(chaine)
Definition: med_utils.h:324
med_int
int med_int
Definition: med.h:342
med.h
ISCRUTE
#define ISCRUTE(entier)
Definition: med_utils.h:313
med_bool
med_bool
Definition: med.h:260
med_float
double med_float
Definition: med.h:336
main
int main(int argc, char **argv)
Definition: Test_MEDinterpBaseFunctionRd.c:31
MED_NONE
#define MED_NONE
Definition: med.h:231
MEDinterpInfo
MEDC_EXPORT med_err MEDinterpInfo(const med_idt fid, const int interpit, char *const interpname, med_geometry_type *const geotype, med_bool *const cellnode, med_int *const nbasisfunc, med_int *const nvariable, med_int *const maxdegree, med_int *const nmaxcoef)
Cette fonction informe des caractéristiques de la fonction d'interpolation n° interpit.
Definition: MEDinterpInfo.c:43
MEDinterpBaseFunctionRd
MEDC_EXPORT 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_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:81
MEDfileClose
MEDC_EXPORT med_err MEDfileClose(med_idt fid)
Fermeture d'un fichier MED.
Definition: MEDfileClose.c:30
med_utils.h