Многотомное издание: Lecture Notes in Computer Science
Страницы: 388-401
Аннотация
In this paper, we present an original preconditioner to solve Poisson equation for strongly heterogeneous media. We suggest using the Conjugate Gradient method with the preconditioner based on the solution of the Poisson equation for homogeneous media. Corresponding operator is easy to invert by spectral method, where spectral decomposition is applied in two spatial directions and the Gauss elimination method is applied to solve a series of 1D problems. We illustrate that use of such precondtioner strongly decreases the number of iterations to solve the original Poisson equation, moreover the number of iterations weakly depends on the problem size. Implementation of the suggested approach using modern GPUs allows solving problems of up to the size of 10003 voxels.