The performance of seasonal ARIMA in PR #1586 is limited by the initialization of the Kalman filter (in particular the matrix P). Solving the Lyapunov/Sylvester equation is very expensive.
We can speed it up with the Bartels-Stewart algorithm if we manage to make a good implementation of the batched real Schur decomposition. Iterative methods to approximate the solution are also a possibility. More info on: https://en.wikipedia.org/wiki/Bartels%E2%80%93Stewart_algorithm
The performance of seasonal ARIMA in PR #1586 is limited by the initialization of the Kalman filter (in particular the matrix
P
). Solving the Lyapunov/Sylvester equation is very expensive.We can speed it up with the Bartels-Stewart algorithm if we manage to make a good implementation of the batched real Schur decomposition. Iterative methods to approximate the solution are also a possibility. More info on: https://en.wikipedia.org/wiki/Bartels%E2%80%93Stewart_algorithm