def METHOD_NAME():
    tape = Tape()
    set_working_tape(tape)
    mesh = IntervalMesh(10, 0, 1)
    V = FunctionSpace(mesh, "CG", 1)
    c = Constant((3.0, 4.0, 7.0))
    u = TrialFunction(V)
    v = TestFunction(V)
    sol = Function(V)
    a = inner(grad(u), grad(v))*dx
    L = inner(c, c)*v*dx
    bc = DirichletBC(V, 0.0, "on_boundary")
    solve(a == L, sol, bc)
    J = assemble(c**2*sol**2*dx)
    Jhat = ReducedFunctional(J, Control(c))
    assert(taylor_test(Jhat, c, Constant((1, 1, 1))) > 1.9)