An efficient parallel solution for Caputo fractional reactiondiffusion equation
The computational complexity of Caputo fractional reactiondiffusion equation is \(O(MN^2)\) compared with \(O(MN)\) of traditional reactiondiffusion equation, where \(M\), \(N\) are the number of time steps and grid points. A efficient parallel solution for Caputo fractional reactiondiffusion equation with explicit difference method is proposed. The parallel solution, which is implemented with MPI parallel 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 sum of constant vector multiplication is optimized. As to the authors knowledge, this is the first parallel solution for Caputo fractional reactiondiffusion equation. The experimental results show that the parallel solution compares well with the analytic solution. The parallel solution on single Intel Xeon X5540 CPU runs more than three times faster than the serial solution on single X5540 CPU core, and scales quite well on a distributed memory cluster system.
fractional differential equation reaction-diffusion equation high performance computing parallel computing finite difference method