QR decomposition: Difference between revisions

Content added Content deleted
m (Changed signValue function)
m (Changed signValue function)
Line 236: Line 236:
R has (sign: R -> Integer) => coerce(sign(r)$R)$R
R has (sign: R -> Integer) => coerce(sign(r)$R)$R
zero? r => r
zero? r => r
-- hack, c.f. complex numbers
if sqrt(r*r) = r then 1 else -1
R has (sqrt: R -> R) =>
if sqrt(r*r) = r then 1 else -1
1
householder(a) ==
householder(a) ==
m := nrows a
m := nrows a