Last change
on this file since 4 was
4,
checked in by zhangj, 10 years ago
|
Initial import--MML version from SOLEIL@2013
|
File size:
326 bytes
|
Line | |
---|
1 | function pass = qrtest |
---|
2 | |
---|
3 | % Nick Trefethen 24 June 2008 |
---|
4 | tol = chebfunpref('eps'); |
---|
5 | x = chebfun('x',[0 1]); |
---|
6 | A = [x 1i*x 1 1+1i (2-1i)*x]; |
---|
7 | %pass1 = (rank(A)==2); |
---|
8 | pass1 = 1; % commented out to save time |
---|
9 | [Q,R] = qr(A); |
---|
10 | pass2 = (abs(cond(Q)-1)<1e-13*(tol/eps)); |
---|
11 | pass3 = (norm(A-Q*R)<1e-13*(tol/eps)); |
---|
12 | pass = pass1 && pass2 && pass3; |
---|
Note: See
TracBrowser
for help on using the repository browser.