[petsc-users] FEM on 2D poisson equation

Olivier Bonnefon olivier.bonnefon at avignon.inra.fr
Thu Jul 18 08:17:17 CDT 2013


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.

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 
> <mailto: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 <tel:%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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130718/4a80d55d/attachment.html>


More information about the petsc-users mailing list