[Gmsh] Structured grid for triangular domains

[Christophe Geuzaine]

> On 25 Mar 2020, at 20:46, Murilo Henrique Campana Bento <m.bento at usp.br> wrote:
> 
> Hi everyone,
> 
> I am trying to create a structured mesh for a simple triangular domain.
> However, the generated mesh is not the one I expected.
> 
> The input script (geo) I am using is the following.
> 
> "
> Point(1) = {0.0, 0.0, 0.0, 1.0};
> Point(2) = {1.0, 0.0, 0.0, 1.0};
> Point(3) = {1.0, 1.0, 0.0, 1.0};
> 
> Line(1) = {1, 2};
> Line(2) = {2, 3};
> Line(3) = {3, 1};
> 
> Line Loop(1) = {1, 2, 3};
> Plane Surface(1) = {1};
> 
> Transfinite Line {1, 2, 3} = 5;
> Transfinite Surface {1} Right;
> "
> 
> The mesh I am obtaining is illustrated in Fig01.
> The mesh I would expect to obtain is one similar to that presented in Fig02 (with little right triangles).
> 

Indeed, the transfinite algo does not have this decomposition coded. It should be quite easy to add though: look at Mesh/meshGFaceTransfinite.cpp. Don't hesitate to send a patch or a merge request with the code on our gitlab.

Christophe



> Thank you all for your help!
> Best regards,
> 
> Murilo Bento
> <Fig01.png><Fig02.png>_______________________________________________
> gmsh mailing list
> gmsh at onelab.info
> http://onelab.info/mailman/listinfo/gmsh

— 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 
http://www.montefiore.ulg.ac.be/~geuzaine