Speeding up forecasting
A(b) is computed two times: for forecasting and analysing. However, forecasting requires A(b) * u rather than A(b).
For now, the same routine is used for both but creating a new subroutine for the forecast would lead to a significant speedup.