File: qr-4.m

package info (click to toggle)
octave 2.0.13.95-1
  • links: PTS
  • area: main
  • in suites: slink
  • size: 26,052 kB
  • ctags: 16,420
  • sloc: cpp: 67,184; fortran: 41,514; ansic: 26,607; sh: 7,291; makefile: 4,089; lex: 1,961; yacc: 1,852; perl: 1,676; lisp: 1,664; exp: 123
file content (10 lines) | stat: -rw-r--r-- 286 bytes parent folder | download | duplicates (9) 
1
2
3
4
5
6
7
8
9
10
 
a = [0, 2; 2, 1; 1, 2];

[q, r, p] = qr (a);

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

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