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: [EMAIL PROTECTED] --------------------------------------------------------------------------------------- _______________________________________________ gmsh mailing list [email protected] http://www.geuz.org/mailman/listinfo/gmsh
