[Gmsh] 2D unstructured mesh problem

[Umberto Paoletti]

Thanks.

I understand, the "line" instruction defines only a segment between two 
points.
After not using the program for a while, I did not remember this 
fundamental "detail", which is also clearly stated in the manual.
It's my fault.

However, I would like to make one comment.
My mistake has been supported by the fact the program seems to accept 
the instruction also with a list of more than two points.
The geometry also is built correctly, and the problem comes out only at 
some corners during meshing.

Thank you again
Umberto



David Colignon wrote:
> Hi,
>
> your Line and Line Loop definitions are wrong.
>
> You should have a look at the tutorial t1.geo
>
> Cheers,
>
> Dave
>
> -- 
> David Colignon, Ph.D.
> Collaborateur Logistique F.R.S.-FNRS
> CÉCI - Consortium des Équipements de Calcul Intensif
> ACE - Applied & Computational Electromagnetics
> Institut Montefiore B28
> Université de Liège
> 4000 Liège - BELGIQUE
> Tél: +32 (0)4 366 37 32
> Fax: +32 (0)4 366 29 10
> WWW:    http://hpc.montefiore.ulg.ac.be/
> Agenda: 
> http://www.google.com/calendar/embed?src=david.colignon%40gmail.com
>
>
> Umberto Paoletti wrote:
>> Hello,
>>
>> I found a problem in this very simple example.
>> Some corners of the internal hole seem to be skipped by the mesh 
>> algorithm.
>> Why?
>>
>> Thanks
>> Umberto
>>
>>
>> unit =  0.001000;
>>
>> lc = 0.1*unit;
>>
>> edge = 2*unit;
>> holedge = 0.5*unit;
>>
>> z = 0*unit;
>>
>> //external edge
>> x1 = 0*unit;
>> y1 = 0*unit;
>> x2 = x1+edge;
>> y2 = y1;
>> x3 = x2;
>> y3 = y2+edge;
>> x4 = x3-edge;
>> y4 = y3;
>> Point(1) = {x1,y1,z, 1*lc};
>> Point(2) = {x2,y2,z, 1*lc};
>> Point(3) = {x3,y3,z, 1*lc};
>> Point(4) = {x4,y4,z, 1*lc};
>> Line(1) = {1:4,1};
>> Line Loop(2) = {1};
>>
>> x5 = 1*unit;
>> y5 = 1*unit;
>> x6 = x5+holedge;
>> y6 = y5;
>> x7 = x6;
>> y7 = y6+holedge;
>> x8 = x7-holedge;
>> y8 = y7;
>>
>> //hole edge
>> Point(5) = {x5,y5,z, 1*lc};
>> Point(6) = {x6,y6,z, 1*lc};
>> Point(7) = {x7,y7,z, 1*lc};
>> Point(8) = {x8,y8,z, 1*lc};
>> Line(3) = {5:8,5};
>> Line Loop(4) = {3};
>>
>>
>> Plane Surface(1) = {2,4};
>>
>>
>>
>
>

-- 
---------------------------------------------------------------------------------------
Umberto Paoletti, Dr. Eng.,
Building A1, Room 402,
Graduate School of Engineering,
Department of Electrical Engineering,
Kyoto University, Nishikyo-ku,
Katsura, 615-8510, Kyoto, Japan
Tel.: +81-(0)75-383-2246
Fax.: +81-(0)75-383-2245
E-mail: paoletti at circuit.kuee.kyoto-u.ac.jp
---------------------------------------------------------------------------------------