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
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
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
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
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
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