
<html>

  <head>


    <meta http-equiv="content-type" content="text/html; charset=UTF-8">

  </head>

  <body text="#000000" bgcolor="#FFFFFF">

    <p>Hello everyone.</p>

    <p>I ran into artifacts in my gradient field when doing parallel

      computations, and to investigate I reduced<br>

      the problem down to an analytical testcase (2D heat equation with

      gaußian heat kernel). My code converges nicely to the<br>

      analytic solution. When I visualize the error field in Paraview

      There is no artifacts visible, but when visualizing the gradient

      of my field I find artifacts which look suspiciously<br>

      like some bug in overlap or parallel implementation. (Vertical

      line error in the middle of the domain<br>

      when I run on two cores).<br>

      <br>

      I write the gradient to file like this:<br>

      <br>

      <i>                using DGFG

        =Dune::PDELab::DiscreteGridFunctionGradient<GFS,RHSType>;</i><i><br>

      </i><i>                using VTKFG =

        Dune::PDELab::VTKGridFunctionAdapter<DGFG>;</i><i><br>

      </i><i>               

vtkSequenceWriter_p->addVertexData(std::make_shared<VTKFG>(std::make_shared<DGFG>(*gfs_p,

        *z_p), "gradient"));</i><i><br>

      </i></p>

    <p>After a quick look into the evaluate() method of said function it

      is not obvious to me what the error might be, as it is simply a<br>

      local evaluation of the gradients in each lfs. Is it necessary to

      do some overlap-communication before running this method?<i><br>

      </i>I work with a yaspgrid with overlap=1 on version 2.6.</p>

    <p>Thanks in advance,</p>

    <p>Michael Wenske<br>

      <i></i></p>

  </body>

</html>