[petsc-users] Parallel TS for ODE

[Zhang, Hong]

> On Mar 31, 2021, at 2:53 AM, Francesco Brarda <brardafrancesco at gmail.com> wrote:
> 
> 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.

Are you considering SIR on a spatial domain? If so, you can parallelize your model in the spatial domain using DM. Splitting the three variables in the ODE among processors would not scale.

Hong (Mr.)  

> When I preallocate the space for the Jacobian matrix, is it better to decide the local or global space?
> 
> Best,
> Francesco