#%%
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)

min_x, max_x = domain.get_mesh_bounding_box()[:, 0]
eps = 1e-8 * (max_x - min_x)
omega = domain.create_region('Omega', 'all')

# %%
from sfepy.discrete import (
    FieldVariable, Material, Integral, Function,
    Equation, Equations, Problem
    )

field = 
Field.from_args('fu', np.float64, 'vector', omega, approx_order=2)