[Gmsh] making an equilateral mesh.
David Colignon
David.Colignon at ulg.ac.be
Wed Oct 7 12:20:22 CEST 2009
Hi Colin,
what about this ?
Point(1) = {0,0,0,0.1};
Extrude {1,0,0} { Point{1}; }
Extrude {0.5,0.866025403784,0} { Line{1}; }
Transfinite Line { 1 , 2 , 3 , 4 } = 5 ;
Transfinite Surface {5} Left ;
Cheers,
Dave
--
David Colignon, Ph.D.
Collaborateur Logistique du F.R.S.-FNRS
CÉCI - Consortium des Équipements de Calcul Intensif
ACE - Applied & Computational Electromagnetics
Sart-Tilman 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
Dr Colin J. Cotter wrote:
> Dear Gmsh list,
> I have been trying to make a mesh of equilateral triangles by
> extruding a horizontal line at a 60 degree angle (later I will use
> periodic boundary conditions so I have an equilateral mesh on the
> torus). I have noticed that if you do this then the extrusion divides
> the quadrilaterals in the wrong direction and you don't get a Delaunay
> mesh (the equilateral one):
>
> Point(1) = {0,0,0,0.1};
> Extrude {1,0,0} {
> Point{1}; Layers{5};
> }
> Extrude {0.5,0.866025403784,0} {
> Line{1}; Layers{5};
> }
>
> but it works if you extrude at a -60 degree angle instead:
>
> Point(1) = {0,0,0,0.1};
> Extrude {1,0,0} {
> Point{1}; Layers{5};
> }
> Extrude {-0.5,0.866025403784,0} {
> Line{1}; Layers{5};
> }
>
> Is there a way to force Gmsh to produce a Delaunay mesh in the first case?
>
> all the best
> --cjc
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh