def METHOD_NAME():
    mesh = UnitSquareMesh(1,1)
    V = FunctionSpace(mesh, "CG", 1)
    u = Function(V)
    c = Constant(1., domain=mesh)
    project(u+c, u)
    project(u+c*u**2, u)
    J = assemble(u**2*dx)
    rf = ReducedFunctional(J, Control(c))
    h = Constant(0.1)
    assert taylor_test(rf, Constant(3., domain=mesh), h)