0

Following is the example from this book. My question is, what value of λ did he put in? He did not explain that, can anybody explain how did he get 1.24?

enter image description here

enter image description here

1 Answers1

1

For the matrix $T_j$, he calculates the eigenvalues. We have the characteristic polynomial:

$$\tag 1 \det(T_j - \lambda I) = -\lambda(\lambda^2 -0.625) \implies \lambda_{1,2,3} = 0, \pm \sqrt{0.625}$$

The spectral radius is given by the maximum magnitude eigenvalue:

$$\rho(A) = \displaystyle \max ~ |\lambda_i| $$

From $(1)$, we have:

$$\rho(T_j) = \displaystyle \max ~ |\lambda_i| = \sqrt{0.625}$$

Now, when $A$ is positive definite and tridiagonal, we have:

$$\omega = \dfrac{2}{1 + \sqrt{1 - [\rho(T_j)]^2}}$$

The calculation follows.

There is a write-up of the theory of this choice in Section 3.

Moo
  • 11,311
  • But what if the spectral radius is greater than 1? The expression inside the square root will become negative and hence $\omega$ becomes imaginary? – Aditi Narware Nov 16 '17 at 13:39
  • @AditiNarware: See https://math.stackexchange.com/questions/2390456/matrix-with-spectral-radius-greater-or-equal-1-such-that-fixed-point-iteration-c – Moo Nov 16 '17 at 15:23
  • @AditiNarware: Also see page $7$ of http://www.it.uom.gr/teaching/linearalgebra/ExamplesToIterativeMethods.pdf – Moo Nov 16 '17 at 15:31