TL;DR: This paper presents a division-free algorithm for evaluating the determinants of a class of quasi-tridiagonal matrices that can be viewed as perturbations of general tridi diagonal matrices and its competitiveness with other related algorithms and MATLAB built-in function.
TL;DR: In this article, a hybrid multigrid-Thomas algorithm is proposed to solve Poisson's equation as part of the spatial discretization of a time evolving PDE system.
TL;DR: An efficient numerical algorithm for evaluating the determinants of general bordered k-tridiagonal matrices in linear time based on a novel incomplete block-diagonalization approach which preserves the low-rank structure and sparsity of the original matrix and its competitiveness with Gaussian elimination algorithm and MATLAB built-in function.
TL;DR: In this paper, it was shown that the closure of the numerical range of an n + 1 -periodic tridiagonal operator is equal to the closed range of a 2(n + 1 ) × 2 (n+ 1 ) complex matrix.
TL;DR: In this paper , it was shown that the cp-rank of any completely positive irreducible tridiagonal doubly stochastic matrix is equal to its rank, and that the same is true for symmetric pentadiagonal matrices.
Abstract: Abstract We provide a decomposition that is sufficient in showing when a symmetric tridiagonal matrix A A is completely positive. Our decomposition can be applied to a wide range of matrices. We give alternate proofs for a number of related results found in the literature in a simple, straightforward manner. We show that the cp-rank of any completely positive irreducible tridiagonal doubly stochastic matrix is equal to its rank. We then consider symmetric pentadiagonal matrices, proving some analogous results and providing two different decompositions sufficient for complete positivity. We illustrate our constructions with a number of examples.
TL;DR: This paper presents an efficient algorithm for numerically computing the inverses of n -square ( p, q )-tridiagonal matrices under a certain condition, based on a bidiagonalization approach which preserves the banded structure and sparsity of the original matrix.
TL;DR: In this paper , an algorithm for inverting complex tridiagonal Hermitian matrices with optimal computational efforts is presented, which leads to closed-form expressions for the elements of inverse matrices.
Abstract: In this paper we give an algorithm for inverting complex tridiagonal Hermitian matrices with optimal computational efforts. For matrices of a special form and, in particular, for Toeplitz matrices, the derived formulas lead to closed-form expressions for the elements of inverse matrices.
TL;DR: This paper first analyzes the potential parallel process of solving tridiagonal sparse linear equations by Gaussian elimination method and matrix splitting method, and designs and implements these two parallel algorithms to solve sparselinear equations.
Abstract: There are many practical problems in real life, which are finally attributed to solving large sparse linear equations. In order to solve sparse linear equations in parallel, this paper first analyzes the potential parallel process of solving tridiagonal sparse linear equations by Gaussian elimination method and matrix splitting method, and designs and implements these two parallel algorithms to solve sparse linear equations, Through different process tests and performance analysis, it shows that the Gaussian elimination method has poor time performance, while the matrix splitting method has good efficiency in both space and time.
TL;DR: A method for computing the positive integer powers of the anti-tridiagonal matrix corresponding to these matrices corresponding to two certain types of tridiagonal matrices of arbitrary order is presented.
Abstract: In this paper, we firstly derive a general expression for the entries of the
m
th (
m
∈
ℕ
) power for two certain types of tridiagonal matrices of arbitrary order. Secondly, we present a method for computing the positive integer powers of the anti-tridiagonal matrix corresponding to these matrices. Also, we give Maple 18 procedures in order to verify our calculations.
TL;DR: In this article , the forward error analysis of the extended cyclic reduction algorithm for solving the block-tridiagonal linear system is studied, and the critical point is to find out that the zeros of matrix polynomial $B_{i}^{(r)}$ are eigenvalues of a principal submatrix of the coefficient matrix.
Abstract: The extended cyclic reduction algorithm developed by Swarztrauber in 1974 was used to solve the block-tridiagonal linear system. The paper fills in the gap of theoretical results concerning the zeros of matrix polynomial $B_{i}^{(r)}$ with respect to a tridiagonal matrix which are computed by Newton's method in the extended cyclic reduction algorithm. Meanwhile, the forward error analysis of the extended cyclic reduction algorithm for solving the block-tridiagonal system is studied. To achieve the two aims, the critical point is to find out that the zeros of matrix polynomial $B_{i}^{(r)}$ are eigenvalues of a principal submatrix of the coefficient matrix.
TL;DR: In this paper , a new efficient and reliable hybrid numerical algorithm for evaluating general n-th order tridiagonal determinants in linear time is presented. But the algorithm is suited for implementation using computer languages such as FORTRAN, PASCAL, ALGOL, MAPLE, MACSYMA and MATHEMATICA.
Abstract: The principal minors of a tridiagonal matrix satisfy two-term and three-term recurrences [1, 2]. Based on these facts, the current article presents a new efficient and reliable hybrid numerical algorithm for evaluating general n-th order tridiagonal determinants in linear time. The hybrid numerical algorithm avoid all symbolic computations. The algorithm is suited for implementation using computer languages such as FORTRAN, PASCAL, ALGOL, MAPLE, MACSYMA and MATHEMATICA. Some illustrative examples are given. Test results indicate the superiority of the hybrid numerical algorithm.
TL;DR: In this article , the inverse eigenvalue problem of a bordered anti-tridiagonal matrix is studied and the uniqueness of its solution is discussed. But the solution is based on the recursive expression of the solution.
Abstract: The inverse matrix problem is a hot and active research topic in computational mathematics[1]. It has broad applications in engineering and scientific calculation, and owns a strong background in physics and practical significance[2]. This paper explores the inverse eigenvalue problem of a bordered anti-tridiagonal matrix. It first illustrates the existence and the uniqueness of its solution, the elaborates on the recursive expression of the solution and uses one numerical example to show the effectiveness of the algorithm, and finally concludes that this work is significant and points out suggestions for further study.
TL;DR: In this paper , the authors show how some standard general results can be used to uncover the spectral theory of tridiagonal and related matrices more elegantly and simply than existing approaches.
Abstract: Abstract This paper aims to show how some standard general results can be used to uncover the spectral theory of tridiagonal and related matrices more elegantly and simply than existing approaches. As a typical example, we apply the theory to the special tridiagonal matrices in recent papers on orthogonal polynomials arising from Jordan blocks. Consequently, we find that the polynomials and spectral theory of the special matrices are expressible in terms of the Chebyshev polynomials of second kind, whose properties yield interesting results. For special cases, we obtain results in terms of the Fibonacci numbers and Legendre polynomials.
TL;DR: In this article , a fast algorithm for solving the special tridiagonal quasi-toeplitz system is presented where the bandwidth of a quasi-Toeplit is larger than the one of Toeplite.
Abstract: Abstract In this paper, a fast algorithm for solving the special tridiagonal quasi-Toeplitz system is presented where the bandwidth of a quasi-Toeplitz is larger than the one of Toeplitz. Our algorithm is quite competitive with the classic LU method. Some examples demonstrate the good efficiency and stability of our algorithm.
TL;DR: A tridiagonal system of N equations with N unknowns is a system of simultaneous algebraic equations with nonzero coefficients only on the main diagonal, the lower diagonal, and the upper diagonal as discussed by the authors .
Abstract: A system of simultaneous algebraic equations with nonzero coefficients only on the main diagonal, the lower diagonal, and the upper diagonal is called a tridiagonal system of equations. Consider a tridiagonal system of N equations with N unknowns, u1, u2, u3,… uN as follows:
TL;DR: In this article , a hybrid multigrid-Thomas algorithm was proposed to solve Poisson's equation as part of the spatial discretization of a time-evolving PDE system.
TL;DR: In this article , it was shown that the closure of the numerical range of an n + 1 -periodic tridiagonal operator is the same as the closed range of a 2 (n + 1 ) × 2 ( n+ 1 ) complex matrix.
Abstract: We present some old and new results in the enumeration of random walks in one dimension, mostly developed in works of enumerative combinatorics. The relation between the trace of the $n$-th power of a tridiagonal matrix and the enumeration of weighted paths of $n$ steps allows an easier combinatorial enumeration of the paths. It also seems promising for the theory of tridiagonal random matrices .
TL;DR: A parallel algorithm for solving tridiagonal systems of linear equations is proposed, utilizing a two-step approach with reduced communication interactions, and demonstrated to be efficient on both common- and distributed-memory supercomputers through theoretical estimates and computational experiments.
Abstract: A parallel algorithm for solving a series of matrix equations with a constant tridiagonal matrix and different right-hand sides is proposed and studied. The process of solving the problem is represented in two steps. The first preliminary step is fixing some rows of the inverse matrix of SLAEs. The second step consists in calculating solutions for all right-hand sides. For reducing the communication interactions, based on the formulated and proved main parallel sweep theorem, we propose an original algorithm for calculating share components of the solution vector. Theoretical estimates validating the efficiency of the approach for both the common- and distributed-memory supercomputers are obtained. Direct and iterative methods of solving a 2D Poisson equation, which include procedures of tridiagonal matrix inversion, are realized using the mpi technology. Results of computational experiments on a multicomputer demonstrate a high efficiency and scalability of the parallel sweep algorithm.