[Gmsh] Tolerance option for the meshing module?

Jö Fahlke jorrit at jorrit.de
Mon Dec 6 16:05:57 CET 2010


Package: gmsh
Version: 2.4.2.dfsg-5
Severity: normal

Hi!

I'm trying to mesh the attached .geo file.  It consists of an inner sphere of
radius 50e-9 and an outer spherical shell around it of radius 500e-9.  The
code to generate the spheres is taken from the tutorial file t6.geo.

When meshing on the command line (log attached) gmsh segfaults when writing
the mesh file and the resulting mesh file has a size of 0 bytes.  Meshing from
the GUI does not result in a segfault, however my grid manager (ALUGrid) barfs
on the resulting .msh file because of some inconsistencies.

When I rescale the spheres to 50e9 and 500e9 respectively, meshing works fine.
So my question is: is there a ways to specify a tolerance for meshing (similar
to the Geometry.Tolerance option)?  Setting Geometry.Tolerance did not seem to
have any effect on meshing.

I guess I can work around the problem by meshing a normal-sized geometry and
then rescale all coordinates in the .msh file using a perl-script or similar.
But it would certainly be nicer to produce the desired mesh file using gmsh
directly.

Bye,
Jö.

-- System Information:
Debian Release: squeeze/sid
  APT prefers testing-proposed-updates
  APT policy: (500, 'testing-proposed-updates'), (500, 'testing')
Architecture: i386 (x86_64)

Kernel: Linux 2.6.32-5-amd64 (SMP w/2 CPU cores)
Locale: LANG=de_DE.UTF-8, LC_CTYPE=de_DE.UTF-8 (charmap=UTF-8)
Shell: /bin/sh linked to /bin/dash

Versions of packages gmsh depends on:
ii  dpkg                    1.15.8.5         Debian package management system
ii  install-info            4.13a.dfsg.1-6   Manage installed documentation in 
ii  libatlas3gf-base [libla 3.8.3-27         Automatically Tuned Linear Algebra
ii  libblas3gf [libblas.so. 1.2-8            Basic Linear Algebra Reference imp
ii  libc6                   2.11.2-7         Embedded GNU C Library: Shared lib
ii  libcgns2                2.5.4-3          CFD General Notation System librar
ii  libfltk1.1              1.1.10-2+b1      Fast Light Toolkit - shared librar
ii  libgcc1                 1:4.4.5-8        GCC support library
ii  libgfortran3            4.4.5-8          Runtime library for GNU Fortran ap
ii  libgl1-mesa-glx [libgl1 7.7.1-4          A free implementation of the OpenG
ii  libglu1-mesa [libglu1]  7.7.1-4          The OpenGL utility library (GLU)
ii  libjpeg62               6b1-1            The Independent JPEG Group's JPEG 
ii  liblapack3gf [liblapack 3.2.1-8          library of linear algebra routines
ii  libopencascade-foundati 6.3.0.dfsg.1-6   OpenCASCADE CAE platform shared li
ii  libopencascade-modeling 6.3.0.dfsg.1-6   OpenCASCADE CAE platform shared li
ii  libpng12-0              1.2.44-1         PNG library - runtime
ii  libstdc++6              4.4.5-8          The GNU Standard C++ Library v3
ii  libx11-6                2:1.3.3-4        X11 client-side library
ii  zlib1g                  1:1.2.3.4.dfsg-3 compression library - runtime

Versions of packages gmsh recommends:
ii  mpi-default-bin               0.6        Standard MPI runtime programs

gmsh suggests no packages.

-- no debconf information

-- 
Interpunktion, Orthographie und Grammatik der Email ist frei erfunden.
Eine Übereinstimmung mit aktuellen oder ehemaligen Regeln wäre rein
zufällig und ist nicht beabsichtigt.
-------------- next part --------------
// radius of scactterer
inner_r = 50e-9;
// radius of domain
outer_r = 10.0 * inner_r;
// edge length
el = inner_r / 2.0;

Geometry.Tolerance = inner_r * 1e-6;

x = 0;
y = 0;
z = 0;

Function CheeseHole 

  // In the following commands we use the reserved variable name
  // `newp', which automatically selects a new point number.  This
  // number is chosen as the highest current point number, plus
  // one.  (Note that, analogously to `newp', the variables `newc',
  // `news', `newv' and `newreg' select the highest number amongst
  // currently defined curves, surfaces, volumes and `any entities
  // other than points', respectively.)

  p1 = newp; Point(p1) = {x,  y,  z  } ;
  p2 = newp; Point(p2) = {x+r,y,  z  } ;
  p3 = newp; Point(p3) = {x,  y+r,z  } ;
  p4 = newp; Point(p4) = {x,  y,  z+r} ;
  p5 = newp; Point(p5) = {x-r,y,  z  } ;
  p6 = newp; Point(p6) = {x,  y-r,z  } ;
  p7 = newp; Point(p7) = {x,  y,  z-r} ;

  c1 = newreg; Circle(c1) = {p2,p1,p7};
  c2 = newreg; Circle(c2) = {p7,p1,p5};
  c3 = newreg; Circle(c3) = {p5,p1,p4};
  c4 = newreg; Circle(c4) = {p4,p1,p2};
  c5 = newreg; Circle(c5) = {p2,p1,p3};
  c6 = newreg; Circle(c6) = {p3,p1,p5};
  c7 = newreg; Circle(c7) = {p5,p1,p6};
  c8 = newreg; Circle(c8) = {p6,p1,p2};
  c9 = newreg; Circle(c9) = {p7,p1,p3};
  c10 = newreg; Circle(c10) = {p3,p1,p4};
  c11 = newreg; Circle(c11) = {p4,p1,p6};
  c12 = newreg; Circle(c12) = {p6,p1,p7};

  // We need non-plane surfaces to define the spherical holes. Here we
  // use ruled surfaces, which can have 3 or 4 sides:

  l1 = newreg; Line Loop(l1) = {c5,c10,c4};   Ruled Surface(newreg) = {l1};
  l2 = newreg; Line Loop(l2) = {c9,-c5,c1};   Ruled Surface(newreg) = {l2};
  l3 = newreg; Line Loop(l3) = {c12,-c8,-c1}; Ruled Surface(newreg) = {l3};
  l4 = newreg; Line Loop(l4) = {c8,-c4,c11};  Ruled Surface(newreg) = {l4};
  l5 = newreg; Line Loop(l5) = {-c10,c6,c3};  Ruled Surface(newreg) = {l5};
  l6 = newreg; Line Loop(l6) = {-c11,-c3,c7}; Ruled Surface(newreg) = {l6};
  l7 = newreg; Line Loop(l7) = {-c2,-c7,-c12};Ruled Surface(newreg) = {l7};
  l8 = newreg; Line Loop(l8) = {-c6,-c9,c2};  Ruled Surface(newreg) = {l8};

  // We then store the surface loops identification numbers in a list
  // for later reference (we will need these to define the final
  // volume):

  theloop = newreg ; 

  Surface Loop(theloop) = {l8+1,l5+1,l1+1,l2+1,l3+1,l7+1,l6+1,l4+1};

Return

r = inner_r ;
Call CheeseHole ;
inner_loop = theloop;

r = outer_r ;
Call CheeseHole ;
outer_loop = theloop;

inner_vol = newreg;
Volume(inner_vol) = { inner_loop };

outer_vol = newreg;
Volume(outer_vol) = { outer_loop, inner_loop };

Physical Volume(1) = { inner_vol };
Physical Volume(2) = { outer_vol };


Field[1] = Box;
Field[1].VIn = el;
Field[1].VOut = el;
Background Field = 1;

-------------- next part --------------
-*- mode: compilation; default-directory: "/home/joe/src/dune-opt/dune-pm/appl/fetd/" -*-
Compilation started at Mon Dec  6 16:02:50

export LANG=C; cd ~/src/dune-opt/dune-pm/appl/fetd && gmsh -3 scattering-3D.geo
Info    : Running 'gmsh -3 scattering-3D.geo'
Info    : Started on Mon Dec  6 16:02:50 2010
Info    : Reading 'scattering-3D.geo'
Info    : Read 'scattering-3D.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1 (Circle)
Info    : Meshing curve 2 (Circle)
Info    : Meshing curve 3 (Circle)
Info    : Meshing curve 4 (Circle)
Info    : Meshing curve 5 (Circle)
Info    : Meshing curve 6 (Circle)
Info    : Meshing curve 7 (Circle)
Info    : Meshing curve 8 (Circle)
Info    : Meshing curve 9 (Circle)
Info    : Meshing curve 10 (Circle)
Info    : Meshing curve 11 (Circle)
Info    : Meshing curve 12 (Circle)
Info    : Meshing curve 30 (Circle)
Info    : Meshing curve 31 (Circle)
Info    : Meshing curve 32 (Circle)
Info    : Meshing curve 33 (Circle)
Info    : Meshing curve 34 (Circle)
Info    : Meshing curve 35 (Circle)
Info    : Meshing curve 36 (Circle)
Info    : Meshing curve 37 (Circle)
Info    : Meshing curve 38 (Circle)
Info    : Meshing curve 39 (Circle)
Info    : Meshing curve 40 (Circle)
Info    : Meshing curve 41 (Circle)
Info    : Mesh 1D complete (0 s)
Info    : Meshing 2D...
Info    : Meshing surface 14 (Ruled surface, MeshAdapt)
Info    : Meshing surface 16 (Ruled surface, MeshAdapt)
Info    : Meshing surface 18 (Ruled surface, MeshAdapt)
Info    : Meshing surface 20 (Ruled surface, MeshAdapt)
Info    : Meshing surface 22 (Ruled surface, MeshAdapt)
Info    : Meshing surface 24 (Ruled surface, MeshAdapt)
Info    : Meshing surface 26 (Ruled surface, MeshAdapt)
Info    : Meshing surface 28 (Ruled surface, MeshAdapt)
Info    : Meshing surface 43 (Ruled surface, MeshAdapt)
Info    : Meshing surface 45 (Ruled surface, MeshAdapt)
Info    : Meshing surface 47 (Ruled surface, MeshAdapt)
Info    : Meshing surface 49 (Ruled surface, MeshAdapt)
Info    : Meshing surface 51 (Ruled surface, MeshAdapt)
Info    : Meshing surface 53 (Ruled surface, MeshAdapt)
Info    : Meshing surface 55 (Ruled surface, MeshAdapt)
Info    : Meshing surface 57 (Ruled surface, MeshAdapt)
Info    : Mesh 2D complete (3.06819 s)
Info    : Meshing 3D...
Info    : Meshing volume 59 (Frontal)
Info    : Region 59 Face 28, 0 intersect
Info    : Region 59 Face 22, 0 intersect
Info    : Region 59 Face 14, 0 intersect
Info    : Region 59 Face 16, 0 intersect
Info    : Region 59 Face 18, 0 intersect
Info    : Region 59 Face 26, 0 intersect
Info    : Region 59 Face 24, 0 intersect
Info    : Region 59 Face 20, 0 intersect
Info    : CalcLocalH: 52 Points 0 Elements 100 Surface Elements 
Info    : Check subdomain 1 / 1 
Info    : 100 open elements 
Info    : Meshing subdomain 1 of 1 
Info    : 100 open elements 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : Delaunay meshing 
Info    : number of points: 52 
Info    : blockfill local h 
Info    : number of points: 52 
Info    : Points: 52 
Info    : Elements: 225 
Info    : 0 open elements 
Info    : Num open: 0 
Info    : free: 0, fixed: 225 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : 0 open elements 
Info    : Num open: 0 
Info    : free: 0, fixed: 225 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : 0 degenerated elements removed 
Info    : Remove intersecting 
Info    : Remove outer 
Info    : tables filled 
Info    : outer removed 
Info    : 100 open elements 
Info    : 52 points, 95 elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : free: 95, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : free: 95, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 95, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : free: 95, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : free: 95, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 95, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : free: 95, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : free: 95, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 95, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : free: 95, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 95, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 95 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 100 open elements 
Info    : 100 open faces found 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : free: 0, fixed: 0 
Info    : CombineImprove 
Info    : 0 elements combined 
Info    : free: 0, fixed: 0 
Info    : ImproveMesh 
Info    : 100 open elements 
Info    : Call remove problem 
Info    : 100 open elements 
Info    : Elements before Remove: 0 
Info    : Elements after Remove: 0 
Info    : 100 open elements 
Info    : 100 open elements 
Info    : 100 open faces 
Info    : Meshing volume 60 (Frontal)
Info    : Region 60 Face 57, 0 intersect
Info    : Region 60 Face 51, 0 intersect
Info    : Region 60 Face 43, 0 intersect
Info    : Region 60 Face 45, 0 intersect
Info    : Region 60 Face 47, 0 intersect
Info    : Region 60 Face 55, 0 intersect
Info    : Region 60 Face 53, 0 intersect
Info    : Region 60 Face 49, 0 intersect
Info    : Region 60 Face 28, 0 intersect
Info    : Region 60 Face 22, 0 intersect
Info    : Region 60 Face 14, 0 intersect
Info    : Region 60 Face 16, 0 intersect
Info    : Region 60 Face 18, 0 intersect
Info    : Region 60 Face 26, 0 intersect
Info    : Region 60 Face 24, 0 intersect
Info    : Region 60 Face 20, 0 intersect
Info    : CalcLocalH: 5698 Points 0 Elements 11388 Surface Elements 
Info    : Check subdomain 1 / 1 
Info    : 11388 open elements 
Info    : Meshing subdomain 1 of 1 
Info    : 11388 open elements 
Info    : Use internal rules 
Info    : 11388 open elements 
Info    : Delaunay meshing 
Info    : number of points: 5698 
Info    : blockfill local h 
Info    : number of points: 18489 
Info    : Points: 18489 
Info    : Elements: 114661 
Info    : 0 open elements 
Info    : Num open: 0 
Info    : free: 0, fixed: 114661 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : 0 open elements 
Info    : Num open: 0 
Info    : free: 0, fixed: 114661 
Info    : SwapImprove  
Info    : 0 swaps performed 
Info    : 0 degenerated elements removed 
Info    : Remove intersecting 
Info    : Remove outer 
Info    : tables filled 
Info    : outer removed 
Info    : 11388 open elements 
Info    : 18344 points, 100257 elements 
Info    : 11288 open elements 
Info    : 11288 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 11288 open elements 
Info    : 11288 open faces found 
Info    : free: 43190, fixed: 57067 
Info    : ImproveMesh 
Info    : free: 43190, fixed: 57067 
Info    : CombineImprove 
Info    : 656 elements combined 
Info    : free: 40203, fixed: 55209 
Info    : ImproveMesh 
Info    : free: 40203, fixed: 55209 
Info    : SwapImprove  
Info    : 4726 swaps performed 
Info    : free: 37067, fixed: 56022 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 37067, fixed: 56022 
Info    : ImproveMesh 
Info    : free: 37067, fixed: 56022 
Info    : CombineImprove 
Info    : 204 elements combined 
Info    : free: 36246, fixed: 55481 
Info    : ImproveMesh 
Info    : free: 36246, fixed: 55481 
Info    : SwapImprove  
Info    : 1177 swaps performed 
Info    : free: 35611, fixed: 55895 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 35611, fixed: 55895 
Info    : ImproveMesh 
Info    : free: 35611, fixed: 55895 
Info    : CombineImprove 
Info    : 52 elements combined 
Info    : free: 35410, fixed: 55752 
Info    : ImproveMesh 
Info    : free: 35410, fixed: 55752 
Info    : SwapImprove  
Info    : 523 swaps performed 
Info    : free: 35276, fixed: 55836 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 35276, fixed: 55836 
Info    : ImproveMesh 
Info    : free: 35276, fixed: 55836 
Info    : CombineImprove 
Info    : 23 elements combined 
Info    : free: 35215, fixed: 55745 
Info    : ImproveMesh 
Info    : 11288 open elements 
Info    : Call remove problem 
Info    : 11288 open elements 
Info    : Elements before Remove: 90960 
Info    : Elements after Remove: 68574 
Info    : 15654 open elements 
Info    : 15654 open elements 
Info    : 15654 open faces 
Info    : start tetmeshing 
Info    : Use internal rules 
Info    : 15654 open elements 
Info    : 15654 open faces found 
Info    : free: 27340, fixed: 41234 
Info    : ImproveMesh 
Info    : free: 27340, fixed: 41234 
Info    : CombineImprove 
Info    : 462 elements combined 
Info    : free: 25924, fixed: 39564 
Info    : ImproveMesh 
Info    : free: 25924, fixed: 39564 
Info    : SwapImprove  
Info    : 1185 swaps performed 
Info    : free: 25383, fixed: 39583 
Info    : SwapImprove2  
Info    : 0 swaps performed 
Info    : free: 25383, fixed: 39583 
Info    : ImproveMesh 
Info    : free: 25383, fixed: 39583 
Info    : CombineImprove 
Info    : 171 elements combined 
Info    : free: 24798, fixed: 39067 
Info    : ImproveMesh 
Info    : free: 24798, fixed: 39067 
Info    : SwapImprove  
Info    : Mesh 3D complete (22.5294 s)
Info    : 16780 vertices 75675 elements
Info    : Writing 'scattering-3D.msh'
/bin/bash: line 1:  5143 Segmentation fault      gmsh -3 scattering-3D.geo

Compilation exited abnormally with code 139 at Mon Dec  6 16:03:16
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 828 bytes
Desc: Digital signature
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20101206/949dae69/attachment.asc>