Systolic Array For Solving Toeplitz Systems Of Equations

Many problems of geophysics, image processing and time series analysis involve the problem of solving Toeplitz systems of equations. We present a fast parallel 0 (mn) algorithm that solves both square and over-determined Toeplitz systems of equations. The solution is obtained directly from the triangular factorization without using back-substitution. This avoids separate factorization and back-substitution sections, which complicate architectural implementation. This also enables us to eliminate intermediate memory to store the triangular factor. The parallel implementation is carried out in two steps. First, Regular Iterative Algorithms (RIAs) for solving Toeplitz system of equations are formulated systematically from the mathematical description of our algorithm. The advantage of having RIAs is that the process of mapping the algorithms on regular processor arrays can be done in a systematic manner.