[petsc-users] FEM on 2D poisson equation

Matthew Knepley knepley at gmail.com
Thu Jul 18 06:12:31 CDT 2013


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


-- 
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/480e70fc/attachment.html>


More information about the petsc-users mailing list