Skip to content
Snippets Groups Projects
Commit 1cc4879b authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

simple mesh read example

parent 2a1271e5
No related branches found
No related tags found
No related merge requests found
......@@ -53,6 +53,9 @@ target_link_libraries(mainPost shared)
add_executable(mainSimple mainSimple.cpp)
target_link_libraries(mainSimple shared)
add_executable(mainReadMesh mainReadMesh.cpp)
target_link_libraries(mainReadMesh shared)
add_executable(mainGeoFactory mainGeoFactory.cpp)
target_link_libraries(mainGeoFactory shared)
......
#include <stdio.h>
#include "Gmsh.h"
#include "GModel.h"
#include "MTetrahedron.h"
#include "MTriangle.h"
int main(int argc, char **argv)
{
GmshInitialize(argc, argv);
GModel *m = new GModel();
m->readMSH("../problem.msh");
for(GModel::riter it = m->firstRegion(); it != m->lastRegion(); ++it){
GRegion *r = *it;
for(unsigned int i = 0; i < r->tetrahedra.size(); i++)
printf(" tet %d: vertices %d %d %d %d\n", r->tetrahedra[i]->getNum(),
r->tetrahedra[i]->getVertex(0)->getNum(),
r->tetrahedra[i]->getVertex(1)->getNum(),
r->tetrahedra[i]->getVertex(2)->getNum(),
r->tetrahedra[i]->getVertex(3)->getNum());
}
for(GModel::fiter it = m->firstFace(); it != m->lastFace(); ++it){
GFace *f = *it;
for(unsigned int i = 0; i < f->triangles.size(); i++)
printf(" tri %d: vertices %d %d %d\n", f->triangles[i]->getNum(),
f->triangles[i]->getVertex(0)->getNum(),
f->triangles[i]->getVertex(1)->getNum(),
f->triangles[i]->getVertex(2)->getNum());
}
delete m;
GmshFinalize();
}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment