bdf2 refactored
control of matrix recalculation implemented.
Works for higher order bdf.
For fixed step 2x speedup (probably should be improved), for adaptive cca 20% (due to keeping one weakform and calling set_ext on it)
times may differ for other examples