Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter March 9, 2016

Solving 3D Time-Fractional Diffusion Equations by High-Performance Parallel Computing

  • Wei Zhang EMAIL logo and Xing Cai


Numerically solving time-fractional diffusion equations, especially in three space dimensions, is a daunting computational task. This is due to the huge requirements of both computation time and memory storage. Compared with solving integer-ordered diffusion equations, the costs for time and storage both increase by a factor that equals the number of time steps involved. Aiming to overcome these two obstacles, we study in this paper three programming techniques: loop unrolling, vectorization and parallelization. For a representative numerical scheme that adopts finite differencing and explicit time integration, the performance-enhancing techniques are indeed shown to dramatically reduce the computation time, while allowing the use of many CPU cores and thereby a large amount of memory storage. Moreover, we have developed simple-to-use performance models that support our empirical findings, which are based on using up to 8192 CPU cores and 12.2 terabytes.


Wei Zhang has been funded by a grant from the China Scholarship Council. A part of the computer time used for the paper has been provided by a grant from the NOTUR project of Norway.


[1] M. Caputo, Linear models of dissipation whose Q is almost frequency independent–II.Geophys. J. Int.13 (1967), 529–539.10.1111/j.1365-246X.1967.tb02303.xSearch in Google Scholar

[2] B. Chapman, G. Jost, R. van der Pas, Using OpenMP: Portable Shared Memory Parallel Programming. MIT Press, Cambridge, MA (2007).Search in Google Scholar

[3] J. Chen, F. Liu, Q. Liu, X. Chen, V. Anh, I. Turner, K. Burrage, Numerical simulation for the three-dimension fractional sub-diffusion equation.Appl. Math. Model.38 (2014), 3695–3705.10.1016/j.apm.2014.03.031Search in Google Scholar

[4] W. Chen, L. Ye, H. Sun, Fractional diffusion equations by the Kansa method.Comput. Math. Appl.59 (2010), 1614–1620.10.1016/j.camwa.2009.08.004Search in Google Scholar

[5] K. Diethelm, An effcient parallel algorithm for the numerical solution of fractional differential equations.Fract. Calc. Appl. Anal.14 , No 3 (2012), 475–490; DOI: 10.2478/s13540-011-0029-1; in Google Scholar

[6] Z. Fu, W. Chen, H. Yang, Boundary particle method for Laplace transformed time fractional diffusion equations.J. Comput. Phys.235 (2013), 52–66.10.1016/ in Google Scholar

[7] G. Gao, Z. Sun, A compact finite difference scheme for the fractional sub-diffusion equations.J. Comput. Phys.230 (2011), 586–595.10.1016/ in Google Scholar

[8] P. Gepner, V. Gamayunov, D.L. Fraser, A compact finite difference scheme for the fractional sub-diffusion equations. Procedia Comput. Sci.4 (2011), 452–460.10.1016/j.procs.2011.04.047Search in Google Scholar

[9] S. Goedecker, A. Hoisie, Performance Optimization of Numerically Intensive Codes . SIAM, Philadelphia (2001).10.1137/1.9780898718218Search in Google Scholar

[10] C. Gong, W. Bao, G. Tang, A parallel algorithm for the Riesz fractional reaction-diffusion equation with explicit finite difference method.Fract. Calc. Appl. Anal.16 , No 3 (2013), 654–669; DOI: 10.2478/s13540-013-0041-8; in Google Scholar

[11] W. Gropp, E. Lusk, A. Skjellum, Using MPI: Portable Parallel Programming with the Message-Passing Interface . MIT Press, Cambridge, MA (1994).10.7551/mitpress/7056.001.0001Search in Google Scholar

[12] R. Hilfer, Applications of Fractional Calculus in Physics . World Scientific, Singapore (2000).10.1142/3779Search in Google Scholar

[13] Intel Intrinsics Guide; in Google Scholar

[14] C. Li, Z. Zhao, Y. Chen, Numerical approximation of nonlinear fractional differential equations with subdiffusion and superdiffusion.Comput. Math. Appl.62 (2011), 855–875.10.1016/j.camwa.2011.02.045Search in Google Scholar

[15] Y. Lin, C. Xu, Finite difference/spectral approximations for the time-fractional diffusion equation.J. Comput. Phys.225 (2007), 1533–1552.10.1016/ in Google Scholar

[16] R. Magin, Fractional calculus models of complex dynamics in biological tissues.Comput. Math. Appl.59 (2010), 1586–1593.10.1016/j.camwa.2009.08.039Search in Google Scholar

[17] R. Metzler, J. Klafter, The random walk’s guide to anomalous diffusion: a fractional dynamics approach.Phys. Rep.339 (2000), 1–77.10.1016/S0370-1573(00)00070-3Search in Google Scholar

[18] I. Podlubny, Fractional Differential Equations . Academic Press, New York (1999).Search in Google Scholar

[19] R. Stern, F. Effenberger, H. Fichtner, T. Sch¨afer, The space-fractional diffusion-advection equation: Analytical solutions and critical assessment of numerical solutions.Fract. Calc. Appl. Anal.17 , No 1 (2014), 171–190; DOI: 10.2478/s13540-014-0161-9; in Google Scholar

[20] The STREAM2 Home Page; in Google Scholar

[21] H. Sun, W. Chen, C. Li, Y. Chen, Finite difference schemes for variable-order time fractional diffusion equation. Int. J. Bifurcat. Chaos22 (2012), 1250085.10.1142/S021812741250085XSearch in Google Scholar

[22] H. Sun, Y. Zhang, W. Chen, D.M. Reeves, Use of a variable-index fractional-derivative model to capture transient dispersion in heterogeneous media.J. Contam. Hydrol.157 (2014), 47–58.10.1016/j.jconhyd.2013.11.002Search in Google Scholar PubMed

[23] D. Valério, J. Sá da Costa, Variable-order fractional derivatives and their numerical approximations.Signal Processing91 (2011), 470–483.10.1016/j.sigpro.2010.04.006Search in Google Scholar

[24] W. Wei, Effective Use of Multicore-based Parallel Computers for Scientific Computing . PhD thesis, University of Oslo (2012).Search in Google Scholar

[25] W. Zhang, W. Wei, X. Cai, Performance modeling of serial and parallel implementations of the fractional Adams-Bashforth-Moulton method. Fract. Calc. Appl. Anal.17 , No 3 (2014), 617–637; DOI: 10.2478/s13540-014-0189-x; fca.2014.17.issue-3/issue-files/fca.2014.17.issue-3.xml.10.2478/s13540-014-0189-xSearch in Google Scholar

[26] Y. Zhang, D.A. Benson, D.M. Reeves, Time and space nonlocalities underlying fractional-derivative models: Distinction and literature review of field applications.Adv. Water Resour.32 (2009), 561–581.10.1016/j.advwatres.2009.01.008Search in Google Scholar

[27] Y. Zhang, Z. Sun, Alternating direction implicit schemes for the two-dimensional fractional sub-diffusion equation.J. Comput. Phys.230 (2011), 8713–8728.10.1016/ in Google Scholar

[28] Y. Zhang, Z. Sun, X. Zhao, Compact alternating direction implicit scheme for the two-dimensional fractional diffusion-wave equation. SIAM J. Numer. Anal.50 (2012), 1535–1555.10.1137/110840959Search in Google Scholar

  1. Please cite to this paper as published in:

    Fract. Calc. Appl. Anal., Vol. 19, No 1 (2016), pp. 140–160, DOI: 10.1515/j_fca-2016-0008

Received: 2014-12-26
Revised: 2015-7-19
Published Online: 2016-3-9
Published in Print: 2016-2-1

© 2016 Diogenes Co., Sofia

Downloaded on 6.12.2023 from
Scroll to top button