[dune-pdelab] Alexander Time stepper in PDELab

[Jared Okiro]

Dear PDELab Developers,

I have implemented a highly stiff Reaction-Diffusion system with four 
components
using the DG Weighted Interior Penalty method given in PDELab.
I have used Alexander3Parameter as the time stepper and it seems to work 
well
with my system. I now need to make this time step adaptive with error 
control.
I have seen that the matrix B does not have embedding. Can someone help 
me do this
(will appreciate some code snippet)so that I can estimate the error 
between two successive time steps.
Any advice will be highly appreciated.

Sincerely,
Jared

-- 
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