def METHOD_NAME():
    tape = Tape()
    set_working_tape(tape)
    mesh = IntervalMesh(10, 0, 1)
    V = FunctionSpace(mesh, "Lagrange", 1)
    c = Constant(1, domain=mesh)
    f = Function(V)
    f.vector()[:] = 1
    u = Function(V)
    v = TestFunction(V)
    bc = DirichletBC(V, c, "on_boundary")
    F = inner(grad(u), grad(v)) * dx - f ** 2 * v * dx
    solve(F == 0, u, bc)
    J = assemble(c ** 2 * u * dx)
    Jhat = ReducedFunctional(J, Control(c))
    c.block_variable.tlm_value = Constant(1)
    tape.evaluate_tlm()
    assert (taylor_test(Jhat, Constant(c), Constant(1), dJdm=J.block_variable.tlm_value) > 1.9)