Thanks alot for all of your suggestions. I think I have a better insight about the direction now.
On Thu, Oct 6, 2016 at 2:03 AM, Dave May <[email protected]> wrote: > > > On 5 October 2016 at 18:49, Matthew Knepley <[email protected]> wrote: > >> On Wed, Oct 5, 2016 at 11:19 AM, E. Tadeu <[email protected]> wrote: >> >>> Matt, >>> >>> Do you know if there is any example of solving Navier Stokes using a >>> staggered approach by using a different DM object such as DMPlex? >>> >> >> SNES ex62 can do P2/P1 Stokes, which is similar. Is that what you want to >> see? >> >> For real structured grid, staggered mesh stuff like MAC, I would just do >> this on a single DMDA, but think of it as being staggered, and expand my >> stencil as necessary. >> > > Following that up, for a DMDA example using a staggered grid, take a look > at snes/ex30.c > > http://www.mcs.anl.gov/petsc/petsc-current/src/snes/ > examples/tutorials/ex30.c.html > > Thanks, > Dave > > >> >> Thanks, >> >> Matt >> >> >>> >>> Thanks, >>> Edson >>> >>> >>> On Tue, Oct 4, 2016 at 11:12 PM, Matthew Knepley <[email protected]> >>> wrote: >>> >>>> On Tue, Oct 4, 2016 at 9:02 PM, Somdeb Bandopadhyay <[email protected] >>>> > wrote: >>>> >>>>> Dear all, >>>>> I want to write a solver for incompressible navier stokes >>>>> using python and I want to use PETsc (particularly dmda & ksp) for this. >>>>> May I know if this type of work is feasible/already done? >>>>> >>>> >>>> How do you plan to discretize your system? DMDA supports only >>>> collocation discretizations, so some sort of penalty for pressure would >>>> have to be employed. >>>> >>>> Thanks, >>>> >>>> Matt >>>> >>>> >>>>> I intend to run my solver in a cluster and so am slightly >>>>> concerned about the performance if I use python with petsc. >>>>> My deepest apologies if this mail of mine caused you any >>>>> inconvenience. >>>>> >>>>> Somdeb >>>>> >>>> >>>> >>>> >>>> -- >>>> 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 >>>> >>> >>> >> >> >> -- >> 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 >> > >
