Movemesh() to generate a circle

General Discussion
1

Hi there, if we use

Th = movemesh(Th0, [xcos(y), xsin(y)])

to generate a circular mesh

image
image671×676 25.3 KB

from a rectangle:

image
image222×663 1.46 KB

How to collapse the overlapped line at y=0 so that it becomes a proper structured mesh?

Seems there is parameter ptmerge, but it is only for 3D?

Thank you for your help in advance!

2

Sorry due to missing code you have to use two 2 mesh add glut this 2 mesh with operator +

mesh Th=square(32,10,[(1+y)*cos(pi*x),(1+y)*sin(x*pi)]) + square(32,10,[(1+y)*cos(pi*(x+1)),(1+y)*sin((x+1)*pi)])  ;
Th =change(Th,rmInternalEdges=1);
plot(Th,wait=1);
1 Like
3

Thank you very much!

4

This bug was fixed in develop branch,

1 Like