[Gmsh] Problem with 2d background mesh

Jose Paulo Moitinho de Almeida moitinho at civil.ist.utl.pt
Wed Aug 13 20:32:43 CEST 2008


Hello

For the attached geometry and background mesh I am getting strange results 
which start at the 1d mesh and produce different results for 2.0.8 and 2.2.3.

As "usual" I have elements of quite different sizes...

-----------------------------
With 2.2.3 (the version I tested first)

$ gmsh aaa.geo -bgm aaa.bgm -noview a.script

Press 1 -> No nodes are generated on the upper left edge.
Press 2 -> The mesh is "cropped" to the left of the singular point.

$ gmsh bbb.geo -bgm aaa.bgm -noview a.script

(File bbb.geo is the same as aaa.geo, just switched lines 1 and 6).

Press 1 -> Way too many points on the upper left edge.
Press 2 -> The mesh is totally irregular near that edge.

-----------------------------
Switched to 2.0.8

$ gmsh aaa.geo -bgm aaa.bgm -noview a.script

Press 1 -> Same as for aaa.geo with 2.2.3.
Press 2 -> Lots of "Error: Colinear points in circum circle computation", but 
the mesh looks almost good. (1316 nodes, 2630 elements)

$ gmsh bbb.geo -bgm aaa.bgm -noview a.script

Press 1 -> Looks good!
Press 2 -> Again lots of "Error: Colinear points in circum circle 
computation", the mesh looks good (1314 nodes, 2626 elements)

-----------------------------

Changing the algorithm does not seem to help. But the problem is really on 
providing info about the background mesh to the 1d generator, correct?

I tried to look at the source code (BGM_MeshSize?), but must admit my failure 
on understanding "what is coming from where".

Any suggestions?

Thanks

ZP
-------------- next part --------------
A non-text attachment was scrubbed...
Name: a.script
Type: text/x-objcsrc
Size: 114 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20080813/f6a4f99c/attachment.script>
-------------- next part --------------
Point(1) = {0,0,0,1};
Point(2) = {1,0,0,1};
Point(3) = {1,0.5,0,1};
Point(4) = {1,0.8,0,1};
Point(5) = {1,1,0,1};
Point(6) = {0.5,1,0,.001};
Point(7) = {0,1,0,.001};
Point(8) = {0,0.8,0,1};
Point(9) = {0.5,0.8,0,1};
Point(10) = {0.5,0.5,0,1};
Point(11) = {0,0.5,0,1};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,5};
Line(5) = {5,6};
Line(6) = {6,7};
Line(7) = {7,8};
Line(8) = {8,9};
Line(9) = {9,10};
Line(10) = {10,11};
Line(11) = {11,1};
Line Loop(12) = {1,2,3,4,5,6,7,8,9,10,11};
Plane Surface(13) = {12};
Physical Surface(14) = {13};
Physical Line(101) = {1};
Physical Line(102) = {2};
Physical Line(103) = {3};
Physical Line(104) = {4};
Physical Line(105) = {5};
Physical Line(106) = {6};
Physical Line(107) = {7};
Physical Line(108) = {8};
Physical Line(109) = {9};
Physical Line(110) = {10};
Physical Line(111) = {11};
-------------- next part --------------
Point(1) = {0,0,0,1};
Point(2) = {1,0,0,1};
Point(3) = {1,0.5,0,1};
Point(4) = {1,0.8,0,1};
Point(5) = {1,1,0,1};
Point(6) = {0.5,1,0,.001};
Point(7) = {0,1,0,.001};
Point(8) = {0,0.8,0,1};
Point(9) = {0.5,0.8,0,1};
Point(10) = {0.5,0.5,0,1};
Point(11) = {0,0.5,0,1};
Line(6) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,5};
Line(5) = {5,6};
Line(1) = {6,7};
Line(7) = {7,8};
Line(8) = {8,9};
Line(9) = {9,10};
Line(10) = {10,11};
Line(11) = {11,1};
Line Loop(12) = {1,2,3,4,5,6,7,8,9,10,11};
Plane Surface(13) = {12};
Physical Surface(14) = {13};
Physical Line(101) = {1};
Physical Line(102) = {2};
Physical Line(103) = {3};
Physical Line(104) = {4};
Physical Line(105) = {5};
Physical Line(106) = {6};
Physical Line(107) = {7};
Physical Line(108) = {8};
Physical Line(109) = {9};
Physical Line(110) = {10};
Physical Line(111) = {11};
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aaa.bgm.bz2
Type: application/x-bzip2
Size: 40031 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20080813/f6a4f99c/attachment.bz2>