<div dir="ltr">On Thu, Jul 18, 2013 at 8:17 AM, Olivier Bonnefon <span dir="ltr"><<a href="mailto:olivier.bonnefon@avignon.inra.fr" target="_blank">olivier.bonnefon@avignon.inra.fr</a>></span> wrote:<br><div class="gmail_extra">
<div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div bgcolor="#FFFFFF" text="#000000">
    It is what I wanted, it works. <br>
    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.<br></div></blockquote><div><br></div><div>The way I would add time dependence is to convert this from a SNES example into a TS example. I can help you</div>
<div>do this since I want to start using TS by default. Does this sound reasonable?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
    Thanks a lot.<br>
    <br>
    Olivier B<br>
    <br>
    On 07/18/2013 01:12 PM, Matthew Knepley wrote:
    <blockquote type="cite">
      <div dir="ltr">On Thu, Jul 18, 2013 at 5:08 AM, Olivier Bonnefon <span dir="ltr"><<a href="mailto:olivier.bonnefon@avignon.inra.fr" target="_blank">olivier.bonnefon@avignon.inra.fr</a>></span>
        wrote:<br>
        <div class="gmail_extra">
          <div class="gmail_quote">
            <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
              <br>
              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 ?<br>
            </blockquote>
            <div><br>
            </div>
            <div>There is, but there it is still somewhat problematic. I
              use FIAT to generate the basis function tabulation,</div>
            <div>so you have to configure with</div>
            <div><br>
            </div>
            <div>  --download-fiat --download-scientificpython
              --download-generator</div>
            <div><br>
            </div>
            <div>and you need mesh generation and partitioning</div>
            <div><br>
            </div>
            <div>  --download-triangle --download-chaco</div>
            <div><br>
            </div>
            <div>and then you can run SNES ex12 using Builder (which
              will make the header file)</div>
            <div><br>
            </div>
            <div>  python2.7 ./config/builder2.py check
              src/snes/examples/tutorials/ex12.c</div>
            <div><br>
            </div>
            <div>Jed and I are working on an all C version of tabulation
              which would mean that you could bypass</div>
            <div>
              the Python code generation step. Once the header is
              generated for the element you want, then</div>
            <div>you can just run the example as normal.</div>
            <div><br>
            </div>
            <div>    Matt</div>
            <div> </div>
            <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
              Thanks a lot.<span><font color="#888888"><br>
                  <br>
                  Olivier Bonnefon<br>
                  <br>
                  -- <br>
                  Olivier Bonnefon<br>
                  INRA PACA-Avignon, Unité BioSP<br>
                  Tel: <a href="tel:%2B33%20%280%294%2032%2072%2021%2058" value="+33432722158" target="_blank">+33 (0)4 32 72
                    21 58</a><br>
                  <br>
                </font></span></blockquote>
          </div>
          <br>
          <br clear="all"><span class="HOEnZb"><font color="#888888">
          <div><br>
          </div>
          -- <br>
          What most experimenters take for granted before they begin
          their experiments is infinitely more interesting than any
          results to which their experiments lead.<br>
          -- Norbert Wiener
        </font></span></div><span class="HOEnZb"><font color="#888888">
      </font></span></div><span class="HOEnZb"><font color="#888888">
    </font></span></blockquote><span class="HOEnZb"><font color="#888888">
    <br>
    <br>
    <pre cols="72">-- 
Olivier Bonnefon
INRA PACA-Avignon, Unité BioSP
Tel: <a href="tel:%2B33%20%280%294%2032%2072%2021%2058" value="+33432722158" target="_blank">+33 (0)4 32 72 21 58</a></pre>
  </font></span></div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>