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 <knep...@gmail.com> wrote: > On Tue, Sep 17, 2019 at 6:15 PM Manuel Valera via petsc-users < > petsc-users@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/> >