|
|
1.1 root 1: // subexec rot(f,g,cs,sn)
2: rho = g; if abs(f) > abs(g), rho = f;
3: cs = 1.0; sn = 0.0; z = 1.0;
4: r = norm(<f g>); if rho < 0, r = -r; r
5: if r <> 0.0, cs = f/r
6: if r <> 0.0, sn = g/r
7: if abs(f) > abs(g), z = sn;
8: if abs(g) >= abs(f), if cs <> 0, z = 1/cs;
9: f = r;
10: g = z;
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.