When using the following method to solve a nonlinear PDE
J = derivative(F, w)
solve(F == 0, w, J=J)
Instead of using div(arg) for a 1D problem I have used (arg).dx(0) as suggested by Miro Kutcha.
Is there another alternative for this (besides writing it term by term) that can work in higher dimensions?