• 中国计算机学会会刊
  • 中国科技核心期刊
  • 中文核心期刊

Computer Engineering & Science

Previous Articles     Next Articles

A parallel 3D poisson equation solver
based on discrete sine transform

LIN Shi-wei1,2,3,ZHANG Wei-min1,2,FANG Min-quan1,2,LI Song1,2   

  1. (1.College of Computer,National University of Defense Technology,Changsha 410073;
    2.Naval Hydrologic Meteorological Center,Beijing 100161;
    3.Academy of Ocean Science and Engineering,National University of Defense Technology,Changsha 410073,China)
     
  • Received:2017-01-11 Revised:2017-03-26 Online:2017-08-25 Published:2017-08-25

Abstract:

Poisson equations are widely applied in physical and engineering problems. Most of numerical methods for solving 3D Poisson equations have no remarkable parallelism, and use the global iteration methods which limit the computational efficiency and stability. We abandon the idea of global iteration and combine the discrete sine transform theory (DST) with 27-point four-order difference scheme to modify and parallelize the 3D Poisson equation solver at the algorithm level, which separates the whole problem to several smaller independent ones and greatly improves the stability and parallel performance. For a given discrete form, common parameters can be used to solve different Possion equations so that the programming efficiency is greatly improved. The algorithm is implemented on the basis of the shared memory parallel model. Experimental results show that  for the given instance, the new algorithm has a good acceleration effect and the final error of the computing result is about 10e-5, within an acceptable rang. And the result accuracy gets improved with the increase of problem scale.
 

Key words: 3D poisson equation, discrete sine transform, parallel, OpenMP