Documentation of TimeSteppingParameters
I just had a look at the time stepping methods in PDELab. Unfortunately they are not documented in onestepparameter.hh.
The parameters are two matrices a, b\in \mathbb{R}^{s\times s+1}
and a vector d\in\mathbb{R}^{s+1}
.
Reverse engineering led me to somewhat like
\begin{aligned}
u^{(0)} &= u_n\\
\sum_{j=0}^s a_{ij}u^{(j)} &= h\sum_{j=0}^s b_{ij}f(t_n+hd_j, u^{(j)}) \qquad \forall i=0,\ldots,s-1\\
u_{n+1} &= u^{(s)}.
\end{aligned}
Is that correct? Do I overlooked some documentation somewhere? If not, I could prepare a MR that adds the documentation in onestepparameter.hh.
Edited by Nils-Arne Dreier