def METHOD_NAME():
    m = UnitSquareMesh(1, 1)
    V = VectorFunctionSpace(m, 'CG', 1)
    f = Function(V)
    c = Constant([10, 11])
    f.assign(c)
    assert np.allclose(f.dat.data_ro[:, 0], 10)
    assert np.allclose(f.dat.data_ro[:, 1], 11)