A New Fast Direct Method For Solving Quasi-Toeplitz System of Equations

The objective of this study is to present a novel, efficient, and fast direct method for solving linear systems of equations whose coefficient matrix is a tridiagonal Quasi-Toeplitz matrix. Such matrices are frequently encountered in the discretization of second-order differential equation problems with Neumann boundary conditions, the discretization of Quasi-Birth-and-Death processes known as QBD matrix equations, and other related applications. The presented method has been demonstrated to produce favorable results in terms of run time for moderately large systems when compared to classical direct methods, such as the LU and PLU methods.