I am really happy to communicate in your forum asking for your valuable help.

In my current studies I try to implement a program in Freefem++ using the discontinuous Galerkin time method. I know that when I have discontinuities in space I work in P1dc or P2dc space adjusting my program.

My main goal is to implement a program in Freefem++ using jumps in time. So I try to implement a program in discontinuous Galerkin using linear polynomials in time. I have implemented theta scheme for a parabolic problem like heat equation with nonlinear terms.

I suppose that what I have to do is to use a scheme like Crank Nicolson method using time jumps.

My question is how I can add time jumps for the weak form implementation.

Any contribution will be appreciated.