[petsc-users] using DMDA with python

[Matthew Knepley]

On Wed, Oct 5, 2016 at 11:19 AM, E. Tadeu <e.tadeu at gmail.com> 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.

  Thanks,

     Matt


>
>   Thanks,
>     Edson
>
>
> On Tue, Oct 4, 2016 at 11:12 PM, Matthew Knepley <knepley at gmail.com>
> wrote:
>
>> On Tue, Oct 4, 2016 at 9:02 PM, Somdeb Bandopadhyay <sb020287 at 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
>>
>
>


-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161005/4bd4e78b/attachment.html>