def METHOD_NAME():
    import sfepy
    from sfepy.base.base import Struct
    from sfepy.discrete.fem import Mesh, FEDomain, Field
    mesh = Mesh.from_file('meshes/2d/rectangle_tri.mesh',
                          prefix_dir=sfepy.data_dir)
    domain = FEDomain('domain', mesh)
    dim = domain.shape.dim
    min_x, max_x = domain.get_mesh_bounding_box()[:,0]
    eps = 1e-8 * (max_x - min_x)
    omega = domain.create_region('Omega', 'all')
    gamma1 = domain.create_region('Gamma1',
                                  'vertices in x < %.10f' % (min_x + eps),
                                  'facet')
    gamma2 = domain.create_region('Gamma2',
                                  'vertices in x > %.10f' % (max_x - eps),
                                  'facet')
    field = Field.from_args('fu', nm.float64, 'vector', omega,
                            approx_order=2)
    return Struct(dim=dim, omega=omega, gamma1=gamma1, gamma2=gamma2,
                  field=field)