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

1.1       root        1: alfa = <>; beta = 0; q = <>; p = p(:,1)/norm(p(:,1));
                      2: t = A'*p(:,1);
                      3: alfa(1) = norm(t);
                      4: q(:,1) = t/alfa(1);
                      5: X = p(:,1)*(alfa(1)*q(:,1))'
                      6: e(1) = norm(A-X,1)
                      7: for j = 2:r, exec('pdqstep',ip); ...
                      8:    X = X + p(:,j)*(alfa(j)*q(:,j)+beta(j)*q(:,j-1))', ...
                      9:    e(j) = norm(A-X,1)

unix.superglobalmegacorp.com

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