Cluster computing for the large scale discrete fractional Cable equation

Abstract This paper presents a numerical simulation technique for the fractional Cable equation in large scale domain. Special attention is given to the parallel execution of the fractional weighted average finite difference method (FWA-FDM) on distributed system with explicit message passing, where the fractional derivative is defined in Riemann–Liouville sense. The resultant huge system of equations is studied using precondition conjugate gradient method (PCG), with the implementation of cluster computing on it. The proposed approach fulfills the suitability for the implementation on Linux PC cluster through the minimization of inter-process communication. To examine the efficiency and accuracy of the proposed method, numerical test experiments using different number of the Linux PC cluster nodes are studied. The performance metrics clearly show the benefit of using the proposed approach on the Linux PC cluster in terms of execution time reduction and speedup with respect to the sequential running in a single PC.

[1]  K. Jacobson,et al.  Structural mosaicism on the submicron scale in the plasma membrane. , 1998, Biophysical journal.

[2]  C. Lubich.  Discretized fractional calculus , 1986 .

[3]  J. Shewchuk.  An Introduction to the Conjugate Gradient Method Without the Agonizing Pain , 1994 .

[4]  M. Benzi.  Preconditioning techniques for large linear systems: a survey , 2002 .

[5]  W. Webb,et al.  Constrained diffusion or immobile fraction on cell surfaces: a new interpretation. , 1996, Biophysical journal.

[6]  Eduardo Cuesta,et al.  A Numerical Method for an Integro-Differential Equation with Memory in Banach Spaces: Qualitative Properties , 2003, SIAM J. Numer. Anal..

[7]  W. Rall.  Core Conductor Theory and Cable Properties of Neurons , 2011 .

[8]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[9]  W. Webb,et al.  Measurement of molecular diffusion in solution by multiphoton fluorescence photobleaching recovery. , 1999, Biophysical journal.

[10]  I. Turner,et al.  Stability and convergence of two new implicit numerical methods for the fractional cable equation , 2009 .

[11]  K. Morton,et al.  Numerical Solution of Partial Differential Equations , 1995 .

[12]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[13]  S. B. Yuste,et al.  An Explicit Numerical Method for the Fractional Cable Equation , 2011 .

[14]  N. H. Sweilam,et al.  On the parallel iterative finite difference algorithm for 2-D Poisson's equation with MPI cluster , 2012, 2012 8th International Conference on Informatics and Systems (INFOS).

[15]  A. M. Nagy,et al.  Numerical solution of two-sided space-fractional wave equation using finite difference method , 2011, J. Comput. Appl. Math..

[16]  S. Wearne,et al.  Fractional cable equation models for anomalous electrodiffusion in nerve cells: infinite domain solutions , 2009, Journal of mathematical biology.

[17]  Thomas Rauber,et al.  Parallel Programming: for Multicore and Cluster Systems , 2010, Parallel Programming, 3rd Ed..

[18]  S. Wearne,et al.  Fractional cable models for spiny neuronal dendrites. , 2008, Physical review letters.

[19]  I. Podlubny.  Fractional differential equations , 1998 .

[20]  I. Turner,et al.  Time fractional advection-dispersion equation , 2003 .

[21]  I. Turner,et al.  Two New Implicit Numerical Methods for the Fractional Cable Equation , 2011 .

[22]  W. Webb,et al.  Automated detection and tracking of individual and clustered cell surface low density lipoprotein receptor molecules. , 1994, Biophysical journal.