int main (int argc, char **argv)
{
med_err ret = 0;
med_idt fid;
/* la dimension du maillage */
med_int mdim = 2;
/* nom du maillage de longueur maxi MED_NAME_SIZE */
char maa[MED_NAME_SIZE+1] = "maa1";
/* le nombre de noeuds */
med_int nnoe = 4;
/* table des coordonnees
(dimension * nombre de noeuds) */
med_float coo[8] = {0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0};
med_float coo_2[8] = {0.0, 10.0, 20.0, 30.0, 40.0, 50.0, 60.0, 70.0};
/* tables des noms et des unites des coordonnees
(dimension*MED_SNAME_SIZE+1) */
/* 12345678901234561234567890123456*/
char nomcoo[2*MED_SNAME_SIZE+1] = "x y ";
char unicoo[2*MED_SNAME_SIZE+1] = "cm cm ";
/* tables des noms, numeros, numeros de familles des noeuds
autant d'elements que de noeuds - les noms ont pout longueur
MED_SNAME_SIZE */
/* 1234567890123456123456789012345612345678901234561234567890123456*/
char nomnoe[4*MED_SNAME_SIZE+1] = "nom1 nom2 nom3 nom4 ";
med_int numnoe[4] = {1,2,3,4};
med_int nufano[4] = {0,1,2,2};
const med_float phi1=0;
const med_float phi2=M_PI_4;
const med_float phi3=M_PI_2;
med_float trsf1[7]= {0,0,0, cos(phi1/2), 0, sin(phi1/2), 0 };
med_float trsf2[7]= {0,0,0, cos(phi2/2), 0, sin(phi2/2), 0 };
med_float trsf3[7]= {0,0,0, cos(phi3/2), 0, sin(phi3/2), 0 };
/* ouverture du fichier */
if ((fid = MEDfileOpen("Test_MEDmeshNodeCoordinateTrsfWr.med",MODE_ACCES)) < 0){
MESSAGE("Erreur à l'ouverture du fichier : ");
return -1;
}
/* Creation du maillage "maa" de type MED_NON_STRUCURE
et de dimension 2 */
if (MEDmeshCr( fid, maa, mdim, mdim, MED_UNSTRUCTURED_MESH,
"un maillage pour Test_MEDmeshNodeCoordinateTrsfWr","s", MED_SORT_DTIT,
MED_CARTESIAN, nomcoo, unicoo) < 0) {
MESSAGE("Erreur a la creation du maillage : "); SSCRUTE(maa);
ret = -1;
}
/* Ecriture des coordonnees des noeuds en mode MED_FULL_INTERLACE :
(X1,Y1, X2,Y2, X3,Y3, ...) dans un repere cartesien */
if (MEDmeshNodeCoordinateWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
MED_FULL_INTERLACE,nnoe, coo) < 0) {
MESSAGE("Erreur a l'ecriture des coordonnees des noeuds");
ret = -1;
}
/* Ecriture des noms des noeuds (optionnel dans un maillage MED) */
if (MEDmeshEntityNameWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_NODE,MED_NONE,nnoe,nomnoe) < 0) {
MESSAGE("Erreur a l'ecriture des noms des noeuds");
ret = -1;
}
/* Ecriture des numeros des noeuds (optionnel dans un maillage MED) */
if (MEDmeshEntityNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_NODE,MED_NONE,nnoe,numnoe) < 0) {
MESSAGE("Erreur a l'ecriture des numeros des noeuds");
ret = -1;
}
/* Ecriture des numeros de famille des noeuds */
if (MEDmeshEntityFamilyNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_NODE,MED_NONE,nnoe,nufano) < 0) {
MESSAGE("Erreur a l'ecriture des numeros de familles des noeuds");
ret = -1;
}
/* Transformation des coordonnees des noeuds coo_1 :
*/
fprintf(stdout,"Un message d'erreur est attendu: \n");
if ( MEDmeshNodeCoordinateTrsfWr(fid,maa,MED_NO_DT,MED_NO_IT,0.4,trsf1) >= 0) {
MESSAGE("Erreur a l'ecriture de la transformation géométrique n°1");
MESSAGE("Aucune transformation géométrique à la séquence de calcul MED_NO_DT,MED_NO_IT n'est possible.");
}
fprintf(stdout,"Fin du message d'erreur attendu\n :");
/* Transformation des coordonnees des noeuds coo_1 :
*/
if ( MEDmeshNodeCoordinateTrsfWr(fid,maa,1,1,0.4,trsf1) < 0) {
MESSAGE("Erreur a l'ecriture de la transformation géométrique n°1");
ret = -1;
}
/* Ecriture des coordonnees des noeuds en mode MED_FULL_INTERLACE :
(X1,Y1, X2,Y2, X3,Y3, ...) dans un repere cartesien */
fprintf(stdout,"Un message d'erreur est attendu : \n");
if (MEDmeshNodeCoordinateWr(fid,maa,1,1, 0.5,
MED_FULL_INTERLACE,nnoe, coo_2) >= 0) {
MESSAGE("Erreur a l'ecriture des coordonnees des noeuds");
MESSAGE("Aucune ecriture des coordonnees des noeuds n'est possible sur une séquence de calcul comportant "\
"une transformation géométrique.");
}
fprintf(stdout,"Fin du message d'erreur attendu :\n");
/* Ecriture des coordonnees des noeuds en mode MED_FULL_INTERLACE :
(X1,Y1, X2,Y2, X3,Y3, ...) dans un repere cartesien */
if (MEDmeshNodeCoordinateWr(fid,maa,2,1, 0.5,
MED_FULL_INTERLACE,nnoe, coo_2) < 0) {
MESSAGE("Erreur a l'ecriture des coordonnees des noeuds");
ret = -1;
}
/* Transformation des coordonnees des noeuds coo_2 :
*/
fprintf(stdout,"Un message d'erreur est attendu\n :");
if ( MEDmeshNodeCoordinateTrsfWr(fid,maa,2,1,0.5,trsf2) >= 0) {
MESSAGE("Erreur a l'ecriture de la transformation géométrique n°2");
MESSAGE("Aucune transformation géométrique à une séquence de calcul comportant de nouvelles coordonnées n'est possible.");
}
fprintf(stdout,"Fin du message d'erreur attendu\n :");
/* Transformation des coordonnees des noeuds coo_2 :
*/
if ( MEDmeshNodeCoordinateTrsfWr(fid,maa,2,2,0.6,trsf2) < 0) {
MESSAGE("Erreur a l'ecriture de la transformation géométrique n°2");
ret = -1;
}
/* Transformation des coordonnees des noeuds coo_3 :
*/
if ( MEDmeshNodeCoordinateTrsfWr(fid,maa,3,2,0.7,trsf3) < 0) {
MESSAGE("Erreur a l'ecriture de la transformation géométrique n°3");
ret = -1;
}
if ( MEDmeshComputationStepCr(fid,maa, 3,2, 3,3, 3.3) < 0) {
MESSAGE("Erreur a la creation d'un pas de temps du maillage maa");
}
/* Modification des numeros des noeuds (optionnel dans un maillage MED) */
numnoe[2]=200;
if (MEDmeshEntityNumberWr(fid,maa,3,3,MED_NODE,MED_NONE,nnoe,numnoe) < 0) {
MESSAGE("Erreur a l'ecriture des numeros des noeuds");
ret = -1;
}
if ( MEDmeshComputationStepCr(fid,maa, 3,3, 3,4, 3.4) < 0) {
MESSAGE("Erreur a la creation d'un pas de temps du maillage maa");
}
/*TODO : Ecrire un test de comparaison interne next prev et ordre de découverte itératif */
/* Fermeture du fichier */
if (MEDfileClose(fid) < 0) {
MESSAGE("Erreur a la fermeture du fichier test4.med");
return -1;
}
return ret;
}
int main (int argc, char **argv) {
med_idt fid;
const char meshname[MED_NAME_SIZE+1] = "2D unstructured mesh";
const med_int spacedim = 2;
const med_int meshdim = 2;
const char axisname[2*MED_SNAME_SIZE+1] = "x y ";
const char unitname[2*MED_SNAME_SIZE+1] = "cm cm ";
const med_float coordinates[30] = { 2.,1., 7.,1., 12.,1., 17.,1., 22.,1.,
2.,6., 7.,6., 12.,6., 17.,6., 22.,6.,
2.,11., 7.,11., 12.,11., 17.,11., 22.,11.};
const med_int nnodes = 15;
const med_int triaconnectivity[24] = { 1,7,6, 2,7,1, 3,7,2, 8,7,3,
13,7,8, 12,7,13, 11,7,12, 6,7,11 };
const med_int ntria3 = 8;
const med_int quadconnectivity[16] = {3,4,9,8, 4,5,10,9,
15,14,9,10, 13,8,9,14};
const med_int nquad4 = 4;
const char familyname[MED_NAME_SIZE+1] = "BOUNDARY_VERTICES";
const char groupname[MED_LNAME_SIZE+1] = "MESH_BOUNDARY_VERTICES";
const med_int familynumbers[15] = { 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1 };
int ret=-1;
/* open MED file */
fid = MEDfileOpen("UsesCase_MEDmesh_10.med",MED_ACC_CREAT);
if (fid < 0) {
MESSAGE("ERROR : file creation ...");
goto ERROR;
}
/* write a comment in the file */
if (MEDfileCommentWr(fid,"A 2D unstructured mesh : 15 nodes, 12 cells") < 0) {
MESSAGE("ERROR : write file description ...");
goto ERROR;
}
/* mesh creation : a 2D unstructured mesh */
if (MEDmeshCr(fid, meshname, spacedim, meshdim, MED_UNSTRUCTURED_MESH,
"A 2D unstructured mesh","",MED_SORT_DTIT,MED_CARTESIAN, axisname, unitname) < 0) {
MESSAGE("ERROR : mesh creation ...");
goto ERROR;
}
/* nodes coordinates in a cartesian axis in full interlace mode
(X1,Y1, X2,Y2, X3,Y3, ...) with no iteration and computation step
*/
if (MEDmeshNodeCoordinateWr(fid, meshname, MED_NO_DT, MED_NO_IT, 0.0,
MED_FULL_INTERLACE, nnodes, coordinates) < 0) {
MESSAGE("ERROR : nodes coordinates ...");
goto ERROR;
}
/* cells connectiviy is defined in nodal mode */
if (MEDmeshElementConnectivityWr(fid, meshname, MED_NO_DT, MED_NO_IT, 0.0, MED_CELL, MED_TRIA3,
MED_NODAL, MED_FULL_INTERLACE, ntria3, triaconnectivity) < 0) {
MESSAGE("ERROR : triangular cells connectivity ...");
goto ERROR;
}
if (MEDmeshElementConnectivityWr(fid, meshname, MED_NO_DT, MED_NO_IT, 0.0, MED_CELL, MED_QUAD4,
MED_NODAL, MED_FULL_INTERLACE, nquad4, quadconnectivity) < 0) {
MESSAGE("ERROR : quadrangular cells connectivity ...");
goto ERROR;
}
/* create family 0 : by default, all mesh entities family number is 0 */
if (MEDfamilyCr(fid, meshname,MED_NO_NAME, 0, 0, MED_NO_GROUP) < 0) {
MESSAGE("ERROR : family creation ...");
goto ERROR;
}
/* create a family for boundary vertices : by convention a nodes family number is > 0,
and an element family number is < 0 */
if (MEDfamilyCr(fid, meshname,familyname, 1, 1, groupname) < 0) {
MESSAGE("ERROR : family creation ...");
goto ERROR;
}
/* write family number for nodes */
if (MEDmeshEntityFamilyNumberWr(fid, meshname, MED_NO_DT, MED_NO_IT,
MED_NODE, MED_NONE, nnodes, familynumbers) < 0) {
MESSAGE("ERROR : nodes family numbers ...");
goto ERROR;
}
ret=0;
ERROR:
/* close MED file */
if (MEDfileClose(fid) < 0) {
MESSAGE("ERROR : close file ...");
ret=-1;
}
return ret;
}
int main (int argc, char **argv)
{
med_idt fid;
char maa[MED_NAME_SIZE+1] = "maa1";
med_int mdim = 3;
med_int n=2;
/* connectivite nodale */
med_int indexf[3] = {1,5,9};
med_int nf = 3;
med_int indexn[9] = {1,4,7,10,13,16,19,22,25};
med_int nn = 9;
med_int conn[24] = {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,
20,21,22,23,24};
/* connectivite descendante */
med_int indexf2[3] = {1,5,9};
med_int nf2 = 3;
med_int conn2[8] = {1,2,3,4,5,6,7,8};
med_int indexn2[8] = {MED_TRIA3,MED_TRIA3,MED_TRIA3,MED_TRIA3,
MED_TRIA3,MED_TRIA3,MED_TRIA3,MED_TRIA3};
med_int nn2=8;
/* 12345678901234561234567890123456*/
char nom[MED_SNAME_SIZE*2+1]="polyedre1 polyedre2 ";
med_int num[2] = {1,2};
med_int fam[2] = {0,-1};
char nomcoo[3*MED_SNAME_SIZE+1] = "x y z ";
char unicoo[3*MED_SNAME_SIZE+1] = "cm cm cm ";
/* Creation du fichier test25.med */
if ((fid = MEDfileOpen("test25.med",MODE_ACCES)) < 0) {
MESSAGE("Erreur a la creation du fichier test25.med");
return -1;
}
printf("Creation du fichier test25.med \n");
/* Creation du maillage */
if (MEDmeshCr( fid, maa, mdim, mdim, MED_UNSTRUCTURED_MESH,
"un maillage pour test25","s", MED_SORT_DTIT,
MED_CARTESIAN, nomcoo, unicoo) < 0) {
MESSAGE("Erreur a la creation du maillage");
return -1;
}
printf("Creation du maillage \n");
/* Ecriture de la connectivite des mailles polyedres en mode nodal */
if (MEDmeshPolyhedronWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,MED_CELL,MED_NODAL,
nf,indexf,nn,indexn,conn) < 0) {
MESSAGE("Erreur a l'ecriture de la connectivite des mailles MED_POLYGONE");
return -1;
}
printf("Ecriture des connectivites de mailles de type MED_POLYEDRE en mode nodal \n");
/* Ecriture des connectivites des mailles polyedres en mode descendant */
if (MEDmeshPolyhedronWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,MED_CELL,MED_DESCENDING,
nf2,indexf2,nn2,indexn2,conn2) < 0) {
MESSAGE("Erreur a l'ecriture des connectivites des mailles MED_POLYEDRE en mode descendant");
return -1;
}
printf("Ecriture des connectivites de mailles de type MED_POLYEDRE en mode descendant \n");
/* Ecriture des noms des polyedres */
/* ecriture (optionnelle) des noms des polyedres */
if (MEDmeshEntityNameWr(fid,maa,MED_NO_DT,MED_NO_IT,
MED_CELL,MED_POLYHEDRON,n,nom) < 0) {
MESSAGE("Erreur a l'ecriture des noms des polyedres");
return -1;
}
printf("Ecriture des noms des polyedres \n");
/* ecriture (optionnelle) des numeros des polyedres */
if (MEDmeshEntityNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,
MED_CELL,MED_POLYHEDRON,n,num) < 0) {
MESSAGE("Erreur a l'ecriture des numeros des polyedres");
return -1;
}
printf("Ecriture des numeros des polyedres \n");
/* ecriture des numeros des familles des polyedres */
if (MEDmeshEntityFamilyNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,
MED_CELL,MED_POLYHEDRON,n,fam) < 0) {
MESSAGE("Erreur a l'ecriture des numeros de familles des polyedres");
return -1;
}
printf("Ecriture des numeros des familles des polyedres \n");
/* Fermeture du fichier */
if (MEDfileClose(fid) < 0) {
MESSAGE("Erreur a la fermeture du fichier");
return -1;
}
printf("Fermeture du fichier test25.med \n");
return 0;
}
int main (int argc, char **argv)
{
med_err ret = 0;
med_idt fid;
/* la dimension du maillage */
med_int mdim = 2;
/* nom du maillage de longueur maxi MED_NAME_SIZE */
char maa[MED_NAME_SIZE+1] = "maa1";
/* le nombre de noeuds */
med_int nnoe = 4;
/* table des coordonnees
(dimension * nombre de noeuds) */
med_float coo[8] = {0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0};
med_float coo_2[8] = {0.0, 10.0, 20.0, 30.0, 40.0, 50.0, 60.0, 70.0};
/* tables des noms et des unites des coordonnees
(dimension*MED_SNAME_SIZE+1) */
/* 12345678901234561234567890123456*/
char nomcoo[2*MED_SNAME_SIZE+1] = "x y ";
char unicoo[2*MED_SNAME_SIZE+1] = "cm cm ";
/* tables des noms, numeros, numeros de familles des noeuds
autant d'elements que de noeuds - les noms ont pout longueur
MED_SNAME_SIZE */
/* 1234567890123456123456789012345612345678901234561234567890123456*/
char nomnoe[4*MED_SNAME_SIZE+1] = "nom1 nom2 nom3 nom4 ";
med_int numnoe[4] = {1,2,3,4};
med_int nufano[4] = {0,1,2,2};
/* Creation du fichier "test4.med" */
/* ouverture du fichier */
if ((fid = MEDfileOpen("test4.med",MODE_ACCES)) < 0){
MESSAGE("Erreur à l'ouverture du fichier : ");
return -1;
}
/* Creation du maillage "maa" de type MED_NON_STRUCURE
et de dimension 2 */
if (MEDmeshCr( fid, maa, mdim, mdim, MED_UNSTRUCTURED_MESH,
"un maillage pour test4","s", MED_SORT_DTIT,
MED_CARTESIAN, nomcoo, unicoo) < 0) {
MESSAGE("Erreur a la creation du maillage : "); SSCRUTE(maa);
ret = -1;
}
/* Ecriture des coordonnees des noeuds en mode MED_FULL_INTERLACE :
(X1,Y1, X2,Y2, X3,Y3, ...) dans un repere cartesien */
if (MEDmeshNodeCoordinateWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
MED_FULL_INTERLACE,nnoe, coo) < 0) {
MESSAGE("Erreur a l'ecriture des coordonnees des noeuds");
ret = -1;
}
/* Ecriture des attributs des noeuds du maillage */
if (MEDmeshAttributeWr(fid,maa, 0, nnoe, 0) < 0 ) {
MESSAGE("Erreur a l'ecriture des attributs des noeuds du maillage");
ret = -1;
}
/* Ecriture des noms des noeuds (optionnel dans un maillage MED) */
if (MEDmeshEntityNameWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_NODE,MED_NONE,nnoe,nomnoe) < 0) {
MESSAGE("Erreur a l'ecriture des noms des noeuds");
ret = -1;
}
/* Ecriture des numeros des noeuds (optionnel dans un maillage MED) */
if (MEDmeshEntityNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_NODE,MED_NONE,nnoe,numnoe) < 0) {
MESSAGE("Erreur a l'ecriture des numeros des noeuds");
ret = -1;
}
/* Ecriture des numeros de famille des noeuds */
if (MEDmeshEntityFamilyNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_NODE,MED_NONE,nnoe,nufano) < 0) {
MESSAGE("Erreur a l'ecriture des numeros de familles des noeuds");
ret = -1;
}
/* Ecriture des coordonnees des noeuds en mode MED_FULL_INTERLACE :
(X1,Y1, X2,Y2, X3,Y3, ...) dans un repere cartesien */
if (MEDmeshNodeCoordinateWr(fid,maa,0,5,0.5,
MED_FULL_INTERLACE,nnoe, coo_2) < 0) {
MESSAGE("Erreur a l'ecriture des coordonnees des noeuds");
ret = -1;
}
/* Ecriture des coordonnees des noeuds en mode MED_FULL_INTERLACE :
(X1,Y1, X2,Y2, X3,Y3, ...) dans un repere cartesien */
if (MEDmeshNodeCoordinateWr(fid,maa,2,0,2.0,
MED_FULL_INTERLACE,nnoe, coo) < 0) {
MESSAGE("Erreur a l'ecriture des coordonnees des noeuds");
ret = -1;
}
/*TODO : Ecrire un test de comparaison next prev et ordre de découverte itératif */
MEDmeshComputationStepCr(fid,maa,-1,-1,0,2,0.25);
MEDmeshComputationStepCr(fid,maa,0,2,0,3,0.35);
/* Fermeture du fichier */
if (MEDfileClose(fid) < 0) {
MESSAGE("Erreur a la fermeture du fichier test4.med");
return -1;
}
return ret;
}
int main (int argc, char **argv)
{
med_idt fid;
char maa[MED_NAME_SIZE+1] = "maa1";
med_int mdim = 3;
med_int index[4] = {1,6,12,17};
med_int con[16] = {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16};
med_int n=3,ni = 4;
/* 123456789012345612345678901234561234567890123456 */
char nom[MED_SNAME_SIZE*3+1]="poly1 poly2 poly3 ";
med_int num[3] = {1,2,3};
med_int fam[3] = {0,-1,-2};
char nomcoo[3*MED_SNAME_SIZE+1] = "x y z ";
char unicoo[3*MED_SNAME_SIZE+1] = "cm cm cm ";
/* Creation du fichier test23.med */
fid = MEDfileOpen("test23.med",MODE_ACCES);
if (fid == -1) {
MESSAGE("Erreur a la creation du fichier test23.med");
return -1;
}
printf("Creation du fichier test23.med \n");
/* Creation du maillage */
if (MEDmeshCr( fid, maa, mdim, mdim, MED_UNSTRUCTURED_MESH,
"un maillage pour test23","s", MED_SORT_DTIT,
MED_CARTESIAN, nomcoo, unicoo) < 0) {
MESSAGE("Erreur a la creation du maillage");
return -1;
}
printf("Creation du maillage \n");
/* Ecriture de la connectivite des mailles polygones en mode nodal */
if (MEDmeshPolygonWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,MED_CELL,MED_NODAL,
ni,index,con) < 0) {
MESSAGE("Erreur a l'ecriture de la connectivite des mailles MED_POLYGONE");
return -1;
}
printf("Ecriture des connectivites de mailles de type MED_POLYGONE en mode nodal \n");
/* Ecriture des noms des polygones */
/* ecriture (optionnelle) des noms des polygones */
if (MEDmeshEntityNameWr(fid,maa,MED_NO_DT,MED_NO_IT,
MED_CELL,MED_POLYGON,n,nom) < 0) {
MESSAGE("Erreur a l'ecriture des noms des polygones");
return -1;
}
printf("Ecriture des noms des polygones \n");
/* ecriture (optionnelle) des numeros des polygones */
if (MEDmeshEntityNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,
MED_CELL,MED_POLYGON,n,num) < 0) {
MESSAGE("Erreur a l'ecriture des numeros des polygones");
return -1;
}
printf("Ecriture des numeros des polygones \n");
/* ecriture des numeros des familles des polygones */
if (MEDmeshEntityFamilyNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,
MED_CELL,MED_POLYGON,n,fam) < 0) {
MESSAGE("Erreur a l'ecriture des numeros de familles des polygones");
return -1;
}
printf("Ecriture des numeros des familles des polygones \n");
/* Fermeture du fichier */
if (MEDfileClose(fid) < 0) {
MESSAGE("Erreur a la fermeture du fichier");
return -1;
}
printf("Fermeture du fichier test23.med \n");
return 0;
}