pysdic.Mesh.from_meshio

classmethod Mesh.from_meshio(mesh, element_type=None, load_properties=True)

Create a Mesh instance from a meshio.Mesh object.

The following fields are extracted:

• mesh.points → vertices

• mesh.cells[0].data → connectivity

• mesh.point_data → vertices properties as arrays of shape (N, A)

• mesh.cell_data → elements properties as arrays of shape (M, B)

See also

• Mesh.to_meshio() for the reverse operation.

• meshio documentation for more information.

• Mesh for more information on the Mesh class and element types.

Parameters:

• mesh (meshio.Mesh) – A meshio Mesh object to extract the first cell block and create the Mesh instance.

• element_type (Optional[str], optional) – The expected type of elements in the mesh, by default None.

• load_properties (bool, optional) – If True, properties are extracted from the meshio.Mesh object, by default True.

Returns:

A Mesh instance created from the meshio.Mesh object.

Return type:

Mesh

Raises:

• TypeError – If the input is not a meshio.Mesh object. If element_type is not a string. If load_properties is not a boolean.

• ValueError – If the mesh structure is invalid. If element_type is invalid or does not match the mesh structure.

Examples

Lets create a mesh using meshio.Mesh and convert it to a Mesh instance.

import numpy as np
import meshio
from pysdic import Mesh

points = np.array([[0, 0, 0], [1, 0, 0], [0, 1, 0], [0, 0, 1]])
cells = [("triangle", np.array([[0, 1, 2], [0, 1, 3], [0, 2, 3], [1, 2, 3]]))]

mesh = meshio.Mesh(points=points, cells=cells)

Create a Mesh instance from the meshio.Mesh object.

mesh3d = Mesh.from_meshio(mesh, element_type="triangle_3")
print(mesh3d.vertices)
# Output: PointCloud with 4 points [[0, 0, 0], [1, 0, 0], [0, 1, 0], [0, 0, 1]]