Hi
This question is more related to Petsc than fenics, but though I might try to post it here.
In petsc it seems like the framework for solving time dependent problems (using the TS objects) expects that the differential equations are only first order in time. What do you do when you for instance want to solve a second order ode like the 1D wave equation (u_tt = c^2*u_xx) using central differences? Can this problem be reformulated to a first order problem, or are petsc only handling second first order problems?
Thanks in advance.