allpole filter…
This runs a block-computing version of the allpole filter. Good luck.
From the source code:
float * ballpole(float *x, long *jcount, long npoles, float *past, float *c, float *out, long nvals)