If a is convergent to g then a is summable to g but

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: ’s equation on a square is not diagonally dominant (e.g., look at the fifth row in the 9 × 9 system on p. 564). But such systems are always positive definite (Example 7.6.2), and there is a classical theorem stating that if A is positive definite, then the Gauss–Seidel iteration converges to the solution of Ax = b for every initial vector x(0). The same cannot be said for Jacobi’s method, but there are matrices (the M-matrices of Example 7.10.7, p. 626) having properties resembling positive definiteness for which Jacobi’s method is guaranteed to converge—see (7.10.29). D E Example 7.10.6 The successive overrelaxation (SOR) method improves on Gauss–Seidel by introducing a real number ω = 0, called a relaxation parameter, to form the splitting A = M − N, where M = ω −1 D − L and N = (ω −1 − 1)D + U. As before, D is the diagonal part of A ( aii = 0 is assumed) and −L and −U contain the entries occurring below and above the diagonal of A, respectively. Since M−1 = ω (D − ω L)−1 = ω (I − ω D−1 L)−1 , the SOR iteration matrix is −1 Hω = M −1 N = (D − ω L) T H −1 (1 − ω )D + ω U = (I − ω D IG R −1 L) (1 − ω )I + ω D−1 U , and the k th SOR iterate emanating from (7.10.16) is x(k ) = Hω x(k − 1) + ω (I − ω D−1 L)−1 D−1 b. (7.10.21) This is the Gauss–Seidel iteration when ω = 1. Using ω > 1 is called overrelaxation, while taking ω < 1 is referred to as underrelaxation. Writing (7.10.21) in the form (I − ω D−1 L)x(k ) = (1 − ω )I + ω D−1 U x(k − 1) + ω D−1 b and considering the ith component on both sides of this equality produces Y P O C xi (k ) = (1 − ω )xi (k − 1) + ω bi − aij xj (k ) − aij xj (k − 1) . (7.10.22) aii j <i j >i The matrix splitting approach is elegant and unifying, but it obscures the simple idea behind SOR. To understand the original motivation, write the Gauss–Seidel iterate in (7.10.20) as xi (k ) = xi (k − 1) + ck , where ck is the “correction term” ck = 1 bi − aii n aij xj (k ) − j <i aij xj...
View Full Document

This document was uploaded on 03/06/2014 for the course MA 5623 at City University of Hong Kong.

Ask a homework question - tutors are online