Create mesh from list of nodes and elements

Question

Hello everybody

Is it possible to obtain a mesh from lists of the node coordinates and the element connections?

A simple example for one Element:

nodes = np.array([[0., 0., 0.],
                  [1., 0., 0.],
                  [0., 1., 0.],
                  [0., 0., 1.]])
elements = np.array([[1,2,3,4]])

And I want to do something like this:

from fenics import *
mesh = Mesh(nodes, elements)

Is this possible and how?

Thanks in advance!

Answer 1

Yes, use the MeshEditor-class:

from dolfin import *
import numpy as np
nodes = np.array([[0., 0., 0.],
                  [1., 0., 0.],
                  [0., 1., 0.],
                  [0., 0., 1.]])

cells = np.array([[0,1,2,3]], dtype=np.uintp)

mesh = Mesh()
editor = MeshEditor()
editor.open(mesh, 3, 3)
editor.init_vertices(4)
editor.init_cells(1)

[editor.add_vertex(i,n) for i,n in enumerate(nodes)]
[editor.add_cell(i,n) for i,n in enumerate(cells)]
editor.close()

plot(mesh)
interactive()

Comments

Thank you for your response!
It works perfectly.