#%%
import numpy as np
from sfepy.discrete.fem import Mesh, FEDomain, Field

mesh_path = 'C:/Users/lzy71/miniconda3/envs/lego/lib/site-packages/sfepy/meshes/2d/rectangle_tri.mesh'

#%%
mesh = Mesh.from_file(mesh_path)
domain = 
FEDomain('domain', mesh)