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