A Fast Direct Solution of Poisson's Equation Using Fourier Analysis
Journal of the ACM (JACM)
A Block Fourier Decomposition Method
PARA '02 Proceedings of the 6th International Conference on Applied Parallel Computing Advanced Scientific Computing
Hi-index | 0.00 |
The block Fourier decomposition method recently proposed by the first author is a special method for decoupling any block tridiagonal matrix of the form K = block-tridiag [B, A, B], where A and B are square submatrices, into diagonal blocks. Unlike the traditional fast Poisson solver, block cyclic reductions, or the FACR algorithm, this approach does not require A and B be symmetric or commute. Presented in this paper is a parallel solver using this block decomposition method to solve linear systems whose coefficient matrices are of the form of K. We describe the computational procedure and implementation for parallel executions on distributed workstations. The performance from our numerical experiments is reported to demonstrate the usefulness of this approach.