Thanks Joe for the response. I want to implement Alexander3Parameter time stepping scheme.<br>Do I need to do jacobian_ volume( ) of the temporal operator then?.<br>Regards,<br>Jared<br><br><div class="gmail_quote">On Tue, Sep 20, 2011 at 12:35 PM, Jö Fahlke <span dir="ltr"><<a href="mailto:jorrit@jorrit.de">jorrit@jorrit.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Am Tue, 20. Sep 2011, 08:44:57 +0200 schrieb Jared Okiro:<br>
<div class="im">> I am a new dune user, and is developing a DG code for a reaction-diffusion<br>
> system.<br>
> When I look at the example code in<br>
> dune/pdelab/localoperator/convectiondiffusiondg.hh,<br>
> I see two ways of doing the residual assembly i.e alpha volume and jacobian<br>
> volume.<br>
> Can someone help me understand the difference between these two functions in<br>
> setting up<br>
> the system?. I am asking this because in standard Galerkin methods, I only<br>
> see alpha volume<br>
> but not jacobian volume also.<br>
<br>
</div>Well, jacobian_volume() is needed when you need to assemble and solve a<br>
matrix, i.e. for stationary or implicit instationary problems.  If you don't<br>
want to implement both of alpha_volume() and jacobian_volume() you can derive<br>
your local operator from NumericalJacobianVolume.  This implements<br>
jacobian_volume() using alpha_volume() an numerical differenciation.<br>
Alternatively, if you have a *linear* problem and don't want to implement<br>
alpha_volume(), you can derive your local operator from<br>
JacobianBasedAlphaVolume.<br>
<div class="im"><br>
> Finally, my problem is of the form:<br>
><br>
> *\frac{\partial u}{\partial t} = D \ Delta u + f(u) = 0  in \Omega,* where<br>
> f(u) is a non-linear function.<br>
><br>
> Do I need to do jacobian volume of the temporal operator also or the alpha<br>
> volume is sufficient?<br>
<br>
</div>That depends on your time-stepping scheme.<br>
<br>
Bye,<br>
Jö.<br>
<font color="#888888"><br>
--<br>
In the beginning the Universe was created.  This has made a lot of<br>
people very angry and been widely regarded as a bad move.<br>
-- Douglas Adams<br>
</font><br>-----BEGIN PGP SIGNATURE-----<br>
Version: GnuPG v1.4.10 (GNU/Linux)<br>
<br>
iQIVAwUBTnhsXumx9CKgc0OWAQgCpg//YHdlNTg511KPQ+fvHFQXAlJOj2oq2PBY<br>
UOMK+6i4erESoEfvoV1Fi0Qp/0heFE5Lp7KGMc7fAXns7i8kSmtmIgY5o2meSb4c<br>
8aZLMOZWHouALzNxLlgN1MVLWC9CFToMxEIK3TNGfIO4CC1vyd02CTDvx2B9mL+Z<br>
2aSlEjW7KTZ/iah6XVei+bLQecSKF09ea97y1ywOD2ruKd1pdJXimihW+ikLoesp<br>
L1mq+1jKrx5Gn+Y+XMwjwWVTifLqVerYAUFOhp/kleQeCcR4zVnZjC2kcsmWD8Xq<br>
eAKVQhOsmAB6IbqUndj0Yo7S/WxMXG9A12nXCkrKsEgUU/grRbesFhM41VL1ZPNb<br>
SeKUoKFcKbqrjACRr9JozyU5kdW0S1C9EIl3CQx2lBr0Mt6IbwCG/QzeOJXP4C5a<br>
gfjoBYcfN31uq2jjJaWESGTcZVzyMD2+XXRz7DPt0fgH8eHwPJovzZEfwTv+7EbW<br>
X2t/3gl18LD7mq5OijmJTtjRnVgEzXXssaZ8fq/fa7u1WXT8RJC3FDuzXXDeC7cp<br>
dmpQssmp55aP8A5DbxHUqvSdwx/QOxps9Cy7HJADRyFe7KyEvdadnby7xk7Wlvp7<br>
PtoLAMmT+0ulG89gBW3xHZTVuyypc/4fOt6RpupLMqeL2p/UGLUUdsnKHnxgV0Pr<br>
vRckb1wTOlM=<br>
=bm7o<br>
-----END PGP SIGNATURE-----<br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div>Jared Ouma Okiro,</div>
<div>Institute for Analysis and Numerics,</div>
<div>Otto-von- Guericke Universität,</div>
<div>Universitätsplatz 02, Gebaude 02 / 115<br>39106 Magdeburg, Germany<br></div>
<div>Office Tel: +493916712660<br>Fax:         +493916718073<br>Mobile:     +4917657479787<br>Home Tel: +4939155917376<br><a href="mailto:Email%3Ajared.okiro@st.ovgu.de" target="_blank">Email:jared.okiro@st.ovgu.de</a><br>
</div><br>