improved speed of ERK methods forward sensitivity propagation
added AD support for forward and backward derivatives with matrix seeds
improved forward sensitivity propagation in ERK methods using the new AD functionality
improved forward-over-backward second order sensitivity propagation in ERK methods using the new AD functionality
added custom functionionality for the new three-sweeps-propagation approach to make it at least competitive with forward-over-backward, currently lacking specific AD functionality for the operators involved