int main(int argc,char *argv[]) {
MMG5_pMesh mmgMesh;
MMG5_pSol mmgSol;
int ier,k;
fprintf(stdout," -- TEST MMG3DLIB \n");
/** ------------------------------ STEP I -------------------------- */
/** 1) Initialisation of mesh and sol structures */
/* args of InitMesh: mesh=&mmgMesh, sol=&mmgSol */
mmgMesh = NULL;
mmgSol = NULL;
MMG5_Init_mesh(&mmgMesh,&mmgSol,NULL);
/** 2) Build mesh in MMG5 format */
/** Two solutions: just use the MMG5_loadMesh function that will read a .mesh(b)
file formatted or manually set your mesh using the MMG5_Set* functions */
/** Manually set of the mesh */
/** a) give the size of the mesh: 12 vertices, 12 tetra, 20 triangles, 0 edges */
if ( !MMG5_Set_meshSize(mmgMesh,12,12,20,0) ) exit(EXIT_FAILURE);
/** b) give the vertices: for each vertex, give the coordinates, the reference
and the position in mesh of the vertex */
mmgMesh->point[1].c[0] = 0.; mmgMesh->point[1].c[1] = 0.; mmgMesh->point[1].c[2] = 0.; mmgMesh->point[1].ref = 0;
/* or with the api function :
if ( !MMG5_Set_vertex(mmgMesh,0 ,0 ,0 ,0, 1) ) exit(EXIT_FAILURE); */
mmgMesh->point[2].c[0] = 0.5; mmgMesh->point[2].c[1] = 0; mmgMesh->point[2].c[2] = 0; mmgMesh->point[2].ref = 0;
mmgMesh->point[3].c[0] = 0.5; mmgMesh->point[3].c[1] = 0; mmgMesh->point[3].c[2] = 1; mmgMesh->point[3].ref = 0;
mmgMesh->point[4].c[0] = 0; mmgMesh->point[4].c[1] = 0; mmgMesh->point[4].c[2] = 1; mmgMesh->point[4].ref = 0;
mmgMesh->point[5].c[0] = 0; mmgMesh->point[5].c[1] = 1; mmgMesh->point[5].c[2] = 0; mmgMesh->point[5].ref = 0;
mmgMesh->point[6].c[0] = 0.5; mmgMesh->point[6].c[1] = 1; mmgMesh->point[6].c[2] = 0; mmgMesh->point[6].ref = 0;
mmgMesh->point[7].c[0] = 0.5; mmgMesh->point[7].c[1] = 1; mmgMesh->point[7].c[2] = 1; mmgMesh->point[7].ref = 0;
mmgMesh->point[8].c[0] = 0; mmgMesh->point[8].c[1] = 1; mmgMesh->point[8].c[2] = 1; mmgMesh->point[8].ref = 0;
mmgMesh->point[9].c[0] = 1; mmgMesh->point[9].c[1] = 0; mmgMesh->point[9].c[2] = 0; mmgMesh->point[9].ref = 0;
mmgMesh->point[10].c[0] = 1; mmgMesh->point[10].c[1] = 1; mmgMesh->point[10].c[2] = 0; mmgMesh->point[10].ref = 0;
mmgMesh->point[11].c[0] = 1; mmgMesh->point[11].c[1] = 0; mmgMesh->point[11].c[2] = 1; mmgMesh->point[11].ref = 0;
mmgMesh->point[12].c[0] = 1; mmgMesh->point[12].c[1] = 1; mmgMesh->point[12].c[2] = 1; mmgMesh->point[12].ref = 0;
/*tetra*/
mmgMesh->tetra[1].v[0] = 1; mmgMesh->tetra[1].v[1] = 2; mmgMesh->tetra[1].v[2] = 4; mmgMesh->tetra[1].v[3] = 8; mmgMesh->tetra[1].ref = 1;
/* or with the api function :
if ( !MMG5_Set_tetrahedra(mmgMesh,1 ,2 ,4 ,8, 1) ) exit(EXIT_FAILURE); */
mmgMesh->tetra[2].v[0] = 8; mmgMesh->tetra[2].v[1] = 3; mmgMesh->tetra[2].v[2] = 2; mmgMesh->tetra[2].v[3] = 7; mmgMesh->tetra[2].ref = 1;
mmgMesh->tetra[3].v[0] = 2; mmgMesh->tetra[3].v[1] = 5; mmgMesh->tetra[3].v[2] = 6; mmgMesh->tetra[3].v[3] = 8; mmgMesh->tetra[3].ref = 1;
mmgMesh->tetra[4].v[0] = 8; mmgMesh->tetra[4].v[1] = 5; mmgMesh->tetra[4].v[2] = 1; mmgMesh->tetra[4].v[3] = 2; mmgMesh->tetra[4].ref = 1;
mmgMesh->tetra[5].v[0] = 2; mmgMesh->tetra[5].v[1] = 7; mmgMesh->tetra[5].v[2] = 8; mmgMesh->tetra[5].v[3] = 6; mmgMesh->tetra[5].ref = 1;
mmgMesh->tetra[6].v[0] = 2; mmgMesh->tetra[6].v[1] = 4; mmgMesh->tetra[6].v[2] = 3; mmgMesh->tetra[6].v[3] = 8; mmgMesh->tetra[6].ref = 1;
mmgMesh->tetra[7].v[0] = 2; mmgMesh->tetra[7].v[1] = 9; mmgMesh->tetra[7].v[2] = 3; mmgMesh->tetra[7].v[3] = 7; mmgMesh->tetra[7].ref = 2;
mmgMesh->tetra[8].v[0] = 7; mmgMesh->tetra[8].v[1] = 11; mmgMesh->tetra[8].v[2] = 9; mmgMesh->tetra[8].v[3] = 12; mmgMesh->tetra[8].ref = 2;
mmgMesh->tetra[9].v[0] = 9; mmgMesh->tetra[9].v[1] = 6; mmgMesh->tetra[9].v[2] = 10; mmgMesh->tetra[9].v[3] = 7; mmgMesh->tetra[9].ref = 2;
mmgMesh->tetra[10].v[0] = 7; mmgMesh->tetra[10].v[1] = 6; mmgMesh->tetra[10].v[2] = 2; mmgMesh->tetra[10].v[3] = 9; mmgMesh->tetra[10].ref = 2;
mmgMesh->tetra[11].v[0] = 9; mmgMesh->tetra[11].v[1] = 12; mmgMesh->tetra[11].v[2] = 7; mmgMesh->tetra[11].v[3] = 10; mmgMesh->tetra[11].ref = 2;
mmgMesh->tetra[12].v[0] = 9; mmgMesh->tetra[12].v[1] = 3; mmgMesh->tetra[12].v[2] = 11; mmgMesh->tetra[12].v[3] = 7; mmgMesh->tetra[12].ref = 2;
/** 3) Build sol in MMG5 format */
/** Two solutions: just use the MMG5_loadMet function that will read a .sol(b)
file formatted or manually set your sol using the MMG5_Set* functions */
/** Manually set of the sol */
/** a) give info for the sol structure: sol applied on vertex entities,
number of vertices=12, the sol is scalar*/
if ( !MMG5_Set_solSize(mmgMesh,mmgSol,MMG5_Vertex,12,MMG5_Scalar) )
exit(EXIT_FAILURE);
/** b) give solutions values and positions */
for(k=1 ; k<=12 ; k++) {
mmgSol->m[k] = 0.5;
/* or with the api function :
if ( !MMG5_Set_scalarSol(mmgSol,0.5,k) ) exit(EXIT_FAILURE); */
}
/** 4) If you don't use the API functions, you MUST call
the MMG5_Set_handGivenMesh() function. Don't call it if you use
the API functions */
MMG5_Set_handGivenMesh(mmgMesh);
/** 5) (not mandatory): check if the number of given entities match with mesh size */
if ( !MMG5_Chk_meshData(mmgMesh,mmgSol) ) exit(EXIT_FAILURE);
/** ------------------------------ STEP II -------------------------- */
/** library call */
ier = MMG5_mmg3dlib(mmgMesh,mmgSol);
if ( ier == MMG5_STRONGFAILURE ) {
fprintf(stdout,"BAD ENDING OF MMG3DLIB: UNABLE TO SAVE MESH\n");
return(ier);
} else if ( ier == MMG5_LOWFAILURE )
fprintf(stdout,"BAD ENDING OF MMG3DLIB\n");
/** ------------------------------ STEP III -------------------------- */
/** get results */
/** Two solutions: just use the MMG5_saveMesh/MMG5_saveMet functions
that will write .mesh(b)/.sol formatted files or manually get your mesh/sol
using the MMG5_getMesh/MMG5_getSol functions */
/** 1) Automatically save the mesh */
/** a) (not mandatory): give the ouptut mesh name using MMG5_Set_outputMeshName
(by default, the mesh is saved in the "mesh.o.mesh" file */
MMG5_Set_outputMeshName(mmgMesh,"result0.mesh");
/** b) function calling */
MMG5_saveMesh(mmgMesh);
/** 2) Automatically save the solution */
/** a) (not mandatory): give the ouptut sol name using MMG5_Set_outputSolName
(by default, the mesh is saved in the "mesh.o.sol" file */
MMG5_Set_outputSolName(mmgMesh,mmgSol,"result0.sol");
/** b) function calling */
MMG5_saveMet(mmgMesh,mmgSol);
/** 3) Free the MMG3D5 structures */
MMG5_Free_all(mmgMesh,mmgSol,NULL);
return(ier);
}
int main(int argc,char *argv[]) {
MMG5_pMesh mesh;
MMG5_pSol met;
int ier;
char stim[32];
atexit(endcod);
tminit(ctim,TIMEMAX);
chrono(ON,&ctim[0]);
/* assign default values */
mesh = NULL;
met = NULL;
MMG5_Init_mesh(&mesh,&met);
/* reset default values for file names */
MMG5_Free_names(mesh,met);
/* command line */
if ( !MMG5_parsar(argc,argv,mesh,met) ) return(1);
/* load data */
fprintf(stdout,"\n -- INPUT DATA\n");
chrono(ON,&ctim[1]);
/* read mesh file */
if ( !MMG5_loadMesh(mesh) ) {
MMG5_Free_all(mesh,met );
return(MMG5_STRONGFAILURE);
}
if ( !MMG5_Set_solSize(mesh,met,MMG5_Vertex,0,MMG5_Scalar) ) {
MMG5_Free_all(mesh,met);
return(MMG5_STRONGFAILURE);
}
/* read metric if any */
ier = MMG5_loadMet(mesh,met);
if ( !ier ) {
MMG5_Free_all(mesh,met);
return(MMG5_STRONGFAILURE);
}
if ( !MMG5_parsop(mesh,met) )
RETURN_AND_FREE(mesh,met,MMG5_LOWFAILURE);
chrono(OFF,&ctim[1]);
printim(ctim[1].gdif,stim);
fprintf(stdout," -- DATA READING COMPLETED. %s\n",stim);
ier = MMG5_mmg3dlib(mesh,met );
if ( ier != MMG5_STRONGFAILURE ) {
chrono(ON,&ctim[1]);
if ( mesh->info.imprim )
fprintf(stdout,"\n -- WRITING DATA FILE %s\n",mesh->nameout);
if ( !MMG5_saveMesh(mesh) ) {
MMG5_Free_all(mesh,met );
return(EXIT_FAILURE);
}
if ( !MMG5_saveMet(mesh,met) ) {
MMG5_Free_all(mesh,met);
return(EXIT_FAILURE);
}
chrono(OFF,&ctim[1]);
if ( mesh->info.imprim )
fprintf(stdout," -- WRITING COMPLETED\n");
}
/* free mem */
chrono(OFF,&ctim[0]);
printim(ctim[0].gdif,stim);
fprintf(stdout,"\n MMG3D: ELAPSED TIME %s\n",stim);
MMG5_Free_all(mesh,met);
return(ier);
}