[petsc-users] Parallel TS for ODE

[Matthew Knepley]

On Wed, Mar 31, 2021 at 3:54 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.
>

You do not need one. We use it in examples because it makes it easy to
create the data.


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

Since you are producing all the Jacobian values, it is whatever is easier
in your code I think.

  THanks,

    Matt


> Best,
> Francesco



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20210331/b52f5568/attachment.html>