Solving the Caputo Fractional Reaction-Diffusion Equation on GPU

We present a parallel GPU solution of the Caputo fractional reaction-diffusion equation in one spatial dimension with explicit finite difference approximation. The parallel solution, which is implemented with CUDA programming model, consists of three procedures: preprocessing, parallel solver, and postprocessing. The parallel solver involves the parallel tridiagonal matrix vector multiplication, vector-vector addition, and constant vector multiplication. The most time consuming loop of vector-vector addition and constant vector multiplication is optimized and impressive performance improvement is got. The experimental results show that the GPU solution compares well with the exact solution. The optimized GPU solution on NVIDIA Quadro FX 5800 is 2.26 times faster than the optimized parallel CPU solution on multicore Intel Xeon E5540 CPU.

[1]  Chunye Gong,et al.  A Parallel Algorithm for the Two-Dimensional Time Fractional Diffusion Equation with Implicit Difference Method , 2014, TheScientificWorldJournal.

[2]  Xiaoyun Jiang,et al.  Analytical solutions to time-fractional partial differential equations in a two-dimensional multilayer annulus , 2012 .

[3]  Rabha W. Ibrahim,et al.  Ulam-Hyers Stability for Cauchy Fractional Differential Equation in the Unit Disk , 2012 .

[4]  A. M. Mathai,et al.  Solution of Generalized Fractional Reaction-Diffusion Equations , 2006 .

[5]  V. Gafiychuk,et al.  Mathematical modeling of time fractional reaction-diffusion systems , 2008 .

[6]  Rui Zhang,et al.  A Three-Stage Optimization Algorithm for the Stochastic Parallel Machine Scheduling Problem with Adjustable Production Rates , 2013 .

[7]  Xiaolin Cao,et al.  JASMIN: a parallel software infrastructure for scientific computing , 2010, Frontiers of Computer Science in China.

[8]  R. Agarwal,et al.  Stability of functional equations in single variable , 2003 .

[9]  Saad Zagloul Rida,et al.  On the solutions of time-fractional reaction–diffusion equations , 2010 .

[10]  Janusz Brzdęk,et al.  Ulam's Type Stability , 2012 .

[11]  Janusz Brzdęk,et al.  On Some Recent Developments in Ulam's Type Stability , 2012 .

[12]  Fawang Liu,et al.  Numerical simulation for the 3D seepage flow with fractional derivatives in porous media , 2008 .

[13]  Allaberen Ashyralyev,et al.  On the Numerical Solution of Fractional Parabolic Partial Differential Equations with the Dirichlet Condition , 2012 .

[14]  Stephen A. Jarvis,et al.  On the Acceleration of Wavefront Applications using Distributed Many-Core Architectures , 2012, Comput. J..

[15]  Chunye Gong,et al.  An efficient parallel solution for Caputo fractional reaction–diffusion equation , 2014, The Journal of Supercomputing.

[16]  S. Wearne,et al.  Fractional Reaction-Diffusion , 2000 .

[17]  Arak M. Mathai,et al.  Further solutions of fractional reaction-diffusion equations in terms of the H-function , 2007, J. Comput. Appl. Math..

[18]  O. P. Singh,et al.  An analytic algorithm for the space-time fractional advection-dispersion equation , 2011, Comput. Phys. Commun..

[19]  Dorian Popa,et al.  Hyers-Ulam stability of the linear recurrence with constant coefficients , 2005 .

[20]  Fawang Liu,et al.  The time fractional diffusion equation and the advection-dispersion equation , 2005, The ANZIAM Journal.

[21]  Haowei Huang,et al.  GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method , 2011, J. Comput. Phys..

[22]  Tao Tang,et al.  Exploiting hierarchy parallelism for molecular dynamics on a petascale heterogeneous system , 2013, J. Parallel Distributed Comput..

[23]  Allaberen Ashyralyev,et al.  Finite Difference and Iteration Methods for Fractional Hyperbolic Partial Differential Equations with the Neumann Condition , 2012 .

[24]  Yufeng Xu,et al.  The short memory principle for solving Abel differential equation of fractional order , 2011, Comput. Math. Appl..

[25]  Fanhai Zeng,et al.  Spectral approximations to the fractional integral and derivative , 2012 .

[26]  Chunye Gong,et al.  A parallel algorithm for the Riesz fractional reaction-diffusion equation with explicit finite difference method , 2013 .

[27]  Wei Liu,et al.  Efficient parallel implementation of large scale 3D structured grid CFD applications on the Tianhe-1A supercomputer , 2013 .

[28]  Haowei Huang,et al.  Particle transport with unstructured grid on GPU , 2012, Comput. Phys. Commun..

[29]  Chunye Gong,et al.  A Domain Decomposition Method for Time Fractional Reaction-Diffusion Equation , 2014, TheScientificWorldJournal.

[30]  Kai Diethelm,et al.  An efficient parallel algorithm for the numerical solution of fractional differential equations , 2011 .

[31]  Ajda Fošner,et al.  On the generalized Hyers–Ulam stability of module left (m, n)-derivations , 2012 .

[32]  Yong-Soo Jung,et al.  On the generalized Hyers-Ulam stability of module left derivations , 2008 .

[33]  Chen Jing-hua,et al.  An Implicit Approximation for the Caputo Fractional Reaction-Dispersion Equation , 2007 .

[34]  Yi Jiang,et al.  Parallelizing a High-Order CFD Software for 3D, Multi-block, Structural Grids on the TianHe-1A Supercomputer , 2013, ISC.

[35]  Samuel Williams,et al.  Optimization of sparse matrix-vector multiplication on emerging multicore platforms , 2009, Parallel Comput..

[36]  Kai Lu,et al.  The TianHe-1A Supercomputer: Its Hardware and Software , 2011, Journal of Computer Science and Technology.