def METHOD_NAME():
    parent_mesh = UnitSquareMesh(100, 100, quadrilateral=True)
    points = [[0.1, 0.1], [0.2, 0.2], [1.1, 1.0]]
    V = FunctionSpace(parent_mesh, 'CG', 2)
    # The point (1.1, 1.0) will still be included in the vertex-only mesh
    vom = VertexOnlyMesh(parent_mesh, points, tolerance=30.0)
    # The tolerance property has been changed - this will print 30.0
    print(parent_mesh.tolerance)
    # This doesn't generate an error
    Function(V).at((1.1, 1.0), tolerance=20.0)
    # The tolerance property has been changed again - this will print 20.0
    print(parent_mesh.tolerance)
    try:
        # This generates an error
        Function(V).at((1.1, 1.0), tolerance=1.0)
    except PointNotInDomainError:
        # But the tolerance property has still been changed - this will print 1.0
        print(parent_mesh.tolerance)
    assert vom