Can freefem++be used to calculate discontinuous finite element implicit schemes? Is it through Newton iteration? Does Freefrem++have built-in Newton iteration?
You could use PETSc SNES , e.g., FreeFem-sources/newton-vi-2d-PETSc.edp at master · FreeFem/FreeFem-sources · GitHub