[petsc-users] FEM on 2D poisson equation

[Olivier Bonnefon]

On 07/18/2013 03:26 PM, Matthew Knepley wrote:
> On Thu, Jul 18, 2013 at 8:17 AM, Olivier Bonnefon 
> <olivier.bonnefon at avignon.inra.fr 
> <mailto: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?
Yes, of course. My goal is to simulate diffusive equation with non 
linear sources, for example Lotka-Voltera competion.

Olivier B
>
>   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
>>     <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  <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/874ec63a/attachment-0001.html>