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