Annotation of researchv10dc/cmd/matlab/exec/lanstep, revision 1.1.1.1

1.1       root        1: alfa(j) = q(:,j)'*r;
                      2: r = r - alfa(j)*q(:,j);
                      3: if ort <> 0, for k = 1:j-1, r = r - r'*q(:,k)*q(:,k);
                      4: beta(j) = norm(r);
                      5: q(:,j+1) = r/beta(j);
                      6: r = A*q(:,j+1) - beta(j)*q(:,j);
                      7: if j > 1, T = diag(beta(1:j-1),1); T = diag(alfa) + T + T'; eig(T)

unix.superglobalmegacorp.com

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.