def METHOD_NAME(mesh):
    x = SpatialCoordinate(mesh)
    # test function space
    v = FunctionSpace(mesh, "DG", 1)
    # Create limiter
    limiter = VertexBasedLimiter(v)
    # advecting velocity
    u0 = conditional(real(x[0]) < 0.5, 1., 0.)
    u = Function(v).interpolate(u0)
    limiter.apply(u)
    assert np.max(u.dat.data_ro) <= 1.0, "Failed by exceeding max values"
    assert np.min(u.dat.data_ro) >= 0.0, "Failed by exceeding min values"