Mercurial > octave-nkf

view test/octave.test/linalg/qr-2.m @ 4720:e759d01692db ss-2-1-53

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2004-01-23 04:13:37 by jwe]
author jwe
date Fri, 23 Jan 2004 04:13:37 +0000
parents 873b1ef30534
children ace8d8d26933
line wrap: on
line source

a = [0, 2, 1; 2, 1, 2];

[q, r, p] = qr (a);  # not giving right dimensions. XXX FIXME XXX

[qe, re, pe] = qr (a, 0);

(size (q) == [2, 2] && size (r) == [2, 3] && size (p) == [3, 3]
 && abs (q * r - a * p) < sqrt (eps)
 && size (qe) == [2, 2] && size (re) == [2, 3] && size (pe) == [1, 3]
 && abs (qe * re - a(:,pe)) < sqrt (eps))