File: qr-2.m

package info (click to toggle)
octave2.1 1%3A2.1.73-13
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 37,028 kB
  • ctags: 20,874
  • sloc: cpp: 106,508; fortran: 46,978; ansic: 5,720; sh: 4,800; makefile: 3,186; yacc: 3,132; lex: 2,892; lisp: 1,715; perl: 778; awk: 174; exp: 134
file content (10 lines) | stat: -rw-r--r-- 332 bytes parent folder | download | duplicates (9) 
1
2
3
4
5
6
7
8
9
10
 
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))