[petsc-users] using DMDA with python

[Dave May]

On 5 October 2016 at 18:49, Matthew Knepley <knepley at gmail.com> wrote:

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

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 <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/036b1238/attachment.html>