[petsc-users] using DMDA with python

[Somdeb Bandopadhyay]

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 <dave.mayhem23 at gmail.com> wrote:

>
>
> 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/20161006/9bb5663e/attachment.html>