Basic features with a real symmetric matrix (and normally huge n>106 and sparse)
ˆA of dimension n×n:
- Lanczos' algorithm generates a sequence of real tridiagonal matrices Tk of dimension k×k with k≤n, with the property that the extremal eigenvalues of Tk are progressively better estimates of ˆA' extremal eigenvalues.* The method converges to the extremal eigenvalues.
- The similarity transformation is
ˆT=ˆQTˆAˆQ,
with the first vector
ˆQˆe1=ˆq1.
We are going to solve iteratively
ˆT=ˆQTˆAˆQ,
with the first vector ˆQˆe1=ˆq1.
We can write out the matrix ˆQ in terms of its column vectors
ˆQ=[ˆq1ˆq2…ˆqn].