GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion

Farquhar, Megan E., Moroney, Timothy J., Yang, Qianqian, & Turner, Ian W. (2016) GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion. SIAM Journal on Scientific Computing, 38(3), C127-C149.

[img] Submitted Version (PDF 1MB)
Administrators only | Request a copy from author

View at publisher


The efficient computation of matrix function vector products has become an important area of research in recent times, driven in particular by two important applications: the numerical solution of fractional partial differential equations and the integration of large systems of ordinary differential equations. In this work we consider a problem that combines these two applications, in the form of a numerical solution algorithm for fractional reaction diffusion equations that after spatial discretisation, is advanced in time using the exponential Euler method. We focus on the efficient implementation of the algorithm on Graphics Processing Units (GPU), as we wish to make use of the increased computational power available with this hardware. We compute the matrix function vector products using the contour integration method in [N. Hale, N. Higham, and L. Trefethen. Computing Aα, log(A), and related matrix functions by contour integrals. SIAM J. Numer. Anal., 46(5):2505–2523, 2008]. Multiple levels of preconditioning are applied to reduce the GPU memory footprint and to further accelerate convergence. We also derive an error bound for the convergence of the contour integral method that allows us to pre-determine the appropriate number of quadrature points. Results are presented that demonstrate the effectiveness of the method for large two-dimensional problems, showing a speedup of more than an order of magnitude compared to a CPU-only implementation.

Impact and interest:

1 citations in Scopus
Search Google Scholar™

Citation counts are sourced monthly from Scopus and Web of Science® citation databases.

These databases contain citations from different subsets of available publications and different time periods and thus the citation count from each is usually different. Some works are not in either database and no count is displayed. Scopus includes citations from articles published in 1996 onwards, and Web of Science® generally from 1980 onwards.

Citations counts from the Google Scholar™ indexing service can be viewed at the linked Google Scholar™ search.

ID Code: 84722
Item Type: Journal Article
Refereed: Yes
Additional URLs:
Keywords: Matrix functions, Contour integral method, fractional reaction-diffusion equation, fractional Laplacian, Exponential Euler Method, Krylov subspace methods
DOI: 10.1137/15M1021672
ISSN: 1095-7197
Divisions: Current > Research Centres > ARC Centre of Excellence for Mathematical & Statistical Frontiers (ACEMS)
Current > Schools > School of Mathematical Sciences
Current > QUT Faculties and Divisions > Science & Engineering Faculty
Copyright Owner: Copyright 2015 Society for Industrial and Applied Mathematics
Deposited On: 14 Jun 2015 23:02
Last Modified: 03 Jun 2016 09:07

Export: EndNote | Dublin Core | BibTeX

Repository Staff Only: item control page