Hi again Sir, Thank you very much for the quick response. I am planning to implement a mustiphase algorithm on collocated grid. I already qrote a C code for 2d case, but it wasn't very generalized . So for the final version, I intend to use python as a script to interact with PETSc kernels.
Somdeb On Wed, Oct 5, 2016 at 10:12 AM, Matthew Knepley <knep...@gmail.com> wrote: > On Tue, Oct 4, 2016 at 9:02 PM, Somdeb Bandopadhyay <sb020...@gmail.com> > 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 >