[Gmsh] nodes don't align with geometry for simple rectangle

Al Danial al.danial at gmail.com
Fri Aug 26 04:55:38 CEST 2005


Thanks for running it on your system.  Wow.  This is a complete
mystery why the same code gives different results for me.  I'll
try removing the package, reinstalling, maybe finding another
machine and running it there....            -- Al


On 8/25/05, CHAMPANEY Laurent <laurent.champaney at meca.uvsq.fr> wrote:
> Hi,
> on my version (the same one : gmsh 1.60.1 via apt-get on debian),
> the mesh (.msh) is correct : see attached file.
> Laurent
> 
> 
> Al Danial a écrit :
> 
> >Here's something curious:  I create a rectangle having dimensions
> >of 6 x 2 units with
> >
> >Point(1) = {0, 0, 0, 0.5};
> >Point(2) = {6, 0, 0, 0.5};
> >Point(3) = {6, 2, 0, 0.5};
> >Point(4) = {0, 2, 0, 0.5};
> >Line (1) = {1, 4};
> >Line (2) = {2, 3};
> >Line (3) = {4, 3};
> >Line (4) = {1, 2};
> >Line Loop (1000006) = {1, 3, -2, -4};
> >Ruled Surface (6) = {1000006};
> >
> >then mesh it with triangular elements.  It gives 66 nodes and
> >130 elements and visually the nodes along the boundary line
> >up perfectly with the geometric Line entities.  So far, so good.
> >
> >However when I save the mesh, the nodal coordinates on the
> >right edge don't match the geometry.  Node 3, which is located
> >where Point 3 is, has coordinates of 5.7, 1.9 instead of
> >6.0, 2.0.  In fact all nodes on the right edge of the rectangle
> >have an X value of 5.7 instead of 6.0, and all nodes on the top
> >edge have a Y value of 1.9 instead of 2.0.
> >
> >I'm using version 1.60-1 (via Debian apt-get).
> >
> >Any hints as to why my boundary nodes don't match the
> >geometric boundary?            -- Al
> >
> >_______________________________________________
> >gmsh mailing list
> >gmsh at geuz.org
> >http://www.geuz.org/mailman/listinfo/gmsh
> >
> >
> >
> 
> --
> Laurent CHAMPANEY - UVSQ Dept de Mecanique
> Laboratoire d'Etudes Mécaniques des Assemblages (CNRS-FRE2481)
> RC 28 - Bat Descartes
> 45, av des Etats Unis - 78035 VERSAILLES Cedex
> Tel 01 39 25 30 26 - Fax 01 39 25 30 15
> http://www.lema.uvsq.fr
> http://www.ipsi.asso.fr (IPSI)
> http://www.librecours.org (Cours en Ligne)
> 
> 
> 
> 
> 
>