Dear PDELab Members,<br><br>I am a new dune user, and is developing a DG code for a reaction-diffusion system.<br>When I look at the example code in dune/pdelab/localoperator/convectiondiffusiondg.hh,<br>I see two ways of doing the residual assembly i.e alpha volume and jacobian volume.<br>
Can someone help me understand the difference between these two functions in setting up<br>the system?. I am asking this because in standard Galerkin methods, I only see alpha volume<br>but not jacobian volume also.<br>Finally, my problem is of the form: <br>
<br clear="all"><b><i>\frac{\partial u}{\partial t} = D \ Delta u + f(u) = 0  in \Omega,</i></b> where f(u) is a non-linear function.<br><br>Do I need to do jacobian volume of the temporal operator also or the alpha volume is sufficient?<br>
Kindly advice so that I don't make a mistake in my code.<br>Kind regards,<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>