[dune-pdelab] Computing Maximum and Average solution with DG

Jared Okiro jared.okiro at st.ovgu.de
Sun May 20 00:28:47 CEST 2012


Hi Stephen,

Thanks for the response and the tip. However, I have tried to increase  
the integration order(upto 4) but still not much improvemenent.I would  
still appreciate further assistance on this.
Thank you.

Kind regards,
-- 
Jared Ouma Okiro
Institute for Analysis and Numerics
Otto von Guericke Universität
Universitätsplatz 02, Gebaude 02 / 115
39106 Magdeburg, Germany
Office Tel: +493916712660
Fax:        +493916718073
Mobile:     +4917657479787
Home Tel:   +4939155917376


Quoting Steffen Müthing <steffen.muething at ipvs.uni-stuttgart.de>:

> Content-Transfer-Encoding: quoted-printable
> Content-Type: text/plain; charset=iso-8859-1
>
> Hi Jared,
>
> I think the problem is due to the different interpolation schemes used by the
> two methods (your function and the VTKWriter, respectively):
>
> You use a first order DG scheme and calculate the integration order in your
> function as
>
> max(2*(dg_order -1),0),
>
> which yields 0 in your case, so you get a zero order integration  
> rule which only
> evaluates your solution at the cell center. For the VTKWriter, you  
> use a vertex
> data set (and the subsampling level is 0), so the VTKWriter evaluates your
> solution at the vertices. And your solution assumes its maximum only  
> at vertex
> positions...
>
> Try increasing the integration order in your function and see if that helps.
>
> Steffen
>
>
>
> Am 19.05.2012 um 18:58 schrieb Jared Okiro:
>
>> <outputtest.cc>
>
> Steffen Müthing
> Universität Stuttgart
> Institut für Parallele und Verteilte Systeme
> Universitätsstr. 38
> 70569 Stuttgart
> Tel: +49 711 685 88429
> Fax: +49 711 685 88340
> Email: steffen.muething at ipvs.uni-stuttgart.de
>
>
>
>
>
>






More information about the dune-pdelab mailing list