How to implement jump(u^3) for nonlinear DG flux in FreeFEM++

General Discussion
1

Hi, I am implementing a DG scheme in FreeFEM++ for a nonlinear diffusion problem.

I need to write the jump of b(u)=u+u^3/3 in intalledges. Writing it directly gives me the compilation error. Can jump() handle nonlinear expressions like u^3? If not, what is the correct way to implement the jump of b(u) across DG faces in FreeFEM++?

2

You can set w=u+u^3/3 (w declared in the same fespace as u), and then put jump(w).

3

Thank you. This is helpful.