Performance Analysis of Effective Methods for Solving Band Matrix SLAEs after Parabolic Nonlinear PDEs
For computational scientists solving parabolic PDEs, this work provides a performance comparison of existing methods, but the improvements are incremental.
This paper experimentally compares direct, symbolic, and iterative (SIP with Hotelling-Bodewig) methods for solving band matrix systems from parabolic nonlinear PDEs, finding that the SIP method with Hotelling-Bodewig acceleration achieves up to 2x speedup over standard SIP on HPC clusters.
This paper presents an experimental performance study of implementations of three different types of algorithms for solving band matrix systems of linear algebraic equations (SLAEs) after parabolic nonlinear partial differential equations -- direct, symbolic, and iterative, the former two of which were introduced in Veneva and Ayriyan (arXiv:1710.00428v2). An iterative algorithm is presented -- the strongly implicit procedure (SIP), also known as the Stone method. This method uses the incomplete LU (ILU(0)) decomposition. An application of the Hotelling-Bodewig iterative algorithm is suggested as a replacement of the standard forward-backward substitutions. The upsides and the downsides of the SIP method are discussed. The complexity of all the investigated methods is presented. Performance analysis of the implementations is done using the high-performance computing (HPC) clusters "HybriLIT" and "Avitohol". To that purpose, the experimental setup and the results from the conducted computations on the individual computer systems are presented and discussed.