S4 (Gu et al., 2022) notices that a linear state space model is equivalent to a convolution. Unrolling the recurrence gives a kernel K whose k-th entry is C A^k B, and that kernel can be computed in parallel.
Implement ssm_kernel(A_bar, B_bar, C, L), the impulse response of the scalar system.
Raise the state transition to each lag, multiply by C and B.
No results yet
Run your code to check it against the examples, or submit to run every test.