A GPU parallelized spectral method for elliptic equations in rectangular domains

We design and implement a polynomial-based spectral method on graphic processing units (GPUs). The key to success lies in the seamless integration of the matrix diagonalization technique and the new generation CUDA tools. The method is applicable to elliptic equations in rectangular domains with general boundary condition. We show remarkable speedups of up to 15 times in the 2-D case and more than 35 times in the 3-D case.

[1]  David Gottlieb,et al.  The Spectrum of the Chebyshev Collocation Operator for the Heat Equation , 1983 .

[2]  Christoph H. Keitel,et al.  A real space split operator method for the Klein-Gordon equation , 2009, J. Comput. Phys..

[3]  Christoph H. Keitel,et al.  Accelerating the Fourier split operator method via graphics processing units , 2010, Comput. Phys. Commun..

[4]  John R. Rice,et al.  Direct solution of partial difference equations by tensor product methods , 1964 .

[5]  Michel Deville,et al.  Chebyshev 3-D spectral and 2-D pseudospectral solvers for the Helmholtz equation , 1984 .

[6]  Xi Shao,et al.  Efficient spectral and pseudospectral algorithms for 3D simulations of whistler-mode waves in a plasma , 2011, J. Comput. Phys..

[7]  Kenneth R. Jackson,et al.  A Parallel Implementation on GPUs of ADI Finite Difference Methods for Parabolic PDEs with Applications in Finance , 2010 .

[8]  Jie Shen,et al.  Numerical simulations of jet pinching-off and drop formation using an energetic variational phase-field method , 2006, J. Comput. Phys..

[9]  Jie Shen,et al.  A phase field model for the mixture of two incompressible fluids and its approximation by a Fourier-spectral method , 2003 .

[10]  Gordon Erlebacher,et al.  High-order finite-element seismic wave propagation modeling with MPI on a large GPU cluster , 2010, J. Comput. Phys..

[11]  Matthew G. Knepley,et al.  Finite Element Integration on GPUs , 2013, TOMS.

[12]  Timothy C. Warburton,et al.  Nodal discontinuous Galerkin methods on graphics processors , 2009, J. Comput. Phys..

[13]  Jie Shen,et al.  Spectral and High-Order Methods with Applications , 2006 .

[14]  Jie Shen,et al.  Spectral Methods: Algorithms, Analysis and Applications , 2011 .

[15]  Jie Shen,et al.  A Phase-Field Model and Its Numerical Approximation for Two-Phase Incompressible Flows with Different Densities and Viscosities , 2010, SIAM J. Sci. Comput..

[16]  T. Takahashi,et al.  GPU‐accelerated boundary element method for Helmholtz' equation in three dimensions , 2009 .

[17]  Jie Shen,et al.  Efficient Spectral-Galerkin Method II. Direct Solvers of Second- and Fourth-Order Equations Using Chebyshev Polynomials , 1995, SIAM J. Sci. Comput..

[18]  Ramani Duraiswami,et al.  Fast multipole methods on graphics processors , 2008, J. Comput. Phys..

[19]  Jie Shen,et al.  An overview of projection methods for incompressible flows , 2006 .

[20]  Xiaofeng Yang,et al.  Numerical approximations of Allen-Cahn and Cahn-Hilliard equations , 2010 .

[21]  Jie Shen,et al.  Efficient Spectral-Galerkin Method I. Direct Solvers of Second- and Fourth-Order Equations Using Legendre Polynomials , 1994, SIAM J. Sci. Comput..