Block pivoting implementation of a symmetric Toeplitz solver ?
We improve the solution of symmetric Toeplitz linear systems in multicore systems. We transform the Toeplitz matrix into a Cauchy-like one to obtain some benefits. The problem is partitioned into two half-sized independent problems. We use partial local pivoting to improve the accuracy of the solution. We propose a special scheme to store data in memory that accelerates the algorithm.
pivoting linear systems multicores symmetric toeplitz matrices displacement structure