|
|
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)
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.