Block Levinson solver

BLOCK_LEVINSON(Y, L) solves the matrix equation T * x = y, where T is a symmetric matrix with block Toeplitz structure, and returns the solution vector x. The matrix T is never stored in full (because it is large and mostly redundant), so the...