fix gls bug: the correct is dx_pnew with new values at mid domain
adapt the code to the new petsc version
RK4 for the mesh velocity when it is a known
notebook last state
fix conflicts