[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>