powered by
Compute IRF from AR coefficients and impulse coefficient
.compute_irf(ar_coefs, ar_orders, max_p, beta_impulse, H)