Closed config-i1 closed 4 years ago
This comes ultimately to the function polynomialiser(), which creates polynomials based on the provided ARIMA orders and parameters.
A possible solution here is to rewrite the function in C++
Having done several tests, the slowest part is the polyprod, which is already extremely fast. The second slowest thing is the eigen(), which cannot be improved. A possible solution for ARIMA is to either use bounds="admissible", which is faster or bounds="none", which is dangerous, but very fast.
Nothing else can be done at the moment.
This involves: