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

1.1       root        1: t = A*q(:,j-1) - alfa(j-1)*p(:,j-1);
                      2:    if ort>0, for i = 1:j-1, t = t - t'*p(:,i)*p(:,i);
                      3: beta(j) = norm(t);
                      4: p(:,j) = t/beta(j);
                      5: t = A'*p(:,j) - beta(j)*q(:,j-1);
                      6:    if ort>0, for i = 1:j-1, t = t - t'*q(:,i)*q(:,i);
                      7: alfa(j) = norm(t);
                      8: q(:,j) = t/alfa(j);

unix.superglobalmegacorp.com

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