|
|
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.