Disk Example: Changing Order of Elements in Mesh?

[Chaztikov]

Hi again,

This is somewhat of an extension of my previous question, but lies at the crux of my problems. If there is a better medium through which I should ask questions, please let me know.

Question 3
For the following simple example of a disk (code that follows works)
"Generate Disk Geometry"

geom = pygmsh.built_in.Geometry()

#add points for circle
x0 = (0.,0., 0.)
radius = 1
lcar = 10
#approximate circle by polygon with #edges = 2*num_sections
num_sections = 40
poly = geom.add_circle(x0, radius, lcar, R=None, compound=False, num_sections=num_sections, holes=None, make_surface=True)

#generate mesh from geom
points, cells, point_data, cell_data, field_data = pygmsh.generate_mesh(geom)

import meshio
meshio.write('disk.vtu', points, cells, cell_data=cell_data)

How do I change the 'order' of the elements (i.e. the number of nodes on a given triangle, in this example).

Thanks, and Happy Thanksgiving
Charlie

[nschloe]

If there is a better medium through which I should ask questions, please let me know.

Here is fine.

For second order I think you need to specify

Mesh.ElementOrder = 2;

in gmsh, hence

geom.add_raw_code('Mesh.ElementOrder = 2;')

[Chaztikov]

Thanks! I'll give it a try and let you know

…

On Thu, Nov 23, 2017 at 4:22 AM, Nico Schlömer ***@***.***> wrote: If there is a better medium through which I should ask questions, please let me know. Here is fine. For second order I think you need to specify Mesh.ElementOrder = 2; in gmsh, hence geom.add_raw_code('Mesh.ElementOrder = 2;') — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#107 (comment)>, or mute the thread <https://github.com/notifications/unsubscribe-auth/AI9ocVWGaZTzuYq77iZhbrel1ZBFUoJQks5s5TlQgaJpZM4QoI1O> .