[petsc-users] TS scheme with different DAs
Manuel Valera
mvalera-w at sdsu.edu
Tue Sep 17 19:04:04 CDT 2019
Thanks Matthew, but my code is too complicated to be redone on DMStag now
after spending a long time using DMDAs,
Is there a way to ensure PETSc distributes several DAs in the same way?
besides manually distributing the points,
Thanks,
On Tue, Sep 17, 2019 at 3:28 PM Matthew Knepley <knepley at gmail.com> wrote:
> 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/4d96eab9/attachment.html>
More information about the petsc-users
mailing list