def METHOD_NAME(V, expr, solution, tolerance):
    """Test projection on P1 x Real space"""
    u = TrialFunction(V)
    v = TestFunction(V)
    out = Function(V)
    solve(inner(u, v)*dx == inner(expr, v)*dx, out)
    l2err = sqrt(assemble(inner(out-solution, out-solution)*dx))
    assert abs(l2err) < tolerance