def METHOD_NAME():
    mesh = UnitIntervalMesh(1)
    V = FunctionSpace(mesh, "CG", 1)
    x = SpatialCoordinate(mesh)
    f = Function(V).interpolate(2.0 * x[0])
    # raise exception without dont_raise
    with pytest.raises(PointNotInDomainError):
        f.at(-1)
    # dont_raise=True
    assert f.at(-1, dont_raise=True) is None
    actual = f.at([-1, 0, 0.5], dont_raise=True)
    assert actual[0] is None
    assert np.allclose([0.0, 1.0], actual[1:])