[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #32035] qr.cc test error on octave-dev source
From: |
Tatsuro MATSUOKA |
Subject: |
[Octave-bug-tracker] [bug #32035] qr.cc test error on octave-dev source (3.3.54) on MInGW |
Date: |
Mon, 03 Jan 2011 08:26:46 +0000 |
User-agent: |
Mozilla/5.0 (Windows; U; Windows NT 5.1; en-US) AppleWebKit/534.10 (KHTML, like Gecko) Chrome/8.0.552.224 Safari/534.10 |
URL:
<http://savannah.gnu.org/bugs/?32035>
Summary: qr.cc test error on octave-dev source (3.3.54) on
MInGW
Project: GNU Octave
Submitted by: tmacchant
Submitted on: 2011年01月03日 08時26分40秒
Category: None
Severity: 3 - Normal
Priority: 5 - Normal
Item Group: Incorrect Result
Status: None
Assigned to: None
Originator Name:
Originator Email:
Open/Closed: Open
Discussion Lock: Any
Release: dev
Operating System: Microsoft Windows
_______________________________________________________
Details:
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) < 1e1*eps)
assert(norm(vec(triu(R)-R),Inf) == 0)
assert(norm(vec(Q*R - [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(Q*R - [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" ???
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?32035>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Octave-bug-tracker] [bug #32035] qr.cc test error on octave-dev source (3.3.54) on MInGW,
Tatsuro MATSUOKA <=