def METHOD_NAME(u, m):
    v = TestFunction(V)
    nu = Constant(0.0001)
    # Here m is the initial condition.
    METHOD_NAME = ((u - m) * v
         + m * u.dx(0) * v
         + nu * u.dx(0) * v.dx(0)) * dx
    return METHOD_NAME