[petsc-users] Parallel TS for ODE

[Francesco Brarda]

Hi everyone!

I am trying to solve a system of 3 ODEs (a basic SIR model) with TS. Sequentially works pretty well, but I need to switch it into a parallel version. 
I started working with TS not very long time ago, there are few questions I’d like to share with you and if you have any advices I’d be happy to hear.
First of all, do I need to use a DM object even if the model is only time dependent? All the examples I found were using that object for the other variable when solving PDEs.
When I preallocate the space for the Jacobian matrix, is it better to decide the local or global space?

Best,
Francesco