Bug in nonlinearconvectiondiffusionfem.hh
@chaiyod.kamthorncharoen found a bug in nonlinearconvectiondiffusionfem.hh:
The documentation talks about solving
q(x,u) - D(x) v(u) \nabla w(u) = f(u)
but the code seems to be doing
q(x, w(u)) - D(x) v(w(u)) \nabla w(u) = f(w(u))
The implementation is quite confusing as the the names for u and w are mixed up but I think this is what happens.
Now the question is: Is this a documentation bug or a bug in the localoperator? What is the equation that this local operator is supposed to solve?
PS: There are some obvious documentation bugs in the description of the Neumann boundary condition that should definitely be fixed. But there it is clear what to do.