def METHOD_NAME():
    gdim, shape, degree = 2, "triangle", 1
    domain = ufl.Mesh(element("Lagrange", shape, degree, gdim=gdim, rank=1))
    x = [[0., 0., 0.], [0., 1., 0.], [1., 1., 0.]]
    cells = [[0, 1, 2]]
    mesh = create_mesh(MPI.COMM_WORLD, cells, x, domain)
    assert squared_distance(mesh, mesh.topology.dim, [0], [-1.0, -1.0, 0.0]) == pytest.approx(2.0)
    assert squared_distance(mesh, mesh.topology.dim, [0], [-1.0, 0.5, 0.0]) == pytest.approx(1.0)
    assert squared_distance(mesh, mesh.topology.dim, [0], [0.5, 0.5, 0.0]) == pytest.approx(0.0)