Mercurial > hg > octave-nkf

comparison src/DLD-FUNCTIONS/qr.cc @ 10340:36317747577a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
qr.cc: increase tolerance for test result
author John W. Eaton <jwe@octave.org>
date Fri, 19 Feb 2010 16:31:31 -0500
parents e317791645c4
children 7ed1f2e831ba
comparison
equal deleted inserted replaced
10339:de2d43bcb083 10340:36317747577a
1329 %! 0.694986 + 0.000571i 0.682327 + 0.841712i 0.807537 + 0.166086i 0.192767 + 0.358098i; 1329 %! 0.694986 + 0.000571i 0.682327 + 0.841712i 0.807537 + 0.166086i 0.192767 + 0.358098i;
1330 %! 0.945002 + 0.066788i 0.350492 + 0.642638i 0.579629 + 0.048102i 0.600170 + 0.636938i ] * I; 1330 %! 0.945002 + 0.066788i 0.350492 + 0.642638i 0.579629 + 0.048102i 0.600170 + 0.636938i ] * I;
1331 %! 1331 %!
1332 %! [Q,R] = qr(AA); 1332 %! [Q,R] = qr(AA);
1333 %! [Q,R] = qrdelete(Q,R,3); 1333 %! [Q,R] = qrdelete(Q,R,3);
1334 %! assert(norm(vec(Q'*Q - eye(5)),Inf) < 1e1*eps) 1334 %! assert(norm(vec(Q'*Q - eye(5)),Inf) < 16*eps)
1335 %! assert(norm(vec(triu(R)-R),Inf) == 0) 1335 %! assert(norm(vec(triu(R)-R),Inf) == 0)
1336 %! assert(norm(vec(Q*R - [AA(:,1:2) AA(:,4)]),Inf) < norm(AA)*1e1*eps) 1336 %! assert(norm(vec(Q*R - [AA(:,1:2) AA(:,4)]),Inf) < norm(AA)*1e1*eps)
1337 %! 1337 %!
1338 %!test 1338 %!test
1339 %! AA = [0.091364 0.613038 0.027504 0.999083; 1339 %! AA = [0.091364 0.613038 0.027504 0.999083;