Parallel implementation of 2-dimensional Toeplitz solver on MasPar with applications to image restoration

Image restoration problems can be transformed into problems of solving a linear system Tx=b where T is a block-Toeplitz or near-block-Toeplitz matrix. However for many of these problems, the size of the matrix T is usually very large. For instance, if one is to manipulate a 1024-by-1024 pixel image, then T will be of the order 1024 ²-by-1024². The authors implement a parallel version of their existing 2-dimensional Toeplitz solver in a data-parallel fashion on MasPar (DECmpp) by fully utilizing its massively parallel processing power. The package they developed is portable and easy-to-use. They demonstrate how to use the package to recover a satellite image which is blurred by atmospheric turbulence. The implementation details and performance results are also presented