Mon 03 Jan 2011 08:26:40 AM UTC, original submission:
Hello
This seem not to be reported.
In the fntests
c:\usr\Tatsu\mingwhome\octaves\hg\octave-work\src/DLD-FUNCTIONS\qr.cc PASS47/48 FAIL 1
From fntests.log
***** test
[Q,R] = qr(Ac);
[Q,R] = qrinsert(Q,R,3,uc);
assert(norm(vec(Q'Q - eye(5)),Inf) < 1e1eps)
assert(norm(vec(triu(R)-R),Inf) == 0)
assert(norm(vec(QR - [Ac(:,1:2) uc Ac(:,3)]),Inf) < norm(Ac)1e1*eps)
!!!!! test failed
assert (norm (vec (Q * R - [Ac(:, 1:2), uc, Ac(:, 3)]), Inf) < norm (Ac) * 1e1 * eps) failed
shared variables {
A =
0.091364 0.613038 0.999083
0.594638 0.425302 0.603537
0.383594 0.291238 0.085574
0.265712 0.268003 0.238409
0.669966 0.743851 0.445057
u =
0.85082
0.76426
0.42883
0.53010
0.80683
v =
0.98810
0.24295
0.43167
Ac =
0.620405 + 0.956953i 0.480013 + 0.048806i 0.402627 + 0.338171i
0.589077 + 0.658457i 0.013205 + 0.279323i 0.229284 + 0.721929i
0.092758 + 0.345687i 0.928679 + 0.241052i 0.764536 + 0.832406i
0.912098 + 0.721024i 0.049018 + 0.269452i 0.730029 + 0.796517i
0.112849 + 0.603871i 0.486352 + 0.142337i 0.355646 + 0.151496i
uc =
0.20351 + 0.05401i
0.13141 + 0.43708i
0.29808 + 0.08789i
0.69821 + 0.38844i
0.74871 + 0.25821i
vc =
0.85839 + 0.29468i
0.20820 + 0.93090i
0.86184 + 0.34689i
#**********************
I have tried the following , after Ac and uc are set properly.
>> [Q,R] = qr(Ac);
>> [Q,R] = qrinsert(Q,R,3,uc);
>> norm(vec(Q*R - [Ac(:,1:2) uc Ac(:,3)]),Inf)
ans = 4.7103e-016
>> norm (Ac) * 1e1 * eps
ans = 6.0510e-015
>> norm (vec (Q * R - [Ac(:, 1:2), uc, Ac(:, 3)]), Inf) < norm (Ac) * 1e1 * eps
ans = 1
>> assert(norm(vec(QR - [Ac(:,1:2) uc Ac(:,3)]),Inf) < norm(Ac)1e1*eps)
>>
???????
Why fntests did report
"assert (norm (vec (Q * R - [Ac(:, 1:2), uc, Ac(:, 3)]), Inf) < norm (Ac) * 1e1 * eps) failed" ???
|