Finding Eigenvectors with Lanczos Algorithm

  • Thread starter Thread starter maverick_starstrider
  • Start date Start date
  • Tags Tags
    Eigenvectors
maverick_starstrider
Messages
1,118
Reaction score
7
Hi, I'm applying the Lanczos algorithm to find the minimal eigenvalue of some huge matrix. Now that I've got it I'm trying to find the eigenvector corresponding to this eigenvalue. Now I have looked through book after book after book and I have yet to find an explanation of how to do this that is even the slightest bit approachable and I don't have the time or patience to read through a whole 300 page textbook to understand the authors notations and terminology well enough to glean a probably 10 line algorithm to generate these things. So I'm wondering if someone could simply tell me a bare-bones algorithm or point me to one that accomplishes this (an algorithm not a linear algebra proof from willoughby and callum or some such).

I have a lanczos algorithm that uses two vectors and the standard lanczos algorithm (http://en.wikipedia.org/wiki/Lanczos_algorithm#The_algorithm). I cannot store more than 2 (maybe 3) vectors (so I can't have something like V=[v1,v2,...,vm]). I generate my tridiagonal lanczos matrix Tm, I can solve for the eigenvalue I want now what do I do? Any help is greatly appreciated, this is driving me mad. It's amazing how many books/internet sources "discuss" this but other then a 5 page linear algebra proof fail to provide any clear explanation, much less an algorithm or an example.
 
Physics news on Phys.org
Did you solve this issue eventually?
 
##\textbf{Exercise 10}:## I came across the following solution online: Questions: 1. When the author states in "that ring (not sure if he is referring to ##R## or ##R/\mathfrak{p}##, but I am guessing the later) ##x_n x_{n+1}=0## for all odd $n$ and ##x_{n+1}## is invertible, so that ##x_n=0##" 2. How does ##x_nx_{n+1}=0## implies that ##x_{n+1}## is invertible and ##x_n=0##. I mean if the quotient ring ##R/\mathfrak{p}## is an integral domain, and ##x_{n+1}## is invertible then...
The following are taken from the two sources, 1) from this online page and the book An Introduction to Module Theory by: Ibrahim Assem, Flavio U. Coelho. In the Abelian Categories chapter in the module theory text on page 157, right after presenting IV.2.21 Definition, the authors states "Image and coimage may or may not exist, but if they do, then they are unique up to isomorphism (because so are kernels and cokernels). Also in the reference url page above, the authors present two...
When decomposing a representation ##\rho## of a finite group ##G## into irreducible representations, we can find the number of times the representation contains a particular irrep ##\rho_0## through the character inner product $$ \langle \chi, \chi_0\rangle = \frac{1}{|G|} \sum_{g\in G} \chi(g) \chi_0(g)^*$$ where ##\chi## and ##\chi_0## are the characters of ##\rho## and ##\rho_0##, respectively. Since all group elements in the same conjugacy class have the same characters, this may be...
Back
Top