From: robert <rob...@un...> - 2012-05-24 14:01:54
|
Hi, thanks for the fast reply. Here the files - it's just a simple cube fun.ele: 6 4 0 1 6 7 1 5 2 1 8 3 4 3 3 7 1 2 4 1 7 6 2 5 1 7 3 8 6 7 8 1 5 # Generated by tetgen -pg fun.poly fun.node: 8 3 0 0 1 0 0 0 2 1 0 0 3 1 1 0 4 0 1 0 5 0 0 1 6 1 0 1 7 1 1 1 8 0 1 1 # Generated by tetgen -pg fun.poly fun.face (I don't think you need this one but anyway): 12 0 1 1 2 3 2 6 5 7 3 3 4 1 4 6 2 1 5 7 5 8 6 3 2 7 7 1 5 6 8 8 4 3 9 7 2 6 10 8 5 1 11 7 8 3 12 1 4 8 # Generated by tetgen -pg fun.poly I'll look at the debugging soon ... Robert On 05/24/2012 03:54 PM, John Peterson wrote: > On Thu, May 24, 2012 at 7:27 AM, robert<rob...@un...> wrote: >> Hello, >> >> has anybody tried to use tetgen_io with one of the recent versions of >> ubuntu (11.10, 12.04). >> >> I just tried to read in a tetgen mesh >> >> >> >> #include<iostream> >> #include "libmesh.h" >> #include "mesh.h" >> #include "mesh_generation.h" >> #include<mesh_data.h> >> >> #include<tetgen_io.h> >> #include<mesh_tetgen_wrapper.h> >> #include<mesh_tetgen_interface.h> >> >> using namespace libMesh; >> >> int main (int argc, char** argv) >> { >> LibMeshInit init (argc, argv); >> >> >> Mesh mesh(3); >> >> MeshData meshdata(mesh); >> meshdata.activate(); >> TetGenIO TG(mesh); >> >> mesh.read("fun.1.ele"); >> >> return 0; >> } >> >> >> >> but get the following message: >> >> / >> TetGenIO found the tetgen files to read >> Segmentation fault (core dumped)/ >> >> >> >> if I start it in gdb I get: >> / >> >> >> Program received signal SIGSEGV, Segmentation fault. >> 0xb7d2040c in libMesh::SerialMesh::insert_elem(libMesh::Elem*) () >> from /home/robert/software/libmesh/lib/i686-pc-linux-gnu_opt/libmesh.so/ >> >> >> I configured libmesh only with >> ./configure >> make >> >> so, I did not specify any options. >> >> Last time when I used tetgen (same code) was about 1 year ago with an >> older version of ubuntu 10.?? and an older version of libmesh. > Can you try compiling and running in debug mode? I'd like to see if > you hit the assert in insert_elem() before the segfault occurs. > > This line of code changed in r4787 (Aug 2011) so that fits with the > last time you attempted to use it. > > Can you also send me (or link me, if the file is large) the fun.1.ele mesh file? > > Looks like a bug in the TetGenIO reader... > |