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