LU decomposition of a tridiagonal matrix

Special matrices can be handled even more efficiently. A tri-diagonal matrix is one with non-zero entries along the main diagonal, and one diagonal above and below the main one (see the figure). The set of equations can then be written as

Structure of Tri-diagonal Matrix

The LU decomposition algorithm for solving this set is

The number of multiplications and divisions for a problem with n unknowns and m right-hand sides is

Operation count = 2 (n - 1) + M (3 n - 2 )

If

no pivoting is necessary. For solving two-point boundary value problems and partial differential equations this is often the case.

Structure of Block Tri-diagonal Matrix

A block tri-diagonal matrix (see the figure) frequently arises when solving multiple differential equations, either two-point boundary value problems or partial differential equations. When solving several simultaneous equations it is important to arrange the unknowns in a particular way to achieve the block tri-diagonal form of the matrix and the resulting efficiencies. Suppose we are solving for two unknowns, say temperature and concentration, and that we have a value for each at each node, i=1 to n. Thus we have unknowns {c1,c2,...cn, T1,T2,...,Tn}. If we order the unknowns in this way, however, the resulting linear equations we need to solve when solving two-point boundary value problems or partial differential equations will be dense, that is the non-zero entries of the matrix A will have no special pattern. However, if we arrange the unknowns in the pattern {c1, T1, c2, T2,...cn,Tn} then special patterns emerge. If the finite difference method is used then the block tri-diagonal matrix arises. For other methods other special patterns arise. Solution algorithms are most efficient if these patterns are taken into account in the LU decomposition. Clearly there is a trade-off between the programming time (needed to exploit any special structure) and the value received from a more efficient solution. If the problem is being solved only a few times or the time is extremely small then the special structure need not be exploited and the usual, available packages are suitable. If the problem is being solved thousands of times or the computer time is large (as it is when n is large, say 1000) then the special structure must be exploited. In this example, suppose the block tri-diagonal matrix is composed of blocks that are nb x nb and there are n columns of blocks (along the diagonal). Then the number of multiplications and divisions is (for large n and nb) [Isaacson and Keller, 1966, p. 61]

If n = 100 and nb = 3 then the block tri-diagonal technique takes only 4500 operations while the dense matrix technique takes 9 x 106 operations. Thus the block tri-diagonal technique is 2000 times faster.