def METHOD_NAME(fes, cf):
    test_function = fes.TestFunction()
    L = LinearForm(fes)
    L += InnerProduct(cf, test_function) * dx
    L += InnerProduct(cf, test_function) * ds(element_boundary=True)
    L.Assemble()
    return L