Return to qr CVS log

Up to [Research Unix] / researchv10dc / cmd / matlab / exec

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

1.1       root        1: scale = s(m);
                      2: sm = s(m)/scale; smm1 = s(m-1)/scale; emm1 = e(m-1)/scale;
                      3: sl = s(l)/scale; el = e(l)/scale;
                      4: b = ((smm1 + sm)*(smm1 - sm) + emm1**2)/2;
                      5: c = (sm*emm1)**2;
                      6: shift = sqrt(b**2+c); if b < 0, shift = -shift;
                      7: shift = c/(b + shift)
                      8: f = (sl + sm)*(sl-sm) - shift
                      9: g = sl*el
                     10: for k = l: m-1, exec('qrstep',ip)
                     11: e(m-1) = f

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