Re: [petsc-users] TS scheme with different DAs

2019-09-26 Thread Manuel Valera via petsc-users
Hi all, Just retouching in case my mail was lost, does my problem make enough sense? I could try to be clearer if you like. I tried using TSGetStages() in the middle of the RHS function but that didn't work. Maybe we can find a different approach for this ? Thanks, On Tue, Sep 24, 2019 at 12:58

Re: [petsc-users] TS scheme with different DAs

2019-09-24 Thread Manuel Valera via petsc-users
Hello all, I finally implemented the TS routine operating in several DAs at the same time, hacking it as you suggested. I still have a problem with my algorithm though. It is not DMDA related so there's that. My algorithm needs to update u,v,w with information from the updated T,S,rho. My problem

Re: [petsc-users] TS scheme with different DAs

2019-09-18 Thread Matthew Knepley via petsc-users
On Tue, Sep 17, 2019 at 8:27 PM Smith, Barry F. wrote: > > Don't be too quick to dismiss switching to the DMStag you may find that > it actually takes little time to convert and then you have a much less > cumbersome process to manage the staggered grid. Take a look at > src/dm/impls/stag/examp

Re: [petsc-users] TS scheme with different DAs

2019-09-17 Thread Smith, Barry F. via petsc-users
Don't be too quick to dismiss switching to the DMStag you may find that it actually takes little time to convert and then you have a much less cumbersome process to manage the staggered grid. Take a look at src/dm/impls/stag/examples/tutorials/ex2.c where const PetscInt dof0 = 0, dof1 = 1,d

Re: [petsc-users] TS scheme with different DAs

2019-09-17 Thread Manuel Valera via petsc-users
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 wrote: > O

[petsc-users] TS scheme with different DAs

2019-09-17 Thread Manuel Valera via petsc-users
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 t