<html>
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
  </head>
  <body>
    Dear Theresa,<br>
    <br>
    The standard explicit RK4 method is implemented in the class <span
      class="keyword"></span><a class="code"
href="https://www.dune-project.org/doxygen/pdelab/master/classDune_1_1PDELab_1_1RK4Parameter.html">RK4Parameter</a>
    <br>
    in PDELab (it merely contains<br>
    the entries of the Butcher Tableau, so you could easily add other
    methods as well).<br>
    Sorry for answering this here, as it is actually the dune list and
    not a pdelab list.<br>
    <br>
    Periodic boundary conditions on YaspGrid need to be implemented "in
    a parallel setting",<br>
    i.e. halo cells are added in the periodic directions and a
    communication steps transports<br>
    the solution from boundary to boundary. But I did not try this for a
    long time ...<br>
    <br>
    Which disretization scheme did you intend to use?<br>
    <br>
    Best regards,<br>
    <br>
    Peter<br>
    <br>
    <div class="moz-cite-prefix">Am 04.12.20 um 16:11 schrieb Theresa
      Pollinger:<br>
    </div>
    <blockquote type="cite"
      cite="mid:8c1b827d-aa33-65da-e8df-2ec5ff594851@ipvs.uni-stuttgart.de">
      <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
      <p>Dear DUNE developers,</p>
      <p>as you may know, we are using DUNE to run FEM with the sparse
        grid combination technique, via coupling to our code DisCoTec
        [0].</p>
      <p>Now for another project, we were simulating a pure advection
        problem with another FEM code, which returned rather unexpected
        results. Now I wanted to turn the exact same advection problem
        into a DUNE problem to validate the results, but it turns out
        that this is quite hard. Two things I am struggling with:</p>
      <p>- how do you set periodic boundary conditions in all
        directions? So far I only found this old answer [1], and a test
        using YaspGrid looks very much like periodicity is only set in
        one dimension, even if I supply an all-true argument `periodic`.</p>
      <p>- What is the proper way to use a Runge-Kutta 4th-order scheme
        in DUNE? AFAICT, pdelab only implements up to order 3 RK methods
        [2]?<br>
      </p>
      <p>Maybe you could let me know if there are people who have
        already implemented those things.</p>
      <p>Best,</p>
      <p>Theresa<br>
      </p>
      <p><br>
      </p>
      <p>[0] <a class="moz-txt-link-freetext"
          href="https://github.com/SGpp/DisCoTec" moz-do-not-send="true">https://github.com/SGpp/DisCoTec</a></p>
      <p>[1] <a class="theme markdown__link"
href="https://lists.dune-project.org/pipermail/dune-pdelab/2011-October/000212.html"
          rel="noreferrer" target="_blank" moz-do-not-send="true">https://lists.dune-project.org/pipermail/dune-pdelab/2011-October/000212.html</a></p>
      <p>[2]
        <a class="moz-txt-link-freetext"
href="https://github.com/dune-project/dune-pdelab/blob/master/dune/pdelab/instationary/onestepparameter.hh"
          moz-do-not-send="true">https://github.com/dune-project/dune-pdelab/blob/master/dune/pdelab/instationary/onestepparameter.hh</a><br>
      </p>
      <pre class="moz-signature" cols="72">-- 

Theresa Pollinger

Universität Stuttgart
Institut für Parallele und Verteilte Systeme
Simulation Software Engineering

Universitätsstr. 38
70569 Stuttgart

Phone: +49 711 685 60921</pre>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <pre class="moz-quote-pre" wrap="">_______________________________________________
Dune mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Dune@lists.dune-project.org">Dune@lists.dune-project.org</a>
<a class="moz-txt-link-freetext" href="https://lists.dune-project.org/mailman/listinfo/dune">https://lists.dune-project.org/mailman/listinfo/dune</a></pre>
    </blockquote>
    <br>
  </body>
</html>