[petsc-users] TS scheme with different DAs

[Matthew Knepley]

On Tue, Sep 17, 2019 at 6:15 PM Manuel Valera via petsc-users <
petsc-users at mcs.anl.gov> wrote:

> Hello, petsc users,
>
> I have integrated the TS routines in my code, but i just noticed i didn't
> do it optimally. I was using 3 different TS objects to integrate
> velocities, temperature and salinity, and it works but only for small DTs.
> I suspect the intermediate Runge-Kutta states are unphased and this creates
> the discrepancy for broader time steps, so I need to integrate the 3
> quantities in the same routine.
>
> I tried to do this by using a 5 DOF distributed array for the RHS, where I
> store the velocities in the first 3 and then Temperature and Salinity in
> the rest. The problem is that I use a staggered grid and T,S are located in
> a different DA layout than the velocities. This is creating problems for me
> since I can't find a way to communicate the information from the result of
> the TS integration back to the respective DAs of each variable.
>
> Is there a way to communicate across DAs? or can you suggest
> an alternative solution to this problem?
>

If you have a staggered discretization on a structured grid, I would
recommend checking out DMStag.

  Thanks,

     MAtt


> Thanks,
>
>

-- 
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/20190917/8f232177/attachment.html>