Hi all,

The follow-up on this has been quite slow due to a server room upgrade but I have managed to rerun the problem again with the latest version of Netgen thanks to Christophe. These runs take three or four days so are quite slow especially with slow cpu's.

The memory error seems to have re-occurred in a different section of code. I attach the gdb output for the following command.

gmsh box3d.geo -3 -bin -algo front3d -optimize_netgen -smooth 5 -clscale 1.0 -rand 1.e-8 -o box_c_gmsh_svn.msh

Additionally I have also run with Delaunay instead of frontal to see if this is different. The output is below as well for this command.

gmsh box3d.geo -3 -bin -optimize_netgen -smooth 5 -clscale 1.0 -rand 1.e-8 -o box_c_gmsh_svn.msh

It seems both methods of meshing have different errors when meshing large meshes. It seems from the gdb output the error occurs when allocating more memory. However I also monitored the memory usage which did not exceed 27% of the server unit therefore I imagine something tried to allocate more than the system has in one go, possibly due to faulty allocation request . I find the value s=18446744072204065817 being passed to DynamicMem a little suspicious.

Could somebody suggest the best way forward as we are very keen to overcome this problem.

Thanks

Andre

On 20/12/2011 20:36, Christophe Geuzaine wrote:
On 20 Dec 2011, at 09:40, Andre Nicolle wrote:

Hi Christophe

Thanks for the update, have recompiled and awaiting the output. Will update 
shortly.

It would be great to hear of anyone who has achieved meshes larger than 80 
million elements. We are currently looking to find a limit for some very large 
CFD calculations.

PS : have you tried the Delaunay algorithm?

Regards

Andre

On 20/12/2011 00:32, Christophe Geuzaine wrote:
Hi Andre,

Can you try recompiling? I have updated the version of Netgen used in Gmsh to 
the latest SVN.

Thanks,

Christophe


On 16 Dec 2011, at 11:33, Andre Nicolle wrote:

Hi all,

I have been using Gmsh for many years and found it to be a great program. We 
are now expanding our calculations and now require bigger meshes but have 
started to find this to be a problem when meshing (using netgen) larger than 
100 million elements. Our system has adequate RAM (256Gb) but gmsh exits with a 
segmentation fault, as shown below. At present we are using Ubuntu 10.04 64Bit 
AMD 6100. The problem can be found with 2.42, 2.5 and svn revision 10744 
versions of the code. Interestingly for some reason the latest version of Gmsh 
crashes considerably sooner than the older versions.

Gmsh 2.50 and repository version were compiled using gcc version 4.4.3, Atlas 
Blas with the following command:
cmake ../CMakeLists.txt -DENABLE_FLTK=0

2.42 was from ubuntu repository

The meshing command is:

gmsh box3d.geo -3 -bin -algo front3d -optimize_netgen -smooth 5 -clscale 1.0 
-rand 1.e-8 -o box_c_new.msh

The box3d.geo file is very simple box with resolutions (0.002) set on the 
points as can be seen below .

To help find the problem I have tried to recompile Gmsh with debugging on. I'm 
not sure which flags got set.

cmake ../CMakeLists.txt -DENABLE_FLTK=0 -DCMAKE_BUILD_TYPE=Debug

and rerun Gmsh using gdb, output below:

The error seems to be generated in Netgen moveablemem.cpp which has disappeared 
in the latest version of Netgen 4.9 anyway. Given that Gmsh is still using 
Netgen 4.5, is it possible to link to the later version of netgen library 
installed locally instead of the built-in version to see if this helps.

Thanks

Andre

///////////////////////////////////////////////////////////////////////////////////////////////////////////////

RES = 0.002;
Point(1) = {0, 0, 0, RES};
Point(2) = {1, 0, 0, RES};
Point(3) = {1, 0, 1, RES};
Point(4) = {0, 0, 1, RES};
Point(5) = {0, 1, 0, RES};
Point(6) = {1, 1, 0, RES};
Point(7) = {1, 1, 1, RES};
Point(8) = {0, 1, 1, RES};
Line (1) = {1, 2};
Line (2) = {2, 3};
Line (3) = {3, 4};
Line (4) = {4, 1};
Line (5) = {1, 5};
Line (6) = {5, 6};
Line (7) = {6, 2};
Line (8) = {5, 8};
Line (9) = {8, 4};
Line (10) = {8, 7};
Line (11) = {7, 6};
Line (12) = {7, 3};
Line Loop (15) = {1, 2, 3, 4};
Plane Surface (16) = {15};
Line Loop (17) = {12, -2, -7, -11};
Plane Surface (18) = {17};
Line Loop (19) = {3, -9, 10, 12};
Plane Surface (20) = {19};
Line Loop (21) = {10, 11, -6, 8};
Plane Surface (22) = {21};
Line Loop (23) = {8, 9, 4, 5};
Plane Surface (24) = {23};
Line Loop (25) = {6, 7, -1, 5};
Plane Surface (26) = {25};
Surface Loop(27) = {20, 16, 26, 22, 18, 24};
Volume(28) = {27};

//////////////////////////////////////////////////////////////////////////////////////////////////////////////////

Info    : Started on Thu Dec 15 20:45:53 2011
Info    : Reading 'box3d.geo'...
Info    : Done reading 'box3d.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1 (Line)
Info    : Meshing curve 2 (Line)
Info    : Meshing curve 3 (Line)
Info    : Meshing curve 4 (Line)
Info    : Meshing curve 5 (Line)
Info    : Meshing curve 6 (Line)
Info    : Meshing curve 7 (Line)
Info    : Meshing curve 8 (Line)
Info    : Meshing curve 9 (Line)
Info    : Meshing curve 10 (Line)
Info    : Meshing curve 11 (Line)
Info    : Meshing curve 12 (Line)
Info    : Done meshing 1D (0.01 s)
Info    : Meshing 2D...
Info    : Meshing surface 16 (Plane, Delaunay)
Info    : Meshing surface 18 (Plane, Delaunay)
Info    : Meshing surface 20 (Plane, Delaunay)
Info    : Meshing surface 22 (Plane, Delaunay)
Info    : Meshing surface 24 (Plane, Delaunay)
Info    : Meshing surface 26 (Plane, Delaunay)
Info    : Done meshing 2D (603.87 s)
Info    : Meshing 3D...
Info    : Meshing volume 28 (Frontal)
Info    : Region 28 Face 20, 0 intersect
Info    : Region 28 Face 16, 1 intersect
Info    : Region 28 Face 26, 1 intersect
Info    : Region 28 Face 22, 1 intersect
Info    : Region 28 Face 18, 1 intersect
Info    : Region 28 Face 24, 1 intersect
Info    : CalcLocalH: 1984450 Points 0 Elements 3968896 Surface Elements
Info    : Check subdomain 1 / 1
Info    : 3968896 open elements
Info    : Meshing subdomain 1 of 1
Info    : 3968896 open elements
Info    : Use internal rules
Info    : 3968896 open elements
Info    : Delaunay meshing
Info    : number of points: 1984450
Info    : blockfill local h
Info    : number of points: 44035686

Program received signal SIGSEGV, Segmentation fault.
0x00007ffff577b1c3 in ?? () from /lib/libc.so.6
(gdb) where
#0  0x00007ffff577b1c3 in ?? () from /lib/libc.so.6
#1  0x00007ffff577ee83 in free () from /lib/libc.so.6
#2  0x0000000000f9c2cf in netgen::BaseMoveableMem::ReAlloc (
    this=0x7fffffffb3b8, s=12612)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/general/moveablemem.cpp:159
#3  0x0000000000f950e7 in netgen::MoveableMem<netgen::INDEX_3>::ReAlloc (
    this=0x7fffffffb3b8, s=1051)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/include/../general/moveablemem.hpp:63
#4  0x0000000000f94f2f in netgen::MoveableArray<netgen::INDEX_3, 0>::SetSize (
    this=0x7fffffffb3b0, nsize=1051)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/include/../general/array.hpp:447
#5  0x0000000000f945c2 in netgen::BASE_INDEX_3_CLOSED_HASHTABLE::BaseSetSize (
    this=0x7fffffffb3b0, size=1051)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/general/hashtabl.cpp:294
#6  0x000000000105a6ba in netgen::INDEX_3_CLOSED_HASHTABLE<int>::SetSize (
    this=0x7fffffffb3b0, size=1051)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/include/../general/hashtabl.hpp:870
#7  0x000000000105a4ae in netgen::MeshNB::ResetFaceHT (this=0x7fffffffb3b0,
    size=1051)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/delaunay.cpp:144
#8  0x0000000001053a54 in netgen::AddDelaunayPoint (newpi=..., newp=...,
    tempels=..., mesh=..., tettree=..., meshnb=..., centers=..., radi2=...,
    connected=..., treesearch=..., freelist=..., list=..., insphere=...,
    closesphere=...)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/delaunay.cpp:496
#9  0x0000000001055625 in netgen::Delaunay1 (mesh=..., mp=...,
    adfront=0x205a4d90, tempels=..., oldnp=1984450, startel=..., pmin=...,
    pmax=...)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/delaunay.cpp:765
#10 0x0000000001055d9f in netgen::Meshing3::Delaunay (this=0x7fffffffd3f0,
    mesh=..., domainnr=1, mp=...)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/delaunay.cpp:816
#11 0x000000000106555b in netgen::MeshVolume (mp=..., mesh3d=...)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/meshfunc.cpp:192
#12 0x0000000000f92263 in nglib::NgAddOn_GenerateVolumeMesh (mesh=0x2078cdc0,
    maxh=1e+22) at /home/andre/gmsh/svn/gmsh/contrib/Netgen/nglib_addon.cpp:80
#13 0x0000000000b2a33e in meshGRegion::operator() (this=0x7fffffffe240,
    gr=0x1844820) at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegion.cpp:883
#14 0x0000000000a79d90 in std::for_each<std::_Rb_tree_const_iterator<GRegion*>,
    at /usr/include/c++/4.4/bits/stl_algo.h:4200
#15 0x0000000000a7768c in Mesh3D (m=0x1833a20)
    at /home/andre/gmsh/svn/gmsh/Mesh/Generator.cpp:533
#16 0x0000000000a77c0d in GenerateMesh (m=0x1833a20, ask=3)
    at /home/andre/gmsh/svn/gmsh/Mesh/Generator.cpp:636
#17 0x000000000093e1b4 in GModel::mesh (this=0x1833a20, dimension=3)
    at /home/andre/gmsh/svn/gmsh/Geo/GModel.cpp:510
#18 0x00000000007c4c75 in GmshBatch ()
    at /home/andre/gmsh/svn/gmsh/Common/Gmsh.cpp:170
#19 0x00000000007c42e3 in main (argc=15, argv=0x7fffffffe5d8)
    at /home/andre/gmsh/svn/gmsh/Common/Main.cpp:27

_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh

This GDB was configured as "x86_64-linux-gnu".
For bug reporting instructions, please see:
<http://www.gnu.org/software/gdb/bugs/>...
Reading symbols from /home/andre/gmsh_svn...done.

(gdb) box3d.geo -3 -bin -algo front3d -optimize_netgen -smooth 5 -clscale 1.0 
-rand 1.e-8 -o box_c_gmsh_svn_vvvvbig.msh
Undefined command: "box3d".  Try "help".
(gdb) run box3d.geo -3 -bin -algo front3d -optimize_netgen -smooth 5 -clscale 
1.0 -rand 1.e-8 -o box_c_gmsh_svn_vvvvbig.msh
Starting program: /home/andre/gmsh_svn box3d.geo -3 -bin -algo front3d 
-optimize_netgen -smooth 5 -clscale 1.0 -rand 1.e-8 -o 
box_c_gmsh_svn_vvvvbig.msh
Info    : Running '/home/andre/gmsh_svn box3d.geo -3 -bin -algo front3d 
-optimize_netgen -smooth 5 -clscale 1.0 -rand 1.e-8 -o 
box_c_gmsh_svn_vvvvbig.msh' [1 node(s), max. 1 thread(s)]
Info    : Started on Thu Jan  5 02:10:56 2012
Info    : Reading 'box3d.geo'...
Info    : Done reading 'box3d.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1 (Line)
Info    : Meshing curve 2 (Line)
Info    : Meshing curve 3 (Line)
Info    : Meshing curve 4 (Line)
Info    : Meshing curve 5 (Line)
Info    : Meshing curve 6 (Line)
Info    : Meshing curve 7 (Line)
Info    : Meshing curve 8 (Line)
Info    : Meshing curve 9 (Line)
Info    : Meshing curve 10 (Line)
Info    : Meshing curve 11 (Line)
Info    : Meshing curve 12 (Line)
Info    : Done meshing 1D (0.01 s)
Info    : Meshing 2D...
Info    : Meshing surface 16 (Plane, Delaunay)
Info    : Meshing surface 18 (Plane, Delaunay)
Info    : Meshing surface 20 (Plane, Delaunay)
Info    : Meshing surface 22 (Plane, Delaunay)
Info    : Meshing surface 24 (Plane, Delaunay)
Info    : Meshing surface 26 (Plane, Delaunay)
Info    : Done meshing 2D (614.45 s)
Info    : Meshing 3D...
Info    : Meshing volume 28 (Frontal)
Info    : Region 28 Face 20, 1 intersect
Info    : Region 28 Face 16, 0 intersect
Info    : Region 28 Face 26, 0 intersect
Info    : Region 28 Face 22, 1 intersect
Info    : Region 28 Face 18, 1 intersect
Info    : Region 28 Face 24, 0 intersect
Info    : CalcLocalH: 1984774 Points 0 Elements 3969544 Surface Elements 
Info    : Check subdomain 1 / 1 
Info    : 3969544 open elements 
Info    : Meshing subdomain 1 of 1 
Info    : Meshing subdomain 1 
Info    : 3969544 open elements 
Info    : Use internal rules 
Info    : 3969544 open elements 
Info    : Delaunay meshing 
Info    : number of points: 1984774 
Info    : blockfill local h 
Info    : rel filldist = 0.1 
Info    : blockfill local h 
Info    : number of points: 36356448 
Info    : Points: 36356448 
Info    : Elements: 232456791 
Info    : Points:           36356452 
Info    : Volume Elements:  232456791 
Info    : Volume elements in qualityclasses: 
Info    : 0 - 0.05: 3990345 
Info    : 0.05 -  0.1: 43684 
Info    : 0.1 - 0.15: 563984 
Info    : 0.15 -  0.2: 302279 
Info    : 0.2 - 0.25: 2621707 
Info    : 0.25 -  0.3: 2888946 
Info    : 0.3 - 0.35: 13035092 
Info    : 0.35 -  0.4: 4390813 
Info    : 0.4 - 0.45: 7709941 
Info    : 0.45 -  0.5: 5499547 
Info    : 0.5 - 0.55: 21729097 
Info    : 0.55 -  0.6: 27700612 
Info    : 0.6 - 0.65: 24520344 
Info    : 0.65 -  0.7: 9046535 
Info    : 0.7 - 0.75: 53367882 
Info    : 0.75 -  0.8: 16343034 
Info    : 0.8 - 0.85: 26693979 
Info    : 0.85 -  0.9: 4493082 
Info    : 0.9 - 0.95: 4695785 
Info    : 0.95 -    1: 2820103 
Info    : total error: 4e+16 
Info    : 1844 open elements 
Info    : Num open: 1844 
Info    : free: 23619, fixed: 232433172 
Info    : free: 23619, fixed: 232433172 
Info    : SwapImprove  
Info    : Start SwapImprove 
terminate called after throwing an instance of 'std::bad_alloc'
  what():  std::bad_alloc

Program received signal SIGABRT, Aborted.
0x00007ffff5734a75 in raise () from /lib/libc.so.6
(gdb) 

(gdb) where
#0  0x00007ffff5734a75 in raise () from /lib/libc.so.6
#1  0x00007ffff57385c0 in abort () from /lib/libc.so.6
#2  0x00007ffff5fea8e5 in __gnu_cxx::__verbose_terminate_handler() ()
   from /usr/lib/libstdc++.so.6
#3  0x00007ffff5fe8d16 in ?? () from /usr/lib/libstdc++.so.6
#4  0x00007ffff5fe8d43 in std::terminate() () from /usr/lib/libstdc++.so.6
#5  0x00007ffff5fe8e3e in __cxa_throw () from /usr/lib/libstdc++.so.6
#6  0x00007ffff5fe92cd in operator new(unsigned long) ()
   from /usr/lib/libstdc++.so.6
#7  0x00007ffff5fe9389 in operator new[](unsigned long) ()
   from /usr/lib/libstdc++.so.6
#8  0x0000000000fa5be4 in netgen::BaseDynamicMem::Alloc (this=0x7fffffffafc8, 
    s=18446744061665665224)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/general/
#9  0x0000000000f9f03b in netgen:: <netgen::INDEX_2>::Alloc (
    this=0x7fffffffafc8, s=18446744072204065817)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/include/../general/dynamicmem.hpp:56
#10 0x0000000000f9ef54 in DynamicMem (this=0x7fffffffafc8, 
    s=18446744072204065817)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/include/../general/dynamicmem.hpp:52
#11 0x0000000000f9ed06 in MoveableArray (this=0x7fffffffafc0, 
    asize=-1505485799)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/include/../general/array.hpp:426
#12 0x0000000000f9dfa9 in BASE_INDEX_2_CLOSED_HASHTABLE (this=0x7fffffffafc0, 
    size=-1505485799)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/general/hashtabl.cpp:217
#13 0x0000000000fdf855 in INDEX_2_CLOSED_HASHTABLE (this=0x7fffffffafc0, 
    size=-1505485799)
    at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/include/../general/hashtabl.hpp:1169
#14 0x00000000010317e4 in netgen::MeshOptimize3d::SwapImprove (
    this=0x7fffffffbcd0, mesh=..., goal=netgen::OPT_CONFORM, 
    working_elements=0x0)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/improve3.cpp:639
#15 0x000000000106620d in netgen::Meshing3::Delaunay (this=0x7fffffffd090, 
    mesh=..., domainnr=1, mp=...)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/delaunay.cpp:895
#16 0x0000000001076af0 in netgen::MeshVolume (mp=..., mesh3d=...)
    at /home/andre/gmsh/svn/gmsh/contrib/Netgen/libsrc/meshing/meshfunc.cpp:187
#17 0x0000000000f9b4b6 in nglib::Ng_GenerateVolumeMesh (mesh=0x207b2310, 
        maxh=1e+22) at 
/home/andre/gmsh/svn/gmsh/contrib/Netgen/nglib_gmsh.cpp:215
#18 0x0000000000b33b1e in meshGRegion::operator() (this=0x7fffffffdef0, 
    gr=0x185e950) at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegion.cpp:882
#19 0x0000000000a823d2 in std::for_each<std::_Rb_tree_const_iterator<GRegion*>, 
meshGRegion> (__first=..., __last=..., __f=...)
    at /usr/include/c++/4.4/bits/stl_algo.h:4200
#20 0x0000000000a7fcd8 in Mesh3D (m=0x184db40)
    at /home/andre/gmsh/svn/gmsh/Mesh/Generator.cpp:533
#21 0x0000000000a80259 in GenerateMesh (m=0x184db40, ask=3)
    at /home/andre/gmsh/svn/gmsh/Mesh/Generator.cpp:636
#22 0x0000000000940f60 in GModel::mesh (this=0x184db40, dimension=3)
    at /home/andre/gmsh/svn/gmsh/Geo/GModel.cpp:513
#23 0x00000000007c7985 in GmshBatch ()
    at /home/andre/gmsh/svn/gmsh/Common/Gmsh.cpp:170
#24 0x00000000007c6ff3 in main (argc=15, argv=0x7fffffffe288)
    at /home/andre/gmsh/svn/gmsh/Common/Main.cpp:27

andre@andre-desktop:~$ gdb ./gmsh_svn 
GNU gdb (GDB) 7.1-ubuntu
Copyright (C) 2010 Free Software Foundation, Inc.
License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>
This is free software: you are free to change and redistribute it.
There is NO WARRANTY, to the extent permitted by law.  Type "show copying"
and "show warranty" for details.
This GDB was configured as "x86_64-linux-gnu".
For bug reporting instructions, please see:
<http://www.gnu.org/software/gdb/bugs/>...
Reading symbols from /home/andre/gmsh_svn...done.
(gdb) run box3d.geo -3 -bin -optimize_netgen -smooth 5 -clscale 1.0 -rand 1.e-8 
-o box_c_gmsh_svn_vvvvbig.msh
Starting program: /home/andre/gmsh_svn box3d.geo -3 -bin -optimize_netgen 
-smooth 5 -clscale 1.0 -rand 1.e-8 -o box_c_gmsh_svn_vvvvbig.msh
Info    : Running '/home/andre/gmsh_svn box3d.geo -3 -bin -optimize_netgen 
-smooth 5 -clscale 1.0 -rand 1.e-8 -o box_c_gmsh_svn_vvvvbig.msh' [1 node(s), 
max. 1 thread(s)]
Info    : Started on Mon Jan  9 18:44:03 2012
Info    : Reading 'box3d.geo'...
Info    : Done reading 'box3d.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1 (Line)
Info    : Meshing curve 2 (Line)
Info    : Meshing curve 3 (Line)
Info    : Meshing curve 4 (Line)
Info    : Meshing curve 5 (Line)
Info    : Meshing curve 6 (Line)
Info    : Meshing curve 7 (Line)
Info    : Meshing curve 8 (Line)
Info    : Meshing curve 9 (Line)
Info    : Meshing curve 10 (Line)
Info    : Meshing curve 11 (Line)
Info    : Meshing curve 12 (Line)
Info    : Done meshing 1D (0.01 s)
Info    : Meshing 2D...
Info    : Meshing surface 16 (Plane, Delaunay)
Info    : Meshing surface 18 (Plane, Delaunay)
Info    : Meshing surface 20 (Plane, Delaunay)
Info    : Meshing surface 22 (Plane, Delaunay)
Info    : Meshing surface 24 (Plane, Delaunay)
Info    : Meshing surface 26 (Plane, Delaunay)
Info    : Done meshing 2D (606.28 s)
Info    : Meshing 3D...
Info    : Meshing volume 28 (Delaunay)
Constructing Delaunay tetrahedralization.
Delaunay seconds:  586.04
Creating surface mesh.
Perturbing vertices.
Delaunizing segments.
Constraining facets.
Segment and facet seconds:  304.07
Removing unwanted tetrahedra.
Hole seconds:  8.16
Repairing mesh.
Repair seconds:  4.63

Writing nodes.
Writing elements.
Writing subfaces.
Writing edges.

Output seconds:  7.37
Total running seconds:  911.21

Statistics:

  Input points: 1984774
  Input facets: 3969544
  Input segments: 5954316
  Input holes: 0
  Input regions: 0

  Mesh points: 2695849
  Mesh tetrahedra: 8667265
  Mesh triangles: 20030377
  Mesh subfaces: 5391694
  Mesh subsegments: 6000

Info    : 2695849 points 6000 edges and 5391694 faces in the final mesh

Program received signal SIGBUS, Bus error.
0x00007ffff577bdf0 in ?? () from /lib/libc.so.6

(gdb) where
#0  0x00007ffff577bdf0 in ?? () from /lib/libc.so.6
#1  0x00007ffff577d58e in malloc () from /lib/libc.so.6
#2  0x00007ffff5fe926d in operator new(unsigned long) () from 
/usr/lib/libstdc++.so.6
#3  0x0000000000b46792 in __gnu_cxx::new_allocator<std::_List_node<MTet4*> 
>::allocate(unsigned long, void const*) ()
#4  0x0000000000b4618a in std::_List_base<MTet4*, std::allocator<MTet4*> 
>::_M_get_node() ()
#5  0x0000000000b44f25 in std::list<MTet4*, std::allocator<MTet4*> 
>::_M_create_node (
    this=0x7fffffffd8f0, __x=@0x7ffff80001a8) at 
/usr/include/c++/4.4/bits/stl_list.h:461
#6  0x0000000000b43c4b in std::list<MTet4*, std::allocator<MTet4*> >::_M_insert 
(
    this=0x7fffffffd8f0, __position=..., __x=@0x7ffff80001a8)
    at /usr/include/c++/4.4/bits/stl_list.h:1407
#7  0x0000000000b42ce4 in std::list<MTet4*, std::allocator<MTet4*> >::push_back 
(
    this=0x7fffffffd8f0, __x=@0x7ffff80001a8) at 
/usr/include/c++/4.4/bits/stl_list.h:920
#8  0x0000000000b3d015 in recur_classify (t=0x71271780, theRegion=..., 
faces_bound=..., 
    bidon=0x7b, model=0x184db40, search=...)
    at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegionDelaunayInsertion.cpp:356
#9  0x0000000000b3d20b in recur_classify (t=0x5f8b19c0, theRegion=..., 
faces_bound=..., 
    bidon=0x7b, model=0x184db40, search=...)
---Type <return> to continue, or q <return> to quit---
    at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegionDelaunayInsertion.cpp:386
#10 0x0000000000b3d20b in recur_classify (t=0x600cb700, theRegion=..., 
faces_bound=..., 
    bidon=0x7b, model=0x184db40, search=...)
    at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegionDelaunayInsertion.cpp:386
#11 0x0000000000b3d20b in recur_classify (t=0x2dddd280, theRegion=..., 
faces_bound=..., 
    bidon=0x7b, model=0x184db40, search=...)
    at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegionDelaunayInsertion.cpp:386
#12 0x0000000000b3d20b in recur_classify (t=0x68ad9880, theRegion=..., 
faces_bound=..., 
    bidon=0x7b, model=0x184db40, search=...)
    at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegionDelaunayInsertion.cpp:386
#13 0x0000000000b3d20b in recur_classify (t=0x68ad9680, theRegion=..., 
faces_bound=..., 
    bidon=0x7b, model=0x184db40, search=...)
    at /home/andre/gmsh/svn/gmsh/Mesh/meshGRegionDelaunayInsertion.cpp:386
#14 0x0000000000b3d20b in recur_classify (t=0x6399c200, theRegion=..., 
faces_bound=..., 
    bidon=0x7b, model=0x184db40, search=...)
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh

Reply via email to