Revision as of 22:35, 11 March 2012 (view source) Mark Clements (talk \| contribs) (Revise sign and AlgebraicNumber output) ← Older edit		Revision as of 22:40, 11 March 2012 (view source) Mark Clements (talk \| contribs) (Remove comments) Newer edit →
Line 242: householder(a:Matrix(R)):Matrix(R) == m := nrows a s := coerce(sign(a(1,1)))$R ~~-- OrderedRing~~ e := coerce(unitVector(m))$Matrix(R) u := a + length(column(a,1))se ~~-- RadicalCategory~~ v := u/u(1,1) ~~-- Field~~ beta := (1+1)/((transpose(v)v)(1,1)) scalarMatrix(m,1) - betav*transpose(v)

QR decomposition: Difference between revisions

QR decomposition (view source)