Test_MEDinterpBaseFunctionWr.c
Aller à la documentation de ce fichier.
/* This file is part of MED.
*
* COPYRIGHT (C) 1999 - 2020 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;
/*Exemple 1 :
- Elément de référence de type géométrique MED_TRIA3
- Point X(X1,X2) quelconque dans le plan de l'élément de référence
- Fonctions de base : P1(X)=1-X1-X2 ; P2(X)=X1; P3(X)=X2;
(issu du choix de la base polynomiale (1,X1,X2)
et des trois noeuds de la maille de référence pour
construire l'interpolation)
*/
const char _interpname1[] ="interpname1";
med_geometry_type _geotype1 =MED_TRIA3;
med_bool _cellnodes1 =MED_TRUE;
med_int _nbasisfunc1 =3;
med_int _nvariable1 =2;
med_int _maxdegree1 =1;
med_int _nmaxcoefficient1=3;
const med_int _ncoefficient1_1 = 3;
const med_int const _power1_1[] = {0,0,1,0,0,1};
const med_float const _coefficient1_1[] = {1,-1,-1};
const med_int _ncoefficient1_2 = 1;
const med_int const _power1_2[] = {0,0,1,0,0,0};
const med_float const _coefficient1_2[] = {0,1,0};
const med_int _ncoefficient1_3 = 1;
const med_int const _power1_3[] = {0,0,0,0,0,1};
const med_float const _coefficient1_3[] = {0,0,1};
/* Ouverture en mode creation du fichier "current.med" */
_fid = MEDfileOpen("current.med",MODE_ACCES);
if (_fid < 0) {
MESSAGE("Erreur a la creation du fichier current.med");
return -1;
}
if ( (_ret = MEDinterpBaseFunctionWr( _fid,
_interpname1,
1,
_ncoefficient1_1,
_power1_1,
_coefficient1_1) <0) ) {
MESSAGE("Erreur à l'écriture de la fonction de base n°1 de la fct. d'intp n°1");
goto ERROR;
}
if ( (_ret = MEDinterpBaseFunctionWr( _fid,
_interpname1,
2,
_ncoefficient1_2,
_power1_2,
_coefficient1_2) <0) ) {
MESSAGE("Erreur à l'écriture de la fonction de base n°2 de la fct. d'intp n°1");
goto ERROR;
}
if ( (_ret = MEDinterpBaseFunctionWr( _fid,
_interpname1,
3,
_ncoefficient1_3,
_power1_3,
_coefficient1_3) <0) ) {
MESSAGE("Erreur à l'écriture de la fonction de base n°3 de la fct. d'intp n°1");
goto ERROR;
}
_ret=0;
ERROR:
if (MEDfileClose(_fid) < 0) {
MESSAGE("ERROR : file closing");
return -1;
}
return _ret;
}
MODE_ACCES
#define MODE_ACCES
Definition 3.0.8/test10.c:35
main
int main(int argc, char **argv)
Definition Test_MEDinterpBaseFunctionWr.c:31
MEDfileClose
MEDC_EXPORT med_err MEDfileClose(med_idt fid)
Fermeture d'un fichier MED.
Definition MEDfileClose.c:30
MEDfileOpen
MEDC_EXPORT med_idt MEDfileOpen(const char *const filename, const med_access_mode accessmode)
Ouverture d'un fichier MED.
Definition MEDfileOpen.c:42
MEDinterpBaseFunctionWr
MEDC_EXPORT 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_geometry_type
int med_geometry_type
Definition med.h:194
med_bool
med_bool
Definition med.h:260
MED_TRUE
@ MED_TRUE
Definition med.h:260
MED_TRIA3
#define MED_TRIA3
Definition med.h:203
med_int
int med_int
Definition med.h:333
med_float
double med_float
Definition med.h:327
med_err
herr_t med_err
Definition med.h:323
med_idt
hid_t med_idt
Definition med.h:322
med.h
MESSAGE
#define MESSAGE(chaine)
Definition med_utils.h:324
med_utils.h