[petsc-users] FEM on 2D poisson equation

[Matthew Knepley]

On Thu, Jul 18, 2013 at 8:17 AM, Olivier Bonnefon <
olivier.bonnefon at avignon.inra.fr> wrote:

>  It is what I wanted, it works.
> If I well understand the code, ex12.h contains the P1 implementation. To
> simulate an other system, with time dependences for examples (du/dt), I
> have to adapt the plugin functions.
>

The way I would add time dependence is to convert this from a SNES example
into a TS example. I can help you
do this since I want to start using TS by default. Does this sound
reasonable?

  Thanks,

     Matt


> Thanks a lot.
>
> Olivier B
>
> On 07/18/2013 01:12 PM, Matthew Knepley wrote:
>
> On Thu, Jul 18, 2013 at 5:08 AM, Olivier Bonnefon <
> olivier.bonnefon at avignon.inra.fr> wrote:
>
>> Hello,
>>
>> I have a 2-d heat equation that I want to simulate with Finit Element
>> Method, to do this, I'm looking for an example solving 2D poisson equation
>> with FEM (DMDA or DMPlex). Is there an example like this ?
>>
>
>  There is, but there it is still somewhat problematic. I use FIAT to
> generate the basis function tabulation,
> so you have to configure with
>
>    --download-fiat --download-scientificpython --download-generator
>
>  and you need mesh generation and partitioning
>
>    --download-triangle --download-chaco
>
>  and then you can run SNES ex12 using Builder (which will make the header
> file)
>
>    python2.7 ./config/builder2.py check src/snes/examples/tutorials/ex12.c
>
>  Jed and I are working on an all C version of tabulation which would mean
> that you could bypass
>  the Python code generation step. Once the header is generated for the
> element you want, then
> you can just run the example as normal.
>
>      Matt
>
>
>> Thanks a lot.
>>
>> Olivier Bonnefon
>>
>> --
>> Olivier Bonnefon
>> INRA PACA-Avignon, Unité BioSP
>> Tel: +33 (0)4 32 72 21 58 <%2B33%20%280%294%2032%2072%2021%2058>
>>
>>
>
>
>  --
> 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
>
>
>
> --
> Olivier Bonnefon
> INRA PACA-Avignon, Unité BioSP
> Tel: +33 (0)4 32 72 21 58
>
>


-- 
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/20130718/82f31317/attachment.html>